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

Таблицы узлов и весов, а также оценки погрешности квадратурных формул приводятся в специальных справочниках.

Квадратурные формулы вычисления кратных интегралов иногда называются кубатурными формулами. Кратные интегралы можно вычислять как повторные интегралы, применяя описанные квадратурные формулы. Т. к. при увеличении кратности существенно возрастает количество узлов, то для вычисления кратных интегралов разработан ряд специальных формул.

Вычисление интегралов на ЭВМ обычно осуществляется с помощью стандартных программ. В случае однократных интегралов наиболее употребительны стандартные программы с автоматическим выбором шага.

Лит.: Крылов В. И., Приближенное вычисление интегралов, 2 изд., М., 1967; Бахвалов Н. С., Численные методы, М., 1973; Никольский С. М., Квадратурные формулы, М., 1958; Березин И. С., Жидков Н. П., Методы вычисления, 3 изд., ч. 1, М., 1966; Соболев С. Л., Введение в теорию кубатурных формул, М., 1974; Коробов Н. М., Теоретикочисловые методы в приближенном анализе, М., 1963.

В. И. Лебедев.

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

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

Если решение ищется в виде бесконечного ряда, то за П. р. принимают конечный отрезок ряда. Например, пусть требуется найти решение дифференциального уравнения y' = f ( x, у ) , удовлетворяющее начальным условиям у ( х 0 ) = y 0, причём известно, что f ( x, у ) — аналитическая функция х, у в некоторой окрестности точки ( х 0, y 0 ) . Тогда решение можно искать в виде степенного ряда:

y ( x ) - y ( x 0 ) = .

Коэффициенты A k ряда могут быть найдены либо по формулам:

A 1 = y’ 0 = f ( x 0, y 0 ) ;

либо с помощью неопределенных коэффициентов метода . Метод рядов позволяет находить решение лишь при малых значениях величины х — х 0.

Часто (например, при изучении периодических движений в небесной механике и теории колебаний) встречается случай, когда уравнение состоит из членов двоякого вида: главных и второстепенных, причём второстепенные члены характеризуются наличием в них малых постоянных множителей. Обычно после отбрасывания второстепенных членов получается уравнение, допускающее точное решение. Тогда решение основного уравнения можно искать в виде ряда, первым членом которого является решение уравнения без второстепенных членов, а остальные члены ряда расположены по степеням малых постоянных величин, входящих во второстепенные члены (малых параметров). При этом уравнения для коэффициентов при степенях малых параметров линейны, что облегчает их решение. В роли малого параметра иногда выступают начальные значения (например, при изучении колебаний около положения равновесия). Метод малого параметра был использован при решении задачи о возмущённом движении в небесной механике Л. Эйлером и П. Лапласом . Теоретическое обоснование этого метода дали А. М. Ляпунов и А. Пуанкаре .

К численным методам относятся методы, позволяющие находить П. р. при некоторых значениях аргумента (т. е. получать таблицу приближённых значений искомого решения), пользуясь известными значениями решения в одной или нескольких точках. Такими методами являются, например, метод Эйлера, метод Рунге и целый ряд разностных методов.

Поясним эти методы на примере уравнения

y’ ’ = f ( x, у )

с начальным условием у ( х 0 ) = y 0 . Пусть точное решение этого уравнения представлено в некоторой окрестности точки х 0 в виде ряда по степеням h = х — х 0 Основной характеристикой точности формул П. р. дифференциальных уравнений является требование, чтобы первые k членов разложения в ряд по степеням h П. р. совпадали с первыми k членами разложения в ряд по степеням h точного решения.

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

,

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

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

,

где

;

;

;

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