rk4.m

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

function y = rk4(ef, y0, h, a, b)
y = y0;
n = (b - a) / h;
x = a : h : b;
for i = 1 : n
  k1 = h * feval(ef, x(i), y(i));
  k2 = h * feval(ef, x(i) + h/2, y(i) + k1/2);
  k3 = h * feval(ef, x(i) + h/2, y(i) + k2/2);
  k4 = h * feval(ef, x(i) + h, y(i) + k3);
  y(i+1) = y(i) + (k1 + k2 + k2 + k3 + k3 + k4) / 6;
end