- #1

- 20

- 0

Please help me on this little problem.

I have the code below.

for k=40:10:80

T(k)=273.15+k;

z=[0.2 0.2 0.2 0.4];

W_PR=0.245;

C=4;

omega=[0.344 0.467 0.578 0.789];

Tc=[600 700 500 570];

Pc=[50 70 58 76];

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);

end

for c=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));

end

for c=1:C

Z(c,c)=A_PR(c,c)./5;

V(c)=B_PR(c).*6;

end

end

QUESTION: Each time I run the code, I want result for Z and V at each T(k). The code as it is only gives result for one T value though I want it to run the loop and give result for Z and V for all T(k).

Thanks

adeeyo