Практическая оценка погрешности при применении квадратурных формул

В полученные нами формулы оценки ошибки квадратурной формулы входят величины Mk, ограничивающие абсолютную величину производной порядка от подынтегральной функции (для формул центральных прямоугольников и трапеций,для формулы Симпсона,идля формулы Уэддля). Если величинаM_k неизвестна (а как правило, в достаточно сложных задачах не вычисление интегралов она неизвестна или получить её весьма нелегко), то пользоваться этими оценками для определения величины ошибки конкретного вычисления невозможно. Так что всё, что дают нам формулы оценки ошибки — это порядок квадратурных формул. Однако на этом основании можно получить следующее практическое правило, которое позволяет получить оценку ошибки конкретного вычисления, если квадратурную формулу применить два раза с разными шагами .

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

Следовательно, если — приближённое значение интеграла, точное значение которого равно, то

Отсюда получаем, что

(5)

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

На такой оценке текущей погрешности, как правило, основаны компьютерные программы, вычисляющие значение определённого интеграла с заданной точностью.

Пример. Вычислить интеграл по формуле Ньютона — Лейбница и по приближенным формулам прямоугольников, трапеций и Симпсона, разбивая интервал на 8 равных частей. Оценить в процентах погрешность результатов, полученных по приближенным формулам.