Запишемо рівняння першого порядку
З початковими умовами (1,2)
Нехай xi(i=0,1,2…)-система рівнозначних значень з кроком h i y(xi). Очевидно маємо
В силу другої інтерполяційної формули Ньютона з точністб до різниць четвертого порядку отримуємо:
де
Підставляю вираз (4а) в формулу (3) і враховуючи те, що
З відси отримуємо формулу експоляриціональну Адамса
Для початкового процессу потрібно чотири початкових значення y0, y1, y2, y3, - початковий відрізок, який приділяє, виходячи із початкових умов (2), яким-небуть чисельним методом. Мажна наприклад використати метод Рунге-Кутта або розкласти в ряд Тейлора
Де i=1,2,3 (або i=-1,1,2) із відповідною зміною нумерування. Знаючи ці значення, із рівнянь (1) можна знайти значення похідних
Подальше значення yi (i=4,5…) шуканого розвязку можна крок за кроком обчислювати за формулою Адамса, поповнюючи по мірі можливості таблицю різниць (6)
Вирахувавши перше наближення для
Визначити
а потім знайти друге наближення для більш точній формулі
На практиці крок h вибирають малим, щоб можна було знехтувати членом
Якщо за розбіжність величин
Звичайно крок h зменшують рівно в 2 рази. Можна показати, як в цьому випадку, маючи до деякого значення і таблицю величин хj, yj, Yj=hy’j (j<=i) з кроком
На основі формули (4) будемо мати
Де
Аналогічно при
Що стосується значень Yi-1 i Yi, то вони знаходяться в старій таблиці. Після цього складаємо початковий відрізок для нової таблиці:
і знаходимо кінцеві різниці:
Далі таблиця будується простим способом, подальшою модифікацією формули (5):
Для роботи на компютерах формулу Адамса (5) вигідно використовувати в розкритому виді. Враховуючи, що
Після цього маємо:
Метод Крилова
Для спрощеня запису обмежимось розглядом диференціальних рівнянь першого порядка
З початковими умовами
Введемо спочатку ряд допоміжних формул
В силу формули Адамса отримаємо
Введемо позначення
Формула (2) називається формулою похилого рядка, так як в ній використовуються різниці, які знаходяться на діагоналі таблиці різниць. Враховуючи, що
Із формули (2) будемо мати
Звідси отримуємо першу допоміжну формулу – яку ще можна назвати перша формула ламаного рядка
Далі враховуючи, що
Якщо
Підмітимо, що формулу (5) можна отримати безпосередньо за допомогою інтегрування, в межах від xi до xi+1 розкладанням
Перейдемо до опису метода Крилова послідовних наближень. Перше наближення полягає у тому, щоб знайти наближене значення
Значення які знайшли заносимо в розділ (І) основного бланку (таблиця 1)