WikiPrépaLivrets

Formule de Stirling

La formule de Stirling approche n!n! par 2πn (n/e)n\sqrt{2\pi n}\,(n/e)^n avec une erreur relative qui tend vers 00. Sous forme logarithmique, elle donne ln⁡(n!)\ln(n!) à o(1)o(1) près. Elle sert à estimer des coefficients binomiaux, des probabilités et des termes de séries.

Au programmeDeuxième année (MP, MPI, PSI) · Séries numériques

Énoncé

Théorème (formule de Stirling)

Lorsque n→+∞n\to+\infty,

n!=2πn (ne)n(1+o(1)).n!=\sqrt{2\pi n}\,\left(\frac ne\right)^n\big(1+o(1)\big).

Forme logarithmique.

ln⁡(n!)=nln⁡n−n+12ln⁡n+ln⁡2π+o(1).\ln(n!)=n\ln n-n+\frac12\ln n+\ln\sqrt{2\pi}+o(1).

En particulier ln⁡(n!)=nln⁡n−n+O(ln⁡n)\ln(n!)=n\ln n-n+O(\ln n) et ln⁡(n!)∼nln⁡n\ln(n!)\sim n\ln n.

Terme correctif (hors programme).

n!=2πn (ne)n(1+112n+O(1n2)).n!=\sqrt{2\pi n}\,\left(\frac ne\right)^n\left(1+\frac1{12n}+O\Big(\frac1{n^2}\Big)\right).

Esquisse de démonstration

Équivalent et forme logarithmique. On pose un=ln⁡(n!)−(n+12)ln⁡n+nu_n=\ln(n!)-\left(n+\frac12\right)\ln n+n. Un développement limité donne

un−un+1=(n+12)ln⁡(1+1n)−1=112n2+O(1n3).u_n-u_{n+1}=\left(n+\tfrac12\right)\ln\left(1+\tfrac1n\right)-1=\frac1{12n^2}+O\Big(\frac1{n^3}\Big).

La série télescopique ∑(un−un+1)\sum(u_n-u_{n+1}) converge, donc un→ℓu_n\to\ell, et les intégrales de Wallis donnent ℓ=ln⁡2π\ell=\ln\sqrt{2\pi} (démonstration détaillée dans la fiche Équivalent de Stirling). La forme logarithmique n'est autre que un=ln⁡2π+o(1)u_n=\ln\sqrt{2\pi}+o(1).

Terme correctif. On somme les restes :

un−ℓ=∑k=n+∞(uk−uk+1)=112∑k=n+∞1k2+O(∑k=n+∞1k3).u_n-\ell=\sum_{k=n}^{+\infty}(u_k-u_{k+1})=\frac1{12}\sum_{k=n}^{+\infty}\frac1{k^2}+O\Big(\sum_{k=n}^{+\infty}\frac1{k^3}\Big).

Pour k≥2k\ge2, l'encadrement 1k−1k+1≤1k2≤1k−1−1k\frac1k-\frac1{k+1}\le\frac1{k^2}\le\frac1{k-1}-\frac1k donne par télescopage 1n≤∑k≥n1k2≤1n−1\frac1n\le\sum_{k\ge n}\frac1{k^2}\le\frac1{n-1}, donc ∑k≥n1k2=1n+O(1n2)\sum_{k\ge n}\frac1{k^2}=\frac1n+O(\frac1{n^2}). De même ∑k≥n1k3≤∫n−1+∞dtt3=O(1n2)\sum_{k\ge n}\frac1{k^3}\le\int_{n-1}^{+\infty}\frac{\dd t}{t^3}=O(\frac1{n^2}) par comparaison avec une intégrale. Ainsi un=ln⁡2π+112n+O(1n2)u_n=\ln\sqrt{2\pi}+\frac1{12n}+O(\frac1{n^2}) et, en passant à l'exponentielle, e112n+O(1/n2)=1+112n+O(1n2)e^{\frac1{12n}+O(1/n^2)}=1+\frac1{12n}+O(\frac1{n^2}).

Méthode : quand et comment l'utiliser

Pour un quotient de factorielles, on remplace chaque factorielle par la formule, on simplifie les puissances nnn^n et e−ne^{-n}, puis on regroupe les facteurs 2πn\sqrt{2\pi n}. Pour une limite du type (⋯ )1/n(\cdots)^{1/n} ou pour comparer des ordres de grandeur, la forme logarithmique est plus commode.

Exemples d'application

Coefficient binomial central

Avec (2n)!∼4πn (2ne)2n=2πn 4nn2ne−2n(2n)!\sim\sqrt{4\pi n}\,\left(\frac{2n}{e}\right)^{2n}=2\sqrt{\pi n}\,4^nn^{2n}e^{-2n} et (n!)2∼2πn n2ne−2n(n!)^2\sim2\pi n\,n^{2n}e^{-2n} :

(2nn)=(2n)!(n!)2∼2πn 4n2πn=4nπn.\binom{2n}{n}=\frac{(2n)!}{(n!)^2}\sim\frac{2\sqrt{\pi n}\,4^n}{2\pi n}=\frac{4^n}{\sqrt{\pi n}}.

La probabilité d'obtenir exactement nn piles en 2n2n lancers d'une pièce équilibrée, (2nn)4−n\binom{2n}{n}4^{-n}, est donc équivalente à 1πn\frac1{\sqrt{\pi n}} et tend vers 00. Pour n=10n=10 : (2010)=184756\binom{20}{10}=184756 et 41010π≈187079\frac{4^{10}}{\sqrt{10\pi}}\approx187079.

Précision numérique

Pour n=10n=10 : 10!=362880010!=3628800 et 20π (10/e)10≈3598696\sqrt{20\pi}\,(10/e)^{10}\approx3598696. Le quotient vaut environ 1,008371{,}00837, très proche de 1+1120≈1,008331+\frac1{120}\approx1{,}00833 : l'erreur relative est de l'ordre de 112n\frac1{12n}, comme l'annonce le terme correctif.

Erreurs classiques

  • Déduire n!∼nnn!\sim n^n ou n!∼(n/e)nn!\sim(n/e)^n de ln⁡(n!)∼nln⁡n\ln(n!)\sim n\ln n : un équivalent entre logarithmes ne donne pas d'équivalent entre les quantités. Pour exponentier, il faut la forme logarithmique à o(1)o(1) près.

  • Utiliser le terme correctif 112n\frac1{12n} sans le justifier : il n'est pas au programme.

  • Remplacer une factorielle par la formule dans une somme ou une différence de termes du même ordre : les équivalents ne s'additionnent pas.

  • Se tromper sur (2n)!(2n)! : il fait apparaître 4πn\sqrt{4\pi n} et (2n)2n=4nn2n(2n)^{2n}=4^nn^{2n} ; une erreur sur ces facteurs change la puissance de nn ou de 44 dans le résultat.

Place dans le programme

Deuxième année (MP, MPI, PSI), chapitre « Séries numériques ».

L'équivalent n!∼2πn (n/e)nn!\sim\sqrt{2\pi n}\,(n/e)^n est au programme de deuxième année et peut être cité directement. Le terme correctif 112n\frac{1}{12n} est hors programme : il faut le redémontrer (voir l'esquisse ci-dessous) avant de l'utiliser dans une copie ou à l'oral.

Dans les sujets de concours

Pour vous entraîner sur des sujets complets : Séries numériques et familles sommables (224 sujets), Probabilités finies, discrètes et dénombrement (438 sujets).

Sujets où cette notion intervient explicitement :

Questions fréquentes

Quelle est la formule de Stirling ?

n!∼2πn (n/e)nn!\sim\sqrt{2\pi n}\,(n/e)^n quand n→+∞n\to+\infty, ou sous forme logarithmique ln⁡(n!)=nln⁡n−n+12ln⁡(2πn)+o(1)\ln(n!)=n\ln n-n+\frac12\ln(2\pi n)+o(1).

Comment démontrer la formule de Stirling ?

On montre que ln⁡(n!)−(n+12)ln⁡n+n\ln(n!)-(n+\frac12)\ln n+n converge (série télescopique de terme O(1/n2)O(1/n^2)), puis on identifie la limite ln⁡2π\ln\sqrt{2\pi} grâce aux intégrales de Wallis.

Faut-il écrire formule de Stirling ou formule de Sterling ?

On écrit formule de Stirling, du nom du mathématicien écossais James Stirling (XVIIIe siècle). « Sterling » est une confusion avec la livre sterling.

Quelle est la précision de la formule de Stirling ?

L'erreur relative est équivalente à 112n\frac1{12n} : elle vaut environ 0,840{,}84 % pour n=10n=10 et passe sous 0,10{,}1 % à partir de n=84n=84 environ.

Exercices corrigés utilisant : Formule de Stirling

12 exercices