R = 8.314;
T0 = 273.16;
Tc = 647;
p0 = 611;
alpha = 59000;
beta = 50;
T = [274:1:647];
p = p0.*(T0./T).^(beta/R).*exp((alpha/R).*(1/T0-1./T));
figure(1);
plot(T,p*1e-5);
%grid on;
xlabel('Temperatur (K)');
ylabel('Damptrykk (bar)');
latentfull = 54000.*(1-T/Tc).^(1/3);
latentlin = alpha - beta.*T;
figure(2);
plot(T,latentlin/1000,'--',T,latentfull/1000,'linewidth',2);
%axis([240 680 0 48.000]);
xlabel('Temperatur (K)');
ylabel('Molar fordampningsvarme (kJ/mol)');
handle=legend('Lineær modell','54(1-T/647)^{1/3}');
set(handle, 'Box', 'off');
%%grid on;