Aller au contenu
Wikipédia l'encyclopédie libre

Fonction gamma

Un article de Wikipédia, l'encyclopédie libre.

Pour les articles homonymes, voir gamma (homonymie).

Tracé de la fonction gamma le long de l'axe des réels.

En mathématiques, la fonction gamma (notée par Γ, la lettre majuscule gamma de l'alphabet grec) est une fonction utilisée communément, définie sur l'ensemble des nombres complexes excepté les entiers négatifs ou nuls. Elle est un exemple typique de fonction étudiée par l'analyse complexe.

On a pour tout entier n {\displaystyle n} {\displaystyle n} strictement positif, Γ ( n ) = ( n − 1 ) ! {\displaystyle \Gamma (n)=(n-1)!} {\displaystyle \Gamma (n)=(n-1)!} (où ( n − 1 ) ! {\displaystyle (n-1)!} {\displaystyle (n-1)!} est la factorielle de n − 1 {\displaystyle n-1} {\displaystyle n-1}). Ainsi Γ prolonge-t-elle la fonction factorielle à son domaine (à un décalage de 1 près).

Cette fonction est l'une des plus importantes fonctions spéciales ; ce qui signifie qu'elle joue un rôle pilier en analyse, et qu'elle n'est pas élémentaire (informellement, on ne peut pas l'écrire de manière élémentaire à l'aide des fonctions élémentaires vues au lycée).

Définition

[modifier | modifier le code ]

Définition de base

[modifier | modifier le code ]

La fonction gamma, notée Γ, est définie pour tout z ∈ C {\displaystyle z\in \mathbb {C} } {\displaystyle z\in \mathbb {C} } tel que R e ( z ) > 0 {\displaystyle \mathrm {Re} (z)>0} {\displaystyle \mathrm {Re} (z)>0}, et elle est à valeurs dans C {\displaystyle \mathbb {C} } {\displaystyle \mathbb {C} }. Elle est définie par :

Γ ( z ) = ∫ 0 + ∞ t z − 1 e − t d t {\displaystyle \Gamma (z)=\int _{0}^{+\infty }t^{z-1},円\mathrm {e} ^{-t},円\mathrm {d} t} {\displaystyle \Gamma (z)=\int _{0}^{+\infty }t^{z-1},円\mathrm {e} ^{-t},円\mathrm {d} t}.

C'est une intégrale paramétrée par z {\displaystyle z} {\displaystyle z}, l'intégration se faisant sur t {\displaystyle t} {\displaystyle t}. Cette intégrale impropre converge absolument sur le demi-plan complexe où la partie réelle est strictement positive[1] , et une intégration par parties [1] montre que sur son domaine Γ vérifie l'équation fonctionnelle :

Γ ( z + 1 ) = z Γ ( z ) {\displaystyle \Gamma (z+1)=z\;\Gamma (z)} {\displaystyle \Gamma (z+1)=z\;\Gamma (z)}.

La fonction ainsi définie est également connue sous le nom d'intégrale d'Euler de seconde espèce .

Définition prolongée

[modifier | modifier le code ]
Tracé du module de la fonction gamma (prolongée) autour de l'origine du plan complexe. Les cinq « cheminées roses » correspondent à des pôles.

La fonction Γ peut être prolongée analytiquement en une fonction méromorphe sur C ∖ { 0 , − 1 , − 2 , … } {\displaystyle \mathbb {C} \setminus \{0,-1,-2,\ldots \}} {\displaystyle \mathbb {C} \setminus \{0,-1,-2,\ldots \}}, admettant comme pôles 0, −1, −2,...

C'est cette définition de la fonction gamma dont nous nous servons dans cet article ; convention partagée par de nombreux auteurs. Nous notons toujours ce prolongement Γ.

L'équation fonctionnelle est toujours valide

[modifier | modifier le code ]

Un résultat essentiel est que l'équation fonctionnelle vue au-dessus est encore vraie pour le prolongement, sur tout son domaine.

Démonstration

La démonstration est malgré la longueur très triviale, à condition de connaître le principe des zéros isolés...

On pose E = C ∖ { 0 , − 1 , − 2 , … } {\displaystyle E=\mathbb {C} \setminus \{0,-1,-2,\ldots \}} {\displaystyle E=\mathbb {C} \setminus \{0,-1,-2,\ldots \}}, et f : E → C , z ↦ Γ ( z + 1 ) − z Γ ( z ) {\displaystyle f:E\to \mathbb {C} ,\;z\mapsto \Gamma (z+1)-z\Gamma (z)} {\displaystyle f:E\to \mathbb {C} ,\;z\mapsto \Gamma (z+1)-z\Gamma (z)}.

L'objectif est de montrer que f {\displaystyle f} {\displaystyle f} est nulle sur E {\displaystyle E} {\displaystyle E} (la conclusion sera triviale).

  • Déjà, on remarque E {\displaystyle E} {\displaystyle E} est ouvert.

Pour ce faire, on montre que F = { 0 , − 1 , − 2 , … } {\displaystyle F=\{0,-1,-2,\ldots \}} {\displaystyle F=\{0,-1,-2,\ldots \}} est fermé. Soit ( z n ) ∈ F N {\displaystyle (z_{n})\in F^{\mathbb {N} }} {\displaystyle (z_{n})\in F^{\mathbb {N} }} qui converge vers l {\displaystyle l} {\displaystyle l}. Alors on a classiquement z n + 1 − z n → 0 {\displaystyle z_{n+1}-z_{n}\to 0} {\displaystyle z_{n+1}-z_{n}\to 0}. On a donc z n + 1 − z n ≤ 1 2 {\displaystyle z_{n+1}-z_{n}\leq {\frac {1}{2}}} {\displaystyle z_{n+1}-z_{n}\leq {\frac {1}{2}}} à partir d'un certain rang N {\displaystyle N} {\displaystyle N}. Or, la distance entre deux éléments distincts de F {\displaystyle F} {\displaystyle F} est trivialement ≥ 1 {\displaystyle \geq 1} {\displaystyle \geq 1}. Il s'ensuit que l'on a nécessairement z n + 1 − z n = 0 {\displaystyle z_{n+1}-z_{n}=0} {\displaystyle z_{n+1}-z_{n}=0} pour n ≥ N {\displaystyle n\geq N} {\displaystyle n\geq N}. Et donc (avec une récurrence immédiate) z n = z N {\displaystyle z_{n}=z_{N}} {\displaystyle z_{n}=z_{N}} pour n ≥ N {\displaystyle n\geq N} {\displaystyle n\geq N}, puis en faisant tendre n {\displaystyle n} {\displaystyle n} vers + ∞ {\displaystyle +\infty } {\displaystyle +\infty }, l = z N {\displaystyle l=z_{N}} {\displaystyle l=z_{N}}, et donc l ∈ F {\displaystyle l\in F} {\displaystyle l\in F}.

  • On montre également que E {\displaystyle E} {\displaystyle E} est connexe.

R 2 {\displaystyle \mathbb {R} ^{2}} {\displaystyle \mathbb {R} ^{2}} et C {\displaystyle \mathbb {C} } {\displaystyle \mathbb {C} } sont classiquement homéomorphes (en tant que R {\displaystyle \mathbb {R} } {\displaystyle \mathbb {R} }-EVN, en considérant les topologies induites) via ϕ ( a , b ) = a + i b {\displaystyle \phi (a,b)=a+ib} {\displaystyle \phi (a,b)=a+ib}. A = R 2 ∖ { ( 0 ; 0 ) , ( − 1 ; 0 ) , ( − 2 ; 0 ) , … } {\displaystyle A=\mathbb {R} ^{2}\setminus \{(0;0),(-1;0),(-2;0),\ldots \}} {\displaystyle A=\mathbb {R} ^{2}\setminus \{(0;0),(-1;0),(-2;0),\ldots \}} est classiquement connexe par arcs (on a en effet privé R 2 {\displaystyle \mathbb {R} ^{2}} {\displaystyle \mathbb {R} ^{2}} d'une partie dénombrable). Dès lors, E = ϕ ( A ) {\displaystyle E=\phi (A)} {\displaystyle E=\phi (A)} est également connexe par arcs ; et donc connexe.

  • Ensuite, on peut voir que f {\displaystyle f} {\displaystyle f} est analytique. Cela découle trivialement du caractère analytique de Γ {\displaystyle \Gamma } {\displaystyle \Gamma } (par construction), puis par conservation du caractère par opérations élémentaires (composition triviale, produit trivial, différence).
  • Les hypothèses nécessaires pour appliquer le principe des zéros isolés sont donc réunies. Or, l'équation fonctionnelle invoquée dans la définition de base montre que f {\displaystyle f} {\displaystyle f} est nulle pour tout z ∈ C {\displaystyle z\in \mathbb {C} } {\displaystyle z\in \mathbb {C} } quand R e ( z ) > 0 {\displaystyle \mathrm {Re} (z)>0} {\displaystyle \mathrm {Re} (z)>0} ; f {\displaystyle f} {\displaystyle f} est donc a fortiori nulle sur la boule ouverte B ( 1 , 1 ) {\displaystyle \mathrm {B} (1,1)} {\displaystyle \mathrm {B} (1,1)}. f {\displaystyle f} {\displaystyle f} est donc nulle sur tout son domaine. CQFD.

Calcul de Γ prolongée

[modifier | modifier le code ]

Soit z ∈ C ∖ { 0 , − 1 , − 2 , … } {\displaystyle z\in \mathbb {C} \setminus \{0,-1,-2,\ldots \}} {\displaystyle z\in \mathbb {C} \setminus \{0,-1,-2,\ldots \}} tel que R e ( z ) ≤ 0 {\displaystyle \mathrm {Re} (z)\leq 0} {\displaystyle \mathrm {Re} (z)\leq 0}.

En itérant l'équation fonctionnelle, on a cette relation :

∀ k ∈ N , Γ ( z ) = Γ ( z + k ) z ( z + 1 ) … ( z + k − 1 ) {\displaystyle \forall k\in \mathbb {N} ,\quad \Gamma (z)={\frac {\Gamma (z+k)}{z(z+1)\ldots (z+k-1)}}} {\displaystyle \forall k\in \mathbb {N} ,\quad \Gamma (z)={\frac {\Gamma (z+k)}{z(z+1)\ldots (z+k-1)}}}.

Il suffit alors de prendre k {\displaystyle k} {\displaystyle k} suffisamment grand (i.e. tel que R e ( z + k ) > 0 {\displaystyle \mathrm {Re} (z+k)>0} {\displaystyle \mathrm {Re} (z+k)>0}) pour se ramener à la définition de base.

Par exemple, en usant de la valeur de Γ ( 1 / 2 ) {\displaystyle \Gamma (1/2)} {\displaystyle \Gamma (1/2)} donnée plus bas, on a Γ ( − 3 / 2 ) = Γ ( 1 / 2 ) ( − 3 / 2 ) × ( − 1 / 2 ) = 4 3 Γ ( 1 / 2 ) = 4 π 3 {\displaystyle \Gamma (-3/2)={\frac {\Gamma (1/2)}{(-3/2)\times (-1/2)}}={\frac {4}{3}}\Gamma (1/2)={\frac {4{\sqrt {\pi }}}{3}}} {\displaystyle \Gamma (-3/2)={\frac {\Gamma (1/2)}{(-3/2)\times (-1/2)}}={\frac {4}{3}}\Gamma (1/2)={\frac {4{\sqrt {\pi }}}{3}}}.

Notation de Gauss

[modifier | modifier le code ]

Une notation alternative est celle introduite par Gauss

Π ( z ) = z Γ ( z ) = Γ ( z + 1 ) = ∫ 0 + ∞ t z e − t d t {\displaystyle \Pi (z)=z\Gamma (z)=\Gamma (z+1)=\int _{0}^{+\infty }t^{z}\mathrm {e} ^{-t}\mathrm {d} t} {\displaystyle \Pi (z)=z\Gamma (z)=\Gamma (z+1)=\int _{0}^{+\infty }t^{z}\mathrm {e} ^{-t}\mathrm {d} t}

qui simplifie certaines formules telles que

∀ n ∈ N , Π ( n ) = n ! {\displaystyle \forall n\in \mathbb {N} ,\Pi (n)=n!} {\displaystyle \forall n\in \mathbb {N} ,\Pi (n)=n!}

Legendre serait à l'origine de la popularisation de la notation Γ, préférant poser le premier pole en 0 plutôt qu'en -1[2]

Autres définitions

[modifier | modifier le code ]

Par changement de variable, l'intégrale précédente (pour Re(z) > 0) s'écrit aussi :

Γ ( z ) = 2 ∫ 0 + ∞ u 2 z − 1 e − u 2 d u et Γ ( z ) = ∫ 0 1 ( − ln ⁡ s ) z − 1 d s {\displaystyle \Gamma (z)=2\int _{0}^{+\infty }u^{2z-1}\mathrm {e} ^{-u^{2}},円\mathrm {d} u\quad {\text{et}}\quad \Gamma (z)=\int _{0}^{1}\left(-\ln s\right)^{z-1},円\mathrm {d} s} {\displaystyle \Gamma (z)=2\int _{0}^{+\infty }u^{2z-1}\mathrm {e} ^{-u^{2}},円\mathrm {d} u\quad {\text{et}}\quad \Gamma (z)=\int _{0}^{1}\left(-\ln s\right)^{z-1},円\mathrm {d} s}.

La définition suivante de la fonction gamma par produits infinis, due à Euler, a un sens pour les nombres complexes z qui ne sont pas des entiers négatifs ou nuls[3]  :

Γ ( z ) = lim n → + ∞ n ! n z z ( z + 1 ) ⋯ ( z + n ) = 1 z ∏ k = 1 + ∞ ( 1 + 1 k ) z 1 + z k {\displaystyle \Gamma (z)=\lim _{n\to {+\infty }}{\frac {n!\;n^{z}}{z\;(z+1)\cdots (z+n)}}={\frac {1}{z}},円\prod _{k=1}^{+\infty }{\frac {\left(1+{\frac {1}{k}}\right)^{z}}{1+{\frac {z}{k}}}}} {\displaystyle \Gamma (z)=\lim _{n\to {+\infty }}{\frac {n!\;n^{z}}{z\;(z+1)\cdots (z+n)}}={\frac {1}{z}},円\prod _{k=1}^{+\infty }{\frac {\left(1+{\frac {1}{k}}\right)^{z}}{1+{\frac {z}{k}}}}}.

Elle est équivalente à celle donnée par Schlömilch [4] ,[5] ,[6]  :

Γ ( z ) = e − γ z z ∏ k = 1 + ∞ e z k 1 + z k {\displaystyle \Gamma (z)={\frac {\operatorname {e} ^{-\gamma z}}{z}}\prod _{k=1}^{+\infty }{\frac {\operatorname {e} ^{\frac {z}{k}}}{1+{\frac {z}{k}}}}} {\displaystyle \Gamma (z)={\frac {\operatorname {e} ^{-\gamma z}}{z}}\prod _{k=1}^{+\infty }{\frac {\operatorname {e} ^{\frac {z}{k}}}{1+{\frac {z}{k}}}}}

où γ = ∑ k = 1 ∞ [ 1 k − ln ⁡ ( 1 + 1 k ) ] {\displaystyle \gamma =\sum _{k=1}^{\infty }\left[{\frac {1}{k}}-\ln \left(1+{\frac {1}{k}}\right)\right]} {\displaystyle \gamma =\sum _{k=1}^{\infty }\left[{\frac {1}{k}}-\ln \left(1+{\frac {1}{k}}\right)\right]} est la constante d'Euler-Mascheroni.

Propriétés

[modifier | modifier le code ]

Lien avec la factorielle

[modifier | modifier le code ]

De Γ(z+1) = zΓ(z) et Γ(1) = 1, il se déduit par récurrence que :

∀ n ∈ N , Γ ( n + 1 ) = n ! {\displaystyle \forall ,円n\in \mathbb {N} ,\;\Gamma (n+1)=n!} {\displaystyle \forall ,円n\in \mathbb {N} ,\;\Gamma (n+1)=n!}

Ce qui se réécrit de manière plus commode : ∀ n ∈ N , Π ( n ) = n ! {\displaystyle \forall ,円n\in \mathbb {N} ,\;\Pi (n)=n!} {\displaystyle \forall ,円n\in \mathbb {N} ,\;\Pi (n)=n!}

On interprète donc la fonction gamma comme un prolongement de la factorielle à l'ensemble des nombres complexes (à l'exception des entiers négatifs ou nuls).

Pas de racine

[modifier | modifier le code ]

La fonction gamma ne s'annule jamais. Ce résultat fondamental n'est pas trivial, à moins de connaitre la formule des compléments.

Caractérisations

[modifier | modifier le code ]

Sur l'ensemble des réels > 0

[modifier | modifier le code ]

La fonction gamma (restreinte à R + ∗ {\displaystyle \mathbb {R} _{+}^{*}} {\displaystyle \mathbb {R} _{+}^{*}}) est la seule fonction de R + ∗ {\displaystyle \mathbb {R} _{+}^{*}} {\displaystyle \mathbb {R} _{+}^{*}} dans R + ∗ {\displaystyle \mathbb {R} _{+}^{*}} {\displaystyle \mathbb {R} _{+}^{*}} vérifiant les trois propriétés suivantes (théorème de Bohr-Mollerup) :

  • Γ ( 1 ) = 1 {\displaystyle \Gamma (1)=1,円} {\displaystyle \Gamma (1)=1,円};
  • Γ ( x + 1 ) = x Γ ( x ) {\displaystyle \Gamma (x+1)=x\;\Gamma (x),円} {\displaystyle \Gamma (x+1)=x\;\Gamma (x),円} (sur tout le domaine) ;
  • Γ {\displaystyle \Gamma } {\displaystyle \Gamma } est logarithmiquement convexe (i.e. ln ∘ Γ {\displaystyle \ln \circ ,円\Gamma } {\displaystyle \ln \circ ,円\Gamma } est convexe).

Sur le demi-plan complexe Re(z)>0

[modifier | modifier le code ]

La fonction gamma (restreinte à { z ∈ C ∣ R e ( z ) > 0 } {\displaystyle \{z\in \mathbb {C} \mid \mathrm {Re} (z)>0\}} {\displaystyle \{z\in \mathbb {C} \mid \mathrm {Re} (z)>0\}}) est la seule fonction de { z ∈ C ∣ R e ( z ) > 0 } {\displaystyle \{z\in \mathbb {C} \mid \mathrm {Re} (z)>0\}} {\displaystyle \{z\in \mathbb {C} \mid \mathrm {Re} (z)>0\}} dans C {\displaystyle \mathbb {C} } {\displaystyle \mathbb {C} } vérifiant les quatre propriétés suivantes (théorème de Wielandt) :

  • Γ ( 1 ) = 1 {\displaystyle \Gamma (1)=1,円} {\displaystyle \Gamma (1)=1,円};
  • Γ ( z + 1 ) = z Γ ( z ) {\displaystyle \Gamma (z+1)=z\;\Gamma (z),円} {\displaystyle \Gamma (z+1)=z\;\Gamma (z),円} (sur tout le domaine) ;
  • Γ {\displaystyle \Gamma } {\displaystyle \Gamma } est holomorphe ;
  • Γ {\displaystyle \Gamma } {\displaystyle \Gamma } est bornée sur l'ensemble { z ∈ C ∣ 1 < R e ( z ) < 2 } {\displaystyle \{z\in \mathbb {C} \mid 1<\mathrm {Re} (z)<2\}} {\displaystyle \{z\in \mathbb {C} \mid 1<\mathrm {Re} (z)<2\}}.

Autres propriétés

[modifier | modifier le code ]

Formule des compléments

[modifier | modifier le code ]

La fonction gamma vérifie la formule de réflexion d'Euler, ou formule des compléments :

∀ z ∈ C ∖ Z Γ ( 1 − z ) Γ ( z ) = π sin ⁡ ( π z ) {\displaystyle \forall z\in \mathbb {C} \setminus \mathbb {Z} \quad \Gamma (1-z)\;\Gamma (z)={\pi \over \sin(\pi z)}} {\displaystyle \forall z\in \mathbb {C} \setminus \mathbb {Z} \quad \Gamma (1-z)\;\Gamma (z)={\pi \over \sin(\pi z)}}.

On la démontre aisément à l'aide de la définition d'Euler en produit infini, en utilisant ensuite la factorisation classique sin ⁡ π z = π z ∏ n = 1 ∞ ( 1 − z 2 n 2 ) {\displaystyle \sin \pi z=\pi z\prod _{n=1}^{\infty }\left(1-{\frac {z^{2}}{n^{2}}}\right)} {\displaystyle \sin \pi z=\pi z\prod _{n=1}^{\infty }\left(1-{\frac {z^{2}}{n^{2}}}\right)}.

Formule de multiplication

[modifier | modifier le code ]

La fonction gamma vérifie également la formule de duplication :

Pour tout z ∈ C {\displaystyle z\in \mathbb {C} } {\displaystyle z\in \mathbb {C} } tel que 2 z ∉ { 0 , − 1 , − 2 , … } {\displaystyle 2z\notin \{0,-1,-2,\ldots \}} {\displaystyle 2z\notin \{0,-1,-2,\ldots \}},

Γ ( z ) Γ ( z + 1 2 ) = 2 1 − 2 z π Γ ( 2 z ) {\displaystyle \Gamma (z)\;\Gamma \left(z+{\frac {1}{2}}\right)=2^{1-2z}\;{\sqrt {\pi }}\;\Gamma (2z)} {\displaystyle \Gamma (z)\;\Gamma \left(z+{\frac {1}{2}}\right)=2^{1-2z}\;{\sqrt {\pi }}\;\Gamma (2z)}.


C'est un cas particulier du théorème de multiplication (à l'ordre m > 0 {\displaystyle m>0} {\displaystyle m>0}) :

Pour tout z ∈ C {\displaystyle z\in \mathbb {C} } {\displaystyle z\in \mathbb {C} } tel que m z ∉ { 0 , − 1 , − 2 , … } {\displaystyle mz\notin \{0,-1,-2,\ldots \}} {\displaystyle mz\notin \{0,-1,-2,\ldots \}},

Γ ( z ) Γ ( z + 1 m ) Γ ( z + 2 m ) ⋯ Γ ( z + m − 1 m ) = ( 2 π ) ( m − 1 ) / 2 m 1 / 2 − m z Γ ( m z ) {\displaystyle \Gamma (z)\;\Gamma \left(z+{\frac {1}{m}}\right)\;\Gamma \left(z+{\frac {2}{m}}\right)\cdots \Gamma \left(z+{\frac {m-1}{m}}\right)=(2\pi )^{(m-1)/2}\;m^{1/2-mz}\;\Gamma (mz)} {\displaystyle \Gamma (z)\;\Gamma \left(z+{\frac {1}{m}}\right)\;\Gamma \left(z+{\frac {2}{m}}\right)\cdots \Gamma \left(z+{\frac {m-1}{m}}\right)=(2\pi )^{(m-1)/2}\;m^{1/2-mz}\;\Gamma (mz)}.


Cette fonction apparaît également dans des formules incluant la fonction zêta de Riemann.

La fonction gamma possède un pôle d'ordre 1 en z = −n pour tout entier naturel n. Le résidu de la fonction en ce pôle est donné par :

Res ⁡ ( Γ , − n ) = ( − 1 ) n n ! {\displaystyle \operatorname {Res} (\Gamma ,-n)={\frac {(-1)^{n}}{n!}}} {\displaystyle \operatorname {Res} (\Gamma ,-n)={\frac {(-1)^{n}}{n!}}}.

Dérivées

[modifier | modifier le code ]

La fonction gamma est infiniment dérivable sur R + ∗ {\displaystyle \mathbb {R} _{+}^{*}} {\displaystyle \mathbb {R} _{+}^{*}}(c’est-à-dire p fois dérivable pour tout entier p). Sa dérivée est exprimée à l'aide de la fonction digamma : Γ ′ ( z ) = Γ ( z ) ψ 0 ( z ) . {\displaystyle \Gamma '(z)=\Gamma (z)\psi _{0}(z).,円} {\displaystyle \Gamma '(z)=\Gamma (z)\psi _{0}(z).,円}

Plus généralement, sa dérivée p-ième possède sur R + ∗ {\displaystyle \mathbb {R} _{+}^{*}} {\displaystyle \mathbb {R} _{+}^{*}} l'expression intégrale suivante :

Γ ( p ) ( x ) = ∫ 0 + ∞ ( ln ⁡ t ) p t x − 1 e − t d t {\displaystyle \Gamma ^{(p)}(x)=\int _{0}^{+\infty }{(\ln t)^{p},円t^{x-1},円\operatorname {e} ^{-t},円\mathrm {d} t}} {\displaystyle \Gamma ^{(p)}(x)=\int _{0}^{+\infty }{(\ln t)^{p},円t^{x-1},円\operatorname {e} ^{-t},円\mathrm {d} t}}.

Lien avec les sommes de Gauss

[modifier | modifier le code ]
Cette section est vide, insuffisamment détaillée ou incomplète. Votre aide est la bienvenue ! Comment faire ?

La définition de la fonction gamma sous forme d'intégrale la fait apparaître comme une convolution entre un caractère additif (l'exponentielle) et un caractère multiplicatif ( x ↦ x s {\displaystyle x\mapsto x^{s}} {\displaystyle x\mapsto x^{s}}).

Lien avec d'autres fonctions

[modifier | modifier le code ]

La fonction gamma est reliée à la fonction ζ de Riemann par :

ζ ( s ) Γ ( s ) = ∫ 0 + ∞ t s − 1 e t − 1 d t {\displaystyle \zeta (s),円\Gamma (s)=\int _{0}^{+\infty }{\frac {t^{s-1}}{\mathrm {e} ^{t}-1}},円\mathrm {d} t} {\displaystyle \zeta (s),円\Gamma (s)=\int _{0}^{+\infty }{\frac {t^{s-1}}{\mathrm {e} ^{t}-1}},円\mathrm {d} t}.

Elle est reliée à la fonction êta de Dirichlet par[7]  :

Γ ( s ) η ( s ) = ∫ 0 ∞ x s − 1 e x + 1 d x = ∫ 0 1 ∫ 0 1 ( − ln ⁡ ( x y ) ) s − 2 1 + x y d x d y {\displaystyle \Gamma (s),円\eta (s)=\int _{0}^{\infty }{\frac {x^{s-1}}{\mathrm {e} ^{x}+1}},円\mathrm {d} x=\int _{0}^{1}\int _{0}^{1}{\frac {(-\ln(xy))^{s-2}}{1+xy}},円\mathrm {d} x,円\mathrm {d} y} {\displaystyle \Gamma (s),円\eta (s)=\int _{0}^{\infty }{\frac {x^{s-1}}{\mathrm {e} ^{x}+1}},円\mathrm {d} x=\int _{0}^{1}\int _{0}^{1}{\frac {(-\ln(xy))^{s-2}}{1+xy}},円\mathrm {d} x,円\mathrm {d} y}.

Dans la définition de la fonction gamma sous forme d'intégrale, les bornes de l'intégrale sont fixées ; la fonction gamma incomplète est la fonction obtenue en modifiant la borne inférieure ou la borne supérieure.

La fonction gamma est reliée à la fonction bêta par la formule :

B ( x , y ) = Γ ( x ) Γ ( y ) Γ ( x + y ) . {\displaystyle \mathrm {B} (x,y)={\frac {\Gamma (x)\;\Gamma (y)}{\Gamma (x+y)}}.} {\displaystyle \mathrm {B} (x,y)={\frac {\Gamma (x)\;\Gamma (y)}{\Gamma (x+y)}}.}

Le logarithme de la fonction gamma est parfois appelé lngamma ou loggamma . Il intervient notamment dans la résolution des problèmes de propagation d’ondes [8]  : l'équation fonctionnelle de la fonction lngamma est :

ln ⁡ Γ ( z ) = ln ⁡ Γ ( z + 1 ) − ln ⁡ ( z ) {\displaystyle \ln \Gamma (z)=\ln \Gamma (z+1)-\ln(z)} {\displaystyle \ln \Gamma (z)=\ln \Gamma (z+1)-\ln(z)}.

Si l’on connaît les valeurs de la fonction sur une bande de largeur 1 en Re(z), on obtient par cette relation les valeurs dans une bande voisine de même largeur, et l’on peut répéter ce procédé. Partant d’un z avec Re(z) >> 1 pour lequel on connaît une bonne approximation, on peut ainsi atteindre la valeur pour un z quelconque.

Rocktaeschel (1922[9] , suivant une indication de Gauss) propose l'approximation pour Re(z) grand :

ln ⁡ Γ ( z ) ≈ ( z − 1 2 ) ln ⁡ z − z + 1 2 ln ⁡ ( 2 π ) {\displaystyle \ln \Gamma (z)\approx (z-{\tfrac {1}{2}})\ln z-z+{\tfrac {1}{2}}\ln(2\pi )} {\displaystyle \ln \Gamma (z)\approx (z-{\tfrac {1}{2}})\ln z-z+{\tfrac {1}{2}}\ln(2\pi )}.

On peut en déduire une approximation de ln Γ(z) pour Re(z) plus petit, en utilisant[10]  :

ln ⁡ Γ ( z − m ) = ln ⁡ Γ ( z ) − ∑ k = 1 m ln ⁡ ( z − k ) {\displaystyle \ln \Gamma (z-m)=\ln \Gamma (z)-\sum _{k=1}^{m}\ln(z-k)} {\displaystyle \ln \Gamma (z-m)=\ln \Gamma (z)-\sum _{k=1}^{m}\ln(z-k)}.

La dérivée logarithmique de la fonction gamma est appelée fonction digamma. Les dérivées d'ordre supérieur sont les fonctions polygamma.

Un analogue de la fonction gamma sur un corps fini ou un anneau fini est fourni par les sommes de Gauss.

En remarquant que la fonction gamma ne s'annule jamais, on peut définir son inverse (en). Or, on peut prolonger cet inverse à C {\displaystyle \mathbb {C} } {\displaystyle \mathbb {C} } tout entier (les inverses des pôles de Γ vaudront 0). On obtient alors une fonction entière (donc, informellement, se comportant un peu tel un polynôme).

Valeurs particulières

[modifier | modifier le code ]

Cette section indique quelques valeurs particulières de la fonction gamma (en) et de ses dérivées.

La valeur de Γ(1/2) = √π est celle de l'intégrale de Gauss ; elle peut aussi se déduire de la formule des compléments. Cette valeur permet, par récurrence, de déterminer les autres valeurs de la fonction gamma pour les demi-entiers positifs :

Γ ( 3 / 2 ) = π 2 , Γ ( 5 / 2 ) = 3 π 4 , … , {\displaystyle \Gamma (3/2)={\frac {\sqrt {\pi }}{2}},\quad \Gamma (5/2)={\frac {3{\sqrt {\pi }}}{4}},\ldots ,} {\displaystyle \Gamma (3/2)={\frac {\sqrt {\pi }}{2}},\quad \Gamma (5/2)={\frac {3{\sqrt {\pi }}}{4}},\ldots ,}
Γ ( n + 1 2 ) = ( n − 1 2 ) Γ ( n − 1 2 ) = ( n − 1 2 ) ( n − 3 2 ) ⋯ 3 2 1 2 Γ ( 1 2 ) = ( 2 n ) ! 2 2 n n ! π {\displaystyle \Gamma \left(n+{\frac {1}{2}}\right)=\left(n-{\frac {1}{2}}\right)\Gamma \left(n-{\frac {1}{2}}\right)=\left(n-{\frac {1}{2}}\right)\left(n-{\frac {3}{2}}\right)\cdots {\frac {3}{2}},円{\frac {1}{2}},円\Gamma \left({\frac {1}{2}}\right)={\frac {(2n)!}{2^{2n}n!}}{\sqrt {\pi }}} {\displaystyle \Gamma \left(n+{\frac {1}{2}}\right)=\left(n-{\frac {1}{2}}\right)\Gamma \left(n-{\frac {1}{2}}\right)=\left(n-{\frac {1}{2}}\right)\left(n-{\frac {3}{2}}\right)\cdots {\frac {3}{2}},円{\frac {1}{2}},円\Gamma \left({\frac {1}{2}}\right)={\frac {(2n)!}{2^{2n}n!}}{\sqrt {\pi }}}

mais aussi négatifs, par exemple :

Γ ( − 1 / 2 ) = − 2 π {\displaystyle \Gamma (-1/2)=-2{\sqrt {\pi }}} {\displaystyle \Gamma (-1/2)=-2{\sqrt {\pi }}}.

En ce qui concerne ses dérivées, avec γ la constante d'Euler-Mascheroni :

Γ ′ ( n + 1 ) = Γ ( n + 1 ) ψ 0 ( n + 1 ) = n ! ( − γ + ∑ 1 ≤ k ≤ n 1 k ) {\displaystyle \Gamma '(n+1)=\Gamma (n+1)\psi _{0}(n+1)=n!\left(-\gamma +\sum _{1\leq k\leq n}{\frac {1}{k}}\right)} {\displaystyle \Gamma '(n+1)=\Gamma (n+1)\psi _{0}(n+1)=n!\left(-\gamma +\sum _{1\leq k\leq n}{\frac {1}{k}}\right)} ;
Γ ′ ( n + 1 2 ) = Γ ( n + 1 2 ) ψ 0 ( n + 1 2 ) = ( 2 n ) ! 2 2 n n ! π ( − γ − 2 ln ⁡ 2 + ∑ 1 ≤ k ≤ n 2 2 k − 1 ) {\displaystyle \Gamma '\left(n+{\frac {1}{2}}\right)=\Gamma \left(n+{\frac {1}{2}}\right)\psi _{0}\left(n+{\frac {1}{2}}\right)={\frac {(2n)!}{2^{2n}n!}}{\sqrt {\pi }}\left(-\gamma -2\ln 2+\sum _{1\leq k\leq n}{\frac {2}{2k-1}}\right)} {\displaystyle \Gamma '\left(n+{\frac {1}{2}}\right)=\Gamma \left(n+{\frac {1}{2}}\right)\psi _{0}\left(n+{\frac {1}{2}}\right)={\frac {(2n)!}{2^{2n}n!}}{\sqrt {\pi }}\left(-\gamma -2\ln 2+\sum _{1\leq k\leq n}{\frac {2}{2k-1}}\right)} ;
Γ ″ ( 1 / 2 ) = π ( γ + 2 ln ⁡ ( 2 ) ) 2 + π 5 / 2 2 , Γ ″ ( 1 ) = γ 2 + π 2 6 , Γ ″ ( 2 ) = ( 1 − γ ) 2 + π 2 6 − 1 {\displaystyle \Gamma ''(1/2)={\sqrt {\pi }}(\gamma +2,円\ln(2))^{2}+{\frac {\pi ^{5/2}}{2}},\quad \Gamma ''(1)=\gamma ^{2}+{\frac {\pi ^{2}}{6}},\quad \Gamma ''(2)=(1-\gamma )^{2}+{\frac {\pi ^{2}}{6}}-1} {\displaystyle \Gamma ''(1/2)={\sqrt {\pi }}(\gamma +2,円\ln(2))^{2}+{\frac {\pi ^{5/2}}{2}},\quad \Gamma ''(1)=\gamma ^{2}+{\frac {\pi ^{2}}{6}},\quad \Gamma ''(2)=(1-\gamma )^{2}+{\frac {\pi ^{2}}{6}}-1}.

On connaît quelques résultats de transcendance et même d'indépendance algébrique sur les valeurs de Γ en certains points rationnels.

Une conjecture de Rohrlich [11] prédit que toute relation multiplicative de la forme

π b / 2 ∏ Γ ( a k ) m k ∈ Q ¯ , ( b ∈ Z , a k ∈ Q , m k ∈ Z ) {\displaystyle \pi ^{b/2}\prod \Gamma (a_{k})^{m_{k}}\in {\overline {\mathbb {Q} }},\quad (b\in \mathbb {Z} ,\;a_{k}\in \mathbb {Q} ,\;m_{k}\in \mathbb {Z} )} {\displaystyle \pi ^{b/2}\prod \Gamma (a_{k})^{m_{k}}\in {\overline {\mathbb {Q} }},\quad (b\in \mathbb {Z} ,\;a_{k}\in \mathbb {Q} ,\;m_{k}\in \mathbb {Z} )}

(où Q désigne le corps des nombres algébriques) se déduit des trois relations standard :

Γ ( z + 1 ) = z Γ ( z ) , Γ ( 1 − z ) Γ ( z ) = π sin ⁡ ( π z ) , ∏ 0 ≤ k < n Γ ( z + k n ) = ( 2 π ) ( n − 1 ) / 2 n − n z + 1 / 2 Γ ( n z ) {\displaystyle \Gamma (z+1)=z\Gamma (z),\quad \Gamma (1-z)\Gamma (z)={\pi \over \sin(\pi z)},\quad \prod _{0\leq k<n}\Gamma \left(z+{\frac {k}{n}}\right)=(2\pi )^{(n-1)/2}n^{-nz+1/2}\Gamma (nz)} {\displaystyle \Gamma (z+1)=z\Gamma (z),\quad \Gamma (1-z)\Gamma (z)={\pi \over \sin(\pi z)},\quad \prod _{0\leq k<n}\Gamma \left(z+{\frac {k}{n}}\right)=(2\pi )^{(n-1)/2}n^{-nz+1/2}\Gamma (nz)}.

Formule asymptotique de Stirling

[modifier | modifier le code ]

De Γ(z) et Γ(z+1)

[modifier | modifier le code ]

La formule de Stirling donne un équivalent au voisinage de l'infini de la factorielle :

n ! = 2 π n n + 1 2 e − n + μ ( n )  pour  n ∈ N   , {\displaystyle n,円!={\sqrt {2\pi }},円n^{n+{\frac {1}{2}}}{\rm {e}}^{-n+\mu (n)}{\text{ pour }}n\in \mathbb {N} \ ,} {\displaystyle n,円!={\sqrt {2\pi }},円n^{n+{\frac {1}{2}}}{\rm {e}}^{-n+\mu (n)}{\text{ pour }}n\in \mathbb {N} \ ,}

avec μ la fonction de Binet :

μ ( z ) = ∑ k = 1 ∞ B 2 k 2 k ( 2 k − 1 ) z 2 k − 1   , {\displaystyle \mu (z)=\sum _{k=1}^{\infty }{\frac {B_{2k}}{2k(2k-1)z^{2k-1}}}\ ,} {\displaystyle \mu (z)=\sum _{k=1}^{\infty }{\frac {B_{2k}}{2k(2k-1)z^{2k-1}}}\ ,}

et Bi les nombres de Bernoulli. Sachant que Γ(n+1) = n! sur N, cet équivalent se généralise à la fonction gamma :

Γ ( z + 1 ) = 2 π z z + 1 2 e − z + μ ( z )  pour  z ∈ C ∖ Z − ∗   , {\displaystyle \Gamma (z+1)={\sqrt {2\pi }},円z^{z+{\frac {1}{2}}}{\rm {e}}^{-z+\mu (z)}{\text{ pour }}z\in \mathbb {C} \setminus \mathbb {Z_{-}^{*}} \ ,} {\displaystyle \Gamma (z+1)={\sqrt {2\pi }},円z^{z+{\frac {1}{2}}}{\rm {e}}^{-z+\mu (z)}{\text{ pour }}z\in \mathbb {C} \setminus \mathbb {Z_{-}^{*}} \ ,}

d’où :

Γ ( z ) = Γ ( z + 1 ) z = 2 π z z − 1 2 e − z + μ ( z )  pour  z ∈ C ∖ Z −   . {\displaystyle \Gamma (z)={\frac {\Gamma (z+1)}{z}}={\sqrt {2\pi }},円z^{z-{\frac {1}{2}}}{\rm {e}}^{-z+\mu (z)}{\text{ pour }}z\in \mathbb {C} \setminus \mathbb {Z_{-}} \ .} {\displaystyle \Gamma (z)={\frac {\Gamma (z+1)}{z}}={\sqrt {2\pi }},円z^{z-{\frac {1}{2}}}{\rm {e}}^{-z+\mu (z)}{\text{ pour }}z\in \mathbb {C} \setminus \mathbb {Z_{-}} \ .}

En calculant les premiers termes de eμ grâce à la formule exponentielle, on obtient le développement asymptotique :

Γ ( z ) = 2 π z z − 1 2 e − z [ 1 + 1 12 z + 1 288 z 2 − 139 51840 z 3 − 571 2488320 z 4 + 163879 209018880 z 5 + O ( 1 z 6 ) ]   . {\displaystyle \Gamma (z)={\sqrt {2\pi }},円z^{z-{\frac {1}{2}}}{\rm {e}}^{-z}\left[1+{\frac {1}{12z}}+{\frac {1}{288z^{2}}}-{\frac {139}{51840z^{3}}}-{\frac {571}{2488320z^{4}}}+{\frac {163879}{209018880z^{5}}}+{\mathcal {O}}\left({\frac {1}{z^{6}}}\right)\right]\ .} {\displaystyle \Gamma (z)={\sqrt {2\pi }},円z^{z-{\frac {1}{2}}}{\rm {e}}^{-z}\left[1+{\frac {1}{12z}}+{\frac {1}{288z^{2}}}-{\frac {139}{51840z^{3}}}-{\frac {571}{2488320z^{4}}}+{\frac {163879}{209018880z^{5}}}+{\mathcal {O}}\left({\frac {1}{z^{6}}}\right)\right]\ .}

De Γ(z+ 1⁄2)

[modifier | modifier le code ]

L’équivalent en z+1⁄2 vaut :

Γ ( z + 1 2 ) = 2 π z z e − z + β ( z )   , {\displaystyle \Gamma \left(z+{\frac {1}{2}}\right)={\sqrt {2\pi }},円z^{z}{\rm {e}}^{-z+\beta (z)}\ ,} {\displaystyle \Gamma \left(z+{\frac {1}{2}}\right)={\sqrt {2\pi }},円z^{z}{\rm {e}}^{-z+\beta (z)}\ ,}

avec :

β ( z ) = ∑ k = 1 ∞ ( 1 2 2 k − 1 − 1 ) B 2 k 2 k ( 2 k − 1 ) z 2 k − 1   , {\displaystyle \beta (z)=\sum _{k=1}^{\infty }{\frac {\left({\frac {1}{2^{2k-1}}}-1\right)B_{2k}}{2k(2k-1)z^{2k-1}}}\ ,} {\displaystyle \beta (z)=\sum _{k=1}^{\infty }{\frac {\left({\frac {1}{2^{2k-1}}}-1\right)B_{2k}}{2k(2k-1)z^{2k-1}}}\ ,}

d’où le développement asymptotique :

Γ ( z + 1 2 ) = 2 π z z e − z [ 1 − 1 24 z + 1 1152 z 2 + 1003 414720 z 3 − 4027 39813120 z 4 − 5128423 6688604160 z 5 + O ( 1 z 6 ) ]   . {\displaystyle \Gamma \left(z+{\frac {1}{2}}\right)={\sqrt {2\pi }},円z^{z}{\rm {e}}^{-z}\left[1-{\frac {1}{24z}}+{\frac {1}{1152z^{2}}}+{\frac {1003}{414720z^{3}}}-{\frac {4027}{39813120z^{4}}}-{\frac {5128423}{6688604160z^{5}}}+{\mathcal {O}}\left({\frac {1}{z^{6}}}\right)\right]\ .} {\displaystyle \Gamma \left(z+{\frac {1}{2}}\right)={\sqrt {2\pi }},円z^{z}{\rm {e}}^{-z}\left[1-{\frac {1}{24z}}+{\frac {1}{1152z^{2}}}+{\frac {1003}{414720z^{3}}}-{\frac {4027}{39813120z^{4}}}-{\frac {5128423}{6688604160z^{5}}}+{\mathcal {O}}\left({\frac {1}{z^{6}}}\right)\right]\ .}

Cas général

[modifier | modifier le code ]

De manière plus générale, pour |a| < |z|, l’équivalent en z + a ∉ Z- vaut :

Γ ( z + a ) = 2 π z z + a − 1 2 exp ⁡ ( − z − ∑ k = 2 ∞ B k ( a ) k ( k − 1 ) ( − z ) k − 1 ) {\displaystyle \Gamma (z+a)={\sqrt {2\pi }},円z^{z+a-{\frac {1}{2}}}\exp \left(-z-\sum _{k=2}^{\infty }{\frac {B_{k}(a)}{k(k-1)(-z)^{k-1}}}\right)} {\displaystyle \Gamma (z+a)={\sqrt {2\pi }},円z^{z+a-{\frac {1}{2}}}\exp \left(-z-\sum _{k=2}^{\infty }{\frac {B_{k}(a)}{k(k-1)(-z)^{k-1}}}\right)}

où Bk sont les polynômes de Bernoulli.

Démonstration

Par généralisation sur les complexes de la formule de Stirling, on sait que, pour z ∉ Z- :

Γ ( z ) = 2 π z z − 1 2 exp ⁡ ( − z + ∑ k = 1 ∞ B 2 k 2 k ( 2 k − 1 ) z 2 k − 1 ) {\displaystyle \Gamma (z)={\sqrt {2\pi }},円z^{z-{\frac {1}{2}}}\exp \left(-z+\sum _{k=1}^{\infty }{\frac {B_{2k}}{2k(2k-1)z^{2k-1}}}\right)} {\displaystyle \Gamma (z)={\sqrt {2\pi }},円z^{z-{\frac {1}{2}}}\exp \left(-z+\sum _{k=1}^{\infty }{\frac {B_{2k}}{2k(2k-1)z^{2k-1}}}\right)}.

Les nombres de Bernoulli de rang impair supérieur ou égal à 3 étant nuls, on peut également écrire, par changement de variable i = 2k et en introduisant les termes (nuls) de rang impair :

Γ ( z ) = 2 π z z − 1 2 exp ⁡ ( − z + ∑ i = 2 ∞ B i i ( i − 1 ) z i − 1 ) {\displaystyle \Gamma (z)={\sqrt {2\pi }},円z^{z-{\frac {1}{2}}}\exp \left(-z+\sum _{i=2}^{\infty }{\frac {B_{i}}{i(i-1)z^{i-1}}}\right)} {\displaystyle \Gamma (z)={\sqrt {2\pi }},円z^{z-{\frac {1}{2}}}\exp \left(-z+\sum _{i=2}^{\infty }{\frac {B_{i}}{i(i-1)z^{i-1}}}\right)},

d’où :

Γ ( z + a ) = 2 π ( z + a ) z + a − 1 2 exp ⁡ ( − ( z + a ) + ∑ i = 2 ∞ B i i ( i − 1 ) ( z + a ) i − 1 ) {\displaystyle \Gamma (z+a)={\sqrt {2\pi }},円(z+a)^{z+a-{\frac {1}{2}}}\exp \left(-(z+a)+\sum _{i=2}^{\infty }{\frac {B_{i}}{i(i-1)(z+a)^{i-1}}}\right)} {\displaystyle \Gamma (z+a)={\sqrt {2\pi }},円(z+a)^{z+a-{\frac {1}{2}}}\exp \left(-(z+a)+\sum _{i=2}^{\infty }{\frac {B_{i}}{i(i-1)(z+a)^{i-1}}}\right)}.

z étant non nul, on peut factoriser z+a en z×ばつ(1+a/z) :

Γ ( z + a ) = 2 π z z + a − 1 2 ( 1 + a z ) z + a − 1 2 exp ⁡ ( − z − a + ∑ i = 2 ∞ B i i ( i − 1 ) z i − 1 ( 1 + a z ) i − 1 ) = 2 π z z + a − 1 2 exp ⁡ ( − z + ( z + a − 1 2 ) ln ⁡ ( 1 + a z ) − a + ∑ i = 2 ∞ B i i ( i − 1 ) z i − 1 ( 1 + a z ) − ( i − 1 ) ) . {\displaystyle {\begin{aligned}\Gamma (z+a)&={\sqrt {2\pi }},円z^{z+a-{\frac {1}{2}}}\left(1+{\tfrac {a}{z}}\right)^{z+a-{\frac {1}{2}}}\exp \left(-z-a+\sum _{i=2}^{\infty }{\frac {B_{i}}{i(i-1)z^{i-1}(1+{\frac {a}{z}})^{i-1}}}\right)\\&={\sqrt {2\pi }},円z^{z+a-{\frac {1}{2}}}\exp \left(-z+\left(z+a-{\frac {1}{2}}\right)\ln \left(1+{\frac {a}{z}}\right)-a+\sum _{i=2}^{\infty }{\frac {B_{i}}{i(i-1)z^{i-1}}}\left(1+{\frac {a}{z}}\right)^{-(i-1)}\right).\end{aligned}}} {\displaystyle {\begin{aligned}\Gamma (z+a)&={\sqrt {2\pi }},円z^{z+a-{\frac {1}{2}}}\left(1+{\tfrac {a}{z}}\right)^{z+a-{\frac {1}{2}}}\exp \left(-z-a+\sum _{i=2}^{\infty }{\frac {B_{i}}{i(i-1)z^{i-1}(1+{\frac {a}{z}})^{i-1}}}\right)\\&={\sqrt {2\pi }},円z^{z+a-{\frac {1}{2}}}\exp \left(-z+\left(z+a-{\frac {1}{2}}\right)\ln \left(1+{\frac {a}{z}}\right)-a+\sum _{i=2}^{\infty }{\frac {B_{i}}{i(i-1)z^{i-1}}}\left(1+{\frac {a}{z}}\right)^{-(i-1)}\right).\end{aligned}}}

Ayant posé |a| < |z|, on a |a/z| < 1, ce qui permet de développer d’une part la série de Taylor du logarithme ln(1 + x) (valable pour |x| < 1) et d’autre part le binôme négatif (1 + x)-n (valable pour |x| < 1 et n ∈ N*) :

ln ⁡ ( 1 + a z ) = − ∑ k = 1 ∞ ( − a z ) k k = − ∑ k = 1 ∞ ( − a ) k k z k = − ∑ k = 2 ∞ ( − a ) k − 1 ( k − 1 ) z k − 1 {\displaystyle \ln \left(1+{\tfrac {a}{z}}\right)=-\sum _{k=1}^{\infty }{\frac {\left(-{\frac {a}{z}}\right)^{k}}{k}}=-\sum _{k=1}^{\infty }{\frac {(-a)^{k}}{k,円z^{k}}}=-\sum _{k=2}^{\infty }{\frac {(-a)^{k-1}}{(k-1),円z^{k-1}}}} {\displaystyle \ln \left(1+{\tfrac {a}{z}}\right)=-\sum _{k=1}^{\infty }{\frac {\left(-{\frac {a}{z}}\right)^{k}}{k}}=-\sum _{k=1}^{\infty }{\frac {(-a)^{k}}{k,円z^{k}}}=-\sum _{k=2}^{\infty }{\frac {(-a)^{k-1}}{(k-1),円z^{k-1}}}},
( 1 + a z ) − ( i − 1 ) = ∑ j = 0 ∞ ( i + j − 2 j ) ( − a z ) j = ∑ j = 0 ∞ ( i + j − 2 ) ! ( − a ) j ( i − 2 ) ! j ! z j . {\displaystyle \left(1+{\tfrac {a}{z}}\right)^{-(i-1)}=\sum _{j=0}^{\infty }{\binom {i+j-2}{j}}\left(-{\frac {a}{z}}\right)^{j}=\sum _{j=0}^{\infty }{\frac {(i+j-2)!,円(-a)^{j}}{(i-2)!,円j!,円z^{j}}}.} {\displaystyle \left(1+{\tfrac {a}{z}}\right)^{-(i-1)}=\sum _{j=0}^{\infty }{\binom {i+j-2}{j}}\left(-{\frac {a}{z}}\right)^{j}=\sum _{j=0}^{\infty }{\frac {(i+j-2)!,円(-a)^{j}}{(i-2)!,円j!,円z^{j}}}.}

On a donc d’une part, par le développement du logarithme :

z ln ⁡ ( 1 + a z ) = − ∑ k = 1 ∞ ( − a ) k k z k − 1 = a − ∑ k = 2 ∞ ( − a ) k k z k − 1 {\displaystyle z\ln \left(1+{\tfrac {a}{z}}\right)=-\sum _{k=1}^{\infty }{\frac {(-a)^{k}}{k,円z^{k-1}}}=a-\sum _{k=2}^{\infty }{\frac {(-a)^{k}}{k,円z^{k-1}}}} {\displaystyle z\ln \left(1+{\tfrac {a}{z}}\right)=-\sum _{k=1}^{\infty }{\frac {(-a)^{k}}{k,円z^{k-1}}}=a-\sum _{k=2}^{\infty }{\frac {(-a)^{k}}{k,円z^{k-1}}}}

et :

a ln ⁡ ( 1 + a z ) = − a ∑ k = 2 ∞ ( − a ) k − 1 ( k − 1 ) z k − 1 = ∑ k = 2 ∞ ( − a ) k ( k − 1 ) z k − 1 {\displaystyle a\ln \left(1+{\tfrac {a}{z}}\right)=-a\sum _{k=2}^{\infty }{\frac {(-a)^{k-1}}{(k-1),円z^{k-1}}}=\sum _{k=2}^{\infty }{\frac {(-a)^{k}}{(k-1),円z^{k-1}}}} {\displaystyle a\ln \left(1+{\tfrac {a}{z}}\right)=-a\sum _{k=2}^{\infty }{\frac {(-a)^{k-1}}{(k-1),円z^{k-1}}}=\sum _{k=2}^{\infty }{\frac {(-a)^{k}}{(k-1),円z^{k-1}}}},

d’où :

( z + a − 1 2 ) ln ⁡ ( 1 + a z ) − a = a + ∑ k = 2 ∞ ( 1 k − 1 − 1 k ) ( − a ) k z k − 1 + 1 2 ∑ k = 2 ∞ ( − a ) k − 1 ( k − 1 ) z k − 1 − a = ∑ k = 2 ∞ ( − a ) k k ( k − 1 ) z k − 1 + 1 2 ∑ k = 2 ∞ ( − a ) k − 1 ( k − 1 ) z k − 1 = ∑ k = 2 ∞ ( − a ) k + k 2 ( − a ) k − 1 k ( k − 1 ) z k − 1 . {\displaystyle {\begin{aligned}\left(z+a-{\tfrac {1}{2}}\right)\ln \left(1+{\tfrac {a}{z}}\right)-a&=a+\sum _{k=2}^{\infty }\left({\frac {1}{k-1}}-{\frac {1}{k}}\right){\frac {(-a)^{k}}{z^{k-1}}}+{\frac {1}{2}}\sum _{k=2}^{\infty }{\frac {(-a)^{k-1}}{(k-1),円z^{k-1}}}-a\\&=\sum _{k=2}^{\infty }{\frac {(-a)^{k}}{k,円(k-1),円z^{k-1}}}+{\frac {1}{2}}\sum _{k=2}^{\infty }{\frac {(-a)^{k-1}}{(k-1),円z^{k-1}}}\\&=\sum _{k=2}^{\infty }{\frac {(-a)^{k}+{\frac {k}{2}},円(-a)^{k-1}}{k,円(k-1),円z^{k-1}}}.\end{aligned}}} {\displaystyle {\begin{aligned}\left(z+a-{\tfrac {1}{2}}\right)\ln \left(1+{\tfrac {a}{z}}\right)-a&=a+\sum _{k=2}^{\infty }\left({\frac {1}{k-1}}-{\frac {1}{k}}\right){\frac {(-a)^{k}}{z^{k-1}}}+{\frac {1}{2}}\sum _{k=2}^{\infty }{\frac {(-a)^{k-1}}{(k-1),円z^{k-1}}}-a\\&=\sum _{k=2}^{\infty }{\frac {(-a)^{k}}{k,円(k-1),円z^{k-1}}}+{\frac {1}{2}}\sum _{k=2}^{\infty }{\frac {(-a)^{k-1}}{(k-1),円z^{k-1}}}\\&=\sum _{k=2}^{\infty }{\frac {(-a)^{k}+{\frac {k}{2}},円(-a)^{k-1}}{k,円(k-1),円z^{k-1}}}.\end{aligned}}}

On a d’autre part, par le développement du binôme négatif, puis en procédant au changement de variable k=i+j :

B i i ( i − 1 ) z i − 1 ( 1 + a z ) − ( i − 1 ) = ∑ j = 0 ∞ B i ( i + j − 2 ) ! ( − a ) j i ! j ! z i − 1 + j = ∑ k = i ∞ B i ( k − 2 ) ! ( − a ) k − i i ! ( k − i ) ! z k − 1 = ∑ k = i ∞ ( k i ) B i ( − a ) k − i k ( k − 1 ) z k − 1   . {\displaystyle {\frac {B_{i}}{i,円(i-1),円z^{i-1}}}\left(1+{\frac {a}{z}}\right)^{-(i-1)}=\sum _{j=0}^{\infty }{\frac {B_{i},円(i+j-2)!,円(-a)^{j}}{i!,円j!,円z^{i-1+j}}}=\sum _{k=i}^{\infty }{\frac {B_{i},円(k-2)!,円(-a)^{k-i}}{i!,円(k-i)!,円z^{k-1}}}=\sum _{k=i}^{\infty }{\binom {k}{i}}{\frac {B_{i},円(-a)^{k-i}}{k,円(k-1),円z^{k-1}}}\ .} {\displaystyle {\frac {B_{i}}{i,円(i-1),円z^{i-1}}}\left(1+{\frac {a}{z}}\right)^{-(i-1)}=\sum _{j=0}^{\infty }{\frac {B_{i},円(i+j-2)!,円(-a)^{j}}{i!,円j!,円z^{i-1+j}}}=\sum _{k=i}^{\infty }{\frac {B_{i},円(k-2)!,円(-a)^{k-i}}{i!,円(k-i)!,円z^{k-1}}}=\sum _{k=i}^{\infty }{\binom {k}{i}}{\frac {B_{i},円(-a)^{k-i}}{k,円(k-1),円z^{k-1}}}\ .}

Puisque ( k i ) = 0 {\displaystyle {\tbinom {k}{i}}=0} {\displaystyle {\tbinom {k}{i}}=0} pour k < i, et i valant au moins 2, on peut étendre la somme ci-dessus pour k allant de 2 (en deçà, on aurait la forme indéterminée 0/0) à i – 1 (somme de i – 2 termes, donc au pire une somme vide, valide, si i = 2) :

B i i ( i − 1 ) z i − 1 ( 1 + a z ) − ( i − 1 ) = ∑ k = 2 ∞ ( k i ) B i ( − a ) k − i k ( k − 1 ) z k − 1 {\displaystyle {\frac {B_{i}}{i,円(i-1),円z^{i-1}}}\left(1+{\frac {a}{z}}\right)^{-(i-1)}=\sum _{k=2}^{\infty }{\binom {k}{i}}{\frac {B_{i},円(-a)^{k-i}}{k,円(k-1),円z^{k-1}}}} {\displaystyle {\frac {B_{i}}{i,円(i-1),円z^{i-1}}}\left(1+{\frac {a}{z}}\right)^{-(i-1)}=\sum _{k=2}^{\infty }{\binom {k}{i}}{\frac {B_{i},円(-a)^{k-i}}{k,円(k-1),円z^{k-1}}}}.

On rappelle que les polynômes de Bernoulli vérifient :

B k ( x ) = ∑ i = 0 k ( k i ) B i x k − i = B 0 x k + k B 1 x k − 1 + ∑ i = 2 k ( k i ) B i x k − i = x k − k 2 x k − 1 + ∑ i = 2 k ( k i ) B i x k − i {\displaystyle B_{k}(x)=\sum _{i=0}^{k}{\binom {k}{i}}B_{i},円x^{k-i}=B_{0},円x^{k}+k,円B_{1},円x^{k-1}+\sum _{i=2}^{k}{\binom {k}{i}}B_{i},円x^{k-i}=x^{k}-{\frac {k}{2}},円x^{k-1}+\sum _{i=2}^{k}{\binom {k}{i}}B_{i},円x^{k-i}} {\displaystyle B_{k}(x)=\sum _{i=0}^{k}{\binom {k}{i}}B_{i},円x^{k-i}=B_{0},円x^{k}+k,円B_{1},円x^{k-1}+\sum _{i=2}^{k}{\binom {k}{i}}B_{i},円x^{k-i}=x^{k}-{\frac {k}{2}},円x^{k-1}+\sum _{i=2}^{k}{\binom {k}{i}}B_{i},円x^{k-i}},

ainsi que :

B k ( − x ) = ( − 1 ) k [ B k ( x ) + k x k − 1 ] = ( − 1 ) k B k ( x ) − k ( − x ) k − 1 {\displaystyle B_{k}(-x)=(-1)^{k}\left[B_{k}(x)+k,円x^{k-1}\right]=(-1)^{k}B_{k}(x)-k,円(-x)^{k-1}} {\displaystyle B_{k}(-x)=(-1)^{k}\left[B_{k}(x)+k,円x^{k-1}\right]=(-1)^{k}B_{k}(x)-k,円(-x)^{k-1}},

d’où :

∑ i = 2 ∞ B i i ( i − 1 ) z i − 1 ( 1 + a z ) − ( i − 1 ) = ∑ k = 2 ∞ ∑ i = 2 ∞ ( k i ) B i ( − a ) k − i k ( k − 1 ) z k − 1 = ∑ k = 2 ∞ B k ( − a ) − ( − a ) k + k 2 ( − a ) k − 1 k ( k − 1 ) z k − 1 = ∑ k = 2 ∞ ( − 1 ) k B k ( a ) − k ( − a ) k − 1 − ( − a ) k + k 2 ( − a ) k − 1 k ( k − 1 ) z k − 1 = ∑ k = 2 ∞ ( − 1 ) k B k ( a ) − ( − a ) k − k 2 ( − a ) k − 1 k ( k − 1 ) z k − 1 = − ∑ k = 2 ∞ B k ( a ) k ( k − 1 ) ( − z ) k − 1 − [ ( z + a − 1 2 ) ln ⁡ ( 1 + a z ) − a ] . {\displaystyle {\begin{aligned}\sum _{i=2}^{\infty }{\frac {B_{i}}{i,円(i-1),円z^{i-1}}}\left(1+{\frac {a}{z}}\right)^{-(i-1)}&=\sum _{k=2}^{\infty }\sum _{i=2}^{\infty }{\binom {k}{i}}{\frac {B_{i},円(-a)^{k-i}}{k,円(k-1),円z^{k-1}}}=\sum _{k=2}^{\infty }{\frac {B_{k}(-a)-(-a)^{k}+{\frac {k}{2}},円(-a)^{k-1}}{k,円(k-1),円z^{k-1}}}\\&=\sum _{k=2}^{\infty }{\frac {(-1)^{k}B_{k}(a)-k,円(-a)^{k-1}-(-a)^{k}+{\frac {k}{2}},円(-a)^{k-1}}{k,円(k-1),円z^{k-1}}}\\&=\sum _{k=2}^{\infty }{\frac {(-1)^{k}B_{k}(a)-(-a)^{k}-{\frac {k}{2}},円(-a)^{k-1}}{k,円(k-1),円z^{k-1}}}\\&=-\sum _{k=2}^{\infty }{\frac {B_{k}(a)}{k,円(k-1),円(-z)^{k-1}}}-\left[\left(z+a-{\tfrac {1}{2}}\right)\ln \left(1+{\tfrac {a}{z}}\right)-a\right].\end{aligned}}} {\displaystyle {\begin{aligned}\sum _{i=2}^{\infty }{\frac {B_{i}}{i,円(i-1),円z^{i-1}}}\left(1+{\frac {a}{z}}\right)^{-(i-1)}&=\sum _{k=2}^{\infty }\sum _{i=2}^{\infty }{\binom {k}{i}}{\frac {B_{i},円(-a)^{k-i}}{k,円(k-1),円z^{k-1}}}=\sum _{k=2}^{\infty }{\frac {B_{k}(-a)-(-a)^{k}+{\frac {k}{2}},円(-a)^{k-1}}{k,円(k-1),円z^{k-1}}}\\&=\sum _{k=2}^{\infty }{\frac {(-1)^{k}B_{k}(a)-k,円(-a)^{k-1}-(-a)^{k}+{\frac {k}{2}},円(-a)^{k-1}}{k,円(k-1),円z^{k-1}}}\\&=\sum _{k=2}^{\infty }{\frac {(-1)^{k}B_{k}(a)-(-a)^{k}-{\frac {k}{2}},円(-a)^{k-1}}{k,円(k-1),円z^{k-1}}}\\&=-\sum _{k=2}^{\infty }{\frac {B_{k}(a)}{k,円(k-1),円(-z)^{k-1}}}-\left[\left(z+a-{\tfrac {1}{2}}\right)\ln \left(1+{\tfrac {a}{z}}\right)-a\right].\end{aligned}}}

Donc, pour |a| < |z| :

Γ ( z + a ) = 2 π z z + a − 1 2 exp ⁡ ( − z − ∑ k = 2 ∞ B k ( a ) k ( k − 1 ) ( − z ) k − 1 ) {\displaystyle \Gamma (z+a)={\sqrt {2\pi }},円z^{z+a-{\frac {1}{2}}}\exp \left(-z-\sum _{k=2}^{\infty }{\frac {B_{k}(a)}{k(k-1)(-z)^{k-1}}}\right)} {\displaystyle \Gamma (z+a)={\sqrt {2\pi }},円z^{z+a-{\frac {1}{2}}}\exp \left(-z-\sum _{k=2}^{\infty }{\frac {B_{k}(a)}{k(k-1)(-z)^{k-1}}}\right)}.

En posant a valant respectivement 0, 1⁄2 et 1, et connaissant les valeurs particulières des polynômes de Bernoulli en ces points, on retrouve immédiatement les équivalents en z, z + 1⁄2 et z + 1 mentionnés plus hauts.

La première occurrence d'un produit qui donnera naissance ultérieurement à la fonction gamma est due à Daniel Bernoulli [12] dans une lettre à Christian Goldbach.

En notation moderne[13]

x ! = lim n → ∞ ( n + 1 + x 2 ) x − 1 ∏ i = 1 n i + 1 i + x {\displaystyle x!=\lim _{n\rightarrow \infty }\left(n+1+{\frac {x}{2}}\right)^{x-1}\prod _{i=1}^{n}{\frac {i+1}{i+x}}} {\displaystyle x!=\lim _{n\rightarrow \infty }\left(n+1+{\frac {x}{2}}\right)^{x-1}\prod _{i=1}^{n}{\frac {i+1}{i+x}}}.

En 1729 également, Euler entreprend l'étude de ce produit et lui donne sa forme intégrale[14] ,[15] .

C'est Legendre qui, en 1811, note cette fonction Γ {\displaystyle \Gamma } {\displaystyle \Gamma }, en apportant de nombreux compléments à son étude[14] ,[16] .

L'article de Borwein et Corless[17] passe en revue trois siècles de travaux mathématiques sur la fonction gamma.

Notes et références

[modifier | modifier le code ]
  1. ↑ a et b Voir par exemple le début de ce devoir corrigé sur Wikiversité.
  2. ↑ (en) Harold M. Edwards, Riemann's Zeta Function, États-Unis , Dover Publications , .
  3. ↑ Pour le cas particulier où z est un réel strictement positif, voir l'article Théorème de Bohr-Mollerup. Pour le cas général, voir cet exercice corrigé sur Wikiversité.
  4. ↑ (de) O. Schlömilch, « Einiges über die Eulerischen Integrale der zweiten Art », Archiv der Mathematik und Physik , vol. 4,‎ , p. 171 (lire en ligne).
  5. ↑ (en) J. L. W. V. Jensen, « An elementary exposition of the theory of the Gamma function », Ann. of Math. , 2e série, vol. 17, no 3,‎ , p. 124-166 (JSTOR 2007272 ) (p. 128).
  6. ↑ « En 1844, 32 ans avant le célèbre travail de Weierstrass sur les fonctions entières » : (en) S. S. Dragomir, Ravi Agarwal (en) et N. S. Barnett, « Inequalities for Beta and Gamma functions via some classical and new integral inequalities », J. Inequal. Appl. (nl) , vol. 5, no 2,‎ , p. 103-165 (lire en ligne) (p. 107).
  7. ↑ (en) Jesús Guillera et Jonathan Sondow, « Double integrals and infinite products for some classical constants via analytic continuations of Lerch's transcendent », The Ramanujan Journal , vol. 16, no 3,‎ , p. 247-270 (DOI 10.1007/s11139-007-9102-0 , arXiv math/0506319 ).
  8. ↑ (en) Karl Rawer, Wave Propagation in the Ionosphere, Dordrecht, Kluwer Academic Publishers, .
  9. ↑ D'après (de) O. R. Rocktäschel, Methoden zur Berechnung der Gammafunktion für komplexes Argument, université technique de Dresde, , thèse de doctorat.
  10. ↑ (de) P. E. Böhmer, Differenzengleichungen und bestimmte Integrale, Leipzig, Köhler Verlag, .
  11. ↑ (en) Serge Lang, Complex Analysis, Springer, coll. « GTM » (no 103), , 489 p. (ISBN 978-0-387-98592-3, lire en ligne), p. 418.
  12. ↑ Paul Heinrich Fuss, Correspondance mathématique et physique de quelques célèbres géomètres du XVIIIe siècle, vol. II, St. Pétersbourg, Académie impériale des sciences, (lire en ligne), p. 324-325.
  13. ↑ (en) Detlef Gronau, « Why is the gamma function so as it is? », Teaching Mathematics and Computer Science, vol. 1, no 1,‎ , p. 43-53.
  14. ↑ a et b G. K. Srinivasan, « The Gamma function: An Eclectic Tour », The American Mathematical Monthly , vol. 114, no 4,‎ , p. 297-315 (DOI 10.1080/00029890.2007.11920418 )
  15. ↑ L. Euler, « De progressionibus transcendentibus seu quarum termini generales algebraice dari nequeunt », sur scholarlycommons.pacific.edu.
  16. ↑ A.-M. Legendre, Exercices de calcul intégral sur divers ordres de transcendantes et sur les quadratures, t. 1, Vve Courcier (Paris), (lire en ligne), p. 221
  17. ↑ (en) Jonathan M. Borwein et Robert M. Corless, Gamma and Factorial in the Monthly, 17 mars 2017 arXiv:1703.05349

Voir aussi

[modifier | modifier le code ]

Sur les autres projets Wikimedia :

Articles connexes

[modifier | modifier le code ]

Bibliographie

[modifier | modifier le code ]

Liens externes

[modifier | modifier le code ]

AltStyle によって変換されたページ (->オリジナル) /