Cd13_1c2 - 2265.62125*t 925.46463 E3=0.0041*t^4-0.07593*t^3 0.52276*t^2-1.59457*t 1.84445 E4=0 ra=-k1*x(1 rc=k2*x(2 rb=k1*x(1-k2*x(2 r if

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

View Full Document Right Arrow Icon
% "cde13_1c2" % m-file to solve CD example 13.1 Part c.2 % % x(1)=ca % x(2)=cb % x(3)=cc % x(4)=cabar % x(5)=cbbar % x(6)=ccbar % xdot(1)=d(ca)/dt, xdot(2)=d(cb)/dt, etc. % function xdot=cde13_1c2(t,x) f cao=1; k1=1; k2=1; E1=0.47219*t^4-1.30733*t^3+0.31723*t^2+0.85688*t+0.20909; E2=3.83999*t^6-58.16185*t^5+366.20970*t^4-1224.66963*t^3+2289.84857*t^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*t+925.46463; E3=0.0041*t^4-0.07593*t^3+0.52276*t^2-1.59457*t+1.84445; E4=0; ra=-k1*x(1); rc=k2*x(2); rb=k1*x(1)-k2*x(2); r if t <=1.82 E=E1; elseif t<=2.8 E=E2; elseif t<6 E=E3; else E=E4; end e xdot(1,:)=ra; xdot(2,:)=rb; xdot(3,:)=rc; xdot(4,:)=x(1)*E; xdot(5,:)=x(2)*E; xdot(6,:)=x(3)*E;...
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