% function es(T)
function y=es(T)
a=17.27;
b=35.86;
y=3.8/0.62197*exp(a*((T-273)/(T-b)));
