Gammafunksjon

Gammafunksjonen er en matematisk funksjon av stor betydning i ren matematikk og med mange anvendelser innen teoretisk fysikk. Den er betegnet med en stor, gresk gamma Γ og har den fundamentale egenskapen

hvor argumentet z  også kan være et komplekst tall. Når det er et naturlig tall n, faller den sammen med fakultetsfunksjonen,   = 1⋅2⋅3⋅....⋅n.

Leonhard Euler viste rundt 1730 at den analytiske formen av funksjonen kan defineres ved integralet



når realdelen av z er større enn null. Denne definisjonen gjør det mulig å vise at funksjonen eksisterer i hele det komplekse planet hvor den er regulær bortsett fra  punktene z = 0, -1, - 2, - 3, etc. hvor den har enkle poler og divergerer.

Denne definisjonen av gammafunksjonen kalles ofte for Eulers andre integral.  Hans første integral brukes til definisjon av betafunksjonen som er nært forbundet med gammafunksjonen.

Egenskaper
Fra Eulers integral følger at gammafunksjonen har samme, fundamentale egenskap som fakultetsfunksjonen har. Det kan vises ved en delvis integrasjon som gir

Her er første ledd på venstre side alltid null så lenge som realdelen til z er et positivt tall, uavhengig hva imaginærdelen er. Det andre leddet inneholder det definerende  integralet for gammafunksjone slik at man dermed har bevist det funksjonelle kravet .

Ut fra dette resultatet kan forskjellige verdier av funksjonen bygges opp hvis en verdi er kjent. Den kan for eksempel velges å være



Herav følger nå at for et  positivt heltall n vil funksjonen ta verdien Γ(n) = 1⋅2⋅3⋅...⋅(n - 1) = (n - 1)!. At , tilsvarer at . Formelt har man da  også at  som betyr at gammafunksjonen for argumentet  må ha en uendelig stor verdi, det vil si at den divergerer i dette punktet.

Halvtallige argument

Verdier for gammafunksjonen når argumentet er halvtallig, kan bygges opp fra Γ(1/2). Denne verdien er gitt som



som fremkommer ved å sette t = x 2 slik at dt = 2xdx og så benytte verdien for det gaussiske integralet. Fra dette resultatet har man derfor at  = √π /2, mens Γ(-1/2) = - 2√π  som følger fra .

Man kan på denne måten utlede verdier for alle halvtallige argument. Resultatet kan sammenfattes i de to formlene

hvor n er et positivt heltall. Legg merke til at Γ(1/2 + n)⋅Γ(1/2 - n) =  π⋅(-1)n som er en enkel utgave av Eulers refleksjonsformel.

Det er ikke mulig å finne slike algebraiske resultat for argument med andre brøker. For slike argument må man benytte numeriske metoder, som for eksempel Stirlings formel. Denne er kun gyldig for store verdier av argumentet. Men sammen med egenskapen  kan man fra slike verdier bygge seg nedover til verdier for mindre argument.

Stirlings formel

Omtrent på samme tid som Euler fant et integral for fakultetetsfunksjonen, fant James Stirling et algebraisk uttrykk for funksjonen som kan benyttes for store argument. Dette kan skrives som



Er n > 10, kommer det ledende bidraget fra faktorene n n e -n som kan mer praktisk beregnes fra logaritmen



Det er i denne grensen Stirlings formel vanligvis  blir anvendt. For eksempel, i statistisk mekanikk er n gitt ved antall partikler i systemet og derfor av orden 1023  for et makroskopisk system.

Ved bruk av sadelpunktsmetoden kan korreksjoner til formelen beregnes slik at den kan brukes med forbausende god nøyaktighet. helt ned til n = 2. Den tar da formen



Høyere ordens korreksjoner vil ikke nødvendigvis fortsatt gi mer nøyaktige resultat da disse ikke utgjør en konvergent rekke, men er kun asymptotisk konvergerende.

En viktig egenskap ved Stirlings formel i sammenligning med Eulers integral, er at den består av algebraiske ledd som kan benyttes også for argument som ikke er heltallige. Den vil derfor kunne brukes for argument som er relle eller til og med komplekse tall så lenge som disse har en positiv realdel. Men igjen kan man benytte  når dette ikke er tilfelle.

Historie

På begynnelsen av 1700-tallet var det av en viss interesse blant flere matematikere å finne ut hvordan et matematisk uttrykk som i utgangspunktet kun er gyldig for  heltallige argument, kan generaliseres til kontinuerlige funksjoner som kan interpolere mellom  diskrete verdier. Et enkelt eksempel er en potens an som er veldefinert for heltallige n, men som det er ønskelig å kunne beregne for brøker og reelle tall.

Den tyske matematiker Christian Goldbach ville gjerne finne ut om dette var mulig for fakultetsfunksjonen. Han henvendte seg til Daniel Bernoulli som arbeidet ved Det russiske vitenskapsakademiet i St. Petersburg og mottok et svar den 6. oktober 1729. Her foreslo Bernoulli i et PS å bruke produktet



hvor A er et vilkårlig stort tall. Uttrykket inneholder derfor også et like stort antall faktorer. Han illustrerte nøyaktigheten av denne formelen ved å regne ut 3! som han med fant lik med 6 og 1/204 med valget  A = 16 . Det må betraktes som et overmåte godt resultat. Derimot fant han  for (3/2)! verdien 1.3005 med A = 8,  mens den riktige er 1.3293. Gronau har påpekt at dette sannsynligvis skyldes en regnefeil da formelen gir 1.3291.D. Gronau, Why is the gamma function so as it is?, Teaching Mathematics and Computer Science, 1 (1), 43-53 (2003).  I brevet ga han ingen forklaring på hvordan han var kommet frem til dette uttrykket. Det var ikke uvanlig på denne tiden å holde sine metoder mer eller mindre hemmelige.

Eulers formel

Euler arbeidet på denne tiden også ved akademiet i St. Petersburg og må ha diskutert disse problemstillingene med Bernoulli. Den 13. oktober skriver han til Goldbach om samme sak og foreslår å bruke formelen



for n-fakultet hvor N er et stort tall som skal tas mot uendelig. På denne formen er det lettere å ha kontroll med konvergensen til produktet. Men ser man bort fra det og multipliserer faktorene sammen, tar den en litt enklere form



Dette uttrykket er likt med Bernoullis resultat i denne grensen, men konvergerer langsommere. Likevel er det herfra at gammafunksjonen tok sin endelige form.

Formelen til Euler oppfyller formelt det viktigste kravet (n + 1)! = (n + 1)⋅n !. Det ser man ved å betrakte



Men her er siste faktor lik med 1 i grensen hvor N blir veldig stor slik at det formelle kravet er oppfylt. Likedan ser man at  og . Dermed vil  og så videre når N går mot uendelig. Likeså viktig var det å ha en formel som også kan benyttes når argumentet n ikke er heltallig.

Euler presenterte disse resultatene den 28. november 1729 i et større arbeid De progressionibus transcendentibus seu quarum termini generales algebraice dari nequeunt (Om transcendentale progresjoner, det vil si slike hvor den generelle termen ikke kan uttrykkes algebraisk) for akademiet i St. Petersburg.L. Euler, On transcendental progressions, that is, those whose general terms cannot be given algebraically, Commentarii academiae scientiarum Petropolitanae 5, 36-57 (1738), engelsk oversettelse av S. Langton, Euler Archive. Der viste han også at formelen ville gi riktig verdi for (1/2)!. Det følger direkte fra produktformlen hvor faktorene ordnes slik at deres produkt gjenkjennes,R. Tambs-Lyche, Lærebok i Matematisk Analyse, Volum III, Gyldendal Norsk Forlag, Oslo (1959).



Den siste kvadratroten her gir 1/√2 i grensen hvor N blir veldig stor, mens det første produktet er ikke noe annet enn kvadratrotroten  av π /2 som følger fra Wallis' formel. Derfor gir formelen det ønskede resultatet f(1/2) = (1/2)⋅√π.

Eulers integral
I dette store arbeidet fra november, 1729 hadde Euler vist at hans produktformel for fakultetsfunksjonen i grensen N → ∞ kunne skrives om til et kompakt integral,



Ved å innføre den nye variable t = - lns i integralet slik at dt = - ds/s, tar det den mer velkjente formen



Dette elegante resultatet beskrev han i et nytt brev til Goldbach i den 8. januar 1730. Den interpolernede formel var ikke lenger algebraisk, men istedet transcendent. I det samme arbeidet hvor han kom frem til denne integralformen, hadde han også undersøkt et annet integral som i ettertid kalles hans betafunksjon.

Navnet gammafunksjon og notasjonen Γ(n + 1) = n ! ble innført rundt 1810 av Legendre. Alternativt foreslo Gauss å betegne funksjonen med symbolet . Denne notasjonen slo ikke an, men kan likevel av og til finnes. En grunn for at det ble slik, kan være at symbolet Π er standard for et matematisk produkt. Men i tillegg kan det også skyldes at forholdet mellom disse to store matematikere ikke alltid var like godt.M. Kline, Mathematical Thought from Ancient to Modern Times, Volume 2, Oxford University Press, Oxford (1972). ISBN 978-0-19-506136-9.

Produktformelen
For et generelt argument z med en positiv realdel, er gammafunksjonen gitt ved integralet



Dette uttrykket er i overensstemmelse med produktformelen. For å se det, kan man bruke uttrykket



for eksponentialfunksjonen i grensen der tallet N blir uendelig stort. Dermed kan man skrive



Ved å utføre N - 1 delvise integrasjoner etter hverandre, gir dette



Da man har at z Γ(z) = Γ(z + 1) = z !, er dette den opprinnelige produktformelen til Euler i grensen hvor N går mot uendelig.

Alternativt kan man se formelen skrevet på formen



etter å ha tatt N mot uendelig og benyttet de gjensidige kansellasjonene mellom faktorene i de to brøkene som inngår i produktet. Denne versjonen tilsvarer Eulers aller første definisjon av gammafunksjonen med raskere konvergens.

Eulers refleksjonsformel

Senere i sitt vitenskapelige arbeid oppdaget Euler en egenskap ved gammafunksjonen som viste seg å være viktig. Den følger fra hans produktformel og kan anvendes på den spesielle kombinasjonen



I det siste leddet i denne utledningen opptrer det uendelige produktet som Euler benyttet i 1735 ved løsningen av Basel-problemet.

Dette elegante resultatet kalles Eulers refleksjonsformel. Navnet er naturlig da den forbinder verdien av gammafunksjonen i et reelt punkt z med verdien i punktet - z som man kan tenke seg fremkommer ved «refleksjon» om den imaginære aksen i det komplekse planet.

En enkel anvendelse av formelen er å ta z = 1/2. Da er 1 - z = 1/2 som sammen med sinπ /2 = 1 gir  i overensstemmelse med andre beregningsmetoder.

Mer generelt kan man sette z = 1/2 + x slik at 1 - z = 1/2 - x. Da tar refleksjonsformelen formen



fordi



Komplekse argument
Refleksjonsformelen har mange anvendelser som at den kan benyttes til å beregne størrelsen av gammafunksjonen for komplekse argument z = x + iy. Dette er basert på egenskapen til funksjonen at dens komplekskonjugerte



hvor det komplekskonjugerte argument z* = x - iy. Dette følger fra dens definisjon.

For et rent, imaginært argument z = iy gir formelen



Men nå kan den vanlige sinusfunksjonen for imaginært argument uttrykkes ved en hyperbolsk sinusfunksjon ved at  Dermed har man resultatet



Et annet, enkelt eksempel gjelder de spesielle, komplekse argumentene langs linjen z = 1/2 + iy. Da blir 1 - z = 1/2 - iy =  z*. Dermed er



da cos iπ y = coshπ y.

For et generelt, komplekst argument z = x + iy kan ikke refleksjonsformelen benyttes på noen enkel måte. Derimot kan man gå tilbake til produktformelen som gir



hvor det uendelige produktet konvergerer raskt.

Digammafunksjonen

Mange av egenskapene til gammafunksjonen ligger skjult i den deriverte Γ'(x) = d Γ/dx. Fra det definerende integralet finner man



Spesielt viktig er verdien



hvor γ = 0.577215... er Euler-Mascheronis konstant.

I stedet for å arbeide med den deriverte av gammafunksjonen, er det mer vanlig å betrakte den relaterte funksjonen



som kalles digammafunksjonen. Da Γ(1) = 1, har man den spesielle verdien ψ(1) = - γ. Generelt kan den beregnes fra



Når argumentet x er et positivt heltall n,  tar den verdier som er harmoniske tall,



Da  ψ(1) = - γ, betyr det at det harmoniske tallet H 0 = 0 som kommer i tillegg til den opprinnelige definisjonen av disse. Med   betyr det at  . Digammafunklsjonen skifter derfor fortegn mellom x = 1 og x = 2 og har et nullpunkt i dette intervallet.

Se også

 Fakultetsfunksjon
 Betafunksjon
 Euler-Mascheronis konstant

 Referanser 

Litteratur

 M.L. Boas, Mathematical Methods in the Physical Sciences, John Wiley & Sons, New York (1983). ISBN 0-471-04409-1.
 J. Mathews and R.L. Walker, Mathematical Methods of Physics,  W.A. Benjamin, New York (1970). ISBN 0-8053-7002-1.
 M. Abramowitz and I.A. Stegun,  Handbook of Mathematical Functions, Dover, New York (1972).
 P.J. Davis, Leonhard Euler's Integral: A Historical Profile of the Gamma Function: In Memoriam: Milton Abramowitz , The American Mathematical Monthly, 66 (10), 849-869 (1959).

Eksterne lenker

 E.W. Weisstein, MathWorld, Gamma Function

Kategori:Funksjoner