Формула парабол (формула Симпсона или квадратурная формула).

Разделим отрезок интегрирования [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,2482725418…