Math & Numbers

Численный решатель Рунге—Кутты 4-го порядка (RK4)

Решает dy/dx = f(x,y), y(x₀)=y₀ классическим методом Рунге—Кутты 4-го порядка: на каждом шаге вычисляются четыре наклона k1–k4 и y_{n+1} = y_n + h/6·(k1+2k2+2k3+k4); глобальная погрешность O(h⁴) — при уменьшении шага вдвое ошибка падает примерно в 16 раз. Выводит пошаговую таблицу, значения k1–k4 первых шагов и конечную оценку; необязательное точное решение y(x) добавляет столбец ошибки. Пример: dy/dx = y, y(0)=1, h=0,1 → y(1) ≈ 2,718282, практически e.