M < 0.1. В этом случае каждое приближение содержит по крайней мере на один верный знак больше, чем предыдущее.

9.1.2. Контрольные вопросы

Что понимается под приближенным решением уравнения? Что такое погрешность?

Вчем состоит процедура отделения корня?

Вчем состоит суть метода деления отрезка пополам?

Какова процедура отыскания корня уравнения методом хорд, каков критерий выбора начального приближения?

В чем состоит суть метода итераций, каково достаточное условие применимости данного метода?

9.2. Примеры численного интегрирования

9.2.1. Вопросы для самостоятельного изучения

Рассмотрим непрерывную на [a,b] функцию f (x) . Пусть требуется найти

b

∫ f (x)dx . Формула Ньютона-Лейбница позволяет свести эту задачу к задаче на-

a

хождения первообразной для функции f (x) . С помощью этой формулы можно

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

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

Поэтому важно знать общие универсальные методы вычисления определенных интегралов. Это приближенные методы, которые позволяют вычислить интеграл по значению подынтегральной функции во множестве точек отрезка [a,b] . Соответствующие формулы называются формулами численного интег-

рирования. Они основываются на определении определенного интеграла, как предела интегральных сумм

126

9.2.1.1. Формулы прямоугольников

b

Пусть требуется найти ∫ f (x)dx . Построим интегральную сумму следую-

a

При достаточно мелком разбиении отрезка [a,b] (т.е. достаточно малом

шаге h) интегральную сумму можно считать приближенным значением интеграла. Таким образом, получаем формулу

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

Название этой формулы связано с ее геометрической интерпретацией

b

(рис. 9.2.1). Интеграл ∫ f (x)dx – это площадь криволинейной трапеции, а при-

a

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

127

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

В данном случае ступенчатая фигура состоит из прямоугольников с основание h и высотой, равной значению подынтегральной функции в левом конце отрезка [xk−1, xk ].

Если при построении интегральной суммы точку ξk [xk−1, xk ], k =1,n ,

выбрать так, чтобы ξk = xk , то аналогично получим еще один вид формулы

прямоугольников (формула правых прямоугольников)

9.2.1.2. Формула трапеций

128

Более точную приближенную формулу для вычисления интеграла можно получить, если заменить график функции f (x) ломаной линией, состоящей из отрезков прямых, проходящих через точки (xk , f (xk )).

В результате получим, что площадь криволинейной трапеции приближенно заменяется суммой площадей трапеций высоты h и с основаниями, равными yk−1 и yk (рис. 9.2.2):

Эта формула называется формулой трапеций. Погрешность формулы (39) можно оценить по формуле

9.2.1.3. Формула Симпсона

Еще более точное значение определенного интеграла можно получить с помощью формулы Симпсона.

При ее построении отрезок интегрирования [a,b] разбивается на четное

число 2n равных частей длины h = b2−na . На каждом k -ом сдвоенном отрезке

[x2(k−1) , x2k ] график функции приближенно заменяется отрезком параболы, про-

ходящей через точки (x2(k−1) , f (x2(k−1) )), (x2k−1, f (x2k−1)), (x2k , f (x2k )).

Можно показать, что уравнение этой параболы будет иметь вид: 129