Формула парабол (формула Симпсона или квадратурная формула) (Томас Симпсон (1710-1761)- английский математик)

Разделим отрезок интегрирования [a, b] на четное число отрезков (2m). Площадь криволинейной трапеции, ограниченной графиком функции f(x) заменим на площадь криволинейной трапеции, ограниченной параболой второй степени с осью симметрии, параллельной оси Оу и проходящей через точки кривой, со значениями f(x₀), f(x₁), f(x₂).

Для каждой пары отрезков построим такую параболу.

у

0 х₀ х₁ х₂ х₃ х₄ х

Уравнения этих парабол имеют вид Ax² + Bx + C, где коэффициенты А, В, С могут быть легко найдены по трем точкам пересечения параболы с исходной кривой.

(1)

Обозначим .

Если принять х₀ = -h, x₁ = 0, x₂ = h, то (2)

Тогда уравнения значений функции (1) имеют вид:

C учетом этого: .

Отсюда уравнение (2) примет вид:

Тогда

Складывая эти выражения, получаем формулу Симпсона:

Чем больше взять число m, тем более точное значение интеграла будет получено.

Пример. Вычислить приближенное значение определенного интеграла

с помощью формулы Симпсона, разбив отрезок интегрирования на 10 частей.

По формуле Симпсона получим:

Точное значение этого интеграла – 91.173.

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

Для сравнения применим к этой же задаче формулу трапеций.

Формула трапеций дала менее точный результат по сравнению с формулой Симпсона.

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

Принцип этого метода состоит в том, чтобы заменить подинтегральную функцию по формуле Тейлора и почленно проинтегрировать полученную сумму.

Пример. С точностью до 0,001 вычислить интеграл

Т.к. интегрирование производится в окрестности точки х=0, то можно воспользоваться для разложения подинтегральной функции формулой Маклорена.

Разложение функции cosx имеет вид:

Зная разложение функции cosх легко найти функцию 1 – cosx:

В этой формуле суммирование производится по п от 1 до бесконечности, а в предыдущей – от 0 до бесконечности. Это – не ошибка, так получается в результате преобразования.

Теперь представим в виде ряда подинтегральное выражение.

Теперь представим наш интеграл в виде:

В следующем действии будет применена теорема о почленном интегрировании ряда. (Т.е. интеграл от суммы будет представлен в виде суммы интегралов членов ряда).

Вообще говоря, со строго теоретической точки зрения для применения этой теоремы надо доказать, что ряд сходится и, более того, сходится равномерно на отрезке интегрирования [0, 0,5]. Эти вопросы будут подробно рассмотрены позже (См. Действия со степенными рядами.) Отметим лишь, что в нашем случае подобное действие справедливо хотя бы по свойствам определенного интеграла (интеграл от суммы равен сумме интегралов).

Итак:

Итого, получаем:

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

Для справки: Точное (вернее – более точное) значение этого интеграла: 0,2482725418…