ex9_8 - rc=3*k2b*x(2); rb=k1a*x(1)/2-k2b*x(2); r

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

View Full Document Right Arrow Icon
% "ex9_8" % m-file to solve example 9.8 % % x(1)=Ca % x(2)=Cb % x(3)=Cc % x(4)=T % xdot(1)=d(Ca)/dt, xdot(2)=d(Cb)/dt, xdot(3)=d(Cc)/dt, etc. % function xdot=ex9_8(t,x) f Ca0=4; v0=240; k1a=1.25*exp((9500/1.987)*(1/320-1/x(4))); k2b=0.08*exp((7000/1.987)*(1/290-1/x(4))); ra=-k1a*x(1); V=100+v0*t;
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: rc=3*k2b*x(2); rb=k1a*x(1)/2-k2b*x(2); r xdot(1,:)=ra+(Ca0-x(1))*v0/V; xdot(2,:)=rb-x(2)*v0/V; xdot(3,:)=rc-x(3)*v0/V; xdot(4,:)=(35000*(298-x(4))-Ca0*v0*30*(x(4)-305)+((-6500)*(-k1a*x(1))+(8000)*(-k2b*x(2)))*V)/((x(1)*30+x(2)*60+x(3)*20)*V+100*35); k...
View Full Document

Ask a homework question - tutors are online