Большая Советская Энциклопедия (ПР)

  Основная идея метода Эйлера заключается в применении метода рядов для вычисления приближённых значений решения у (х ) в точках x₁ , x₂ ,..., x_n некоторого фиксированного отрезка [х₀ , b ] Так, для того чтобы вычислить у (х₁ ), где х₁ = х₀ + h, h = (b — x₀ )/n, представляют у (х₁ ) в виде конечного числа членов ряда по степеням h = х₁ — х₀ . Например, ограничиваясь первыми двумя членами ряда, получают для вычисления у (xk ) формулы:

,

Это т. н. метод ломаных Эйлера (на каждом отрезке [xk_, x_k+1 ] интегральная кривая заменяется прямолинейным отрезком — звеном ломаной Эйлера). Погрешность метода пропорциональна h² .

  В методе Рунге вместо того, чтобы отыскивать производные, находят такую комбинацию значений f (x, у ) в некоторых точках, которая даёт с определённой точностью несколько первых членов степенного ряда для точного решения уравнения. Например, правая часть формулы Рунге:

,

где

;

;

;

дает первые пять членов степенного ряда с точностью до величин порядка h⁵ .

  В разностных формулах П. р. удаётся несколько раз использовать уже вычисленные значения правой части. Решение ищется в виде линейной комбинации у (xi ), h_i и разностей Dⁱ h_j , где

h_j = hf (x_j , y_j ); Dh_j = h_j+1 - h_j ;

Dⁱ h_j = D^i-1 h_j+1 - D^i-1 h_j .

  Примером разностной формулы П. р. является экстраполяционная формула Адамса. Так, формула Адамса, учитывающая «разности» 3-го порядка:

даёт решение у (х ) в точке x_k с точностью до величин порядка h⁴ .

  Для уравнений 2-го порядка можно получить формулы численного интегрирования путём двукратного применения

Формула	k = 2	k = 3	k = 4
(1 + x )3 » 1 + 3x	0,04	0,012	0,004
	0,06	0,022	0,007
	0,19	0,062	0,020
	0,20	0,065	0,021
	0,31 (17°48')	0,144 (8°15')	0,067 (3°50')
	0,10 (5°43')	0,031 (l'48')	0,010 (0°34')
	0,25 (14°8')	0,112 (6°25')	0,053 (3°2')
	0,14	0,47	0,015
	0,04	0,014	0,004
	0,25	0,119	0,055

формулы Адамса. Норвежский математик К. Стёрмер получил формулу:

особенно удобную для решения уравнений вида у'' = f (x, у ). По этой формуле находят D² y_n-1 , а затем y_n+1 = y_n +Dy_{n+ 1} + D² y_n-1 . Найдя y_n+1 , вычисляют y’’_n+1 = f (x_n+1 , y_n+1 ), находят разности и повторяют процесс далее.

  Указанные выше численные методы распространяются и на системы дифференциальных уравнений.

  Значение численных методов решения дифференциальных уравнений особенно возросло с распространением ЭВМ.

  Кроме аналитических и численных методов, для П. р. дифференциальных уравнений применяются графические методы. В простейшем из них строят поле направлений, определяемое дифференциальным уравнением, т. е. в некоторых точках рисуют направления касательной к интегральной кривой, проходящей через эту точку. Затем проводят кривую так, чтобы касательные к ней имели направления поля (см. Графические вычисления ).

  Лит.: Березин И. С., Жидков Н. П., Методы вычислений, 2 изд., т. 2, М.. 1962; Бахвалов Н. С., Численные методы, М., 1973: Коллатц Л., Численные методы решения дифференциальных уравнений, пер. с нем., М., 1953; Милн В. Э., Численное решение дифференциальных уравнений, пер, с англ., М., 1955.

Приближённые вычисления

Приближённые вычисле'ния, вычисления, в которых данные и результат (или по крайней мере только результат) являются числами, лишь приближённо представляющими истинные значения соответствующих величин. П. в. возникают в связи с численным решением задач и обусловлены неточностями, которые присущи формулировке задачи и способам её решения. Общие правила и теорию методов П. в. принято называть численными методами .

Приближённые формулы

Приближённые фо'рмулы, математические формулы, получаемые из формул вида f (x ) = f* (x ) + e(х ), где e(х ) рассматривается как погрешность и после оценки отбрасывается. Таким образом, П. ф. имеет вид f (x ) » f* (x ).

  Например, П. ф. (1 + х )² » 1 + 2x получается из точной формулы для (1 + х )² при малых |x |; этой формулой можно пользоваться при вычислении с точностью до сотых, тысячных, десятитысячных, если |x | соответственно не больше 0,0707..., 0,0223..., 0,00707... Эта П. ф. даёт результат тем более точный, чем х ближе к 0. Но так бывает не всегда. Например, точность П. ф. tg тем больше, чем х ближе к p/2.

  Выше (стр. 555) приведено несколько наиболее употребительных П. ф., причём показано, какого числа не должно превосходить |x |, чтобы формула давала k точных десятичных знаков.