euler.m

MATLAB科学计算与工程分析源代码源程序3

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