我们说的都对+池羽凡12.01+DAY11

[池羽凡]

求常微分方程（ODE）的数值解

[t,y]=ode23(@C7Ex9,[0 10],[2 0]);
plot(t,y);

function xdot=C7Ex9(t,x)
%用于第 7 章练习 9 的函数
%设 x1=y,x2=y’,则 x1’=x2,x2’=y’’=exp(-t)-x1+2x2;
xdot=zeros(2,1);
xdot(1)=x(2);
xdot(2)=exp(-t)-x(1)+2*x(2);
end