Оценка погрешности численного интегрирования

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

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

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

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

Например, выражение для остаточного члена формулы трапеций имеет вид:

Для остаточного члена обобщенной формулы трапеций получаем:

, где

- это интегральная сумма для , т.е.

, где - это величина, стремящаяся к нулю при (). Таким образом, для погрешности обобщенной формулы трапеций имеем формулу вида:

, где (75)

Аналогично, для погрешности обобщенной формулы Симпсона выводится формула , где выражается через интеграл от производной четвертого порядка подинтегральной функции. Тогда первая формула Рунге (31) дает оценку погрешности значения интеграла, вычисленного о обобщенной квадратурной формуле с шагом () через это приближение и приближение, вычисленное с шагом ():