function h=latent_heat(t)
a=[0.014615072315;4.20542536016;-7.03894624e-4;7.94732605e-6;-3.1034444e-8;7.9997e-11];
x=rot90(a,2);
h1=polyval(x,t);
b=[2500.904774201;1.83443741725;-3.47293054e-4;-2.0969031e-6;-3.5911751e-8;3.8508e-11];
y=rot90(b,2);
h2=polyval(y,t);
h=h2-h1;
