for k=40:10:80 T(k)=273.15+k; for c=1:C x_PR(1,c)=z(c)/(1+W_PR*(K_PR(c)-1)); x_PR(2,c)=K_PR(c)*x_PR(1,c); kappa_PR=0.37464+1.54226.*omega(c)-0.26992.*omega(c).^2; alpha_PR=(1+kappa_PR.*(1-sqrt(T(k)./Tc(c)))).^2; a_PR(c,c)=0.45724.*R.^2.*Tc(c).^2./Pc(c).*alpha_PR; b_PR(c)=0.07780*R.*Tc(c)./Pc(c); end for c=2:C for n=1:(c-1) a_PR(c,n)=sqrt(a_PR(c,c).*a_PR(n,n)); a_PR(n,c)=a_PR(c,n); end end for c=1:C A_PR(c,c)=a_PR(c,c).*P./(R.*T(k)).^2; B_PR(c)=b_PR(c).*P./(R.*T(k)); Z(c,c)=A_PR(c,c)./5; V(c)=B_PR(c).*6; end for c=1:C Aa_PR(c,c)=a_PR(c,c).*P./(R.*T(k)).^3; Ba_PR(c)=b_PR(c).*P./(R.*T(k)); Za(c,c)=A_PR(c,c)./8; Va(c)=B_PR(c).*9; end