스털링 공식

수학노트
http://bomber0.myid.net/ (토론)님의 2009년 8월 17일 (월) 15:48 판
둘러보기로 가기 검색하러 가기
간단한 소개
  • 팩토리얼의 근사식

\( n! \approx \sqrt{2\pi n}\, \left(\frac{n}{e}\right)^{n}\)

  • 좀더 정확히는 다음과 같이 주어짐.
    \( n!=\sqrt{2\pi n}\left({n\over e}\right)^n \left( 1 +{1\over12n} +{1\over288n^2} -{139\over51840n^3} -{571\over2488320n^4} + \cdots \right)\)

 

증명

 

\(\sum_{i=0}^{n-1} f(i) = \sum_{k=0}^p\frac{B_k}{k!}\left(f^{(k-1)}(n)-f^{(k-1)}(0)\right)+R\)

\(\sum_{i=0}^{n-1} f(i) = \int^n_0f(x)\,dx-\frac{1}{2}(f(n)-f(0))+\frac{1}{12}(f'(n)-f'(0))-\frac{1}{720}(f^{(3)}(n)-f^{(3)}(0))+\frac{1}{30240}(f^{(5)}(n)-f^{(5)}(0))-\frac{1}{1209600}(f^{(7)}(n)-f^{(7)}(0))+\cdots\)

오차항은  \(\left|R\right|\leq\frac{2}{(2\pi)^{2(p+1)}}\int_0^n\left|f^{(p)}(x)\right|\,dx\) 로 주어짐.

  • \(f(x) = \ln(1+x)\) 에서 유도. \(R\) 은 0 으로 수렴한다.
    \(f^{(k-1)}(x)=(-1)^{k-1}\frac{(k-1)!}{x^{k}}\)

\(\frac{B_k}{k!}\left(f^{(k-1)}(n)-f^{(k-1)}(1)\right) =(-1)^{k-1}\frac{B_k}{k}(\frac{1}{n^{k}}-1)\)

\(\int f(x)\,dx=(1+x)\ln (1+x)-x\), \(f(x)=\ln (1+x)\), \(f'(x)=\frac{1}{1+x}\),  \(f^{(2)}(x)=\frac{-1}{(1+x)^2}\), \(f^{(3)}(x)=\frac{2}{(1+x)^3}\),  \(f^{(4)}(x)=\frac{-6}{(1+x)^4}\), \(f^{(k-1)}(x)=(-1)^{k}\frac{(k-1)!}{x^{k-1}}\)

  • 위의 결과를 쓰면

\(\ln(n!)= \sum_{k = 0}^{n-1}\ln(1+k)=(n+1)\ln (n+1) - n-\frac{1}{2}\ln(n+1)+\frac{1}{12}(\frac{1}{1+n}-1)+\cdots\)

그러므로,

\( n! \approx B\sqrt{n}\, \left(\frac{n}{e}\right)^{n}\)

  • 계수 \(\sqrt{2 \pi}\) 는 어떻게 얻을까?

\(\frac{\pi}{2}=\lim_{n\to\infty}{1\over{2n}}\cdot{{2^{4n}\,(n!)^4}\over{((2n)!)^2}}\)

\(\sqrt{\frac{\pi}{2}}=\lim_{n\to\infty}\sqrt{{1\over{2n}}}\cdot{{2^{2n}\,(n!)^2}\over{(2n)!}}=\lim_{n\to\infty}\sqrt{{1\over{2n}}}\cdot{{2^{2n}\,(B\sqrt{n}\, \left(\frac{n}{e}\right)^{n})^2}\over{B\sqrt{2n}\, \left(\frac{2n}{e}\right)^{2n}}}=\frac{B}{2}\)

\(B=\sqrt{2\pi}\)

 

재미있는 사실
  • 드무아브르가 스털링에 앞서서 팩토리얼의 근사식을 유도함.
  • 다만 \(\sqrt{2\pi}\)이라는 상수를 구하지 않고 적당한 상수 \(B\)에 대하여 다음과 수준의 표현을 남김.
    \( n! \approx B\sqrt{n}\, \left(\frac{n}{e}\right)^{n}\)
  • 나중에 \(B=\sqrt{2\pi}\) 는 스털링이 해결함.

In Miscellanea Analytica (1730) appears Stirling’s formula (wrongly attributed to Stirling) which de Moivre used in 1733 to derive the normal curve as an approximation to the binomial. In the second edition of the book in 1738 de Moivre gives credit to Stirling for an improvement to the formula. De Moivre wrote:-

I desisted in proceeding farther till my worthy and learned friend Mr James Stirling, who had applied after me to that inquiry, [discovered that c = √(2 π)].

  • 오늘날 팩토리얼의 근사식은 선구적인 업적을 남긴 드무아브르의 이름은 온데간데 없이 스털링의 공식이라고 불리게 됨.
  • 드무아브르는 조금 섭섭하지 않을런지?
  • 팩토리얼은 정의는 간단할지라도 n이 조금만 커지기 시작하면 계산하기가 그리 만만치 않음. 실용적인 측면에서도 유용한 근사식.

 

관련된 단원

 

 

많이 나오는 질문

 

관련된 고교수학 또는 대학수학

 

관련된 다른 주제들

 

관련도서 및 추천도서

 

참고할만한 자료

 

관련기사

 

 

블로그

 

이미지 검색

 

동영상
원본 주소 "https://wiki.mathnt.net/index.php?title=스털링_공식&oldid=14213"