cd13_1c1 -...

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

View Full Document Right Arrow Icon
% "cde13_1c1" % m-file to solve CD example 13.1 Part c.1 % % 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=ex13_9c1(t,x) f cao=1; k1=1; k2=1; Conv=(cao-x(1))/cao; E1=-2.104*t^4+4.167*t^3-1.596*t^2+0.353*t-0.004;
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: E2=-2.104*t^4+17.037*t^3-50.247*t^2+62.964*t-27.402; ra=-k1*x(1); rc=k2*x(2); rb=k1*x(1)-k2*x(2); r if t <=1.26 E=E1; else E=E2; 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