function T=Tsat(P)
x0 = 273;
sc.x = x0;
sc.res = Psat(x0);
x0 = x0./sc.x;
options = optimset('TolFun', 1e-10);
x = fsolve(@resPsat, x0, options, P, sc);
T = x.*sc.x;
