Используется схема
. | (5.9) |
где k1, k2, k3, k4 - поправки, вычисляемые по формулам
. | (5.10) |
При определении yi+1 по заданному yi необходимо четыре раза вычислять правую часть (5.9) в следующей последовательности: k1, k2, k3, k4. Если предположить достаточную гладкость u(t) (непрерывную дифференцируемость вплоть до производных 4-го порядка) и разложить ui+1, k1, k2, k3, k4 в окрестности t=ti, нетрудно показать, что невязка y = 0(t4), т.е. разностная схема (5.9) имеет 4-й порядок аппроксимации.