Gammafunktion

Die Eulersche Gammafunktion, auch kurz Gammafunktion oder Eulersches Integral zweiter Gattung, ist eine der wichtigsten speziellen Funktionen und wird in den mathematischen Teilgebieten der Analysis und der Funktionentheorie untersucht. Sie wird heute durch ein $\Gamma$ , den griechischen Großbuchstaben Gamma, bezeichnet und ist eine transzendente meromorphe Funktion mit der Eigenschaft

\Gamma (n)=(n-1)!

für jede natürliche Zahl $n$ , wobei mit ${\displaystyle$ die Fakultät bezeichnet wird. Die Motivation zur Definition der Gammafunktion war, die Fakultätsfunktion auf reelle und komplexe Argumente erweitern zu können. Der Schweizer Mathematiker Leonhard Euler löste im Jahr 1729 diese Fragestellung und definierte die Gammafunktion mittels eines unendlichen Produktes. Heute wird die Gammafunktion oft mittels einer Integraldarstellung definiert, die ebenfalls auf Euler zurückgeht.

Die Gammafunktion liegt der Gamma-Wahrscheinlichkeitsverteilung zugrunde.

Einordnung ohne mathematisches Vorwissen

Eine mathematische Funktion funktioniert im Grunde wie eine Rechenmaschine. Man gibt einen Wert in die Funktion ein, und diese liefert dann ein Ergebnis in Abhängigkeit vom Eingabewert, zumindest theoretisch. Damit ist gemeint, dass die Funktion an sich nicht rechnet, sondern meist nur eine Rechenvorschrift formelhaft festhält. Einfaches Beispiel für eine Funktion ist die quadratische Funktion, welche die Eingabe mit sich selbst multipliziert. Formelhaft schreibt man dies als $f(x)=x^{2}$ . Somit ordnet die quadratische Funktion beispielsweise der Zahl $-2$ den Wert $(-2)^{2}$ zu. Rechnet man dies aus, ergibt sich $4$ , also $f(-2)=4$ .

Die Gammafunktion fußt auf einer Vorschrift, die auch als Fakultät bekannt ist. Diese ordnet einer natürlichen Zahl das Produkt aller natürlichen Zahlen bis zu dieser Zahl zu. Bezeichnet wird die Fakultät mit dem Symbol des Ausrufezeichens. Also gilt zum Beispiel

4!=1\cdot 2\cdot 3\cdot 4=24.

Es galt innerhalb der Mathematik als Problem, ob sich diese Vorschrift auch auf Zahlen anderer Art erweitern ließe. Konkret bedeutet das:

Lassen sich Fakultäten auch für beliebige rationale, reelle, komplexe Zahlen berechnen? Wie in etwa könnte man sich $({\tfrac {1}{2}})!$ vorstellen?
Falls solche „universelle“ Vorschriften gefunden werden, welche mathematischen Eigenschaften können ihnen gegeben werden? Zeichnet sich eine dieser Vorschriften als ganz besonders natürlich und strukturell aus? Ist diese besondere Vorschrift eindeutig bestimmt, also „die eine“ verallgemeinerte Fakultät?

Die Antwort auf diese Fragen liefert die Gammafunktion. Für beliebige Werte $z$ liefert $\Gamma (z+1)=z!$ , also gilt zum Beispiel $\Gamma (5)=24.$ Die Verschiebung um 1 von der oben erwähnten Fakultät ist auf eine Konvention aus dem 19. Jahrhundert zurückzuführen. Die Strategie der Verallgemeinerung basiert auf der Beobachtung, dass aus einer vorherigen Fakultät durch Hinzunahme eines weiteren Faktors eine weitere Fakultät gewonnen wird. So gilt etwa $4!\cdot 5=5!$ und ganz allgemein $n!\cdot (n+1)=(n+1)!$ . Demnach sollte die Gammafunktion die Funktionalgleichung $\Gamma (z)\cdot z=\Gamma (z+1)$ erfüllen. Stellt man weitere wichtige Bedingungen, wie Differenzierbarkeit, an $\Gamma (z)$ , so kann diese schließlich eindeutig charakterisiert werden, womit „die“ verallgemeinerte Fakultät gefunden ist. Zwar gibt es zahlreiche andere komplexwertige Funktionen $f$ , die $f(n)=(n-1)!$ für alle natürlichen $n$ erfüllen, durch das Fordern zusätzlicher (als natürlich angesehener) Eigenschaften an $f$ bleibt jedoch nur noch $f=\Gamma$ übrig.

Es gilt dann $({\tfrac {1}{2}})!=\Gamma ({\tfrac {3}{2}})={\tfrac {\sqrt {\pi }}{2}}\approx 0{,}88622$ mit der Kreiszahl $\pi$ . Dieser Zusammenhang lässt sich über das Gaußsche Fehlerintegral erklären, das u. a. mit der Normalverteilung aus der Wahrscheinlichkeitstheorie zusammenhängt.

Geschichte

Als früheste Definition der Gammafunktion gilt die in einem Brief von Daniel Bernoulli an Christian Goldbach vom 6. Oktober 1729 gegebene:

{\Bigl (}A+{\frac {x}{2}}{\Bigr )}^{x-1}{\Bigl (}{\frac {2}{1+x}}\cdot {\frac {3}{2+x}}\cdot {\frac {4}{3+x}}\cdots {\frac {A}{A-1+x}}{\Bigr )}

für unendlich große $A$ , entsprechend heutiger Notation $x!$ oder $\Gamma (x+1)$ . Wenige Tage später, am 13. Oktober^jul. / 24. Oktober 1729^greg., beschrieb Euler ebenfalls in einem Brief an Goldbach die ähnliche, etwas einfachere Formel

{\frac {1\cdot 2\cdot 3\dotsm n}{(1+m)(2+m)\dotsm (n+m)}}\,(n+1)^{m},

Diese von Leonhard Euler entdeckte Formel für die Gammafunktion kann direkt als unendliche Produktreihe in ihrer standardisierten Form so dargestellt werden:

\Gamma (x)={\frac {1}{x}}\prod _{n=1}^{\infty }{\bigl (}1+{\frac {1}{n}}{\bigr )}^{x}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}

\Gamma (x)=\prod _{n=1}^{\infty }{\bigl (}1+{\frac {1}{n}}{\bigr )}^{x-1}{\bigl [}1+(x-1){\frac {1}{n}}{\bigr ]}^{-1}

Sie ist für alle komplexen Zahlen $x\in \mathbb {C}$ gültig und wurde 1812 von Gauß (für den Fall komplexer Zahlen) wiederentdeckt (die genannten Briefe wurden erst 1843 herausgegeben). Am 8. Januar 1730 beschrieb Euler in einem Brief an Goldbach folgendes Integral zur Interpolation der Fakultätsfunktion, das er am 28. November 1729 der St. Petersburger Akademie vorgestellt hatte:

\int \!\mathrm {d} x(-lx)^{n},

in heutiger Notation:

\displaystyle \Gamma (n+1)=\int _{0}^{1}(-\log x)^{n}\mathrm {d} x

Diese Definition wurde von Euler später bevorzugt verwendet und geht durch die Substitution $t=-\log x$ in die Form

\Gamma (n+1)=\int _{0}^{\infty }t^{n}\mathrm {e} ^{-t}\mathrm {d} t

über. Euler entdeckte dieses Integral bei der Untersuchung eines Problems aus der Mechanik, bei dem die Beschleunigung eines Teilchens betrachtet wird.

Adrien-Marie Legendre führte 1809 die griechische Majuskel $\Gamma$ (Gamma) als Funktionssymbol ein. Gauß verwendete 1812 das Funktionssymbol $\Pi$ (Pi) so, dass $\Pi (x)=\Gamma (x+1)$ und somit auch $\Pi (n)=n!$ für nichtnegative ganzzahlige $n$ gilt. Es setzte sich jedoch nicht durch; heute wird $\Pi$ als Symbol für ein Produkt benutzt (analog zu $\Sigma$ für eine Summe).

Definition und elementare Darstellungsformen

Es gibt in der Literatur keine einheitliche Definition für die Gammafunktion.

Häufig wird das Eulersche Integral zweiter Gattung gegeben. Ein Nachteil ist, dass dieses Integral nicht überall konvergiert. Somit ist eine globale Berechnung mittels dieser Definition nur indirekt möglich. Für komplexe Zahlen $z$ mit positivem Realteil ist die Gammafunktion damit das uneigentliche Integral

\Gamma (z)=\int _{0}^{\infty }t^{z-1}{\mathrm {e} }^{-t}\mathrm {d} t.

Die dadurch definierte Funktion ist holomorph, da das Integral (wegen des schnellen Abfallens der Exponentialfunktion) auf kompakten Mengen gleichmäßig konvergiert. Dies ermöglicht den Einsatz des Weierstraßschen Konvergenzsatzes. Mittels meromorpher Fortsetzung lässt sich $\Gamma (z)$ schließlich für alle Werte $\mathbb {C} \setminus \{0,-1,-2,\dotsc \}$ berechnen.

Eine andere Darstellung mittels eines Produktes motiviert die Verallgemeinerung der Fakultät auf direkte Weise. Sie ist gegeben durch:

\Gamma (z)=\lim _{n\to \infty }{\frac {n!\,n^{z}}{z(z+1)(z+2)\dotsm (z+n)}}.

In seinem Buch Number Theory. Analytic and modern tools gibt Henri Cohen eine Definition mittels der Hurwitzschen Zeta-Funktion $\zeta (s,z)$ . Als Begründung hierfür wird eine „einfache Möglichkeit der Verallgemeinerung“ und die „Betonung wichtiger Formeln“ angegeben. Es gilt demnach für komplexe Zahlen $z$ mit positivem Realteil

\Gamma (z)=\exp(\zeta '(0,\ z)-\zeta '(0,\ 1)),

wobei die Ableitung bezüglich der ersten Variablen gebildet ist.

Globale Eigenschaften

Funktionalgleichung und Meromorphie

Die Gammafunktion erfüllt in ihrem Definitionsbereich für alle $z$ die Funktionalgleichung

z\;\Gamma (z)=\Gamma (z+1).

Mittels dieser Relation ist eine induktive Fortsetzung (beispielsweise des Eulerschen Integrals) möglich. Es gilt für alle $n=0,1,2,\dotsc$

\Gamma (z)={\frac {\Gamma (z+n+1)}{z(z+1)(z+2)\cdots (z+n)}}.

Nullstellen und Polstellen

Aus der vorherigen Darstellung kann gefolgert werden, dass $\Gamma (z)$ zu einer auf $\mathbb {C}$ meromorphen Funktion fortgesetzt werden kann, die Pole an den Stellen $z=0,-1,-2,\dotsc$ besitzt. Alle Pole sind einfach und besitzen das Residuum

\operatorname {Res} _{z=-n}\Gamma (z)={\frac {(-1)^{n}}{n!}}

,

hierbei ist $n\in \mathbb {N} _{0}$ . Nullstellen besitzt $\Gamma$ keine. Das macht $\Gamma ^{-1}$ zu einer ganzen Funktion mit ausschließlich einfachen Nullstellen.

Der Satz von Hölder

Der Satz von Hölder (Otto Hölder 1886) ist ein Negativresultat und besagt, dass die Gammafunktion keine algebraische Differentialgleichung erfüllt, deren Koeffizienten rationale Funktionen sind. Das heißt, es gibt keine Differentialgleichung der Form $f(z,y(z),y'(z),\dotsc ,y^{(n)}(z))=0$ mit einer nichtnegativen ganzen Zahl $n$ und einem Polynom $f\neq 0$ in $y,y',\dotsc ,y^{(n)}$ , dessen Koeffizienten rationale Funktionen von $z$ sind, und der Lösung $y=\Gamma$ .

Axiomatische Charakterisierung

Fortsetzung der Fakultät

Die Bedingungen $G(1)=1$ und $G(x+1)=x\cdot G(x)$ , die die Fakultät für natürliche Zahlen eindeutig beschreiben, werden auch von anderen analytischen Funktionen als der Gammafunktion erfüllt. Für positive $x$ erfüllt beispielsweise die Funktion

G(x)=\Gamma (x)\cdot {\bigl (}1+c\,\sin(2\pi x){\bigr )}

für $0<c<1$ die charakteristischen Bedingungen der Gammafunktion. Weierstraß fügte 1854 daher die notwendige und hinreichende Bedingung

\lim _{n\to \infty }{\frac {G(x+n)}{G(n)\,n^{x}}}=1

hinzu, womit aber die Suche nach einer möglichst elementaren oder natürlichen charakterisierenden Eigenschaft nicht beendet war. Emil Artin diskutierte 1931 die mögliche Kennzeichnung durch Funktionalgleichungen.

Der Satz von Bohr-Mollerup

Der Satz von Bohr-Mollerup (Harald Bohr und Johannes Mollerup 1922) erlaubt eine einfache Charakterisierung der Gammafunktion:

Eine Funktion

G\colon \mathbb {R} _{>0}\to \mathbb {R} _{>0}

ist in diesem Bereich genau dann gleich der Gammafunktion, wenn gilt:

$G(1)=1,$
$G(x+1)=x\cdot G(x),$
$G$ ist logarithmisch konvex, das heißt, $x\mapsto \log G(x)$ ist eine konvexe Funktion.

Diese Axiome sind bei Nicolas Bourbaki der Ausgangspunkt für die Darstellung der Theorie der Gammafunktion.

Der Satz von Wielandt

Der Satz von Wielandt über die Gammafunktion (Helmut Wielandt 1939) charakterisiert die Gammafunktion als holomorphe Funktion und besagt:

Eine holomorphe Funktion

G

, definiert auf einem Gebiet

D

, das den Streifen

S=\{x\in \mathbb {C} \mid 1\leq \operatorname {Re} (x)<2\}

enthält, ist genau dann gleich der Gammafunktion auf

D

, wenn gilt:

$G(1)=1,$
$G(x+1)=x\cdot G(x),$
$|G|$ ist auf dem Streifen $S$ beschränkt, das heißt, es existiert ein $c>0$ , sodass $|G(x)|<c$ für alle $x$ aus $S$ .

Genauer gilt $|\Gamma (x)|\leq \Gamma (\operatorname {Re} (x))$ für alle $x$ mit $\operatorname {Re} (x)>0$ .

Weitere Darstellungsformen

Gaußsche und Weierstraßsche Darstellung

Neben der Darstellung der Gammafunktion aus der Definition gibt es noch andere äquivalente Darstellungen. Eine direkte Definition von $\Gamma (x)$ für alle $x\in \mathbb {C} \setminus \{0,-1,-2,\dotsc \}$ gibt die Produktdarstellung der Gammafunktion nach Gauß,

\Gamma (x)=\lim _{n\to \infty }{\frac {n!\,n^{x}}{x(x+1)(x+2)\dotsm (x+n)}},

die für positive reelle Zahlen bereits von Euler 1729 angegeben wurde. Daraus abgeleitet ist die Darstellung von $1/\Gamma$ als Weierstraß-Produkt:

1/\Gamma (x)=x\cdot \prod _{n=1}^{\infty }\left(1+{\frac {x}{n}}\right)\mathrm {e} ^{-x\log({\frac {n+1}{n}})}=x\cdot \mathrm {e} ^{\gamma \,x}\cdot \prod _{n=1}^{\infty }\left(1+{\frac {x}{n}}\right)\mathrm {e} ^{-x/n}

mit der Euler-Mascheroni-Konstanten $\gamma =\lim _{n\to \infty }{\bigl (}({\tfrac {1}{1}}+{\tfrac {1}{2}}+{\tfrac {1}{3}}+\dotsb +{\tfrac {1}{n}})-\log n{\bigr )}$ . Das zweite Produkt wird üblicherweise als Weierstraßsche Darstellung bezeichnet, Karl Weierstraß verwendete jedoch nur das erste.

Eulersche Darstellung

Die Integraldarstellung aus der Definition geht ebenfalls auf Euler 1729 zurück, sie gilt allgemeiner für komplexe Zahlen mit positivem Realteil:

\Gamma (x)=\int _{0}^{\infty }t^{x-1}\mathrm {e} ^{-t}\,\mathrm {d} t,

wenn

\operatorname {Re} (x)>0.

Durch die Zerlegung dieses Integrals folgerte E. F. Prym 1876 eine in ganz $\mathbb {C} \setminus \{0,-1,-2,-3,\dotsc \}$ gültige Darstellung:

\Gamma (x)=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{n!(n+x)}}+\int _{1}^{\infty }t^{x-1}e^{-t}\mathrm {d} t

Eine andere Variante der Eulerschen Integraldarstellung gibt es für $x\in \mathbb {C}$ mit $0<\operatorname {Re} (x)<1$ :

\Gamma (x)=\mathrm {e} ^{\pi \mathrm {i} x/2}\int _{0}^{\infty }t^{x-1}\mathrm {e} ^{-\mathrm {i} t}\,\mathrm {d} t

Aus dieser Darstellung lassen sich zum Beispiel auf elegante Weise die Fresnelschen Integralformeln ableiten.

Hankelsche Darstellung

Der deutsche Mathematiker Hermann Hankel entdeckte eine Integralrepresentation für den Kehrwert der Gammafunktion:

{\frac {1}{\Gamma (x)}}={\frac {1}{2\pi }}\,x^{1-x}\int _{-\pi }^{\pi }\exp {\bigl \{}xy\cot(y)-x\ln {\bigl [}y\csc(y){\bigr ]}{\bigr \}}\,\mathrm {d} y

Jedoch ist die nun gezeigte Definition nur für positive Werte x gültig. Dieses Integral behandelten ebenso die Mathematiker Thomas Schmelzer und Lloyd Trefethen in ihrem Aufsatz Computing the Gamma Function using contour integrals and rational approximations aus dem Jahre 2007.

Diese Definition kann durch innere Substitution nach dem Muster $y=2\arctan(z)$ so umgeformt werden:

{\frac {1}{\Gamma (x)}}={\frac {1}{2\pi }}\,x^{1-x}\int _{-\infty }^{\infty }{\frac {2}{1+z^{2}}}\exp {\bigl \{}{\frac {x}{z}}(1-z^{2})\arctan(z)-x\ln {\bigl [}{\frac {1}{z}}(1+z^{2})\arctan(z){\bigr ]}{\bigr \}}\,\mathrm {d} z

{\frac {1}{\Gamma (x)}}={\frac {1}{2\pi }}\,x^{1-x}\int _{-\infty }^{\infty }{\frac {2\exp {\bigl [}x(1-z^{2})\arctan(z)/z{\bigr ]}}{(1+z^{2})^{x+1}}}{\biggl [}{\frac {z}{\arctan(z)}}{\biggr ]}^{x}\,\mathrm {d} z

Darstellung nach Whittaker und Watson

Für den natürlichen Logarithmus aus der Gammafunktion existieren auch einige Integralrepresentationen für die Gammafunktion. Eine solche Integralrepresentation wurde durch die britischen Mathematiker Edmund Taylor Whittaker und George Neville Watson entdeckt:

\ln {\bigl [}\Gamma (x){\bigr ]}={\bigl (}x-{\frac {1}{2}}{\bigr )}\ln(x)-x+{\frac {1}{2}}\ln(2\pi )+2x\int _{0}^{\infty }{\frac {\arctan(y)}{\exp(2\pi xy)-1}}\,\mathrm {d} y

Diese Formel kann ebenso mit Hilfe der Abel-Plana-Summenformel hergeleitet werden.

Kummersche Reihen

Ernst Eduard Kummer gab 1847 die Fourierentwicklung der logarithmischen Gammafunktion an:

\log \Gamma (x)=\left({\tfrac {1}{2}}-x\right){\bigl (}\gamma +\log(2\pi ){\bigr )}+{\frac {1}{2}}\log {\frac {\pi }{\sin(\pi x)}}+{\frac {1}{\pi }}\sum _{k=2}^{\infty }{\frac {\log k}{k}}\sin(2\pi kx)

für

0<x<1

Sie heißt auch Kummersche Reihe. Bereits 1846 fand Carl Johan Malmstén eine ähnliche Reihe:

\log {\frac {\Gamma ({\tfrac {1}{2}}+x)}{\Gamma ({\tfrac {1}{2}}-x)}}=-2x\,{\bigl (}\gamma +\log(2\pi ){\bigr )}+{\frac {2}{\pi }}\sum _{k=2}^{\infty }(-1)^{k}{\frac {\log k}{k}}\sin(2\pi kx)

für

-{\tfrac {1}{2}}<x<{\tfrac {1}{2}}

Harmonische Reihe

Gegeben ist diese Identität für die Harmonische Reihenfunktion:

\operatorname {H} (x)=\int _{0}^{\infty }{\frac {1-\exp(-xy)}{\exp(y)-1}}\,\mathrm {d} y

Deswegen ist folgende Integralidentität für den Logarithmus naturalis der Fakultätsfunktion gültig:

\gamma \,x+\ln {\bigl [}\Gamma (x+1){\bigr ]}=\int _{0}^{\infty }{\frac {\exp(-xy)+xy-1}{y{\bigl [}\exp(y)-1{\bigr ]}}}\,\mathrm {d} y

Aus der gezeigten Formel kann das Element der Mascheroni-Konstante so entfernt werden:

\ln {\bigl [}\Gamma (x+1){\bigr ]}=\int _{0}^{\infty }{\frac {1}{y}}{\biggl \{}x\exp(-y)-{\frac {1-\exp(-xy)}{\exp(y)-1}}{\biggr \}}\,\mathrm {d} y

Für nähere Herleitungen siehe den Artikel Euler-Mascheroni-Konstante!

Für die Debyeschen Funktionen gilt:

\int _{0}^{\infty }{\frac {x^{z}}{\exp(x)-1}}\,\mathrm {d} x=\Gamma (z+1)\,\zeta (z+1)

Die zuvor genannte Integralidentität für die Harmonische Reihenfunktion kann so dargestellt werden:

\operatorname {H} (x)=\int _{0}^{\infty }{\frac {1-\exp(-xy)}{\exp(y)-1}}\,\mathrm {d} y=\int _{0}^{\infty }{\frac {1}{\exp(y)-1}}{\biggl [}\sum _{m=1}^{\infty }{\frac {(xy)^{2m-1}}{\Gamma (2m)}}-{\frac {(xy)^{2m}}{\Gamma (2m+1)}}{\biggr ]}\mathrm {d} y

=\sum _{m=1}^{\infty }\int _{0}^{\infty }{\frac {1}{\exp(y)-1}}{\biggl [}{\frac {(xy)^{2m-1}}{\Gamma (2m)}}-{\frac {(xy)^{2m}}{\Gamma (2m+1)}}{\biggr ]}\mathrm {d} y=

=\sum _{m=1}^{\infty }{\biggl [}{\frac {x^{2m-1}}{\Gamma (2m)}}\int _{0}^{\infty }{\frac {y^{2m-1}}{\exp(y)-1}}\,\mathrm {d} y-{\frac {x^{2m}}{\Gamma (2m+1)}}\int _{0}^{\infty }{\frac {y^{2m}}{\exp(y)-1}}\,\mathrm {d} y{\biggr ]}=

=\sum _{m=1}^{\infty }{\biggl [}{\frac {\Gamma (2m)\zeta (2m)x^{2m-1}}{\Gamma (2m)}}-{\frac {\Gamma (2m+1)\zeta (2m+1)x^{2m}}{\Gamma (2m+1)}}{\biggr ]}=

=\sum _{m=1}^{\infty }{\bigl [}\zeta (2m)x^{2m-1}-\zeta (2m+1)x^{2m}{\bigr ]}

Die folgende Formel kann darauf aufgestellt werden:

\gamma \,x+\ln {\bigl [}\Gamma (x+1){\bigr ]}=\sum _{m=1}^{\infty }{\biggl [}{\frac {\zeta (2m)}{2m}}\,x^{2m}-{\frac {\zeta (2m+1)}{2m+1}}\,x^{2m+1}{\biggr ]}

Jedoch ist diese Formel nur für Werte $|x|\leq 1$ gültig beziehungsweise konvergent.

Außerdem gilt folgende verallgemeinerte Identität für die Mascheronische Konstante:

\gamma \,x+\ln {\bigl [}\Gamma (x+1){\bigr ]}=\sum _{n=1}^{\infty }{\biggl [}{\frac {x}{n}}-\ln {\biggl (}1+{\frac {x}{n}}{\biggr )}{\biggr ]}

Die soeben genannte Formel mit der Riemannschen Zetafunktion geht dann durch Darstellung der soeben gezeigten Formel mittels Stammfunktion der geometrischen Reihe und anschließenden Einsatz der Definition der Riemannschen Zetafunktion hervor:

\zeta (s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}

Grundlegende Funktionalgleichungen

Die Gammafunktion genügt der Funktionalgleichung

\Gamma (x+1)=x\cdot \Gamma (x)

mit

\Gamma (1)=1.

Mit dem Ergänzungssatz der Gammafunktion (Euler 1749)

\Gamma (x)\cdot \Gamma (1-x)={\frac {\pi }{\sin(\pi x)}}

für

x\in \mathbb {C} \setminus \mathbb {Z}

erhält man $\Gamma ({\tfrac {1}{2}})={\sqrt {\pi }}=1{,}77245\,38509\,05516\,02729\dotso$ (Folge A002161 in OEIS) sowie

\Gamma (-n+{\tfrac {1}{2}})={\frac {n!\,(-4)^{n}}{(2n)!}}\,{\sqrt {\pi }}

und

\Gamma (n+{\tfrac {1}{2}})={\frac {(2n)!}{n!\,4^{n}}}\,{\sqrt {\pi }}

für

n=0,1,2,\dotsc

Mit allgemeiner gewähltem $n$ wird aus der letzten Formel die Legendresche Verdopplungsformel (Legendre 1809)

\Gamma \left({\frac {x}{2}}\right)\cdot \Gamma \left({\frac {x+1}{2}}\right)={\frac {\sqrt {\pi }}{2^{x-1}}}\cdot \Gamma (x)

für

x\in \mathbb {C} \setminus \{0,-1,-2,\dotsc \}.

Diese ist ein Spezialfall der Gaußschen Multiplikationsformel (Gauß 1812)

\Gamma \left({\frac {x}{n}}\right)\cdot \Gamma \left({\frac {x+1}{n}}\right)\cdots \Gamma \left({\frac {x+n-1}{n}}\right)={\frac {(2\pi )^{(n-1)/2}}{n^{\,x-1/2}}}\cdot \Gamma (x)

für

n=1,\,2,\,3,\,\ldots

und

x\in \mathbb {C} \setminus \{0,-1,-2,\dotsc \}.

Gammafunktionswerte von Brüchen

Informationen über elliptische Gammafunktionswerte von Brüchen

Gregory Chudnovsky zeigte 1975, dass jede der Zahlen $\Gamma (1/6)$ , $\Gamma (1/4)$ , $\Gamma (1/3)$ , $\Gamma (2/3)$ , $\Gamma (3/4)$ und $\Gamma (5/6)$ transzendent und algebraisch unabhängig von $\pi$ ist. Sie sind nicht elementar darstellbar, aber können sehr wohl über algebraische Kombinationen von vollständigen elliptischen Integralen erster und zweiter Art dargestellt werden. Hingegen ist beispielsweise von dem Funktionswert $\Gamma (1/5)=4{,}59084\,37119\,98803\,05320\,\dotso$ (Folge A175380 in OEIS) nicht einmal bekannt, ob er irrational ist. Und bei diesem Wert ist eine Darstellung aus einer algebraischen Kombination von vollständigen elliptischen Integralen erster und zweiter Art und aus algebraischen Vorfaktoren als einzige Komponenten in der betroffenen Darstellung nicht möglich. Wenn aber vollständige elliptische Integrale erster Art oder zweiter Art selbst durch eine algebraische Kombination von Gammafunktionswerten rationaler Zahlen dargestellt werden können, dann ist der elliptische Modul von den betroffenen vollständigen elliptischen Integralen komplett immer ein Lambda-Stern-Funktionswert von einer rationalen Zahl. Solche elliptischen Integrale werden im deutschen Sprachraum als Singuläre Elliptische Integralwerte und im englischen Sprachraum als Elliptic Integral Singular Values bezeichnet.

Satz von Fubini

In den anschließenden Abschnitten dieses Artikels wird der Satz von Fubini in einer abgewandelten Form eingesetzt, die nun im Folgenden hergeleitet werden soll:

${\frac {\mathrm {d} }{\mathrm {d} x}}\int _{0}^{1}x\,f(xy)\,\mathrm {d} y=f(x)$ ${\frac {\mathrm {d} }{\mathrm {d} x}}\int _{0}^{1}x\,g(xy)\,\mathrm {d} y=g(x)$ $\lim _{x\rightarrow \infty }{\biggl [}\int _{0}^{1}x\,f(xy)\,\mathrm {d} y{\biggr ]}=\int _{0}^{\infty }f(x)\,\mathrm {d} x$ $\lim _{x\rightarrow \infty }{\biggl [}\int _{0}^{1}x\,g(xy)\,\mathrm {d} y{\biggr ]}=\int _{0}^{\infty }g(x)\,\mathrm {d} x$

Unter Anwendung der Produktregel kann dieser Ausdruck formuliert werden:

{\frac {\mathrm {d} }{\mathrm {d} x}}{\biggl [}\int _{0}^{1}x\,f(xy)\,\mathrm {d} y{\biggr ]}{\biggl [}\int _{0}^{1}x\,g(xy)\,\mathrm {d} y{\biggr ]}=f(x){\biggl [}\int _{0}^{1}x\,g(xy)\,\mathrm {d} y{\biggr ]}+g(x){\biggl [}\int _{0}^{1}x\,f(xy)\,\mathrm {d} y{\biggr ]}

Mit dem Hauptsatz der Differential- und Integralrechnung kann folgende Gleichungskette aufgestellt werden:

\int _{0}^{\infty }{\biggl \{}{\frac {\mathrm {d} }{\mathrm {d} x}}{\biggl [}\int _{0}^{1}x\,f(xy)\,\mathrm {d} y{\biggr ]}{\biggl [}\int _{0}^{1}x\,g(xy)\,\mathrm {d} y{\biggr ]}{\biggr \}}\mathrm {d} x={\biggl \{}{\biggl [}\int _{0}^{1}x\,f(xy)\,\mathrm {d} y{\biggr ]}{\biggl [}\int _{0}^{1}x\,g(xy)\,\mathrm {d} y{\biggr ]}{\biggr \}}_{x=0}^{x=\infty }=

={\biggl [}\int _{0}^{\infty }f(x)\,\mathrm {d} x{\biggr ]}{\biggl [}\int _{0}^{\infty }g(x)\,\mathrm {d} x{\biggr ]}

Und mit dem Satz von Fubini kann diese Gleichungskette aufgestellt werden:

\int _{0}^{\infty }{\biggl \{}{\frac {\mathrm {d} }{\mathrm {d} x}}{\biggl [}\int _{0}^{1}x\,f(xy)\,\mathrm {d} y{\biggr ]}{\biggl [}\int _{0}^{1}x\,g(xy)\,\mathrm {d} y{\biggr ]}{\biggr \}}\mathrm {d} x=\int _{0}^{\infty }{\biggl \{}f(x){\biggl [}\int _{0}^{1}x\,g(xy)\,\mathrm {d} y{\biggr ]}+g(x){\biggl [}\int _{0}^{1}x\,f(xy)\,\mathrm {d} y{\biggr ]}{\biggr \}}\mathrm {d} x=

=\int _{0}^{\infty }\int _{0}^{1}x\,f(x)\,g(xy)+x\,g(x)\,f(xy)\,\mathrm {d} y\,\mathrm {d} x=\int _{0}^{1}\int _{0}^{\infty }x\,f(x)\,g(xy)+x\,g(x)\,f(xy)\,\mathrm {d} x\,\mathrm {d} y

Aus den beiden soeben genannten Gleichungsketten folgt diese Gleichung:

${\color {blue}{\biggl [}\int _{0}^{\infty }f(x)\,\mathrm {d} x{\biggr ]}{\biggl [}\int _{0}^{\infty }g(x)\,\mathrm {d} x{\biggr ]}=\int _{0}^{1}\int _{0}^{\infty }x\,f(x)\,g(xy)+x\,g(x)\,f(xy)\,\mathrm {d} x\,\mathrm {d} y}$

Wenn die Funktion $g(x)$ mit der Funktion $f(x)$ gleichgesetzt wird, dann entsteht dieser Spezialfall:

${\color {blueviolet}{\biggl [}\int _{0}^{\infty }f(x)\,\mathrm {d} x{\biggr ]}^{2}=\int _{0}^{1}\int _{0}^{\infty }2\,x\,f(x)\,f(xy)\,\mathrm {d} x\,\mathrm {d} y}$

Beweise für den elementaren Wert Gamma(1/2)

Der erste nun folgende Beweis für Gamma(1/2) wird über das Wallissche Produkt absolviert:

Das Wallissche Produkt lässt sich auf folgende Weise darstellen: $\prod _{k=1}^{\infty }{\frac {4k(k+1)}{(2k+1)^{2}}}={\frac {\pi }{4}}$ Folgender Bruch hat folgenden Grenzwert: $\lim _{n\rightarrow \infty }{\frac {\Gamma (n+1){\sqrt {n+1}}}{\Gamma (n+3/2)}}=\lim _{n\rightarrow \infty }{\frac {\sqrt {\Gamma (n+1)\Gamma (n+2)}}{\Gamma (n+3/2)}}=1$ Für alle n ∈ ℕ gelten folgende Ausdrücke: $\Gamma (n+1){\sqrt {n+1}}=\prod _{k=1}^{n}{\sqrt {k(k+1)}}$ ${\frac {1}{\Gamma (n+3/2)}}={\frac {1}{\Gamma (3/2)}}\prod _{k=1}^{n}{\frac {2}{2k+1}}$ Folglich gilt diese Formel: $\lim _{n\rightarrow \infty }\left[\prod _{k=1}^{n}{\sqrt {k(k+1)}}{\frac {1}{\Gamma (3/2)}}\prod _{k=1}^{n}{\frac {2}{2k+1}}\right]=1$ Die Formel wird nach Γ(3/2) aufgelöst: $\Gamma (3/2)=\lim _{n\rightarrow \infty }\left[\prod _{k=1}^{n}{\sqrt {k(k+1)}}\prod _{k=1}^{n}{\frac {2}{2k+1}}\right]=$ $=\lim _{n\rightarrow \infty }\left[\prod _{k=1}^{n}{\sqrt {k(k+1)}}{\frac {2}{2k+1}}\right]=\prod _{k=1}^{\infty }{\sqrt {k(k+1)}}{\frac {2}{2k+1}}=$ $={\sqrt {\prod _{k=1}^{\infty }{\frac {4k(k+1)}{(2k+1)^{2}}}}}={\sqrt {\frac {\pi }{4}}}={\frac {1}{2}}{\sqrt {\pi }}$ Daraus folgt: $\Gamma {\bigl (}{\tfrac {1}{2}}{\bigr )}=2\,\Gamma {\bigl (}{\tfrac {3}{2}}{\bigr )}={\sqrt {\pi }}$

Der zweite Beweis für Gamma(1/2) wird über den Satz von Fubini mit der zuvor gezeigten Formel bewerkstelligt:

Der Funktionswert Gamma(1/2) taucht als Integral der Gaußschen Glockenkurve auf: $\Gamma \left({\tfrac {1}{2}}\right)=\int _{0}^{\infty }x^{-1/2}\exp(-x)\,\mathrm {d} x=2\int _{0}^{\infty }\exp \left(-x^{2}\right)\,\mathrm {d} x$ Die oben genannte Formel lautet: ${\color {blueviolet}{\biggl [}\int _{0}^{\infty }f(x)\,\mathrm {d} x{\biggr ]}^{2}=\int _{0}^{1}\int _{0}^{\infty }2\,x\,f(x)\,f(xy)\,\mathrm {d} x\,\mathrm {d} y}$ Eingesetzt entsteht dann folgendes Resultat: ${\color {blueviolet}\left[\int _{0}^{\infty }\exp \left(-x^{2}\right)\,\mathrm {d} x\right]^{2}=\int _{0}^{1}\int _{0}^{\infty }2x\exp \left(-x^{2}\right)\exp \left(-x^{2}y^{2}\right)\,\mathrm {d} x\,\mathrm {d} y=}$ ${\color {blueviolet}=\int _{0}^{1}\int _{0}^{\infty }2x\exp \left[-x^{2}(y^{2}+1)\right]\,\mathrm {d} x\,\mathrm {d} y=\int _{0}^{1}{\frac {1}{y^{2}+1}}\mathrm {d} y={\frac {\pi }{4}}}$ Daraus folgt ebenso: $\Gamma \left({\tfrac {1}{2}}\right)={\sqrt {\pi }}$

Ein dritter Beweis beinhaltet die Dawsonsche Funktion und ähnelt dem zweiten Beweis:

Allgemein gilt für alle natürlichen Zahlen $n\in \mathbb {N}$ diese Formel: $\int _{0}^{\infty }\exp(-x^{n})\,\mathrm {d} x={\frac {1}{n}}\,\Gamma {\bigl (}{\frac {1}{n}}{\bigr )}$ Für $n=2$ gilt dann das: $\int _{0}^{\infty }\exp(-x^{2})\,\mathrm {d} x={\frac {1}{2}}\,\Gamma {\bigl (}{\frac {1}{2}}{\bigr )}$ Mit der Dawson-Minus-Funktion lässt sich das nun gezeigte Integral beschreiben: $\int _{0}^{\infty }\exp(-x^{2})\,\mathrm {d} x=\lim _{x\rightarrow \infty }\exp(-x^{2})\operatorname {D} _{-}(x)$ Die Dawson-Minus-Funktion hat diese Integralidentität: $2\exp(-2\,x^{2})\operatorname {D} _{-}(x)=\int _{0}^{1}2\,x\exp {\bigl [}-x^{2}(y^{2}+1){\bigr ]}\,\mathrm {d} y$ Durch die Bildung der Ursprungsstammfunktion von der nun genannten Formel bezüglich x entsteht diese Formel: $\exp(-2\,x^{2})\operatorname {D} _{-}(x)^{2}=\int _{0}^{1}{\frac {1-\exp {\bigl [}-x^{2}(y^{2}+1){\bigr ]}}{y^{2}+1}}\,\mathrm {d} y$ Über die Grenzwertbildung entsteht dann die anschließende Gleichung: $\lim _{x\rightarrow \infty }\exp(-2\,x^{2})\operatorname {D} _{-}(x)^{2}=\lim _{x\rightarrow \infty }\int _{0}^{1}{\frac {1-\exp {\bigl [}-x^{2}(y^{2}+1){\bigr ]}}{y^{2}+1}}\,\mathrm {d} y=\int _{0}^{1}{\frac {1}{y^{2}+1}}\,\mathrm {d} y={\frac {\pi }{4}}$ Daraus folgt dieses Resultat: $\int _{0}^{\infty }\exp(-x^{2})\,\mathrm {d} x={\frac {1}{2}}{\sqrt {\pi }}$ Und direkt daraus ergibt sich das Endresultat: $\Gamma {\bigl (}{\frac {1}{2}}{\bigr )}={\sqrt {\pi }}$

Beweis für die lemniskatischen Werte Gamma(1/4) und Gamma(3/4)

Mit der lemniskatischen Konstante $\varpi$ gilt diese Formel:

\Gamma \left({\frac {1}{4}}\right)={\sqrt {2\varpi \,{\sqrt {2\pi }}}}=2\,\pi ^{1/4}K({\tfrac {1}{2}}{\sqrt {2}})^{1/2}=3{,}62560\,99082\,21908\,31193\dotso

(Folge A068466 in OEIS).

Und wegen des Ergänzungssatzes und der Legendreschen Identität gilt:

\Gamma \left({\frac {3}{4}}\right)={\sqrt {2}}\pi /\Gamma \left({\frac {1}{4}}\right)={\frac {\sqrt[{4}]{\pi ^{3}}}{{\sqrt[{4}]{2}}{\sqrt {\varpi }}}}=\pi ^{1/4}{\bigl [}2\,E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}{\bigr ]}^{1/2}=1{,}22541\,67024\,65177\,64512\,90983\dotso

Hierbei ist K das vollständige elliptische Integral erster Ordnung:

K(\varepsilon )=\int _{0}^{\pi /2}{\frac {1}{\sqrt {1-\varepsilon ^{2}\sin(\varphi )^{2}}}}\,\mathrm {d} \varphi =\int _{0}^{1}{\frac {2}{\sqrt {(x^{2}+1)^{2}-4\varepsilon ^{2}x^{2}}}}\,\mathrm {d} x

Und E ist das vollständige elliptische Integral zweiter Ordnung:

E(\varepsilon )=\int _{0}^{\pi /2}{\sqrt {1-\varepsilon ^{2}\sin(\varphi )^{2}}}\,\mathrm {d} \varphi =\int _{0}^{1}{\frac {2{\sqrt {(x^{2}+1)^{2}-4\varepsilon ^{2}x^{2}}}}{(x^{2}+1)^{2}}}\,\mathrm {d} x

Diese Formeln werden mit dem oben beschriebenen Integrationsmechanismus in Folgenden bewiesen:

Basierend auf der Definition der Gammafunktion gilt:

\Gamma {\biggl (}{\frac {1}{4}}{\biggr )}=4\int _{0}^{\infty }\exp(-x^{4})\,\mathrm {d} x

\Gamma {\biggl (}{\frac {3}{4}}{\biggr )}=4\int _{0}^{\infty }x^{2}\exp(-x^{4})\,\mathrm {d} x

Die Kombination aus Satz von Fubini und Produktregel wird als Schlüssel verwendet:

${\color {blueviolet}{\biggl [}\int _{0}^{\infty }f(x)\,\mathrm {d} x{\biggr ]}^{2}=\int _{0}^{1}\int _{0}^{\infty }2\,x\,f(x)\,f(xy)\,\mathrm {d} x\,\mathrm {d} y}$

Dann gilt folgende Gleichungskette:

{\color {blueviolet}{\frac {1}{16}}\,\Gamma {\biggl (}{\frac {1}{4}}{\biggr )}^{2}={\biggl [}\int _{0}^{\infty }\exp(-x^{4})\,\mathrm {d} x{\biggr ]}^{2}=\int _{0}^{1}\int _{0}^{\infty }2x\exp(-x^{4})\exp(-x^{4}y^{4})\,\mathrm {d} x\,\mathrm {d} y=}

{\color {blueviolet}=\int _{0}^{1}\int _{0}^{\infty }2x\exp[-x^{4}(y^{4}+1)]\,\mathrm {d} x\,\mathrm {d} y=\int _{0}^{1}\int _{0}^{\infty }\exp[-x^{2}(y^{4}+1)]\,\mathrm {d} x\,\mathrm {d} y=}

{\color {blueviolet}=\int _{0}^{1}{\biggl [}{\frac {1}{\sqrt {y^{4}+1}}}\int _{0}^{\infty }\exp(-x^{2})\,\mathrm {d} x{\biggr ]}\,\mathrm {d} y={\color {ForestGreen}{\biggl [}\int _{0}^{1}{\frac {1}{\sqrt {y^{4}+1}}}\mathrm {d} y{\biggr ]}}{\biggl [}\int _{0}^{\infty }\exp(-x^{2})\,\mathrm {d} x{\biggr ]}=}

{\color {blueviolet}={\color {ForestGreen}{\frac {\varpi }{2{\sqrt {2}}}}}\int _{0}^{\infty }\exp(-x^{2})\,\mathrm {d} x={\frac {\varpi {\sqrt {\pi }}}{4{\sqrt {2}}}}}

Hierbei gelten folgende Formeln über den Arcussinus lemniscatus:

{\frac {1}{\sqrt {y^{4}+1}}}={\frac {\mathrm {d} }{\mathrm {d} y}}{\sqrt {2}}\,\mathrm {arcsl} {\bigl [}y{\bigl (}{\sqrt {y^{4}+1}}+1{\bigr )}^{-1/2}{\bigr ]}

\mathrm {arcsl} {\bigl (}{\sqrt {{\sqrt {2}}-1}}{\bigr )}={\frac {\varpi }{4}}

Auch mit folgender Stammfunktion kann das gezeigte Integral nachgewiesen werden:

{\color {ForestGreen}\int _{0}^{1}{\frac {1}{\sqrt {y^{4}+1}}}\,\mathrm {d} y={\biggl \{}{\frac {1}{2}}F{\biggl [}2\arctan(y);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}{\biggr \}}_{y=0}^{y=1}={\frac {1}{2}}K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}={\frac {\varpi }{2{\sqrt {2}}}}}

Die genannte Kombination aus Satz von Fubini und Produktregel stellt die Verbindung zwischen beiden Integralen her:

${\color {blue}{\biggl [}\int _{0}^{\infty }f(x)\,\mathrm {d} x{\biggr ]}{\biggl [}\int _{0}^{\infty }g(x)\,\mathrm {d} x{\biggr ]}=\int _{0}^{1}\int _{0}^{\infty }x\,f(x)\,g(xy)+x\,g(x)\,f(xy)\,\mathrm {d} x\,\mathrm {d} y}$

In Bezug auf den Eulerschen Ergänzungssatz gilt somit:

{\color {blue}{\frac {1}{16}}\,\Gamma {\biggl (}{\frac {1}{4}}{\biggr )}\,\Gamma {\biggl (}{\frac {3}{4}}{\biggr )}={\biggl [}\int _{0}^{\infty }\exp(-x^{4})\,\mathrm {d} x{\biggr ]}{\biggl [}\int _{0}^{\infty }x^{2}\exp(-x^{4})\,\mathrm {d} x{\biggr ]}=}

{\color {blue}=\int _{0}^{1}\int _{0}^{\infty }x{\bigl [}\exp(-x^{4}){\bigr ]}{\bigl [}x^{2}y^{2}\exp(-x^{4}y^{4}){\bigr ]}+x{\bigl [}\exp(-x^{4}y^{4}){\bigr ]}{\bigl [}x^{2}\exp(-x^{4}){\bigr ]}\,\mathrm {d} x\,\mathrm {d} y=}

{\color {blue}=\int _{0}^{1}\int _{0}^{\infty }x^{3}(y^{2}+1)\exp[-x^{4}(y^{4}+1)]\,\mathrm {d} x\,\mathrm {d} y=\int _{0}^{1}{\frac {y^{2}+1}{4(y^{4}+1)}}\,\mathrm {d} y={\frac {\pi }{8{\sqrt {2}}}}}

Aus diesen beiden Formelketten entsteht durch Zusammenfügung:

\Gamma {\biggl (}{\frac {1}{4}}{\biggr )}=2^{3/4}\pi ^{1/4}\varpi ^{1/2}

\Gamma {\biggl (}{\frac {1}{4}}{\biggr )}\,\Gamma {\biggl (}{\frac {3}{4}}{\biggr )}={\sqrt {2}}\,\pi

Und dies folgt darauf:

\Gamma {\biggl (}{\frac {3}{4}}{\biggr )}=2^{-1/4}\pi ^{3/4}\varpi ^{-1/2}

Beweis für die äquianharmonischen Werte Gamma(1/3) und Gamma(2/3)

Die Gammafunktionswerte der Drittel können ebenso mit Hilfe elliptischer Integrale erster und zweiter Ordnung dargestellt werden:

\Gamma \left({\frac {1}{3}}\right)={\frac {2^{7/9}}{3^{1/12}}}{\pi }^{1/3}{K[\sin({\frac {\pi }{12}})]}^{1/3}

\Gamma \left({\frac {2}{3}}\right)={\frac {2^{2/9}}{3^{5/12}}}{\pi }^{2/3}{K[\sin({\frac {\pi }{12}})]}^{-1/3}={\frac {2^{5/9}}{3^{5/12}}}\,\pi ^{1/3}{\bigl \{}{2{\sqrt {3}}E{\bigl [}\sin {\bigl (}{\frac {\pi }{12}}{\bigr )}{\bigr ]}-({\sqrt {3}}+1)K{\bigl [}\sin {\bigl (}{\frac {\pi }{12}}{\bigr )}{\bigr ]}}{\bigr \}}^{1/3}

Diese Formeln werden mit dem oben beschriebenen Integrationsmechanismus in Folgenden bewiesen:

Basierend auf der Definition der Gammafunktion gilt:

\Gamma {\biggl (}{\frac {1}{3}}{\biggr )}=3\int _{0}^{\infty }\exp(-x^{3})\,\mathrm {d} x

\Gamma {\biggl (}{\frac {2}{3}}{\biggr )}=3\int _{0}^{\infty }x\exp(-x^{3})\,\mathrm {d} x

Die Kombination aus Satz von Fubini und Produktregel wird auch hier als Schlüssel verwendet:

${\color {blueviolet}{\biggl [}\int _{0}^{\infty }f(x)\,\mathrm {d} x{\biggr ]}^{2}=\int _{0}^{1}\int _{0}^{\infty }2\,x\,f(x)\,f(xy)\,\mathrm {d} x\,\mathrm {d} y}$

Dann gilt Folgendes:

{\color {blueviolet}{\frac {1}{9}}\,\Gamma {\biggl (}{\frac {1}{3}}{\biggr )}^{2}={\biggl [}\int _{0}^{\infty }\exp(-x^{3})\,\mathrm {d} x{\biggr ]}^{2}=\int _{0}^{1}\int _{0}^{\infty }2x\exp(-x^{3})\exp(-x^{3}y^{3})\,\mathrm {d} x\,\mathrm {d} y=}

{\color {blueviolet}=\int _{0}^{1}\int _{0}^{\infty }2x\exp[-x^{3}(y^{3}+1)]\,\mathrm {d} x\,\mathrm {d} y=\int _{0}^{1}{\biggl [}{\frac {1}{(y^{3}+1)^{2/3}}}\int _{0}^{\infty }2x\exp(-x^{3})\,\mathrm {d} x{\biggr ]}\mathrm {d} y=}

{\color {blueviolet}={\color {ForestGreen}{\biggl [}\int _{0}^{1}{\frac {1}{(y^{3}+1)^{2/3}}}\mathrm {d} y{\biggr ]}}{\biggl [}\int _{0}^{\infty }2x\exp(-x^{3})\,\mathrm {d} x{\biggr ]}={\color {ForestGreen}{\frac {2}{{\sqrt[{3}]{4}}\,{\sqrt[{4}]{27}}}}K{\bigl [}\sin {\bigl (}{\frac {\pi }{12}}{\bigr )}{\bigr ]}}{\biggl [}\int _{0}^{\infty }2x\exp(-x^{3})\,\mathrm {d} x{\biggr ]}={\frac {2^{4/3}}{3^{7/4}}}K{\bigl [}\sin {\bigl (}{\frac {\pi }{12}}{\bigr )}{\bigr ]}\Gamma \left({\frac {2}{3}}\right)}

Mit folgender Kettenregelumformung und anschließender Stammfunktion kann das gezeigte Integral nachgewiesen werden:

{\color {ForestGreen}\int _{0}^{1}{\frac {1}{(y^{3}+1)^{2/3}}}\,\mathrm {d} y=\int _{0}^{1}{\biggl \{}{\biggl [}{\frac {\mathrm {d} }{\mathrm {d} z}}\,{\frac {z}{(1+{\sqrt {1-z^{3}}})^{2/3}}}{\biggr ]}{\frac {1}{(y^{3}+1)^{2/3}}}{\biggl [}y={\frac {z}{(1+{\sqrt {1-z^{3}}})^{2/3}}}{\biggr ]}{\biggr \}}\,\mathrm {d} z=}

{\color {ForestGreen}=\int _{0}^{1}{\frac {1}{2^{2/3}{\sqrt {1-z^{3}}}}}\,\mathrm {d} z={\biggl \{}{\frac {1}{{\sqrt[{3}]{4}}\,{\sqrt[{4}]{3}}}}F{\biggl [}2\arctan {\biggl (}{\frac {{\sqrt[{4}]{3}}z}{{\sqrt {1+z+z^{2}}}+{\sqrt {1-z}}}}{\biggr )};\cos {\bigl (}{\frac {\pi }{12}}{\bigr )}{\biggr ]}{\biggr \}}_{z=0}^{z=1}=}

{\color {ForestGreen}={\frac {2}{3{\sqrt[{3}]{4}}\,{\sqrt[{4}]{3}}}}K{\bigl [}\cos {\bigl (}{\frac {\pi }{12}}{\bigr )}{\bigr ]}={\frac {2}{{\sqrt[{3}]{4}}\,{\sqrt[{4}]{27}}}}K{\bigl [}\sin {\bigl (}{\frac {\pi }{12}}{\bigr )}{\bigr ]}}

Die genannte Kombination aus Satz von Fubini und Produktregel stellt ebenso an dieser Stelle die Verbindung zwischen beiden Integralen her:

${\color {blue}{\biggl [}\int _{0}^{\infty }f(x)\,\mathrm {d} x{\biggr ]}{\biggl [}\int _{0}^{\infty }g(x)\,\mathrm {d} x{\biggr ]}=\int _{0}^{1}\int _{0}^{\infty }x\,f(x)\,g(xy)+x\,g(x)\,f(xy)\,\mathrm {d} x\,\mathrm {d} y}$

In Bezug auf den Eulerschen Ergänzungssatz gilt deswegen analog zum vorherigen Beispiel:

{\color {blue}{\frac {1}{9}}\,\Gamma {\biggl (}{\frac {1}{3}}{\biggr )}\,\Gamma {\biggl (}{\frac {2}{3}}{\biggr )}={\biggl [}\int _{0}^{\infty }\exp(-x^{3})\,\mathrm {d} x{\biggr ]}{\biggl [}\int _{0}^{\infty }x\exp(-x^{3})\,\mathrm {d} x{\biggr ]}=}

{\color {blue}=\int _{0}^{1}\int _{0}^{\infty }x{\bigl [}\exp(-x^{3}){\bigr ]}{\bigl [}xy\exp(-x^{3}y^{3}){\bigr ]}+x{\bigl [}\exp(-x^{3}y^{3}){\bigr ]}{\bigl [}x\exp(-x^{3}){\bigr ]}\,\mathrm {d} x\,\mathrm {d} y=}

{\color {blue}=\int _{0}^{1}\int _{0}^{\infty }x^{2}(y+1)\exp[-x^{3}(y^{3}+1)]\,\mathrm {d} x\,\mathrm {d} y=\int _{0}^{1}{\frac {1}{3(y^{2}-y+1)}}\,\mathrm {d} y={\frac {2\,\pi }{9{\sqrt {3}}}}}

Aus diesen beiden Formelketten entsteht durch Zusammenfügung:

\Gamma {\biggl (}{\frac {1}{3}}{\biggr )}^{2}=2^{4/3}3^{1/4}K{\bigl [}\sin {\bigl (}{\frac {\pi }{12}}{\bigr )}{\bigr ]}\Gamma \left({\frac {2}{3}}\right)

\Gamma \left({\frac {1}{3}}\right)\Gamma \left({\frac {2}{3}}\right)={\frac {2\pi }{\sqrt {3}}}

Durch Kombination dieser beiden Formeln folgt darauf:

\Gamma \left({\frac {1}{3}}\right)={\frac {2^{7/9}}{3^{1/12}}}{\pi }^{1/3}{K[\sin({\frac {\pi }{12}})]}^{1/3}

\Gamma \left({\frac {2}{3}}\right)={\frac {2^{2/9}}{3^{5/12}}}{\pi }^{2/3}{K[\sin({\frac {\pi }{12}})]}^{-1/3}

Liste elliptischer Gammafunktionswerte

Generell gilt folgende Formel für alle $n\in \mathbb {N}$ :

{\frac {\Gamma (1+1/n)^{2}}{\Gamma (1+2/n)}}=\int _{0}^{1}{\frac {1}{(x^{n}+1)^{2/n}}}\mathrm {d} x=\int _{0}^{1}{\frac {1}{2^{2/n}{\sqrt {1-x^{n}}}}}\mathrm {d} x

Im letzten Schritt wird auf folgende Weise substituiert: $x\mapsto {\frac {x}{(1+{\sqrt {1-x^{n}}})^{2/n}}}$

Auf diese Weise lassen sich alle Gamma-Funktionswerte von rationalen Zahlen ermitteln.

Folgende weitere Funktionswerte der Gammafunktion lassen sich mit elliptischen Integralen erster und zweiter Ordnung darstellen:

$\Gamma (1/6)=2^{11/9}3^{1/3}\pi ^{1/6}K[\sin({\frac {\pi }{12}})]^{2/3}$ $\Gamma (5/6)=2^{4/9}3^{-1/3}\pi ^{1/6}\{{2{\sqrt {3}}E[\sin({\frac {\pi }{12}})]-({\sqrt {3}}+1)K[\sin({\frac {\pi }{12}})]}\}^{2/3}$

$\Gamma (1/8)=2^{17/8}\pi ^{1/8}K({\frac {1}{\sqrt {2}}})^{1/4}K({\sqrt {2}}-1)^{1/2}$ $\Gamma (3/8)=2^{11/8}({\sqrt {2}}-1)^{1/2}\,\pi ^{1/8}[2E({\frac {1}{\sqrt {2}}})-K({\frac {1}{\sqrt {2}}})]^{1/4}K({\sqrt {2}}-1)^{1/2}$ $\Gamma (5/8)=2^{5/8}\pi ^{1/8}K({\frac {1}{\sqrt {2}}})^{1/4}[{\sqrt {2}}E({\sqrt {2}}-1)-K({\sqrt {2}}-1)]^{1/2}$ $\Gamma (7/8)=2^{-1/8}({\sqrt {2}}+1)^{1/2}\pi ^{1/8}[2E({\frac {1}{\sqrt {2}}})-K({\frac {1}{\sqrt {2}}})]^{1/4}[{\sqrt {2}}E({\sqrt {2}}-1)-K({\sqrt {2}}-1)]^{1/2}$

$\Gamma (1/12)=2^{55/36}3^{7/24}({\sqrt {3}}+1)^{1/2}\pi ^{1/12}K({\frac {1}{\sqrt {2}}})^{1/2}K[\sin({\frac {\pi }{12}})]^{1/3}$ $\Gamma (5/12)=2^{29/36}3^{-1/24}({\sqrt {3}}-1)^{1/2}\pi ^{1/12}K({\frac {1}{\sqrt {2}}})^{1/2}\{{2{\sqrt {3}}E[\sin({\frac {\pi }{12}})]-({\sqrt {3}}+1)K[\sin({\frac {\pi }{12}})]}\}^{1/3}$ $\Gamma (7/12)=2^{19/36}3^{1/24}({\sqrt {3}}-1)^{1/2}\pi ^{1/12}[2E({\frac {1}{\sqrt {2}}})-K({\frac {1}{\sqrt {2}}})]^{1/2}K[\sin({\frac {\pi }{12}})]^{1/3}$ $\Gamma (11/12)=2^{-7/36}3^{-7/24}({\sqrt {3}}+1)^{1/2}\pi ^{1/12}[2E({\frac {1}{\sqrt {2}}})-K({\frac {1}{\sqrt {2}}})]^{1/2}\{{2{\sqrt {3}}E[\sin({\frac {\pi }{12}})]-({\sqrt {3}}+1)K[\sin({\frac {\pi }{12}})]}\}^{1/3}$

Basierend auf den weiter oben gezeigten Beweisformeln können die Werte für die Sechstel und Zwölftel mit Hilfe der Legendreschen Verdopplungsformel nachgewiesen werden. Weitere Beziehungen zwischen der Gammafunktion und den elliptischen Integralen:

{\frac {\Gamma (21/20)^{2}}{\Gamma (11/10)}}=\int _{0}^{1}{\frac {1}{\sqrt[{10}]{x^{20}+1}}}\,\mathrm {d} x=2^{9/10}5^{-1/2}\sin({\frac {\pi }{5}})^{1/2}\cos({\frac {\pi }{20}})K\{\sin[{\frac {1}{2}}\arcsin({\sqrt {5}}-2)]\}

{\frac {\Gamma (23/20)^{2}}{\Gamma (13/10)}}=2^{27/10}5^{-5/4}3\sin({\frac {\pi }{5}})^{5/2}\cos({\frac {3\pi }{20}})K\{\sin[{\frac {1}{2}}\arcsin({\sqrt {5}}-2)]\}

{\frac {\Gamma (25/24)^{2}}{\Gamma (13/12)}}=\int _{0}^{1}{\frac {1}{\sqrt[{12}]{x^{24}+1}}}\,\mathrm {d} x=2^{-13/12}3^{-1/4}({\sqrt {3}}-1)({\sqrt {2}}+1)K\{\tan[{\frac {1}{4}}\arcsin({\frac {1}{3}})]\}

{\frac {\Gamma (29/24)^{2}}{\Gamma (17/12)}}=2^{-17/12}3^{-3/4}5({\sqrt {3}}-1)K\{\tan[{\frac {1}{4}}\arcsin({\frac {1}{3}})]\}

Kurvendiskussion

Polstellen und Residuen

Die Gammafunktion hat an den Stellen $x=-n\ \left(n=0,1,2,\dotsc \right)$ Pole erster Ordnung. Aus der Funktionalgleichung erhält man für die Residuen

\operatorname {Res} _{x=-n}\Gamma (x)={\frac {(-1)^{n}}{n!}}.

Alternativ lassen sie sich direkt an der Formel

\Gamma (z)=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{n!}}{\frac {1}{z+n}}+\int _{1}^{\infty }t^{z-1}e^{-t}\mathrm {d} t

ablesen. Da $\Gamma$ keine Nullstellen hat, ist $1/\Gamma$ eine ganze Funktion.

Ableitung und Digammafunktion

Die Ableitung der Gammafunktion stimmt mit dem Produkt aus Gammafunktion und Digammafunktion überein:

{\frac {\mathrm {d} }{\mathrm {d} x}}\,\Gamma (x)=\Gamma '(x)=\psi (x)\,\Gamma (x)

Die Digammafunktion erhalt man, wenn man die Harmonische Reihenfunktion im Koordinatensystem um Eins nach rechts und um die Euler-Mascheroni-Konstante nach unten verschiebt:

\psi (x)=\mathrm {H} (x-1)-\gamma

Für die Gaußsche Pifunktion beziehungsweise Fakultätsfunktion und ihren Kehrwert sind nach Weierstraß diese Definitionen gültig:

\Pi (x)=\Gamma (x+1)=\exp(-\,\gamma \,x)\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}

\pi (x)={\frac {1}{\Pi (x)}}=\exp(\gamma \,x)\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}\exp {\bigl (}-{\frac {x}{n}}{\bigr )}{\bigr ]}

Mit diesen Definitionen kann die erste Ableitung der Gaußschen Pifunktion kenntlich gemacht werden. Die Ableitung eines unendlichen Produkts ist das Produkt aus diesem unendlichen Produkt multipliziert mit der unendlichen Summe aus den Quotienten der Ableitungen der Faktoren dividiert durch die Faktoren selbst:

{\frac {\mathrm {d} }{\mathrm {d} x}}\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}=\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}{\biggl \langle }\sum _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}^{-1}{\biggl \{}{\frac {\mathrm {d} }{\mathrm {d} x}}{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}{\biggr \}}{\biggr \rangle }=

=\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}{\biggl \{}\sum _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}^{-1}\,{\frac {x}{(n+x)^{2}}}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\biggr \}}=

Deswegen gilt in Kombination mit der Produktregel:

{\frac {\mathrm {d} }{\mathrm {d} x}}\,\Pi (x)={\frac {\mathrm {d} }{\mathrm {d} x}}{\biggl \{}\exp(-\,\gamma \,x)\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}{\biggr \}}=

={\bigl [}{\frac {\mathrm {d} }{\mathrm {d} x}}\exp(-\,\gamma \,x){\bigr ]}\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}+\exp(-\,\gamma \,x){\biggl \{}{\frac {\mathrm {d} }{\mathrm {d} x}}\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}{\biggr \}}=

=-\gamma \exp(-\,\gamma \,x)\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}+\exp(-\,\gamma \,x)\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}{\biggl \{}\sum _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}^{-1}\,{\frac {x}{(n+x)^{2}}}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\biggr \}}=

=\exp(-\,\gamma \,x)\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}{\biggl \langle }{\biggl \{}\sum _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}^{-1}\,{\frac {x}{(n+x)^{2}}}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\biggr \}}-\gamma {\biggr \rangle }=

=\exp(-\,\gamma \,x)\prod _{n=1}^{\infty }{\bigl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-1}\exp {\bigl (}{\frac {x}{n}}{\bigr )}{\bigr ]}{\biggl \{}{\biggl [}\sum _{n=1}^{\infty }{\bigl (}{\frac {1}{n}}-{\frac {1}{n+x}}{\bigr )}{\biggr ]}-\gamma {\biggr \}}=\Pi (x){\bigl [}\mathrm {H} (x)-\gamma {\bigr ]}

Denn die Harmonische Reihenfunktion hat diese Definition:

\mathrm {H} (x)=\sum _{n=1}^{\infty }{\bigl (}{\frac {1}{n}}-{\frac {1}{n+x}}{\bigr )}

Somit gilt analog für die Ableitung der beiden Gaußschen Pifunktionen:

{\frac {\mathrm {d} }{\mathrm {d} x}}\,\Pi (x)=\Pi (x){\bigl [}\mathrm {H} (x)-\gamma {\bigr ]}

{\frac {\mathrm {d} }{\mathrm {d} x}}\,\pi (x)=\pi (x){\bigl [}\gamma -\mathrm {H} (x){\bigr ]}

Beweisführung für einen Ableitungswert

Die Steigung der Gammafunktion an der Stelle 1 ist gleich dem Negativen der Euler-Mascheroni-Konstante $\gamma$ :

\Gamma '(1)=-\,\gamma

Diese Tatsache kann so mit der von Euler gebrauchten Formel bewiesen werden:

\Gamma '(1)=\Pi '(0)={\biggl [}{\frac {\mathrm {d} }{\mathrm {d} x}}\int _{0}^{\infty }y^{x}\exp(-y)\,\mathrm {d} y{\biggr ]}(x=0)={\biggl [}\int _{0}^{\infty }{\frac {\mathrm {d} }{\mathrm {d} x}}\,y^{x}\exp(-y)\,\mathrm {d} y{\biggr ]}(x=0)=

={\biggl [}\int _{0}^{\infty }\ln(y)\,y^{x}\exp(-y)\,\mathrm {d} y{\biggr ]}(x=0)=\int _{0}^{\infty }\ln(y)\exp(-y)\,\mathrm {d} y=-\,\gamma

Denn folgende Integralidentität gilt:

\int _{0}^{\infty }\ln(y)\exp(-y)\,\mathrm {d} y=-\int _{0}^{\infty }{\frac {1}{\exp(x)-1}}-{\frac {1}{x\exp(x)}}\,\mathrm {d} x

Die nun gezeigte Integralidentität basiert auf dieser Ableitung:

{\frac {\mathrm {d} }{\mathrm {d} x}}{\biggl \{}\ln[1-\exp(-x)]-{\frac {\ln(x)}{\exp(x)}}{\biggr \}}={\frac {1}{\exp(x)-1}}-{\frac {1}{x\exp(x)}}+{\frac {\ln(x)}{\exp(x)}}

Das drittoberste Integral in der Auflistung kann so bewiesen werden:

\int _{0}^{\infty }{\frac {1}{\exp(x)-1}}-{\frac {1}{x\exp(x)}}\,\mathrm {d} x=\int _{0}^{\infty }{\frac {\exp(-x)+x-1}{x[\exp(x)-1]}}\,\mathrm {d} x=\int _{0}^{\infty }{\frac {1}{x[\exp(x)-1]}}\sum _{n=1}^{\infty }{\biggl [}{\frac {x^{2n}}{(2n)!}}-{\frac {x^{2n+1}}{(2n+1)!}}{\biggr ]}\mathrm {d} x=

=\int _{0}^{\infty }\sum _{n=1}^{\infty }{\biggl \{}{\frac {1}{\exp(x)-1}}{\biggl [}{\frac {x^{2n-1}}{(2n)!}}-{\frac {x^{2n}}{(2n+1)!}}{\biggr ]}{\biggr \}}\mathrm {d} x=\sum _{n=1}^{\infty }\int _{0}^{\infty }{\biggl \{}{\frac {1}{\exp(x)-1}}{\biggl [}{\frac {x^{2n-1}}{(2n)!}}-{\frac {x^{2n}}{(2n+1)!}}{\biggr ]}{\biggr \}}\mathrm {d} x=

=\sum _{n=1}^{\infty }{\biggl [}{\frac {1}{(2n)!}}\int _{0}^{\infty }{\frac {x^{2n-1}}{\exp(x)-1}}\,\mathrm {d} x-{\frac {1}{(2n+1)!}}\int _{0}^{\infty }{\frac {x^{2n}}{\exp(x)-1}}\,\mathrm {d} x{\biggr ]}=\sum _{n=1}^{\infty }{\biggl [}{\frac {(2n-1)!\,\zeta (2n)}{(2n)!}}-{\frac {(2n)!\,\zeta (2n+1)}{(2n+1)!}}{\biggr ]}=

=\sum _{n=1}^{\infty }{\biggl [}{\frac {\zeta (2n)}{2n}}-{\frac {\zeta (2n+1)}{2n+1}}{\biggr ]}=\gamma

Für die Erklärung dieser Zetafunktionsformel siehe Artikel Euler-Mascheroni-Konstante!

Alternativ hierzu kann die Ableitung an der genannten Stelle auch direkt mit der Ableitungsformel für die Gaußsche Pifunktion ermittelt werden:

\Gamma '(1)=\Pi '(0)=\Pi (x){\bigl [}\mathrm {H} (x)-\gamma {\bigr ]}(x=0)=-\,\gamma

MacLaurinsche Reihe für den Gamma-Kehrwert

Die MacLaurinsche Reihe beziehungsweise Taylorsche Reihe für die Gammafunktion und ihren Kehrwert wurde insbesondere durch Leonhard Euler und Lorenzo Mascheroni erforscht. Der Kehrwert der Gammafunktion hat folgenden ersten Ableitungswert und folgenden zweiten Ableitungswert am Koordinatenursprung:

{\frac {\mathrm {d} }{\mathrm {d} x}}\,{\frac {1}{\Gamma (x)}}(x=0)=1

{\frac {\mathrm {d} ^{2}}{\mathrm {d} x^{2}}}\,{\frac {1}{\Gamma (x)}}(x=0)=2\,\gamma

{\frac {\mathrm {d} ^{3}}{\mathrm {d} x^{3}}}\,{\frac {1}{\Gamma (x)}}(x=0)=-{\frac {1}{2}}(\pi ^{2}-6\,\gamma ^{2})

Anders als der Ableitungswert der reziproken Gammafunktion an der Stelle $x=0$ nimmt der Ableitungswert an der Stelle $x=1$ einen nicht elementaren Wert an:

{\frac {\mathrm {d} }{\mathrm {d} x}}\,{\frac {1}{\Gamma (x)}}(x=1)=\gamma

Auch in den drei zuletzt genannten Formeln wird mit dem Kürzel $\gamma$ die Mascheronische Konstante repräsentiert. Der Graph vom Kehrwert der Gammafunktion nimmt im Intervall von $x=0$ bis $x=1$ einen sigmoiden Verlauf an. Die MacLaurinsche Reihe für den Gamma-Kehrwert wurde insbesondere durch Wrench in seinem Werk Concerning Two Series for the Gamma Function aus dem Jahr 1968 akkurat beschrieben. Ebenso erforschten die Mathematiker Bourguet (1883), Davis (1933), Isaacson und Salzer (beide 1943) diese Reihenentwicklung.

Integration

Stammfunktion der Gammafunktion

Das Integral der Gammafunktion selbst ist nicht als elementare Kombination von der Gammafunktion und anderen elementaren Funktionen darstellbar. Diese Tatsache wurde vom Mathematiker Otto Hölder gezeigt. Aber folgende Integraldarstellung existiert für die Stammfunktion der Gammafunktion:

\int _{1}^{x}\Gamma (w)\,\mathrm {d} w=\int _{0}^{\infty }{\frac {y^{x-1}-1}{\ln(y)\exp(y)}}\,\mathrm {d} y

Als Stammfunktion der von Euler verwendeten Integralformel für die Gammafunktion geht diese Formel hervor. Denn das durch den Punkt P(0|1) verlaufende Integral einer verallgemeinerten Exponentialfunktion bezüglich des Ausdrucks im Exponenten ergibt immer das Produkt dieser Exponentialfunktion dividiert durch den Logarithmus Naturalis von der betroffenen konstanten Basis. Beispielsweise gilt:

\int _{1}^{2}\Gamma (w)\,\mathrm {d} w=\int _{0}^{\infty }{\frac {y-1}{\ln(y)\exp(y)}}\,\mathrm {d} y=0{,}92274595068063060514388\ldots

Fransén-Robinson-Konstante

Das uneigentliche Integral von Null bis Unendlich beim Kehrwert der Gammafunktion nimmt den Wert der Fransén-Robinson-Konstante an:

\mathrm {F} =\int _{0}^{\infty }{\frac {1}{\Gamma (x)}}\,\mathrm {d} x

Diese Konstante hat folgende Integralidentität bezüglich der elementaren Funktionen und Werte:

\mathrm {F} =\mathrm {e} +\int _{0}^{\infty }{\frac {\exp(-x)}{\ln(x)^{2}+\pi ^{2}}}\,\mathrm {d} x

\mathrm {F} \approx 2{,}807770242028519365225\ldots

Mit dem Buchstaben e wird an dieser Stelle die Eulersche Zahl ausgedrückt.

Hyperfakultät und Superfakultät

Der Logarithmus Naturalis aus der Eulerschen Gammafunktion beziehungsweise Gaußschen Pifunktion wird mit der Hyperfakultät integriert:

\int _{0}^{x}\ln {\bigl [}\Gamma (y+1){\bigr ]}\,\mathrm {d} y=\ln {\bigl [}\operatorname {hf} (x){\bigr ]}-{\frac {x}{2}}{\bigl [}x+1-\ln(2\,\pi ){\bigr ]}=

=(x+1)\ln {\bigl [}\Gamma (x+1){\bigr ]}-\ln {\bigl [}\operatorname {sf} (x){\bigr ]}-{\frac {x}{2}}{\bigl [}x+1-\ln(2\,\pi ){\bigr ]}

Das Kürzel $\operatorname {hf}$ stellt die Hyperfakultät und das Kürzel $\operatorname {sf}$ stellt die Superfakultät dar.

Auf folgende Weise ist die Superfakultät für alle reellen Werte $x\in \mathbb {R}$ definiert:

\operatorname {sf} (x)=\Gamma (x+1)\exp {\biggl \{}{\frac {x}{2}}{\bigl [}\ln(2\,\pi )-\gamma \,x-x-1{\bigr ]}{\biggr \}}\prod _{n=1}^{\infty }{\biggl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{n}\exp {\bigl (}{\frac {x^{2}}{2n}}-x{\bigr )}{\biggr ]}

Und auf diese Weise kann die Hyperfakultät definiert werden:

\operatorname {hf} (x)=\exp {\biggl \{}{\frac {x}{2}}{\bigl [}x+1-\ln(2\,\pi )-\gamma \,x{\bigr ]}{\biggr \}}\prod _{n=1}^{\infty }{\biggl [}{\bigl (}1+{\frac {x}{n}}{\bigr )}^{-n-x}\exp {\bigl (}{\frac {x^{2}}{2n}}+x{\bigr )}{\biggr ]}

Sukzessiv hierzu kann die Hyperfakultät auch so definiert werden:

\operatorname {hf} (x)=\operatorname {sf} (x)^{-1}\Pi (x)^{x+1}=\operatorname {sf} (x)^{-1}\Gamma (x+1)^{x+1}

Für Hyperfakultät und Superfakultät gelten diese Induktionsformeln, die zur sukzessiven Ermittlung der Werte dieser Funktionen für natürlichzahlige Abszissenwerte dienen:

\operatorname {hf} (0)=1

\operatorname {hf} (x)=x^{x}\operatorname {hf} (x-1)

\operatorname {sf} (0)=1

\operatorname {sf} (x)=\Pi (x)\operatorname {sf} (x-1)

Für alle natürlichen Zahlen $n\in \mathbb {N}$ gelten somit diese Formeln:

\operatorname {hf} (n)=\prod _{m=1}^{n}m^{m}

\operatorname {sf} (n)=\prod _{m=1}^{n}m!

Für Hyperfakultät und Superfakultät werden im nun Folgenden die ersten Zahlen aufgezählt:

$\operatorname {hf} (0)=1$	$\operatorname {hf} (1)=1$	$\operatorname {hf} (2)=4$	$\operatorname {hf} (3)=108$	$\operatorname {hf} (4)=27648$	$\operatorname {hf} (5)=86400000$

$\operatorname {sf} (0)=1$	$\operatorname {sf} (1)=1$	$\operatorname {sf} (2)=2$	$\operatorname {sf} (3)=12$	$\operatorname {sf} (4)=288$	$\operatorname {sf} (5)=34560$

Somit gelten diese Rechenbeispiele für die Integration vom Logarithmus Naturalis der Gammafunktion:

$\int _{0}^{1}\ln {\bigl [}\Gamma (x+1){\bigr ]}\,\mathrm {d} x=-1+{\frac {1}{2}}\ln(2\,\pi )$ $\int _{0}^{2}\ln {\bigl [}\Gamma (x+1){\bigr ]}\,\mathrm {d} x=\ln(4)-3+\ln(2\,\pi )$ $\int _{0}^{3}\ln {\bigl [}\Gamma (x+1){\bigr ]}\,\mathrm {d} x=\ln(108)-6+{\frac {3}{2}}\ln(2\,\pi )$ $\int _{0}^{4}\ln {\bigl [}\Gamma (x+1){\bigr ]}\,\mathrm {d} x=\ln(27648)-10+2\ln(2\,\pi )$ $\int _{0}^{5}\ln {\bigl [}\Gamma (x+1){\bigr ]}\,\mathrm {d} x=\ln(86400000)-15+{\frac {5}{2}}\ln(2\,\pi )$

Zusammenhang mit der Riemannschen ζ-Funktion

Bernhard Riemann brachte 1859 die Gammafunktion mit der Riemannschen ζ-Funktion über die Formel

\Gamma (s)\,\zeta (s)=\int _{0}^{\infty }{\frac {x^{s-1}}{\mathrm {e} ^{x}-1}}\,\mathrm {d} x

und die folgende Feststellung in Beziehung: Der Ausdruck $\Gamma (s/2)\,\pi ^{-s/2}\,\zeta (s)$ „bleibt ungeändert, wenn $s$ in $1-s$ verwandelt wird“, also

\Gamma (s/2)\,\pi ^{-s/2}\,\zeta (s)=\Gamma {\bigl (}(1-s)/2{\bigr )}\,\pi ^{-(1-s)/2}\,\zeta (1-s).

Näherungsweise Berechnung

Stirlingsche Formel

Näherungswerte der Gammafunktion für $x>0$ liefert unter anderem die Stirlingsche Formel, es gilt

\Gamma (x)={\sqrt {2\pi }}\,x^{x-1/2}\,\mathrm {e} ^{-x}\,\mathrm {e} ^{\mu (x)}

mit

0<\mu (x)<1/(12x).

Rekursive Näherung

Aus der Funktionalgleichung

\Gamma (z+1)=z\cdot \Gamma (z),

die eine Art Periodizität beinhaltet, können aus bekannten Funktionswerten in einem Streifen der Breite 1 in $\operatorname {Re} (z)$ die Werte in jedem anderen entsprechenden Streifen rekursiv berechnet werden. Mit

\log \Gamma (z)=\log \Gamma (z+1)-\log z

kann man von einem Streifen auf den benachbarten mit kleinerem Realteil gelangen, und das $m$ -fach. Da es für großes $|z|$ sehr gute Näherungen für $\log \Gamma (z)$ gibt, kann deren Genauigkeit in Bereiche übertragen werden, in denen direkte Anwendung der betreffenden Näherung nicht anzuraten wäre. Nach Rocktäschel empfiehlt sich, wie schon von Carl Friedrich Gauß bemerkt, die aus der Stirling-Formel abgeleitete asymptotische Entwicklung in $z$

\operatorname {Ro} (z):={\frac {1}{2}}\log(2\pi )+\left(z-{\frac {1}{2}}\right)\left(\log \left(z-{\frac {1}{2}}\right)-1\right)\sim \log \Gamma (z)

.

Diese hat zwar im Nahbereich bei $z={\tfrac {1}{2}}$ eine Irregularität, ist aber schon für $|z|>10$ brauchbar. Mit dem Korrekturterm $-{\tfrac {1}{24}}\left(z-{\tfrac {1}{2}}\right)^{-1}$ wird ihr Fehler auf die Größenordnung ${\mathcal {O}}(z^{-3})$ für unbeschränkt wachsendes $|z|$ verringert.

Die $m$ -fache Anwendung dieser Näherung führt auf

\log \Gamma (z)\approx \operatorname {Ro} (z+m)-\sum _{k=0}^{m-1}\log(z+k).

Den komplexen Logarithmus berechnet man über die Polardarstellung von $z$ . Für die meisten Anwendungen, etwa in der Wellenausbreitung, sollte $m=100$ ausreichen.

Unvollständige Gammafunktion

In der Literatur wird dieser Begriff, im Hinblick auf Integrationsgrenzen und Normierung (Regularisierung), nicht einheitlich verwendet.

Häufige Notationen sind:

\gamma (a,x)=\int _{0}^{x}t^{a-1}\mathrm {e} ^{-t}\,\mathrm {d} t

unvollständige Gammafunktion der oberen Grenze

\Gamma (a,x)=\int _{x}^{\infty }t^{a-1}\mathrm {e} ^{-t}\,\mathrm {d} t

unvollständige Gammafunktion der unteren Grenze

\operatorname {P} (a,x)={\frac {\gamma (a,x)}{\Gamma (a)}}

regularisierte (unvollständige) Gammafunktion der oberen Grenze

\operatorname {Q} (a,x)={\frac {\Gamma (a,x)}{\Gamma (a)}}

regularisierte (unvollständige) Gammafunktion der unteren Grenze

Spricht man von einer regularisierten Gammafunktion, so impliziert dies schon, dass sie unvollständig ist.

\Gamma (a,x,y)=\int _{x}^{y}t^{a-1}\mathrm {e} ^{-t}\,\mathrm {d} t

oder

\Gamma (a,x,y)={\frac {1}{\Gamma (a)}}\int _{x}^{y}t^{a-1}\mathrm {e} ^{-t}\,\mathrm {d} t

steht für die verallgemeinerte unvollständige Gammafunktion.

Verallgemeinerung

Eine Verallgemeinerung ist die multivariate Gammafunktion, die in der Wishart-Verteilung anzutreffen ist.

Siehe auch

Literatur

Niels Nielsen: Handbuch der Theorie der Gammafunktion. B. G. Teubner, Leipzig 1906 (im Internetarchiv, dito, dito).
E. T. Whittaker, G. N. Watson: The Gamma function. Kapitel 12 in A course of modern analysis. Cambridge University Press, 4. Ausgabe 1927; Neuauflage 1996, ISBN 0-521-58807-3, S. 235–264 (englisch; im Internetarchiv).
Emil Artin: Einführung in die Theorie der Gammafunktion. B. G. Teubner, Leipzig 1931; The Gamma function. Holt, Rinehart and Winston, New York 1964 (englische Übersetzung von Michael Butler).
Friedrich Lösch, Fritz Schoblik: Die Fakultät (Gammafunktion) und verwandte Funktionen. Mit besonderer Berücksichtigung ihrer Anwendungen. B. G. Teubner, Leipzig 1951.
Philip J. Davis: Leonhard Euler’s integral: A historical profile of the gamma function. The American Mathematical Monthly 66, 1959, S. 849–869 (englisch; 1963 mit dem Chauvenet-Preis ausgezeichnet; bei MathDL).
Konrad Königsberger: Die Gammafunktion. Kapitel 17 in Analysis 1. Springer, Berlin 1990; 6. Auflage 2003, ISBN 3-540-40371-X, S. 351–360.
Reinhold Remmert: Die Gammafunktion. Kapitel 2 in Funktionentheorie 2. Springer, Berlin 1991; mit Georg Schumacher: 3. Auflage 2007, ISBN 978-3-540-40432-3, S. 31–73.
Eberhard Freitag, Rolf Busam: Die Gammafunktion. Kapitel 4.1 in Funktionentheorie 1. Springer, Berlin 1993; 4. Auflage 2006, ISBN 3-540-31764-3, S. 194–212.
Jörg Arndt: Matters Computational, Ideas, Algorithms, Source Code. Springer-Verlag Berlin Heidelberg 2011, ISBN 978-3-642-14763-0, Seite 610
Hermann Hankel: Die Eulerschen Integrale bei unbeschränkter Variabilität des Arguments. Z. Math. Phys., 9 (1864), Seiten 1–21.
Edmund T. Whittaker, George Neville Watson: A Course in Modern Analysis, 4th ed. Cambridge, England: Cambridge University Press, 1990.
Matthias Hirschmanner, Doktor Stefan Krause: Die Gammafunktion. Institut für Analysis und Scientific Computing der Technischen Universität Wien, 2014.
Bourguet, L. "Sur les intégrales Eulériennes et quelques autres fonctions uniformes." Acta Math. 2, 261-295, 1883.
Davis, H. T. Tables of the Higher Mathematical Functions. Bloomington, IN: Principia Press, 1933.
Isaacson, E. and Salzer, H. E. "Mathematical Tables--Errata: 19. J. P. L. Bourget, 'Sur les intégrales Eulériennes et quelques autres fonctions uniformes,' Acta Mathematica, v. 2, 1883, pp. 261-295.' " Math. Tab. Aids Comput. 1, 124, 1943.
Wrench, J. W. Jr. "Concerning Two Series for the Gamma Function." Math. Comput. 22, 617-626, 1968.

Weblinks

Eric W. Weisstein: Gamma Function. In: MathWorld (englisch).
Aaron Krowne u. a.: Gamma Function. In: PlanetMath. (englisch)
Gamma Function. Bei The Wolfram Functions Site (englisch; mit Berechnungsmöglichkeiten).
Videos über die Gamma-Funktion für Anfänger.
https://mathworld.wolfram.com/EllipticIntegralSingularValue.html

Einzelnachweise

↑ Brief (JPG-Datei, 136 kB) von Daniel Bernoulli an Christian Goldbach vom 6. Oktober 1729, abgedruckt in Paul Heinrich Fuss (Hrsg.): Correspondance mathématique et physique de quelques célèbres géomètres du XVIIIème siècle. (Band 2), St.-Pétersbourg 1843, S. 324–325 (französisch).
↑ Peter Luschny: Interpolating the natural factorial n! or The birth of the real factorial function (1729–1826). (Englisch).
1 2 Brief (PDF-Datei, 118 kB) von Leonhard Euler an Christian Goldbach vom 13. Oktober 1729, abgedruckt in Paul Heinrich Fuss (Hrsg.): Correspondance mathématique et physique de quelques célèbres géomètres du XVIIIème siècle. (Band 1), St.-Pétersbourg 1843, S. 3–7 (lateinisch).
1 2 Carl Friedrich Gauß: Disquisitiones generales circa seriem infinitam 1+… Pars I. (30. Januar 1812), Commentationes Societatis Regiae Scientiarum Gottingensis recentiores 2 (classis mathematicae), 1813, S. 26 (lateinisch; auch in Gauß: Werke. Band 3. S. 145).
↑ Brief (PDF-Datei, 211 kB) von Leonhard Euler an Christian Goldbach vom 8. Januar 1730, abgedruckt in Paul Heinrich Fuss (Hrsg.): Correspondance mathématique et physique de quelques célèbres géomètres du XVIIIème siècle. Band 1, St.-Pétersbourg 1843, S. 11–18 (lateinisch).
1 2 Leonhard Euler: De progressionibus transcendentibus, seu quarum termini generales algebraice dari nequeunt. (28. November 1729), Commentarii academiae scientiarum imperialis Petropolitanae 5, 1738, S. 36–57 (lateinisch).
↑ Leonhard Euler: De evolutione integralium per producta infinita. (PDF-Datei, 1,2 MB), Kapitel 9 in Teil 1 des ersten Bandes von Euler: Institutionum calculi integralis. 1768, S. 225–250 (lateinisch).
↑ Adrien-Marie Legendre: Recherches sur diverses sortes d’intégrales définies. (13. November 1809), Mémoires de la classe des sciences mathématiques et physiques de l’Institut de France 10, 1809, S. 477 (französisch).
↑ Adrien-Marie Legendre: Traité des fonctions elliptiques et des intégrales Eulériennes. (Band 2), Huzard-Courcier, Paris 1826, S. 365 (französisch).
↑ O. Hölder: Ueber die Eigenschaft der Gammafunction keiner algebraischen Differentialgleichung zu genügen. 26. Juni 1886, Mathematische Annalen 28, 1887, S. 1–13.
↑ Steven B. Bank, Robert P. Kaufman: A note on Hölder’s theorem concerning the Gamma function. Mathematische Annalen 232, 1978, S. 115–120 (englisch).
↑ Karl Weierstraß: Über die Theorie der analytischen Facultäten. (20. Mai 1854), Journal für die reine und angewandte Mathematik 51, 1856, S. 36.
↑ Nielsen: Handbuch der Theorie der Gammafunktion. 1906, S. 3.
↑ Davis: Leonhard Euler’s integral: A historical profile of the gamma function. 1959, S. 867.
↑ Artin: Einführung in die Theorie der Gammafunktion. 1931, S. 31–35.
↑ Harald Bohr, Johannes Mollerup: Lærebog i matematisk Analyse III. (Lehrbuch der mathematischen Analysis III), Jul. Gjellerups Forlag, København (Kopenhagen) 1922 (dänisch).
↑ Artin: Einführung in die Theorie der Gammafunktion. 1931, S. 12–13.
↑ N. Bourbaki: Éléments de mathématique IV. Fonctions d’une variable réelle. Hermann, Paris 1951 (französisch).
↑ Konrad Knopp: Funktionentheorie II. (5. Auflage), de Gruyter, Berlin 1941, S. 47–49.
↑ Reinhold Remmert: Wielandt’s theorem about the Γ-function. The American Mathematical Monthly 103, 1996, S. 214–220 (englisch).
↑ Brief von Carl Friedrich Gauß an Friedrich Wilhelm Bessel vom 21. November 1811, abgedruckt in Arthur Auwers (Hrsg.): Briefwechsel zwischen Gauss und Bessel, Wilhelm Engelmann, Leipzig 1880, S. 151–155 (Auszug in Gauß: Werke. Band 10.1. S. 362–365).
↑ O. Schlömilch: Einiges über die Eulerischen Integrale der zweiten Art. Archiv der Mathematik und Physik 4, 1844, S. 171.
↑ Remmert: Die Gammafunktion. Kapitel 2 in Funktionentheorie 2. 2007, S. 39.
↑ E. Freitag, R. Busam: Funktionentheorie 1. Springer-Verlag, ISBN 3-540-31764-3, Seite 225.
↑ Siehe Remmert: Funktionentheorie 2. Kapitel 2, S. 51.
↑ DLMF: §5.9 Integral Representations ‣ Properties ‣ Chapter 5 Gamma Function. Abgerufen am 3. Februar 2023.
↑ T. Schmelzer, L. Trefethen: Computing the Gamma Function Using Contour Integrals and Rational Approximations. In: SIAM J. Numer. Anal. 2007 (semanticscholar.org [abgerufen am 3. Februar 2023]).
↑ E. E. Kummer: Beitrag zur Theorie der Function $\textstyle \Gamma (x)=\int _{0}^{\infty }\!\mathrm {e} ^{-v}v^{x-1}dv$ . Journal für die reine und angewandte Mathematik 35, 1847, S. 4.
↑ C. J. Malmstén: De integralibus quibusdam definitis, seriebusque infinitis. (1. Mai 1846), Journal für die reine und angewandte Mathematik 38, 1849, S. 25 (lateinisch).
↑ Ia. V. Blagouchine: Rediscovery of Malmsten’s integrals, their evaluation by contour integration methods and some related results. 2014, Ramanujan J., 35(1), 21–110, doi:10.1007/s11139-013-9528-5. Erratum-Addendum doi:10.1007/s11139-015-9763-z
↑ L. Euler: Remarques sur un beau rapport entre les séries des puissances tant directes que réciproques. (1749), Histoire de l’Académie Royale des Sciences et Belles-Lettres 17 (1761), 1768, S. 96/97 (französisch).
↑ L. Euler: Evolutio formulae integralis $\textstyle \int \!x^{f-1}\mathrm {d} x(lx)^{\frac {m}{n}}$ integratione a valore x=0 ad x=1 extensa. 4. Juli 1771, Novi commentarii academiae scientiarum imperialis Petropolitanae 16, 1772, S. 121 (lateinisch).
↑ Adrien-Marie Legendre: Recherches sur diverses sortes d’intégrales définies. (13. November 1809), Mémoires de la classe des sciences mathématiques et physiques de l’Institut de France 10, 1809, S. 485 (französisch).
↑ Carl Friedrich Gauß: Disquisitiones generales circa seriem infinitam 1+… Pars I. 30. Januar 1812, Commentationes Societatis Regiae Scientiarum Gottingensis recentiores 2 (classis mathematicae), 1813, S. 30 (lateinisch; auch in Gauß: Werke. Band 3, S. 150).
↑ Steven R. Finch: Mathematical constants. Cambridge University Press, Cambridge 2003, ISBN 0-521-81805-2, S. 33 (englisch).
↑ Gregory V. Chudnovsky: Contributions to the theory of transcendental numbers. American Mathematical Society, 1984, ISBN 0-8218-1500-8, S. 8 (englisch).
↑ Eric W. Weisstein: Elliptic Integral Singular Value. Abgerufen am 3. Juni 2022 (englisch).
↑ calculus - Taylor Series of Gamma Function. Abgerufen am 27. April 2023 (englisch).
↑ Eric W. Weisstein: Gamma Function. Abgerufen am 27. April 2023 (englisch).
↑ integration - how to take integral of the Gamma Function. Abgerufen am 9. November 2022 (englisch).
↑ Bernhard Riemann: Über die Anzahl der Primzahlen unter einer gegebenen Größe. (19. Oktober 1859), Monatsberichte der Königlichen Preuß. Akademie der Wissenschaften zu Berlin, 1860, S. 671–680.
↑ Paul Eugen Böhmer: Differenzengleichungen und bestimmte Integrale. K. F. Koehler, Leipzig 1939, S. 108.
↑ Otto Rudolf Rocktäschel: Methoden zur Berechnung der Gammafunktion für komplexes Argument. Dissertation, Dresden 1922, S. 14.
↑ Karl Rawer: Wave Propagation in the Ionosphere. Kluwer Academic Publishers, Dordrecht 1993, ISBN 0-7923-0775-5 (englisch).

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[1] Brief (JPG-Datei, 136 kB) von Daniel Bernoulli an Christian Goldbach vom 6. Oktober 1729, abgedruckt in Paul Heinrich Fuss (Hrsg.): Correspondance mathématique et physique de quelques célèbres géomètres du XVIIIème siècle. (Band 2), St.-Pétersbourg 1843, S. 324–325 (französisch).

[2] Peter Luschny: Interpolating the natural factorial n! or The birth of the real factorial function (1729–1826). (Englisch).

[euler1729-10-3] 1 2 Brief (PDF-Datei, 118 kB) von Leonhard Euler an Christian Goldbach vom 13. Oktober 1729, abgedruckt in Paul Heinrich Fuss (Hrsg.): Correspondance mathématique et physique de quelques célèbres géomètres du XVIIIème siècle. (Band 1), St.-Pétersbourg 1843, S. 3–7 (lateinisch).

[gauss1812-4] 1 2 Carl Friedrich Gauß: Disquisitiones generales circa seriem infinitam 1+… Pars I. (30. Januar 1812), Commentationes Societatis Regiae Scientiarum Gottingensis recentiores 2 (classis mathematicae), 1813, S. 26 (lateinisch; auch in Gauß: Werke. Band 3. S. 145).

[5] Brief (PDF-Datei, 211 kB) von Leonhard Euler an Christian Goldbach vom 8. Januar 1730, abgedruckt in Paul Heinrich Fuss (Hrsg.): Correspondance mathématique et physique de quelques célèbres géomètres du XVIIIème siècle. Band 1, St.-Pétersbourg 1843, S. 11–18 (lateinisch).

[euler1729-11-6] 1 2 Leonhard Euler: De progressionibus transcendentibus, seu quarum termini generales algebraice dari nequeunt. (28. November 1729), Commentarii academiae scientiarum imperialis Petropolitanae 5, 1738, S. 36–57 (lateinisch).

[7] Leonhard Euler: De evolutione integralium per producta infinita. (PDF-Datei, 1,2 MB), Kapitel 9 in Teil 1 des ersten Bandes von Euler: Institutionum calculi integralis. 1768, S. 225–250 (lateinisch).

[8] Adrien-Marie Legendre: Recherches sur diverses sortes d’intégrales définies. (13. November 1809), Mémoires de la classe des sciences mathématiques et physiques de l’Institut de France 10, 1809, S. 477 (französisch).

[9] Adrien-Marie Legendre: Traité des fonctions elliptiques et des intégrales Eulériennes. (Band 2), Huzard-Courcier, Paris 1826, S. 365 (französisch).

[10] O. Hölder: Ueber die Eigenschaft der Gammafunction keiner algebraischen Differentialgleichung zu genügen. 26. Juni 1886, Mathematische Annalen 28, 1887, S. 1–13.

[11] Steven B. Bank, Robert P. Kaufman: A note on Hölder’s theorem concerning the Gamma function. Mathematische Annalen 232, 1978, S. 115–120 (englisch).

[12] Karl Weierstraß: Über die Theorie der analytischen Facultäten. (20. Mai 1854), Journal für die reine und angewandte Mathematik 51, 1856, S. 36.

[13] Nielsen: Handbuch der Theorie der Gammafunktion. 1906, S. 3.

[14] Davis: Leonhard Euler’s integral: A historical profile of the gamma function. 1959, S. 867.

[15] Artin: Einführung in die Theorie der Gammafunktion. 1931, S. 31–35.

[16] Harald Bohr, Johannes Mollerup: Lærebog i matematisk Analyse III. (Lehrbuch der mathematischen Analysis III), Jul. Gjellerups Forlag, København (Kopenhagen) 1922 (dänisch).

[17] Artin: Einführung in die Theorie der Gammafunktion. 1931, S. 12–13.

[18] N. Bourbaki: Éléments de mathématique IV. Fonctions d’une variable réelle. Hermann, Paris 1951 (französisch).

[19] Konrad Knopp: Funktionentheorie II. (5. Auflage), de Gruyter, Berlin 1941, S. 47–49.

[20] Reinhold Remmert: Wielandt’s theorem about the Γ-function. The American Mathematical Monthly 103, 1996, S. 214–220 (englisch).

[21] Brief von Carl Friedrich Gauß an Friedrich Wilhelm Bessel vom 21. November 1811, abgedruckt in Arthur Auwers (Hrsg.): Briefwechsel zwischen Gauss und Bessel, Wilhelm Engelmann, Leipzig 1880, S. 151–155 (Auszug in Gauß: Werke. Band 10.1. S. 362–365).

[22] O. Schlömilch: Einiges über die Eulerischen Integrale der zweiten Art. Archiv der Mathematik und Physik 4, 1844, S. 171.

[23] Remmert: Die Gammafunktion. Kapitel 2 in Funktionentheorie 2. 2007, S. 39.

[24] E. Freitag, R. Busam: Funktionentheorie 1. Springer-Verlag, ISBN 3-540-31764-3, Seite 225.

[25] Siehe Remmert: Funktionentheorie 2. Kapitel 2, S. 51.

[26] DLMF: §5.9 Integral Representations ‣ Properties ‣ Chapter 5 Gamma Function. Abgerufen am 3. Februar 2023.

[27] T. Schmelzer, L. Trefethen: Computing the Gamma Function Using Contour Integrals and Rational Approximations. In: SIAM J. Numer. Anal. 2007 (semanticscholar.org [abgerufen am 3. Februar 2023]).

[28] E. E. Kummer: Beitrag zur Theorie der Function $\textstyle \Gamma (x)=\int _{0}^{\infty }\!\mathrm {e} ^{-v}v^{x-1}dv$ . Journal für die reine und angewandte Mathematik 35, 1847, S. 4.

[29] C. J. Malmstén: De integralibus quibusdam definitis, seriebusque infinitis. (1. Mai 1846), Journal für die reine und angewandte Mathematik 38, 1849, S. 25 (lateinisch).

[30] Ia. V. Blagouchine: Rediscovery of Malmsten’s integrals, their evaluation by contour integration methods and some related results. 2014, Ramanujan J., 35(1), 21–110, doi:10.1007/s11139-013-9528-5. Erratum-Addendum doi:10.1007/s11139-015-9763-z

[31] L. Euler: Remarques sur un beau rapport entre les séries des puissances tant directes que réciproques. (1749), Histoire de l’Académie Royale des Sciences et Belles-Lettres 17 (1761), 1768, S. 96/97 (französisch).

[32] L. Euler: Evolutio formulae integralis $\textstyle \int \!x^{f-1}\mathrm {d} x(lx)^{\frac {m}{n}}$ integratione a valore x=0 ad x=1 extensa. 4. Juli 1771, Novi commentarii academiae scientiarum imperialis Petropolitanae 16, 1772, S. 121 (lateinisch).

[33] Adrien-Marie Legendre: Recherches sur diverses sortes d’intégrales définies. (13. November 1809), Mémoires de la classe des sciences mathématiques et physiques de l’Institut de France 10, 1809, S. 485 (französisch).

[34] Carl Friedrich Gauß: Disquisitiones generales circa seriem infinitam 1+… Pars I. 30. Januar 1812, Commentationes Societatis Regiae Scientiarum Gottingensis recentiores 2 (classis mathematicae), 1813, S. 30 (lateinisch; auch in Gauß: Werke. Band 3, S. 150).

[35] Steven R. Finch: Mathematical constants. Cambridge University Press, Cambridge 2003, ISBN 0-521-81805-2, S. 33 (englisch).

[36] Gregory V. Chudnovsky: Contributions to the theory of transcendental numbers. American Mathematical Society, 1984, ISBN 0-8218-1500-8, S. 8 (englisch).

[37] Eric W. Weisstein: Elliptic Integral Singular Value. Abgerufen am 3. Juni 2022 (englisch).

[38] calculus - Taylor Series of Gamma Function. Abgerufen am 27. April 2023 (englisch).

[39] Eric W. Weisstein: Gamma Function. Abgerufen am 27. April 2023 (englisch).

[40] integration - how to take integral of the Gamma Function. Abgerufen am 9. November 2022 (englisch).

[41] Bernhard Riemann: Über die Anzahl der Primzahlen unter einer gegebenen Größe. (19. Oktober 1859), Monatsberichte der Königlichen Preuß. Akademie der Wissenschaften zu Berlin, 1860, S. 671–680.

[42] Paul Eugen Böhmer: Differenzengleichungen und bestimmte Integrale. K. F. Koehler, Leipzig 1939, S. 108.

[43] Otto Rudolf Rocktäschel: Methoden zur Berechnung der Gammafunktion für komplexes Argument. Dissertation, Dresden 1922, S. 14.

[44] Karl Rawer: Wave Propagation in the Ionosphere. Kluwer Academic Publishers, Dordrecht 1993, ISBN 0-7923-0775-5 (englisch).