常微分方程作业欧拉法与改进欧拉法.pdfVIP

P77 31.利用改进欧拉方法计算下列初值问题，并画出近似解的草图： dy (1) y 1, y(0) 3,0 t 2, t 0.5; dt 代码： %改进欧拉法 function Euler(t0,y0,inv,h) n=round(inv(2)-inv(1))/h; t(1)=t0; y(1)=y0; for i=1:n y1(i+1)=y(i)+h*fun(t(i),y(i)); t(i+1)=t(i)+h; y(i+1)=y(i)+1/2*h*(fun(t(i),y(i))+ fun(t(i+1),y1(i+1))) end plot(t,y,*r) function y=fun(t,y); y=y+1; 调用： Euler(0,3,[0,2], 得到解析解： hold on; y=dsolve(Dy=y+1,(y(0)=3),t); ezplot(y,[0,2]) 图像： dy 2 (2) dt y 4t , y (0) 0.5,0 t 2, t