Exponentielle intégrale

En mathématiques, la fonction exponentielle intégrale, habituellement notée $Ei$ , est définie par :

{\mbox{Ei}}:{\begin{cases}\mathbb {R} ^{*}\to \mathbb {R} \\x\mapsto {\mbox{Ei}}(x)=-\int _{-x}^{\infty }{\frac {{\rm {e}}^{-t}}{t}}\,\mathrm {d} t\,=\int _{-\infty }^{x}{\frac {{\rm {e}}^{t}}{t}}\,\mathrm {d} t.\end{cases}}

Comme l'intégrale de la fonction inverse ( $t\mapsto {\frac {1}{t}}$ ) diverge en 0, cette définition doit être comprise, si $x > 0$ , comme une valeur principale de Cauchy.

Représentation graphique de la fonction exponentielle intégrale.

Lien avec le logarithme intégral modifier

La fonction $Ei$ est liée à la fonction $li$ (logarithme intégral) par :

{\rm {Ei}}(x)=\operatorname {li} \left(\mathrm {e} ^{x}\right).

Développement en série de $Ei$ modifier

Représentation graphique des fonctions

E 1

(en haut) et

Ei

(en bas), pour x > 0.

L'exponentielle intégrale a pour développement en série :

{\mbox{Ei}}(x)=\gamma +\ln |x|+\sum _{k=1}^{\infty }{\frac {x^{k}}{k\cdot k!}},

où $γ$ est la constante d'Euler-Mascheroni.

Les fonctions $E n$ modifier

L'exponentielle intégrale est reliée à une autre fonction, notée $E 1$ définie, pour x > 0, par :

{\rm {E}}_{1}(x)=\int _{x}^{\infty }{\frac {{\rm {e}}^{-t}}{t}}\,\mathrm {d} t=-\int _{-\infty }^{-x}{\frac {{\rm {e}}^{t}}{t}}\,\mathrm {d} t.

On dispose alors de la relation, pour $x > 0$ :

{\rm {Ei}}(-x)=-{\rm {E}}_{1}(x).

Les deux fonctions s'expriment en fonction de la fonction entière définie par :

{\rm {Ein}}(x)=\int _{0}^{x}(1-{\rm {e}}^{-t})\,{\frac {\mathrm {d} t}{t}}=\sum _{k=1}^{\infty }{\frac {(-1)^{k+1}x^{k}}{k\;k!}}.

En effet, on peut montrer que, pour x > 0 :

{\rm {E}}_{1}(x)=-\gamma -\ln(x)+{\rm {Ein}}(x)

et

{\rm {Ei}}(-x)=\gamma +\ln x-{\rm {Ein}}(x).

La relation donnée pour $E 1$ permet d'étendre cette fonction sur tout ouvert simplement connexe du plan complexe ne contenant pas 0, en prenant une détermination du logarithme sur ce plan. On prend généralement comme ouvert le plan complexe privé des réels strictement négatifs.

Plus généralement, on définit, pour tout entier n strictement positif, la fonction $E n$ par :

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

Ces fonctions sont reliées par la relation :

\mathrm {E} _{n}(x)={\frac {{\rm {e}}^{-x}}{x}}-{\frac {n}{x}}\mathrm {E} _{n+1}(x)

Calcul de $E 1$ modifier

La fonction $E 1$ ne possède pas d’expression à l’aide des fonctions élémentaires usuelles, d’après un théorème dû à Liouville. Différentes méthodes peuvent être utilisées afin de calculer $E 1 (x)$ en double précision.

Pour x compris entre 0 et 2,5 modifier

On a :

\mathrm {E_{1}} \left(z\right)=-\gamma -\ln z+\sum _{k=1}^{\infty }{\frac {(-1)^{k+1}z^{k}}{k\;k!}}\qquad \left(\left|\arg z\right|<\pi \right).

Cette série convergente peut théoriquement être utilisée pour calculer $E 1 (x)$ pour tout réel $x > 0$ mais avec les opérations à virgule flottante, le résultat est inexact pour $x > 2,5$ à cause de la perte de précision relative quand on soustrait des nombres d'ordres de grandeur différents.

Pour x > 40 modifier

Erreur relative de l'approximation asymtotique pour diverses valeurs du nombre N de termes de la somme partielle : N = 1 (rouge), 2 (vert), 3 (jaune), 4 (bleu) et 5 (rose)

Il existe une série divergente permettant d'approcher $E 1$ pour les grandes valeurs de $Re(z)$ , obtenue par intégration par parties^[1], qui donne le développement asymptotique suivant :

\mathrm {E_{1}} (z)={\frac {\exp(-z)}{z}}\sum _{n=0}^{N-1}{\frac {n!}{(-z)^{n}}}+{\frac {N!}{(-z)^{N}}}\mathrm {E} _{N+1}(z)

avec $\mathrm {E} _{N+1}(z)=o({\rm {e}}^{-z})$ quand z tend vers $+\infty$ .

Afin d'avoir une précision de 64 bit (double précision), il faut utiliser la valeur $N = 40$ ^[2].

Références modifier

↑ (en) Norman Bleistein et Richard A. Handelsman, Asymptotic Expansions of Integrals, Dover, 1975, 425 p. (ISBN 978-0-486-65082-1, lire en ligne), p. 3.
↑ (en) P. Pecina, « On the Function Inverse to the Exponential Integral Function », Bulletin of the Astronomical Institutes of Czechoslovakia, vol. 37,‎ 1986, p. 8 (lire en ligne)

(en) Milton Abramowitz et Irene Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables [détail de l’édition] (lire en ligne), p. 228-230

Voir aussi modifier

Portail de l'analyse

[1] (en) Norman Bleistein et Richard A. Handelsman, Asymptotic Expansions of Integrals, Dover, 1975, 425 p. (ISBN 978-0-486-65082-1, lire en ligne), p. 3.

[Pecina-2] (en) P. Pecina, « On the Function Inverse to the Exponential Integral Function », Bulletin of the Astronomical Institutes of Czechoslovakia, vol. 37,‎ 1986, p. 8 (lire en ligne)

[1]

[2]