function f = pb(T)
f = 10^(-1686.8/((T-32)/1.8+273)+7.6546)/760;
