% Plotter her kompressibilitet fra 1 til 300 bar beregnet med Peng-Robinson
% og Van der Waals tilstandsligning.

% Gassblandingen som plottes har sammensetting som gitt av vektor x
% x = [0.15; 0.55; 0.30] betyr 15% O2, 55% He, 30% N2.

% x = [0.21; 0.35; 1-0.21-0.35]; % Trimix 15/55
x = [0.21; 0; 1-0.21]; % Luft
% x = [0.32; 0; 1-0.32]; % Nitrox 32%
% x = [0; 1; 0]; % ren He
%  x = [1; 0; 0]; % ren O2
% x = [0; 0; 1]; % ren N2
% x = [0.21; 0.35; 1-0.21-0.35];

Trange = [10;20;30;40]+273;
V = 24e-3;

prange = linspace(1e5, 300e5,50);
ni = length(prange);
Ti = length(Trange);

pv = zeros(Ti,ni);
pr = pv;
pideal = pv;

for i = 1:ni
		for j = 1:Ti
			nideal = prange(i)*V/8.31/Trange(j);
			pv(j,i) = pvdw(Trange(j), V, x*nideal);
			pr(j,i) = prfun(Trange(j), V, x*nideal);
	end
end

plot( ...
prange,pv(1,:)./prange, 'og', ...
prange,pr(1,:)./prange, 'xg', ...
prange,pv(2,:)./prange, 'ob',...
prange,pr(2,:)./prange, 'xb',...
prange,pv(3,:)./prange, 'or', ...
prange,pr(3,:)./prange, 'xr', ...
prange,pv(4,:)./prange, 'om', ...
prange,pr(4,:)./prange, 'xm' ...
)
titlestr = sprintf('Avvik fra idealitet for en trimix %.0f/%.0f', ...
     round(x(1:2)*100));
title(titlestr);
xlabel('Tykk [Pa]');
ylabel('Kompressibilitet');
legend('10°C VdW', '10°C PR', '20°C VdW', '20°C PR', ...
	'30°C VdW', '30°C PR', '40°C VdW', '40°C PR');
