仿真大实验作业.docVIP

大实验 实验一、 （1）运用欧拉法： 主程序段（M函数）： function[x,y]=naler(dyfun,xspan,y0,h) x=xspan(1):h:xspan(2); y(1)=y0; for n=1:length(x)-1; y(n+1)=y(n)+h*feval(dyfun,x(n),y(n)); end x=x;y=y; 运行过程及中间结果： clear; dyfun=inline(y-2*x/y); [x,y]=naler(dyfun,[0,1],1,0.2); p=[x,y] p = 0 1.0000 0.2000 1.2000 0.4000 1.3733 0.6000 1.5315 0.8000 1.6811 1.0000 1.8269 x x = 0 0.2000 0.4000 0.6000 0.8000 1.0000 y y = 1.0000 1.2000 1.3733