Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 2)

Несобственный интеграл равномерно сходится. Его можно рассматривать, например, как интеграл по треугольнику, ограниченному осями координат и прямой x+y = R при R . В двойном интеграле сделаем замену переменных:

Якобиан этой замены

Пределы интегрирования: u меняется от 0 до ∞, v при этом меняется от 0 до 1. В результате получим:

Перепишем опять этот интеграл как повторный, в результате получим:

где Rp > 0, Rv > 0.

1. Производная гамма функции

Интеграл

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

В области , где - произвольное положительное число, этот интеграл сходится равномерно, так как и можно применить признак Вейрштраса. Сходящимся при всех значениях является и весь интеграл так как и второе слагаемое правой части является интегралом, заведомо сходящимся при любом .Легко видеть что интеграл сходится по в любой области где произвольно. Действительно для всех указанных значений и для всех ,и так как сходится, то выполнены условия признака Вейерштрасса. Таким образом , в области интеграл сходится равномерно.

Отсюда вытекает непрерывность гамма функции при .Докажем дифференцируемость этой функции при .Заметим что функция непрерывна при и , и покажем ,что интеграл :

сходится равномерно на каждом сегменте , . Выберем число так , чтобы ; тогда при .Поэтому существует число такое , что и на .Но тогда на справедливо неравенство

и так как интеграл сходится, то интеграл сходится равномерно относительно на . Аналогично для существует такое число , что для всех выполняется неравенство . При таких и всех получим , откуда в силу признака сравнения следует , что интеграл сходится равномерно относительно на . Наконец , интеграл

в котором подынтегральная функция непрерывна в области

, очевидно, сходится равномерно относительно на . Таким образом , на интеграл

сходится равномерно , а, следовательно , гамма-функция бесконечно дифференцируема при любом и справедливо равенство

.

Относительно интеграла можно повторить те же рассуждения и заключить, что

По индукции доказывается , что Г-функция бесконечно дифференцируема при и для ее я -ой производной справедливо равенство

Изучим теперь поведение - функции и построим эскиз ее графика. (см. Приложение 1)

Из выражения для второй производной -функции видно, что для всех . Следовательно, возрастает. Поскольку , то по теореме Роля на сегменте [1,2]производная при и при , т. е. Монотонно убывает на и монотонно возрастает на . Далее , поскольку , то при . При из формулы следует , что при .

Равенство , справедливое при , можно использовать при распространении - функции на отрицательное значение .

Положим для , что . Правая часть этого равенства определена для из (-1,0). Получаем, что так продолженная функция принимает на (-1,0) отрицательные значения и при , а также при функция .

Определив таким образом на , мы можем по той же формуле продолжить ее на интервал (-2,-1). На этом интервале продолжением окажется функция, принимающая положительные значения и такая, что при и . Продолжая этот процесс, определим функцию , имеющею разрывы в целочисленных точках (см. Приложение 1.)

Отметим еще раз, что интеграл

определяет Г-функцию только при положительных значениях , продолжение на отрицательные значения осуществлено нами формально с помощью формулы приведения .

4. Вычисление некоторых интегралов.

Формула Стирлинга

Применим гамма функцию к вычислению интеграла:

где m > -1,n > -1.Полагая , что ,имеем

и на основании (2.8) имеем

(4.1)

В интеграле

Где k > -1,n > 0,достаточно положить

Интеграл

Где s > 0,разложить в ряд

=

где дзетта функция Римана

Рассмотрим неполные гамма функции (функции Прима)

связанные неравенством

Разлагая, в ряд имеем

Переходя к выводу формулы Стирлинга , дающей в частности приближенное значение n! при больших значениях n ,рассмотрим предварительно вспомогательную функцию

(4.2)

Непрерывна на интервале (-1, ) монотонно возрастает от до при изменении от до и обращаются в 0 при u = 0.Так как

то при u > 0 и при u < 0 , далее имеем

И так производная непрерывна и положительна во всем интервале ,удовлетворяет условию

Из предыдущего следует, что существует обратная функция, определенная на интервале непрерывная и монотонно возрастающая в этом интервале,

Обращающаяся в 0 при v=0 и удовлетворяющая условие

(4.3)

Формулу Стирлинга выведем из равенства

полагая ,имеем

Положим далее введенная выше обратная функция, удовлетворяющая условиям u = -1при ,и при .Замечая что(см.4.2)

имеем

,

полагая на конец , ,получим

или

в пределе при т.е. при (см 4.3)

откуда вытекает формула Стирлинга

которую можно взять в виде

(4.4)

где ,при

для достаточно больших полагают

(4.5)

вычисление же производится при помощи логарифмов

если целое положительное число, то и (4.5) превращается в приближенную формулу вычисления факториалов при больших значениях n

приведем без вывода более точную формулу

где в скобках стоит не сходящийся ряд.

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

Для вычисления необходимы формулы:

Г( )

Вычислить интегралы

ПРАКТИЧЕСКАЯ ЧАСТЬ

Для вычисления гамма-функции используется аппроксимация её логарифма. Для аппроксимации гамма-функции на интервале x>0 используется следующая формула (для комплексных z):

Г(z+1)=(z+g+0.5)^z+0.5exp(-(z+g+0.5)) [a₀+a₁/(z+1)+a₂/(z+2)+...+a_n/(z+n)+eps]

Эта формула похожа на аппроксимацию Стирлинга, но в ней имеется корректирующая серия. Для значений g=5 и n=6, проверено, что величина погрешности ε не превышает 2*10^-10. Более того, погрешность не превышает этой величины на всей правой половине комплексной плоскости: z > 0.

Для получения (действительной) гамма-функции на интервале x>0 используется рекуррентная формула Г(z+1)=zГ(z) и вышеприведенная аппроксимация Г(z+1). Кроме того, можно заметить, что удобнее аппроксимировать логарифм гамма-функции, чем ее саму. Во-первых, при этом потребуется вызов только одной математической функции - логарифма, а не двух - экспоненты и степени (последняя все равно использует вызов логарифма), во-вторых, гамма-функция - быстро растущая для больших x, и аппроксимация ее логарифмом снимает вопросы переполнения.

Для аппроксимации Ln(Г(х) - логарифма гамма-функции - получается формула:

log(Г(x))=(x+0.5)log(x+5.5)-(x+5.5)+

log(C₀(C₁+C₂/(x+1)+C₃/(x+2)+...+C₇/(x+8))/x)

Значения коэффициентов C_k - табличные данные (см. в программе).

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

Заключение

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

Благодаря этому они широко применяются в математике и ее приложениях, в механике, термодинамике и в других отраслях современной науки.

Список литературы

1. Специальные функции и их приложения:

Лебедев И.И.,М.,Гостехтериоиздат,1953

2. Математический анализ часть 2:

Ильин О.А., Садовничий В.А., Сендов Бл.Х.,М.,”Московский университет”,1987

3. Сборник задач по математическому анализу:

Демидович Б.П.,М.,Наука,1966

4. Интегралы и ряды специальные функции:

Прудников А.П., Брычков Ю.А.,М.,Наука,1983

5. Специальные функции:

Кузнецов , М.,”Высшая школа”,1965

6.Асимптотика и специальные функции

Ф.Олвер, М.,Наука,1990.

7.Зоопарк чудовищ или знакомство со спецмальными функциями

О.М.Киселёв,

ПРИЛОЖЕНИЯ

Приложение 1 - График гамма-функции действительного переменного

Приложение 2 – График Гамма-функции

Таблица – таблица значений гамма-функции при некоторых значениях аргумента.

Приложение 3 – листинг программы, рисующий таблицу значений гамма-функции при некоторых значениях аргумента.

Приложение 4 – листинг программы, рисующей график гамма-функции

СОДЕРЖАНИЕ

Реферат ...................................3

Введение ...................................4

Теоретическая часть…………………………………………………….5

Бета функция Эйлера…………………………………………….5

Гамма функция. ...................................8

2.1. Определение………………………………………………...8

2.2. Интегральное представление………………………………8

2.3. Область определения и полюсы…………………………..10

2.4. Представление Ганкеля через интеграл по петле………..10

2.5. Предельная форма Эйлера………………………………...12

Характеристики

Список файлов курсовой работы