# Constante d'Euler-Mascheroni

En mathématiques, la constante d'Euler-Mascheroni, ou constante d'Euler, est une constante mathématique définie comme la limite de la différence entre la série harmonique et le logarithme népérien. 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 γ 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 (suite A001620 de l'OEIS) sont : γ ≈ 0,5772156649.
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 4e décimale, 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 2e décimale dans le sens naturel, et la 4e décimale 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 γ. Puis Lorenzo Mascheroni en proposa 32 en 1790, mais avec une erreur à partir de la 20e, erreur corrigée en 1809 par Johann Georg von Soldner. Donald Knuth donne 1 271 décimales 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 2020, le record vérifié semble être détenu par Seungmin Kim et Ian Cutress, avec plus de 600 milliards de décimales (600 000 000 100 pour être précis) en utilisant y-cruncher.
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 242 080 chiffres (Havil 2003, p. 97).
En 2009, Alexander Aptekarev a prouvé qu'au moins l'une des deux constantes d'Euler–Mascheroni et d'Euler–Gompertz est irrationnelle. Ce résultat a été amélioré en 2012 par Tanguy Rivoal en prouvant qu'au moins l'une des deux est transcendante$^,$.

## Formules diverses
### Formules intégrales

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 (Sondow 2003, Sondow 2005) d'exprimer γ 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 (Sondow 2005) :
    $\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 (Sondow 2010) :
    $\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.
* 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.
* 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.

## Généralisation

Il est possible de généraliser le sujet en définissant les constantes suivantes, appelées constantes de Stieltjes :
    $\gamma (m)=\lim _{n\to \infty }\left(\sum _{k=1}^{n}{\frac {(\ln k)^{m}}{k}}-{\frac {(\ln n)^{m+1}}{m+1}}\right)$.
On constate que $\gamma (0)=\gamma$, la constante d'Euler.