ex7_3 - k xdot(1,:)=-k1*x(1)-k2*x(1)*x(2)-k4*x(1)*x(6);

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

View Full Document Right Arrow Icon
% "ex7_3" % m-file to solve example 7.3 % % x(1)=C1 % x(2)=C2 % x(3)=C3 % x(4)=C4 % x(5)=C5 % x(6)=C6 % x(7)=C7 % x(8)=C8 % x(9)=CP5 % x(10)=CP1 % xdot(1)=dC1/dt, xdot(2)=dC2/dt, etc. % function xdot=ex7_3(t,x) f T=1000; k1=10*exp((87500/1.987)*(1/1250-1/T)); k2=8450000*exp((13000/1.987)*(1/1250-1/T)); k3=3200000*exp((40000/1.987)*(1/1250-1/T)); k4=2530000000*exp((9700/1.987)*(1/1250-1/T)); k5=3980000000;
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 xdot(1,:)=-k1*x(1)-k2*x(1)*x(2)-k4*x(1)*x(6); xdot(2,:)=2*k1*x(1)-k2*x(1)*x(2); xdot(3,:)=k2*x(1)*x(2); xdot(4,:)=k2*x(1)*x(2)-k3*x(4)+k4*x(6)*x(1)-k5*x(4)^2; xdot(5,:)=k3*x(4); xdot(6,:)=k3*x(4)-k4*x(6)*x(1); xdot(7,:)=k4*x(1)*x(6); xdot(8,:)=0.5*k5*x(4)^2; xdot(9,:)=k3*(2*k1/k5)^0.5*x(10)^0.5; xdot(10,:)=-k1*x(10)-2*k1*x(10)-(k3*(2*k1/k5)^0.5)*(x(10)^0.5);...
View Full Document

Ask a homework question - tutors are online