function y=qwat(y)
x=y(1);
p=y(2);
es=6.112.*exp(17.502.*x./(x+241.0))*(1.0007+3.46e-6*p);
y=es*622./(1010.0-.378*es);