%%FY1005/TFY4165, Oving 9, Oppgave 5
%%
%%Velg utskrift med mange siffer
format long;
%%R = gasskonstanten = 8.314 J/mol*K
R=8.314;
%%T = absolutt temperatur (K)
T = 118;
%%Verdier for a og b for luft: 1.368  0.0367
%%Enheter: [a] = bar*(L/mol)^2 og [b] = L/mol
a=1.368;
b=0.0367;
%%Gjetter en startverdi for vaeskevolumet V1=0.065 L/mol
V1 = 0.06494246177;
%%Tilhorende trykkverdi:
p1 = (R*T/(V1-b))/100 - a/(V1*V1);
%%Bestem gassvolum V2 som gir samme trykkverdi p1 ved aa gaa i loop
%%fra en volumverdi som helt sikkert er storre enn V2
V2 = 0.2768;
dV = 0.00000000001;
p2 = (R*T/(V2-b))/100 - a/(V2*V2);
while(p2<p1)
    V2 = V2 - dV;
    p2 = (R*T/(V2-b))/100 - a/(V2*V2);                         
end;
p1,p2
V1,V2
%%Beregn areal under p(V) mellom V1 og V2
p = @(v)(R*T./(v-b))/100 - a./(v.*v);
areal = quad(p,V1,V2)
%%Beregn p1*(V2-V1)
rektangel = p1*(V2-V1)
%%Beregn differansen mellom areal og rektangel
differanse = areal - rektangel