cd13_1c3 - ra=-k1*x(1); rc=k2*x(2); rb=k1*x(1)-k2*x(2); r...

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

View Full Document Right Arrow Icon
% "cde13_1c3" % m-file to solve CD example 13.1 Part c.3 % % x(1)=ca % x(2)=cb % x(3)=cc % x(4)=F % xdot(1)=d(ca)/dz, xdot(2)=d(cb)/dz, etc. % function xdot=cde13_1c3(z,x) f cao=1; cbo=0; cco=0; lam=2.52-z; k1=1; k2=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;
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: ra=-k1*x(1); rc=k2*x(2); rb=k1*x(1)-k2*x(2); r if lam <=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; x...
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