euler.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 17 行

function E=euler(f1521,a,b,ya,N)
%  f1521为该问题中的函数f(x,y)。
%  a,b分别为取值范围的左右端点。
%  ya为给定初值y(a)。
%  N为迭代步数。
%  h为步长。
%  输出值为对应每个节点的近似值。
h=(b-a)/N;
T=zeros(1,N+1);
Y=zeros(1,N+1);
T=a:h:b;
Y(1)=ya;
for j=1:N
    Y(j+1)=Y(j)+h*feval('f1521',T(j),Y(j));
end
T=[T' Y']