Lanczos-Approximation für die Gamma-Funktion

Online Rechner und Formeln zur Berechnung der Lanczos-Approximation


Argument
Dezimalstellen
Resultat

Zur Berechnung tragen Sie das Argument x ein. Dann klicken Sie den Button berechnen.

gamma Funktion

Die Lanczos-Approximation ist eine praxisnahe Methode zur numerischen Berechnung der Gamma-Funktion \(\Gamma(z)\) mit konstanter Genauigkeit. Entwickelt von Cornelius Lanczos im Jahr 1964, übertrifft sie für moderate Argumente die Stirling-Formel in Genauigkeit und wird in vielen wissenschaftlichen Bibliotheken eingesetzt (z. B. GNU Scientific Library, Boost, CPython).

Wenn x null oder eine negative ganze Zahl ist, ist das Ergebnis ±∞ und es wird eine Fehlermeldung angezeigt.


Formeln


Für \(\Re(z)>-\bigl(g+\tfrac12\bigr)\) lautet eine Variante der Lanczos-Formel: \[ \Gamma(z+1) = \sqrt{2\pi}\,\bigl(z + g + \tfrac12\bigr)^{\,z+\tfrac12}\, \exp\bigl(-\bigl(z+g+\tfrac12\bigr)\bigr)\;A_g(z), \] wobei \[ A_g(z) = \tfrac12\,p_0(g) + \sum_{k=1}^{N} p_k(g)\,\frac{z}{\,z+k\,} \] schnell konvergiert, sobald die Parameter \(g\) und die Koeffizienten \(p_k(g)\) passend gewählt sind.

Partialbruch-Darstellung

Durch Vorberechnen geeigneter Koeffizienten \(\{c_k\}\) lässt sich \(A_g(z)\) in die Form \[ A_g(z) = c_0 \;+\;\sum_{k=1}^N \frac{c_k}{\,z + k\,} \] umschreiben. Diese Darstellungsweise minimiert Rundungsfehler und erlaubt die Auswertung von \(\Gamma(z)\) mit nur wenigen Multiplikationen, Potenzen, Exponentialfunktionen und einem kleinen Lookup-Tabellenzugriff.