En mathématiques, la fonction gamma inverse est la fonction

Tracé de1/Γ(x) le long de l'axe réel
Fonction gamma inverse sur le plan complexe
z↦1Γ(z),

où Γ(z) désigne la fonction gamma. Comme la fonction gamma est méromorphe et non nulle partout dans le plan complexe, son inverse est une fonction entière. En tant que fonction entière, elle est d'ordre 1 (ce qui signifie que ln ln |1/Γ(z)| croit moins lentement que ln |z| mais est de type infini (ce qui signifie que ln |1/Γ(z)| croit plus vite que tout multiple de |z|, puisque sa croissance est approximativement proportionnel à |z| ln |z| dans le demi-plan gauche).

L'inverse est parfois utilisée comme point de départ pour le calcul numérique de la fonction gamma, et quelques bibliothèques logicielles la fournissent séparément de la fonction gamma standard.

Karl Weierstrass a appelé la fonction gamma inverse la « factorielle » et l'a utilisée dans son développement du théorème de factorisation de Weierstrass.

Développement en produit infini

modifier

Suite aux définitions du produit infini pour la fonction gamma, dues respectivement à Euler et à Weierstrass, on obtient le développement du produit infini suivant pour la fonction gamma réciproque :

1Γ(z)=z∏n=1∞1+zn(1+1n)z1Γ(z)=zeγz∏n=1∞(1+zn)e−zn

où γ = 0,577216... est la constante d'Euler-Mascheroni. Ces développements sont valables pour tous les nombres complexes z .

Développement en série de Taylor

modifier

Le développement en série de Taylor autour de 0 donne[1]:

1 Γ(z) =z+γ z2+(γ22−π212) z3+(γ36−γπ212+ζ(3)3 )z4+⋯ 

où γ est la constante d'Euler-Mascheroni. Pour n > 2, le coefficient an pour le terme zn peut être calculé récursivement comme[2],[3],[4]

an= a2 an−1+∑j=2n−1(−1)j+1 ζ(j) an−j n−1= γ an−1−ζ(2) an−2+ζ(3) an−3−⋯ n−1

où ζ est la fonction zêta de Riemann. Une représentation intégrale de ces coefficients a été récemment trouvée par Fekih-Ahmed (2014)[4] :

an=(−1)nπn!∫0∞e−t Im⁡[ (ln⁡(t)−iπ)n ] dt .

Pour les petites valeurs, on obtient les valeurs suivantes :

Fekih-Ahmed (2014) [4] donne également une approximation pour an :

an≈(−1)n (n−1)!  2 πn   Im⁡( z0(1/2−n) e−nz0 1+z0 ) ,

où z0=−1nexp⁡(W−1(−n)) , et W−1 est la branche négative de la fonction W de Lambert.

Le développement de Taylor autour de 1 a les mêmes coefficients (mais décalés), c'est-à-dire :

1Γ(1+z)=1zΓ(z)=1+γ z+(γ22−π212) z2+(γ36−γπ212+ζ(3)3 )z3+⋯ 

(l'inverse de la fonction pi de Gauss).

Développement asymptotique

modifier

Lorsque |z| tend vers l'infini à arg(z) constant, on a le développement suivant :

ln⁡(1/Γ(z))∼−zln⁡(z)+z+12ln⁡(z2π)−112z+1360z3−11260z5pour |arg⁡(z)|<π

Représentation intégrale

modifier

La fonction gamma inverse est un exemple classique d'application de la méthode du point col[5],[6]. On a par exemple :

1Γ(z)=ezz1−z2iπ∫−ππexp⁡[−z(1−θcot⁡θ+ln⁡θsin⁡θ)]dθ,

Représentation intégrale de contour

modifier

Une représentation intégrale due à Hermann Hankel est

1Γ(z)=i2π∮H(−t)−ze−tdt,

où H est le contour de Hankel, c'est-à-dire le chemin encerclant 0 dans le sens positif, commençant et revenant à l'infini positif par rapport à la branche coupée le long de l'axe réel positif. Selon Schmelzer et Trefethen, l'évaluation numérique de l'intégrale de Hankel est la base de certaines des meilleures méthodes de calcul de la fonction gamma[5].

Représentations intégrales aux entiers positifs

modifier

Pour les entiers positifs n≥1, il existe une intégrale pour la fonction factorielle inverse donnée par [7]

1n!=12π∫−ππe−inteeit dt.

De même, pour tout réel c>0 et z∈ℂ tel que ℜe(z)>0 on a l'intégrale suivante pour la fonction gamma inverse le long de l'axe réel sous la forme de[8] :

1Γ(z)=12π∫−∞∞(c+it)−zec+itdt,

où le cas particulier où z=n+1/2 fournit une relation correspondante pour la fonction factorielle double inverse, 1(2n−1)!!=π2n⋅Γ(n+12).

Intégrale selon l'axe réel

modifier

L'intégration de la fonction gamma inverse le long de l'axe réel positif donne la valeur

∫0+∞1Γ(x)dx≈2,80777024,

qui est connue sous le nom de constante de Fransén-Robinson.

On a également la formule suivante ( [9] chapitre 9, exercice 100)

∫0+∞axΓ(x)dx=aea+a∫0+∞e−axln2(x)+π2dx

Voir aussi

modifier

Références

modifier
(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Reciprocal gamma function » (voir la liste des auteurs).
  1. ↑ (en) Eric W. Weisstein, « Gamma function », sur MathWorld
  2. ↑ (en) J.W. Wrench, « Concerning two series for the gamma function », Mathematics of Computation, vol. 22, no 103,‎ , p. 617–626 (DOI 10.1090/S0025-5718-1968-0237078-4, S2CID 121472614)
  3. ↑ (en) J.W. Wrench, « Erratum: Concerning two series for the gamma function », Mathematics of Computation, vol. 27, no 123,‎ , p. 681–682 (DOI 10.1090/S0025-5718-1973-0319344-9)
  4. 1 2 3 L. Fekih-Ahmed, « On the power series expansion of the reciprocal gamma function », HAL archives,
  5. 1 2 (en) Thomas Schmelzer et Lloyd N. Trefethen, « Computing the Gamma function using contour integrals and rational approximations », SIAM Journal on Numerical Analysis, Society for Industrial and Applied Mathematics, vol. 45, no 2,‎ , p. 558–571 (DOI 10.1137/050646342, lire en ligne Accès payant)
  6. ↑ (en) Nico M. Temme, Special Functions : An Introduction to the Classical Functions of Mathematical Physics, John Wiley & Sons, , 392 p.
  7. ↑ Graham, Knuth, and Patashnik, Concrete Mathematics, Addison-Wesley, , p. 566
  8. ↑ Schmidt, « A Short Note on Integral Transformations and Conversion Formulas for Sequence Generating Functions », Axioms, vol. 8, no 2,‎ , p. 62 (DOI 10.3390/axioms8020062, arXiv 1809.03933)
  9. ↑ (en) Henri Cohen, Number Theory Volume II: Analytic and Modern Tools, vol. 240, coll. « Graduate Texts in Mathematics », (ISBN 978-0-387-49893-5, ISSN 0072-5285, DOI 10.1007/978-0-387-49894-2, lire en ligne)