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

    Разложение функции ln (1 + x) в степенной ряд дает возможность легко получить асимптотическую формулу для факториала n! при nбесконечность. Эта формула называется  формулой Стирлинга; она имеет вид

n! ~ ,     nбесконечность,

(32.76)

т. е. отношение n! к выражению, стоящему в правой части этой формулы, стремится к 1 при nбесконечность.
начало    Действительно, если |x| < 1, то

ln = ln (1 + x) - ln (1 - x) =  (-1)n+1xn/n - (xn/n) = 2xsumk0inf.gif (946 bytes)x2n/(2k + 1).

    Положив x = 1/(2n + 1), n = 1, 2, ... получим

ln (1 + 1/n) = =
= 1 + + + ... > = ,

откуда

(n + 1/2)ln (n + 1/2) > 1.

Пропотенцировав и заметив, что функция ln x возрастает, а 1 = ln e, получим

(n + 1/2)(n+1/2) > e.

(32.77)

Положим

xn = .

(32.78)

Поскольку, согласно неравенству (32.77),

>1,

то последовательность {xn} убывает. Кроме того, она ограничена снизу: xn > 0. Следовательно, она имеет конечный предел. Обозначим его a:

xn = a.

(32.79)

Покажем, что aне равно0. Так как

 + + ... <  + + ... = = ,

то

(n + 1/2)ln (n + 1/2) < 1 +

и, следовательно,

(1 + 1/n)n+1/2 < e1+1/(12n(n+1)).

Поэтому

< e1+1/(12n(n+1)) = .

т. е.

xne-1/(12n) < xn+1e-1/(12n(n+1)).

(32.80)

Последовательность yn = xne-1/(12n), n = 1, 2, ..., будучи произведением двух последовательностей {xn} и {e-1/(12n)}, имеющих предел, также имеет предел, причем

yn = xne-1/(12n)  a.

    Неравенство (32.80) означает, что последовательность {yn} возрастает, поэтому yn, а так как yn > 0, то доказано, что a > 0.
    Из равенства (32.79) следует, что

xn = a(1 + эпсилонn),

(32.81)

где эпсилонn = 0. Подставив (32.81) в (32.78), получим

(32.82)

    Для того чтобы найти значение числа a, вспомним, что по формуле Валлиса (см. (26.9) в п. 26.2)

.

(32.83)

    Из формулы (32.82) следует, что

Подставив это выражение в формулу Валлиса (32.83), получим

(32.84)

откуда a = . Следовательно,

n! = (1 + эпсилонn),     

т. е. формула Стирлинга (32.76) доказана. конец


Разложение элементарных функций в ряд Тейлора  Оглавление