ex4_11 - cb0=0.025 rate=k*x(1*x(2 ca0=0.05 c...

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

View Full Document Right Arrow Icon
% "ex4_11" % m-file to solve example 4.11 % % x(1)=ca % x(2)=cb % x(3)=cc % x(4)=cd % xdot(1)=d(ca)/dt, xdot(2)=d(cb)/dt, xdot(3)=d(cc)/dt, xdot(4)=d(cd)/dt % function xdot=ex4_11(t,x) f v0=5; v00=0.05; v=v0+v00*t; k=2.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: cb0=0.025; rate=k*x(1)*x(2); ca0=0.05; c xdot(1,:)=-k*x(1)*x(2)-v00*x(1)/v; xdot(2,:)=-k*x(1)*x(2)+v00*(cb0-x(2))/v; xdot(3,:)=k*x(1)*x(2)-v00*x(3)/v; xdot(4)=k*x(1)*x(2)-v00*x(4)/v;...
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