文档介绍：微分方程数值解dy/dx=y/2x+x2/2yy(1)=1a)欧拉格式代码:functionf=ww(x,y)f=y/(2*x)+x^2/(2*y);h=;y(1)=1;x(1)=1;n=1;whilex(n)<=2n=n+1;x(n)=x(n-1)+h;y(n)=y(n-1)+h*ww(x(n-1),y(n-1));x(n)y(n)end显示结果:==;h==)中点格式代码:h=;y(1)=1;x(1)=1;n=1;whilex(n)<=2n=n+1;x(n)=x(n-1)+h;y(n)=y(n-1)+h*ww(x(n-1)+h/2,y(n-1)+h*ww(x(n-1),y(n-1))/2);x(n)y(n))预报校正格式代码:h=;y(1)=1;x(1)=1;n=1;whilex(n)<=2n=n+1;x(n)=x(n-1)+h;yy(n)=y(n-1)+h*ww(x(n-1),y(n-1));y(n)=y(n-1)+h*[ww(x(n-1),y(n-1))+ww(x(n),yy(n))]/2;x(n)y(n)end显示结果:)R