function f = pz(T)
f = 10^(-3667.5/(T+459.69)+4.858);
