ex3530.m

来自「数字信号处理的MATLAB实现」· M 代码 · 共 6 行

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