{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

cd13_1c2 - 2265.62125*t 925.46463 E3=0.0041*t^4-0.07593*t^3...

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

View Full Document Right Arrow Icon
% "cde13_1c2" % m-file to solve CD example 13.1 Part c.2 % % x(1)=ca % x(2)=cb % x(3)=cc % x(4)=cabar % x(5)=cbbar % x(6)=ccbar % xdot(1)=d(ca)/dt, xdot(2)=d(cb)/dt, etc. % function xdot=cde13_1c2(t,x) f cao=1; k1=1; k2=1; E1=0.47219*t^4-1.30733*t^3+0.31723*t^2+0.85688*t+0.20909; E2=3.83999*t^6-58.16185*t^5+366.20970*t^4-1224.66963*t^3+2289.84857*t^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: 2265.62125*t+925.46463; E3=0.0041*t^4-0.07593*t^3+0.52276*t^2-1.59457*t+1.84445; E4=0; ra=-k1*x(1); rc=k2*x(2); rb=k1*x(1)-k2*x(2); r if t <=1.82 E=E1; elseif t<=2.8 E=E2; elseif t<6 E=E3; else E=E4; end e xdot(1,:)=ra; xdot(2,:)=rb; xdot(3,:)=rc; xdot(4,:)=x(1)*E; xdot(5,:)=x(2)*E; xdot(6,:)=x(3)*E;...
View Full Document

{[ snackBarMessage ]}