0001 function s=qsee(y) 0002 x=y(:,1); 0003 p=y(:,2); 0004 es=6.112.*exp(17.502.*x./(x+240.97))*.98.*(1.0007+3.46e-6*p); 0005 s=es*622./(p-.378*es);