{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

cde9_2 - tau=.06 V=tau*vol cao=Fao/vol cbo=Fbo/vol...

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

View Full Document Right Arrow Icon
% "cde9_2" % m-file to solve CD example 9.2 % % x(1)=ca % x(2)=cb % x(3)=cc % x(4)=cm % x(5)=T % xdot(1)=d(ca)/dt, xdot(2)=d(cb)/dt, xdot(3)=d(cc)/dt, etc. % function xdot=cde9_1(t,x) f UA=55000; Ta=85; Fao=300; To=68; k=16.96e12*exp(-32400/1.987/(x(5)+460)); Fbo=2000; Fco=0; Fmo=20; ra=-k*x(1); rb=-k*x(1); rc=k*x(1); Thetacp=35+Fbo/Fao*18+Fmo/Fao*19.5; vol=Fao*58.08/(.859/454*2.54^3*12^3)+Fbo*18.02/(.9941/454*2.54^3*12^3)+Fco*76.11/ (1.036/454*2.54^3*12^3)+Fmo*32.04/(.7914/454*2.54^3*12^3);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: tau=.06; V=tau*vol; cao=Fao/vol; cbo=Fbo/vol; cco=Fco/vol; cmo=Fmo/vol; Na=x(1)*V; Nb=x(2)*V; Nc=x(3)*V; Nm=x(4)*V; NCp=Na*35+Nb*18+Nc*46+Nm*19.5; N xdot(1,:)=1/tau*(cao-x(1))+ra; xdot(2,:)=1/tau*(cbo-x(2))+rb; xdot(3,:)=1/tau*(0-x(3))+rc; xdot(4,:)=1/tau*(cmo-x(4)); xdot(5,:)=(UA*(Ta-x(5))-Fao*Thetacp*(x(5)-To)+(-36000)*ra*V)/NCp;...
View Full Document

{[ snackBarMessage ]}