clc
y=[0 .05 .1 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6 .65 .7 .75 .8 .85 .9 .95 1];
a1=.0592; b1=0.0181;
a2=.1491; b2=0.0347;
P=20;
T=400;
R=.08206;
for i=1:21
A(i)=y(i)^2*a1+(1-y(i))^2*a2+2*y(i)*(1-y(i))*(a1*a2)^.5;
B(i)=y(i)*b1+(1-y(i))*b2;
Aprime(i)=2*y(i)*a1-2*(1-y(i))*a2+2*(1-2*y(i))*(a1*a2)^.5;
Bprime=b1-b2;
end
for i=1:21
b=A(i)-B(i)-B(i)^2;
c=A(i)*B(i);
K=[1 -1 b -c];
m=roots(K);
Z(i)=m(1);
end
for i=1:21
Zprime(i)=((B(i)-Z(i))*Aprime(i)+(A(i)+Z(i)+2*B(i)*Z(i))*Bprime)/(3*Z(i)^2-2*Z(i)+A(i)-B(i)-B(i)^2);
V(i)=Z(i)*R*T/P;
Vbar1(i)=R*T*(Z(i)+(1-y(i))*Zprime(i))/P;
Vbar2(i)=R*T*(Z(i)-y(i)*Zprime(i))/P;
end
table1(:,1)=y';
table1(:,2)=Z';
table1(:,3)=V';
table1(:,4)=Vbar1';
table1(:,5)=Vbar2';
disp(' Y CO2 Z V V CO2 bar V C3H6 bar')
disp(table1)
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA
%BY HRA
%BY HRA
%BY HRA
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA
%BY HRA
%BY HRA
%BY HRA
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA
%BY HRA
%BY HRA
%BY HRA
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA
%BY HRA
%BY HRA
%BY HRA
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA
%BY HRA
%BY HRA
%BY HRA
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA
%BY HRA
%BY HRA
%BY HRA
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA
%BY HRA
%BY HRA
%BY HRA
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%BY HRA IRAN PERSIAN GROUP
%B