gusucode.com > 《精通matlab 7》包括配套光盘 > 配套光盘里的程序附件/Ex-35/Ex3530.m

    %例35-30
t0=0;tf=5;
y0=2;
[t,y]=ode23('func',[t0,tf],y0);   %求数值解
y1=sqrt(t+1)+1;             %求精确解