ex13_8 - 0.000865652/2*lam^2+0.028004*lam;

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

View Full Document Right Arrow Icon
% "ex13_8" % m-file to solve example 13.8 % % x(1)=X % xdot(1)=d(X)/dz % function xdot=ex13_8(z,x) f cao=8; k=0.01; lam=200-z; ca=cao*(1-x(1)); E1=4.44658e-10*lam.^4-1.1802e-7*lam^3+1.35358e-5*lam^2-0.000865652*lam+0.028004; E2=-2.64e-9.*lam.^3+1.3618e-6*lam^2-0.00024069*lam+0.015011; F1=4.44658e-10/5*lam^5-1.1802e-7/4*lam^4+1.35358e-5/3*lam^3-
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 0.000865652/2*lam^2+0.028004*lam; F2=-(-9.30769e-8*lam^3+5.02846e-5*lam^2-0.00941*lam+0.618231-1); ra=-k*ca^2; r if lam <=70 E=E1; else E=E2; end e if lam <=70 F=F1; else F=F2; end EF=E/(1-F); E xdot(1,:)=-(ra/cao+E/(1-F)*x(1)); x...
View Full Document

Ask a homework question - tutors are online