euler.m

来自「MATLAB科学计算与工程分析源代码源程序3」· M 代码 · 共 9 行

function [x,y]=euler(f,x0,y0,xf,h)
n=fix((xf-x0)/h);
y(1)=y0;
x(1)=x0;
for i=1:n
    x(i+1)=x0+i*h;
    y(i+1)=y(i)+h*feval(f,x(i),y(i));
end