6.2.9 Оценка точности квадратурных формул

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

Искомый интеграл вычисляется дважды: при делении отрезка напчастей и на2пчастей (при интегрировании по формуле Симпсонапдолжно быть четным). Вслед за этим полученные значения интегралаI_nиI_2nсравниваются, и совпадающие первые десятичные знаки считаются верными.

Пусть R_nиR_2n– погрешности интегрирования по формуле Симпсона, соответственно припи2п отрезках разбиения. Учитывая оценку (6.43), можно составить равенство:

(6.44)

где h_nиh_2n– длина отрезков разбиения (шаг интегрирования) в первом и втором случае. Т.к.Тогда из (6.44) получаем:

(6.45)

Если I – истинное значение интеграла, тоI=I_n+R_n и I=I_2n+R_{2n ,} откуда

I_n+16R_2n= I_2n+R_2n, т.е.

(6.46)

Формула (6.46) удобна для практической оценки погрешности метода Симпсона, но требует двойного счета.

Из оценочных формул (6.39) и (6.42) следует, что ошибка интегрирования по методу трапеций и методу Симпсона уменьшается с уменьшением шага интегрирования. При последовательном увеличении числа отрезков разбиения будем получать значение интеграла, все более и более близкое к истинному. Этот вывод имеет теоретическое значение. В процессе практических вычислений при последовательном удвоении числа отрезков разбиения начинает сильно прогрессировать удельный вес ошибки округления, значение которой с некоторого момента ставит предел достижимой точности результата интегрирования.