eult.m

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

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