R_n

= I

 - S_n

.

  Квадратурная формула содержит 2n + 1 не зависящих от функции f (x ) параметров: n, x_k , A_k (k = 1, 2,..., n ), которые выбирают так, чтобы при f Î W погрешность её была допустимо малой. Точность квадратурной формулы для f Î W характеризует величина rn (W) — точная верхняя грань ½Rn

½ на множестве W:

.

  Пусть

  Квадратурная формула, для которой W_n (W) = r_n (W), называется оптимальной на классе П. Веса и узлы в оптимальной квадратурной формуле могут быть произвольными или подчинёнными определённым связям.

  Различают два класса квадратурных формул: элементарные и составные. Разработано несколько методов построения элементарных квадратурных формул. Пусть w_q (x ), q = 0, 1,..., — полная система функций в классе W, и любая f (x ) Î Q достаточно хорошо приближается линейными комбинациями первых функций w_q (x ). Пусть l (w_q ), q = 0, 1, 2,..., можно вычислить точно. Для каждого n параметры квадратурной формулы можно определить из требования, чтобы

I (w_q ) = S_n (w_q ), q = 0, 1,..., m,

для возможно большего значения m. В методе Ньютона — Котеса в квадратурной формуле выбираются узлы x_k , а определению подлежат веса A_k . В методе Чебышева на веса A_k заранее накладываются некоторые связи [например, A_k = (b -а )/n ], а определению подлежат узлы x_k . В методе Гаусса определяются и веса A_k и узлы x_k . В методе Маркова j узлов (j < n ) считают заранее известными, а определяют веса и оставшиеся узлы. Точность полученных такими методами квадратурных формул существенно повышается при удачном выборе функций w_q (x ).

  Формулы Ньютона — Котеса строятся на основе системы функций w_q = x^q , q = , 1,...; узлы x_k разбивают отрезок интегрирования на равные части. Примерами таких формул являются прямоугольников формула , трапеций формула и Симпсона формула .

  Поскольку заменой переменной интегрирование по [а, b ] сводится к интегрированию по отрезку [-1, 1], то для определения весов и узлов элементарных формул на [а, b ] достаточно знать их для отрезка [-1, 1]. В случае составных формул исходный интеграл представляется в виде:

и для вычисления интегралов по отрезкам [a_i , a_i+1 ] применяются элементарные квадратурные формулы.

  В формулах Гаусса m = 2n — 1, а при а = — 1, b = 1 узлы x_k являются корнями Лежандра многочленаP_n (x ) степени n, а

A_k = 2(1 - x²_k )^-1 (P’_n (x_k ))^-2

  Квадратурная формула Чебышева существует при A_k = l/n, l = b - а и x_k Î [a, b ] лишь для n = 1,..., 7, 9; в ней m = n - 1. Применение равных весов минимизирует вероятностную ошибку, если значения f (x ) содержат независимые случайные ошибки с одинаковой дисперсией.

  При вычислении интегралов от функций с периодом l наиболее употребительны квадратурные формулы типа Гаусса:

.

  Существуют квадратурные формулы для вычисления интегралов вида

где р (х ) — фиксированная, т. н. весовая функция. Её подбирают так, чтобы для всех f Î W функции f (x ) хорошо приближалась линейными комбинациями функций w_q (x ).

  Для приближённого вычисления неопределённых интегралов их представляют как определённые интегралы с переменным верхним пределом и далее применяют указанные выше формулы.

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

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

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

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

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

Приближённое решение

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

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

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

y (x ) - y (x ) =

.

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

A₁ = y’ = f (x , y );