ex8_7 - delcp=6.8-0.0115*x(2)-3.81*10^(-6)*x(2)^2;

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
% "ex8_7" % m-file to solve example 8.7 % % x(1)=X % x(2)=T % xdot(1)=d(X)/dv, xdot(2)=d(T)/dv % function xdot=ex8_7(v,x) f cao=18.8; To=1035; Tr=298; k=8.2*10^14*exp(-34222/x(2)); ra=-k*cao*(1-x(1))/(1+x(1))*To/x(2); cpa=26.63+0.183*x(2)-0.00004586*x(2)^2;
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: delcp=6.8-0.0115*x(2)-3.81*10^(-6)*x(2)^2; deltaH=80770+6.8*(x(2)-Tr)-0.00575*(x(2)^2-Tr^2)-1.27*10^(-6)*(x(2)^3-Tr^3); fao=38.3; f xdot(1,:)=-ra/fao; xdot(2,:)=-ra*(-deltaH)/(fao*cpa+x(1)*delcp);...
View Full Document

Ask a homework question - tutors are online