Вычисление определенных интегралов методом Симпсона
Возьмем для аппроксимации функции f (x) многочлен второй степени . Чтобы определить коэффициенты a,b,c, необходимо на интервале аппроксимации задать три точки. Выберем шаг интегрирования D x и вычислим значения функции на концах интервала и в его середине. Обозначим, как это показано на рис.23.5, значения аргумента через x0, x1, x2 и, соответственно, значения функции через f0, f1, f2. Половину шага интегрирования обозначим .
Уравнение параболы, проходящей через три точки, можно записать в виде
.
Учитывая, что , при интегрировании этой функции в пределах от x0 до x0+2h после преобразований получим:
.
Для всего интервала интегрирования будем иметь:
, (23.5)
где .
Это формула Симпсона.
Как и для рассмотренных ранее методов, вычисления интегралов по формуле Симпсона могут быть по блок-схемам рис.23.2 или рис.23.3.
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | Поиск по сайту:
|