eulg.m

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

function eulg(x0,y0,b,h)
n=fix((b-x0)/h);
for i=1:n
    x(i)=x0+h;
    yi=y0+h*(fa(x0,y0)+fa(x0+h,y0+h*fa(x0,y0)))/2;
    y(i)=y0;
    x0=x0+h;
    y0=yi;
    hold on
end
[x',y']
function fa=fa(x,y)
fa=y-2*x/y;