euler.m

来自「数值分析, 同济大学教材<<现代数值数学和计算>>数值分」· M 代码 · 共 8 行

function y = euler(ef, y0, h, a, b)
y(1) = y0;
n = (b - a) / h;
x = a : h : b;
for i = 1 : n
  y(i+1) = y(i) + h * feval(ef, x(i), y(i)); 
end