% Entalpi.m
% Finner entalpien i et punkt

p0 = 20e5     % Pa
T0 = 251;       % K

p = p0;
T = T0;
V = 5e-5;
n = [1];
dV = V;

% Newtoniterasjoner til V ved ønsket trykk finnes:
while abs(dV)>1e-6*V
    A = co2([T,V,n]');
    dV = (p0 + A.g(2))/(-A.H(2,2))
    V = V + dV;
    pause
end


h =-A.g(1)*T + A.g(3)*n   % Entalpi

