ex13_9d - rc=k1*x(2)*x(1); rb=-k1*x(1)*x(2)-k3*x(2)*x(5);

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

% "ex13_9d" % m-file to solve example 13.9 Part 4 % % x(1)=ca % x(2)=cb % x(3)=cc % x(4)=F % x(5)=cd % x(6)=ce % xdot(1)=d(ca)/dz, xdot(2)=d(cb)/dz, etc. % function xdot=ex13_9d(z,x) f cao=1; cbo=1; cco=0; cdo=0; ceo=0; lam=2.52-z; k1=1; k2=1; k3=1; E1=-2.104*lam^4+4.167*lam^3-1.596*lam^2+0.353*lam-0.004; E2=-2.104*lam^4+17.037*lam^3-50.247*lam^2+62.964*lam-27.402; ra=-k1*x(1)*x(2)-k2*x(1);
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: rc=k1*x(2)*x(1); rb=-k1*x(1)*x(2)-k3*x(2)*x(5); rd=k2*x(1)-k3*x(2)*x(5); re=k3*x(2)*x(5); r if lam &lt;=1.26 E=E1; else E=E2; end e EF=E/(1-x(4)); E xdot(1,:)=-(-ra+(x(1)-cao)*EF); xdot(2,:)=-(-rb+(x(2)-cbo)*EF); xdot(3,:)=-(-rc+(x(3)-cco)*EF); xdot(4,:)=-E; xdot(5,:)=-(-rd+(x(5)-cdo)*EF); xdot(6,:)=-(-re+(x(6)-ceo)*EF);...
View Full Document

This note was uploaded on 09/02/2011 for the course CHEM 102 taught by Professor Kmn during the Spring '11 term at Babson College.

Ask a homework question - tutors are online