euler.m

来自「数值分析各种公式的解法」· M 代码 · 共 23 行

% euler.m 
% MATLAB上提供的“改进的欧拉方法”
function s=euler(fun,x0,xn,y0,n)
% 用euler法计算常微分方程初值问题
% x0和y0分别为初值条件y(x0)=y0中的x0,y0
% xn为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
s=y0;
return

% 首先建立链接向fun的函数
function y=ff2(x,y)
y=y-2*x/y;