ex8_11 - Cto=.1; Ca=Cto*x(1)/ft; Cb=Cto*x(2)/ft;...

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

View Full Document Right Arrow Icon
% "ex8_11" % m-file to solve example 8.11 % % x(1)=fa % x(2)=fb % x(3)=fc % x(4)=T % xdot(1)=d(fa)/dv, xdot(2)=d(fb)/dv, xdot(3)=d(fc)/dv, etc. % function xdot=ex8_11(v,x) f ft=x(1)+x(2)+x(3); k1a=10*exp(4000*(1/300-1/x(4))); k2a=.09*exp(9000*(1/300-1/x(4)));
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Cto=.1; Ca=Cto*x(1)/ft; Cb=Cto*x(2)/ft; Cc=Cto*x(3)/ft; r1a=-k1a*Ca; r2a=-k2a*Ca^2; r xdot(1,:)=r1a+r2a; xdot(2,:)=-r1a; xdot(3,:)=-r2a/2; xdot(4,:)=(4000*(373-x(4))+(-r1a)*20000+(-r2a)*60000)/(90*x(1)+90*x(2)+180*x(3));...
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