euler2.m

来自「有关matlab r2007的实用教程」· M 代码 · 共 19 行

function [xout,yout]=Euler2(fun,x0,xn,y0,n)
%用改进欧拉方法计算常微分方程的初值问题
%x0,y0为初值条件
%xn,yn为x取值区间的最后一个节点的横坐标和纵坐标
%n为区间的等分数目
if nargin<5
    error;
    return
end
h=(xn-x0)/n;
for i=1:n
    yp=y0+h*feval(fun,x0,y0);
    x0=x0+h;
    yc=y0+h*feval(fun,x0,yp);
    y0=(yp+yc)/2;
end
xout=x0;
yout=y0;