ff6.m

来自「研究生期间的数值计算的书配套程序源代码 奉献了」· M 代码 · 共 9 行

x0=0;xf=0.99999;
[x,y]=ode23('eq1',[x0 xf],[0 0]);
 Y=0:0.01:2;
   plot(1,Y,'g.')
   hold on
plot(x,y(:,1),'b*')