4.3. Формула трапеций
Если аппроксимировать подынтегральную функцию в (1) многочленом Лагранжа:
(10)
то, с учетом того, что, согласно (3.1),
(11)
получим квадратурную формулу трапеций
(12)
Основная ее часть имеет смысл площади трапеции. По главному члену погрешности виден второй порядок точности формулы.
Обобщенная формула трапеций имеет вид:
(13)
На равномерной сетке
(14)
По последнему слагаемому можно априорно оценить погрешность основной части этой формулы, имеющей второй порядок точности. Соответствующие формулы средних (6),(7) имеют погрешность вдвое меньшую, чем формулы трапеций и поэтому более предпочтительны.
Если в формуле (14) последнее слагаемое включить в основной результат, то главным членом погрешности становится слагаемое четвертого порядка малости:
(15)
В такой интерпретации формула (14) называется квадратурной формулой Эйлера-Маклорена.
4.4. Формула Симпсона
Еще одна возможность получения квадратурных формул - использование метода Рунге-Ромберга (подраздел 3.3). Для вычисления интеграла
(16)
применим обобщенную формулу трапеций (14) с n=2
(17)
на сетке с шагом h=(b-a)/2. Здесь, в соответствии с обозначениями формулы (3.13),
(18)
На сетке с вдвое более крупным шагом 2h=b-a главная часть формулы трапеций имеет вид
(19)
Учитывая, что r=2, после применения формулы Рунге-Ромберга (3.15) получим уточненную квадратурную формулу
(20)
называемую формулой Симпсона.
Более детальный анализ позволяет определить главный член погрешности
(21)
Таким образом, формула Симпсона (20) имеет четвертый порядок точности. Она точна для многочлена третьей степени.
Обобщенная формула Симпсона для равномерной сетки и четного числа шагов n имеет вид
(22)
4.5. Формулы Гаусса и Маркова
Квадратурная формула Ньютона-Котеса (3) имеет n+1 параметр сi. Поэтому она точна для полинома n-й степени. Гаусс рассматривал узловые точки xi также как варьируемые параметры, общее число которых при этом становится равным 2(n+1). Таким образом, рациональным выбором узловых точек можно получить квадратурные формулы Гаусса, точные для полиномов (2n+1)-й степени. Эти формулы записываются для a=-1, b=1 в виде
(23)
Ниже приведены узлы и веса для формул с различным числом узлов m:
(24)
Если
интеграл задан в произвольном виде,
то к стандартному виду (23) он приводится
заменой переменных =kx+d,
где k=(b-a)/2;
d=(a+b)/2.
Границы области интегрирования часто вызывают особый интерес и поэтому бывает необходимо расположить в них узлы. Если остальные узлы выбирать из соображений наивысшей точности квадратурной формулы, то последняя, называемая формулой Маркова, будет иметь точность на два порядка меньшую, чем соответствующая формула Гаусса. Узлы и веса для квадратурных формул вида (23) имеют следующие значения:
(25)