ex8_7_b - term=-1.27*10^(-6)*(x(2)^3-298^3);

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

View Full Document Right Arrow Icon
% "ex8_7_b" % m-file to solve example 8.7 part b % % x(1)=X % x(2)=T % xdot(1)=d(X)/dv, xdot(2)=d(T)/dv % function xdot=ex8_7_b(v,x) f cao=18.8; To=1035; ra=-cao*3.58*exp(34222*(1/To-1/x(2)))*(1-x(1))*(To/x(2))/(1+x(1)); cpa=26.6+0.183*x(2)-0.0000459*x(2)^2; dcp=6.8-0.0115*x(2)-3.01*10^(-6)*x(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: term=-1.27*10^(-6)*(x(2)^3-298^3); dh=80770+6.8*(x(2)-298)-0.00575*(x(2)^2-298^2)+term; fao=0.0376; ua=16500; Ta=1150; T xdot(1,:)=-ra./fao; xdot(2,:)=(ua.*(Ta-x(2))+ra.*dh)./(fao*(cpa+x(1).*dcp));...
View Full Document

Ask a homework question - tutors are online