function [f]=dcA(cA,T)
global cAf k0 F Vf=F/V*(cAf-cA)-k0*cA*(1+0.5*cA)*exp(27.631-9260/T);
function [g]=dT(cA,T)
global k0 Tf Cps H F Vg=F/V*(Tf-T)-H/Cps*k0*cA*(1+0.5*cA)*exp(27.631-9260/T);
function[cA,T]=Runge(cA,T,h)
k1=dcA(cA,T);
k2=dcA(cA+h/2*k1,T);
k3=dcA(cA+h/2*k2,T);
k4=dcA(cA+h*k3,T);
f1=dT(cA,T);
f2=dT(cA,T+h/2*f1);
f3=dT(cA,T+h/2*f2);
f4=dT(cA,T+h*f3);
cA=cA+h/6*(k1+2*k2+2*k3+k4);
T=T+h/6*(f1+2*f2+2*f3+f4);
clearclch=0.001;
tf=1000;
t0=0;
cAi=0.56;
Ti=305;
for t=t0:h:tfcA=cAi;
T=Ti;[cAi,Ti]=Rung(cA,T,h);
[cA,T]
if(abs(cA-cAi)<1e-6)&&(abs(T-Ti)<1e-6)
break;
end
end
global cAf k0 F Vf=F/V*(cAf-cA)-k0*cA*(1+0.5*cA)*exp(27.631-9260/T);
function [g]=dT(cA,T)
global k0 Tf Cps H F Vg=F/V*(Tf-T)-H/Cps*k0*cA*(1+0.5*cA)*exp(27.631-9260/T);
function[cA,T]=Runge(cA,T,h)
k1=dcA(cA,T);
k2=dcA(cA+h/2*k1,T);
k3=dcA(cA+h/2*k2,T);
k4=dcA(cA+h*k3,T);
f1=dT(cA,T);
f2=dT(cA,T+h/2*f1);
f3=dT(cA,T+h/2*f2);
f4=dT(cA,T+h*f3);
cA=cA+h/6*(k1+2*k2+2*k3+k4);
T=T+h/6*(f1+2*f2+2*f3+f4);
clearclch=0.001;
tf=1000;
t0=0;
cAi=0.56;
Ti=305;
for t=t0:h:tfcA=cAi;
T=Ti;[cAi,Ti]=Rung(cA,T,h);
[cA,T]
if(abs(cA-cAi)<1e-6)&&(abs(T-Ti)<1e-6)
break;
end
end