Constante d'Euler-Mascheroni

Modèle:Voir homonymes Modèle:Sources à lier En mathématiques, la constante d'Euler-Mascheroni, ou constante d'Euler, est une constante mathématique, utilisée principalement en théorie des nombres, définie comme la limite de la différence entre la série harmonique et le logarithme naturel. On la note usuellement $\gamma$ (gamma minuscule).

Liste des nombres γ - ζ(3) - √2 - √3 - √5 - φ - α - e - π - δ
Binaire	0,100 100 111 100 010 001 1…
Décimal	0,577 215 664 901 532 860 6…
Hexadécimal	0,93C 467 E37 DB0 C7A 4D1 B…
Fraction continue	$0+{\cfrac {1}{1+{\cfrac {1}{1+{\cfrac {1}{2+{\cfrac {1}{1+{\cfrac {1}{\ \ddots \ {}}}}}}}}}}}$ (on ignore encore si cette fraction continue se termine ou non).

Définition

La constante d'Euler-Mascheroni $\gamma$ est définie de la manière suivante :

\gamma =\lim _{n\to \infty }\left(1+{\frac {1}{2}}+{\frac {1}{3}}+{\frac {1}{4}}+\dots +{\frac {1}{n}}-\ln(n)\right)

.

De façon condensée, on obtient :

\gamma =\lim _{n\to \infty }\left(\sum _{k=1}^{n}{\frac {1}{k}}-\ln(n)\right)

.

La constante peut également être définie sous la forme explicite d'une série (telle qu'elle fut d'ailleurs introduite par Euler) :

\gamma =\sum _{k=1}^{\infty }\left[{\frac {1}{k}}-\ln \left(1+{\frac {1}{k}}\right)\right]

.

La série harmonique diverge, tout comme la suite de terme général $\ln(n)$ ; l'existence de cette constante indique que les deux expressions sont asymptotiquement liées.

Valeur approchée et propriétés

Les 10 premières décimales de la constante d'Euler-Mascheroni (Modèle:OEIS) sont : Modèle:Formule ≈ 0,577 215 664 9.

Le calcul au moyen de la suite $\sum _{k=1}^{n}{\frac {1}{k}}-\ln(n)$ est extrêmement lent et imprécis. Il présente néanmoins un intérêt pédagogique pour se sensibiliser aux problèmes de propagation d'erreurs d'arrondi. En simple précision, pour 100 000 termes, en sommant dans l'ordre naturel, il y a une erreur sur la Modèle:4e, erreur beaucoup plus faible si la somme est effectuée dans l'ordre inverse (du plus petit au plus grand), ou si on utilise l'algorithme de Kahan (voir somme (algorithmique)). Pour un million de termes, l'erreur atteint la Modèle:2e dans le sens naturel, et la Modèle:4e dans le sens inverse ; par contre, par la méthode de Kahan, on a atteint les 6 décimales exactes.

Des méthodes plus efficaces doivent être mises en œuvre pour obtenir une précision suffisante. Par exemple, l'utilisation de la formule d'Euler-Maclaurin permet d'obtenir des développements asymptotiques tels que :

\gamma =\sum _{k=1}^{n}{\frac {1}{k}}-\ln(n)-{\frac {1}{2n}}+{\frac {1}{12n^{2}}}-{\frac {1}{120n^{4}}}+\dots

.

Cela permit à Euler d'obtenir 16 décimales de $\gamma$ . Puis Lorenzo Mascheroni en proposa 32 en 1790, mais avec une erreur à partir de la Modèle:20e, erreur corrigée en 1809 par Johann Georg von Soldner. Donald Knuth donne Modèle:Unité en 1962, Thomas Papanikolaou donne un million de décimales en 1997, P. Dechimel et X. Gourdon en donnent cent millions deux ans plus tard. En 2008, le record est de dix milliards de décimales, établi par Shigeru Kondo et Steve Pagliarulo.

On ignore toujours si la constante d'Euler-Mascheroni est ou non un nombre rationnel. Cependant, l'analyse en fraction continue de la constante indique que si elle est rationnelle, le dénominateur de sa fraction irréductible possède plus de Modèle:Unité Modèle:Harv.

Formules diverses

Formules intégrales

Constante gamma (somme des aires des surfaces bleues) comme intégrale de 1/E(x) – 1/x de 1 à ∞.

La constante d'Euler-Mascheroni intervient dans plusieurs intégrales :

\gamma =\int _{1}^{\infty }\left({1 \over E(x)}-{1 \over x}\right)\,{\rm {d}}x

(où

E

est la fonction partie entière)

=1-\int _{1}^{\infty }\ {\frac {x-E(x)}{x^{2}}}\,{\rm {d}}x

=-\int _{0}^{1}{\ln \ln \left({\frac {1}{x}}\right)}\,{\rm {d}}x

=\int _{0}^{1}\left({\frac {1}{\ln(x)}}+{\frac {1}{1-x}}\right)\,{\rm {d}}x

=\int _{0}^{\infty }{\left({\frac {1}{1-\mathrm {e} ^{-x}}}-{\frac {1}{x}}\right)\mathrm {e} ^{-x}}\,{\rm {d}}x

=\int _{0}^{\infty }{{\frac {1}{x}}\left({\frac {1}{1+x}}-\mathrm {e} ^{-x}\right)}\,{\rm {d}}x

.

Il est possible (Modèle:Harvsp, Modèle:Harvsp) d'exprimer $\gamma$ sous forme d'une intégrale double (avec ici la série équivalente) :

\gamma =\int _{0}^{1}\int _{0}^{1}{\frac {x-1}{(1-x\,y)\ln(x\,y)}}\,{\rm {d}}x\,{\rm {d}}y=\sum _{n=1}^{\infty }\left({\frac {1}{n}}-\ln \left({\frac {n+1}{n}}\right)\right)

.

Une autre constante s'exprime de façon analogue Modèle:Harv :

\ln \left({\frac {4}{\pi }}\right)=\int _{0}^{1}\int _{0}^{1}{\frac {x-1}{(1+x\,y)\ln(x\,y)}}\,{\rm {d}}x\,{\rm {d}}y=\sum _{n=1}^{\infty }(-1)^{n-1}\left({\frac {1}{n}}-\ln \left({\frac {n+1}{n}}\right)\right)

.

Ces deux constantes sont également liées par deux séries Modèle:Harv :

\gamma =\sum _{n=1}^{\infty }{\frac {N_{1}(n)+N_{0}(n)}{2n(2n+1)}}

\ln \left({\frac {4}{\pi }}\right)=\sum _{n=1}^{\infty }{\frac {N_{1}(n)-N_{0}(n)}{2n(2n+1)}}

où $N_{1}(n)$ et $N_{0}(n)$ sont le nombre de 1 et de 0 dans l'écriture de $n$ en base 2.

On trouvera d'autres expressions non classiques de la constante d'Euler dans l'article « Mesure secondaire ».

Formules en relation avec certaines fonctions analytiques

La constante d'Euler-Mascheroni possède des liens avec d'autres fonctions analytiques particulières :

Fonction gamma :
- $\Gamma (z)=\int _{0}^{\infty }\mathrm {e} ^{-t}t^{z-1}\,{\rm {d}}t={\frac {\mathrm {e} ^{-\gamma z}}{z}}\prod _{n=1}^{\infty }{\frac {n\mathrm {e} ^{z/n}}{n+z}}$ ,
- $\Gamma (x)={\frac {1}{x}}-\gamma +o(1)$ quand x tend vers 0,
- $\Gamma '(1)=\int _{0}^{\infty }{\mathrm {e} ^{-x}\ln x}\,{\rm {d}}x=-\gamma$ ,
- $\Gamma ''(1)=\int _{0}^{\infty }{\mathrm {e} ^{-x}(\ln(x))^{2}}\,{\rm {d}}x=\gamma ^{2}+{\frac {\pi ^{2}}{6}}$ ,
- $\Gamma '(1/2)=4\int _{0}^{\infty }{\mathrm {e} ^{-x^{2}}\ln(x)}\,{\rm {d}}x=-(\gamma +2\ln 2){\sqrt {\pi }}$ ;
Fonction exponentielle intégrale :
- $E_{1}(z)=\int _{z}^{\infty }{\mathrm {e} ^{-t} \over t}\,{\rm {d}}t=\int _{1}^{\infty }{\mathrm {e} ^{-zt} \over t}\,{\rm {d}}t=\mathrm {e} ^{-z}\int _{0}^{\infty }{\mathrm {e} ^{-zt} \over {1+t}}\,{\rm {d}}t$
  $={\mathrm {e} ^{-z} \over z}\int _{0}^{\infty }{\mathrm {e} ^{-t} \over {1+t/z}}\,{\rm {d}}t=-\ln z-\gamma +\sum _{n=1}^{\infty }{(-1)^{n-1}z^{n} \over n\cdot n!}$ ;
Fonction logarithme intégral :
- ${\rm {pour}}\;x>1,\ \mathrm {li} (x)=\gamma +\ln(\ln(x))+\sum _{n=1}^{\infty }{\frac {\ln(x)^{n}}{n\cdot n!}}$ ;
Fonction cosinus intégral :
- $\forall x>0,\ \mathrm {Ci} (x)=\gamma +\ln(x)+\sum _{n=1}^{+\infty }(-1)^{n}{\frac {x^{2n}}{(2n)!(2n)}}$ ;
Fonction psi :
- $\psi (z)={\Gamma '(z) \over \Gamma (z)}=-\gamma -{1 \over z}+\sum _{n=1}^{\infty }{1 \over n}-{1 \over n+z}$ ,
  en particulier, $\psi (1)=\Gamma '(1)=-\gamma$ et $\sum _{k=1}^{n}{1 \over k}=\psi (n+1)+\gamma$ ;
Fonction zêta de Riemann :
- $\zeta (1+x)={\frac {1}{x}}+\gamma +o(1)$ quand x tend vers 0,
- $\sum _{n=2}^{\infty }{\frac {1}{n}}(\zeta (n)-1)=1-\gamma$ ,
- $\sum _{n=1}^{\infty }{\frac {1}{(2n+1)2^{2n}}}(\zeta (2n+1)-1)=1-\gamma -\ln {\frac {3}{2}}$ ,
- $\sum _{n=2}^{\infty }{\frac {(-1)^{n}}{n}}\zeta (n)=\gamma$ .

Formules en relation avec certaines fonctions arithmétiques

Dans ce paragraphe, p désigne un nombre premier.

$\lim _{n\to \infty }{\frac {1}{\ln(n)}}\prod _{p\leq n}\left(1-{\frac {1}{p}}\right)^{-1}=\mathrm {e} ^{\gamma }$ (théorème de Mertens).
$\lim _{n\to \infty }{\frac {1}{\ln(n)}}\prod _{p\leq n}\left(1+{\frac {1}{p}}\right)={\frac {6\mathrm {e} ^{\gamma }}{\pi ^{2}}}$ .
Soit $\Lambda$ la fonction de von Mangoldt, définie sur les entiers par $\Lambda (n)=\ln(p)$ si n est une puissance du nombre premier p et $\Lambda (n)=0$ sinon. Alors $\sum _{n=2}^{\infty }{\frac {\Lambda (n)-1}{n}}=-2\gamma$ .
Soit $d(n)$ le nombre de diviseurs de n (y compris 1 et n lui-même). Alors ${\frac {1}{n}}\sum _{k=1}^{n}d(k)=\ln(n)+2\gamma -1+O\left({\frac {1}{\sqrt {n}}}\right)$ quand n tend vers l'infini^[1].
Soit $\sigma (n)$ la somme des diviseurs de l'entier n. Alors $\limsup {\frac {\sigma (n)}{n\ln(\ln(n))}}=\mathrm {e} ^{\gamma }$ , où lim sup désigne la limite supérieure de la suite^[1].
Soit $\varphi$ la fonction indicatrice d'Euler. Alors $\liminf {\frac {\varphi (n)\ln(\ln(n))}{n}}=\mathrm {e} ^{-\gamma }$ , où lim inf désigne la limite inférieure de la suite^[1].