eulg.m

数值分析各种公式的解法

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;