ex8_6 - k=31.1*exp(7906*(t1-360)/(t1*360));...

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

View Full Document Right Arrow Icon
% "ex8_6" % m-file to solve example 8.6 % % x(1)=X % xdot(1)=d(X)/dv % function xdot=ex8_6(v,x) f global kc cao g cao=9.3; t1=330+43.3*x(1); kc=3.03*exp(-830.3*(t1-333)/(333*t1));
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: k=31.1*exp(7906*(t1-360)/(t1*360)); ra=-k*cao*(1-(1+1/kc)*x(1)); fao=0.9*163; xe=kc/(1+kc); x xdot(1,:)=-ra/fao; 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