eulh.m

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

function eulh(x0,y0,b,h,eps)
n=fix(b-x0)/h;
for i=1:n
    x(i)=x0+h;
    yi=y0+h*fa(x0,y0);
    yj=y0+h*fa(x0+h,yi);
    while abs(yj-yi)>eps
        p=yj;
        yj=y0+h*fa(x0+h,yi);
        yi=p;
end
x0=x0+h;
y0=yj;
y(i)=y0;
hold on
end
[x',y']
function fa=fa(x,y)
fa=y-2*x/y;