clear
clc
//Given
//Volumetric flow rates(litre/hr)
vo=[10;3;1.2;0.5];
//Concentrations (millimol/litre)
CA=[85.7;66.7;50;33.4];
CAo=100;
//Volume(litre)
V=0.1;
//For the stoichiometry 2A-->R
//Expansion factor is
e=(1-2)/2;
//Initialization
XA=zeros(4,1);
rA=zeros(4,1);
//Relation between concentration and conversion
for i=1:4
XA(i)=(1-CA(i)/CAo)/(1+e*CA(i)/CAo);
//Rate of reaction is given by
rA(i)=vo(i)*CAo*XA(i)/V;
//Testing nth order kinetics
//-rA=k*CA^n
//log(-rA)=logk+nlog(CA)
m(i)=log10(CA(i));
n(i)=log10(rA(i));
end
//For nth order plot between n & m should give a straight line
plot(m,n)
coefs=regress(m,n);
printf("Intercept of the graph is %f\n",coefs(1))
printf("Slope of the graph is %f\n",coefs(2))
k=10^coefs(1)
n=coefs(2)
printf("\n Taking n=2,rate of equation(millimol/litre.hr) is %f",k)
printf("CA^2 \n")
disp('The sol slightly differ from that given in book because regress fn is used to calculate the slope')