Метод Симпсона

Значительное повышение точности приближенных формул может быть достигнуто за счет повышения порядка интерполяции. Одним из таких методов приближенного интегрирования является метод парабол. Идея метода исходит из того, что на частичном промежутке дуга некоторой параболы в общем случае теснее прилегает к кривой y=f(x), чем хорда, соединяющая концы дуги этой кривой, и поэтому значения площадей соответствующих элементарных трапеций, ограниченных “сверху” дугами парабол, являются более близкими к значениям площадей соответствующих частичных криволинейных трапеций, ограниченных сверху дугой кривой y=f(x), чем значения площадей соответствующих прямолинейных трапеций. Сущность метода заключается в следующем. Отрезок [a,b] делится на 2n равных частей. Пусть точки деления будут

х0=а, x1, x2, …x2n-2, x2n-1, x2n=b,

а y0, y1, …y2n – соответствующие значения подинтегральной функции на отрезке [a,b]. Произведем квадратичную интерполяцию данной подинтегральной функции на каждом из отрезков разбиения (заменим дугу графика подинтегральной функции дугой параболы с вертикальной осью) (рис.11).

Приведем без вывода формулу парабол в окончательном виде:

(7)

Если подинтегральная функция f(x) имеет на отрезке [a,b] непрерывную четвертую производную, то для поправочного члена формулы (7) имеет место оценка (8)

где М4- максимум модуля четвертой производной подинтегральной функции на отрезке [a,b].

Cравнивая между собой оценки (6) и (8), замечаем, что с увеличением n поправочный член формулы трапеций уменьшается пропорционально

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