Семейство методов Адамса

Известны методы Адамса k-го порядка. Простейший из них при k = 1 повторяет метод Эйлера первого порядка в точности. Метод четвертого порядка на практике принято называть методом Адамса. Рабочую формулу для него получают следующим образом.

Пусть известные в четырех последовательных узлах (k = 4) значение сеточной функции yi3, yi2, yi1, yi и вычисленные первоначально значения правой части (4) fi3, fi2, fi1, fi. В качестве интерполяционного многочлена P3(x) возьмем многочлен Ньютона. В случае h = const конечные разности для правой части в узле xi будут иметь вид

Dfi = fifi–1;

D2fi = fi – 2 fi–1 + fi–2;

D3fi = fi – 3 fi–1 + 3 fi–2fi–3.

Тогда разностная схема метода Адамса запишется в виде

. (26)

По сравнению с методом Рунге-Кутта той же точности можно отметить его экономичность, так как (26) предусматривает на каждом шаге только один раз вычисление правой части в соотношении (4). Однако расчет здесь можно начать только с узла x4. Значения y1, y2, y3 необходимые для вычисления y4 нужно определять одношаговым методом, что несколько усложняет алгоритм вычисления. Кроме того, метод Адамса не позволяет изменять шаг h в процессе счета, что доступно для одношаговых методов.