10.5.3. Формула Эйлера.
Рассмотрим тождество
Проинтегрируем последовательно по частям. Получаем
Применим формулу (64) к последовательным производным
ограничив разложения членами, содержащими производную
Умножив соответственно на
получаем
Умножим равенство (64) на единицу, первое равенство
второе равенство (65) на
третье равенство (65) на нуль, четвертое на
на
и результаты сложим. В правой части полученной таким образом суммы коэффициент при
будет равен нулю, коэффициент при
равный
также нулем. Легко показать, что коэффициент при
равный
есть нуль. Действительно, вынося множитель
получаем
т. е. в силу формулы (56), где взято
Это выражение равно нулю, как показывает формула (59). Поэтому результат сложения равен
По формуле (56) полином в квадратных скобках под знаком интеграла может быть написан в виде
Положим
Тогда равенство (66) можно записать в виде
где
Существенное замечание. Так как числа Бернулли очень быстро возрастают вместе с
то весьма важно заметить, что остаточный член
вообще говоря, бесконечно возрастает вместе с
если
фиксировано. Но если
стремится к нулю, а
фиксировано, то
быстро стремится к нулю. Формула Эйлера (70) представляет собой асимптотический ряд.
Пример. Требуется найти
Результат известен заранее:
Делим промежуток интегрирования на шесть частичных интервалов
и применяем формулу Эйлера при
равном 3.
Вычислим сначала верхнюю границу совершенной ошибки. Так как последовательные производные сохраняют постоянный знак между 1 и 1,6, то
следует взять в виде (72):
Применение формулы Эйлера дает
откуда
Замечание. Из предыдущего примера видно, что применение формулы Эйлера позволяет получить очень большую точность при приближенном вычислении определенного интеграла. Это, однако, верно при условии, что имеется возможность вычислить последовательные значения производных интегрируемой функции на концах промежутка интегрирования и, сверх того, для вычисления совершенной ошибки найти оценку производной высокого порядка внутри этого промежутка. Это, конечно, невозможно, если функция дана в эмпирической форме, и часто бывает затруднительно, если аналитическая форма интегрируемой функции сложна. Отсюда возникает необходимость установить некоторое количество формул достаточной точности, вводящих производные интегрируемой функции только в необязательный первый поправочный член.