ex14_3 - t1=t-80 t if t<80 CTe=CTe1 else CTe=CTe2 end e...

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

View Full Document Right Arrow Icon
% "ex14_3" % m-file to solve example 14.3 % % x(1)=CT1 % x(2)=CT2 % xdot(1)=CT1/dt, xdot(2)=CT2/dt % function xdot=ex14_3(t,x) f beta=0.1; alpha=0.8; tau=40; CTe1=2000-59.6*t+0.64*t^2-0.00146*t^3-1.047*10^(-5)*t^4; CTe2=921-17.3*t+0.129*t^2-0.000438*t^3+5.6*10^(-7)*t^4;
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: t1=t-80; t if t<80 CTe=CTe1; else CTe=CTe2; end e xdot(1,:)=(beta*x(2)-(1+beta)*x(1))/alpha/tau; xdot(2,:)=(beta*x(1)-beta*x(2))/(1-alpha)/tau;...
View Full Document

{[ snackBarMessage ]}

Ask a homework question - tutors are online