eul.m

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

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