Stirlings formel

Stirlings formel er en approksimasjon til fakultetsfunksjonen for store verdier av argumentet  oppkalt etter den skotske matematiker James Stirling. Mest vanlig er det å skrive formelen uttrykt ved logaritmefunksjonen på den korte formen

 

hvor man utelater termer av størrelsesorden ln n. Ekvivalent kan man derfor skrive at n ! = nn e-n.

Navnet på denne formelen er til en viss grad misvisende da den først ble funnet av den franske matematiker Abraham de Moivre som viste at

  

hvor C er en numerisk konstant med en verdi av størrelsesorden en. Det var Stirling som klarte å beregne denne konstanten og fant C = √(2π ) = 2,5066... I de aller fleste tilfeller av praktisk betydning er resultatet til de Moivre mer enn godt nok.

Enkel utledning

Fra definisjonen for  n fakultet har man at

 

Når argumentet n er stort, kan man nå erstatte summasjonen med en integrasjon,

 

Det siste ett-tallet her kan neglisjeres i denne approksimasjonen hvor n er mye større enn 1.

Systematisk utledning
En mer nøyaktig utledning av formelen kan finnes ved bruk av gammafunksjonen til Euler. Den tilsier at

 

Nå kan man skrive at  xn = en lnx slik at integralet er på formen

 

med f(x) = x - n lnx. Denne funksjonen har et minimum for x0 = n og antar positive verdier for større verdier av argumentet og når dette nærmer seg null. Verdien av integralet er derfor med stor nøyaktighet bestemt av hvordan funksjonen  f(x) varierer rundt x0. Dette kommer frem fra en Taylor-utvikling av funksjonen rundt dette punktet,

 

Da  f(x) har et minimum for x =  x0 , er den førstederiverte f' (x0) = 0. Første, ikke-trivielle ledd i Taylor-utviklingen er da gitt verd den andrederiverte f" (x0). Ser vi i denne omgang bort fra høyere ordens ledd, er derved integralet gitt ved

 

I denne tilnærmelsen kan vi nå med tilsvarende nøyaktighet ta den nedre integrasjonsgrensen til - ∞  slik at integralet antar formen

 

som bærer navnet Gauss-integralet. 

For funksjonen f(x) som opptrer her, er f(x0) = n - n lnn og f" (x0) =  1/n. Dermed oppnås resultatet

 

og kan med rette kalles Stirlings formel. Konstanten C som han regnet ut, er ikke noe annet enn verdien av integralet til Gauss.

Denne metoden tillater også beregning av høyere ordens korreksjoner til formelen. Tar man med neste ledd i Taylor-utviklingen som ble benyttet, finner man

 

Her kan man nå benytte at 1 + 1/12n ≈ √(1 + 1/6n) som tillater å skrive formelen på formen

 

Den er nå nøyaktig selv for små verdier av argumentet n. Ved å sette inn numeriske verdier, finner man at for n > 2 gir den et resultat som har en nøyaktighet bedre enn 0,1 %. For større verdier av n, blir nøyaktigheten tilsvarende enda bedre.

Eksterne lenker
 Stirlings Approximation in MathWorld

Funksjoner