function f = pt(T)
f = 10^(-3179.0/(T+459.69)+4.60);
