%trykk.m
%Beregner metningstrykk for H2O ved en gitt temperatur T[=]K
%Utrykket er hentet fra [2] som har hentet det fra [10]

function P=trykk(T)

Tc_H2O=647.3;           %K, kritisk temp. 
Pc_H2O=221.2*10^5;      %Pa, kritisk trykk 
x=1-T/Tc_H2O;

VPA=-7.76451;
VPB=1.45838;
VPC=-2.77580;
VPD=-1.23303;

P1=(1-x)^-1*(VPA*x+VPB*x^1.5+VPC*x^3+VPD*x^6);

P=Pc_H2O*exp(P1);

return