Hi, I want to plot on the interval the solution of the following system of ODEs
with the initial conditions
,
where
.
I wrote the function, by defining the coefficient f, using "piecewise" sintax:
function dzdt=odefunw1(t,z) syms t f=piecewise(0<=t<1, -2*t^2+3*t, t>=1, 1/t); dzdt=zeros(2,1); dzdt(1)=z(2); dzdt(2)=z(1)-f*z(2); end
In the Command Window, I wrote
tspan = [0 1000]; z0 = [0.1 0.1]; [t,z] = ode45(@(t,z) odefunw1(t,z), tspan, z0); plot(t,z(:,1),'g',t,z(:,2),'b')
But several errors occured. Where do I wrong and how could I fix it? Thanks in advance.
Best Answer