cd13_1c4 - 2265.62125*lam+925.46463;

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

View Full Document Right Arrow Icon
% "cde13_1c4" % m-file to solve CD example 13.1 Part c.4 % % 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_1c4(z,x) f cao=1; cbo=0; cco=0; lam=6-z; k1=1; k2=1; E1=0.47219*lam^4-1.30733*lam^3+0.31723*lam^2+0.85688*lam+0.20909; E2=3.83999*lam^6-58.16185*lam^5+366.20970*lam^4-1224.66963*lam^3+2289.84857*lam^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: 2265.62125*lam+925.46463; E3=0.0041*lam^4-0.07593*lam^3+0.52276*lam^2-1.59457*lam+1.84445; ra=-k1*x(1); rc=k2*x(2); rb=k1*x(1)-k2*x(2); r if lam <=1.82 E=E1; elseif lam <=2.8 E=E2; else E=E3; 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