chase2.m

来自「Matlab教程（含M文件）」· M 代码 · 共 9 行

t0=0,tf=2;
[t,y]=ode45('eq2',[t0 tf],[0 0]);
X=1;
Y=0:0.01:2;
plot(X,Y,'-')
hold on
plot(y(:,1),y(:,2),'*')