0001 function y=qsat(y) 0002 x=y(:,1);%temp 0003 p=y(:,2);%pressure 0004 es=6.112.*exp(17.502.*x./(x+241.0)).*(1.0007+3.46e-6*p); 0005 y=es*622./(p-.378*es);