Exponential integral

Exponential integral is a mathematical function of a single argument which has no established notation.

Definitions

For real values of "x", the exponential integral Ei("x") can be defined as :$mbox\{Ei\}(x)=int\_\{-infty\}^xfrac\{e^t\}t,dt.,$ The definition above can be used for positive values of "x", but integral has to be understood in terms of the Cauchy principal value.

For complex values of the argument, the definition becomes ambiguous cite book
last = Abramovitz
first = Milton
others= Abramowitz and Stegun
coauthors =Irene Stegun
title =Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables
publisher = Dover
year = 1964
location = New York
pages =
url = http://www.math.sfu.ca/~cbm/aands
doi =
id =
isbn = 0-486-61272-4 ] (see http://www.math.sfu.ca/~cbm/aands/page_228.htm, formula 5.1.1); in order to avoid confusion, the following notation is used:

:

At positive values of real part of argument, this presentation can be converted to

:$\{\; m\; E\}\_1(z)\; =\; int\_1^infty\; frac\{e^\{-tz\{t\},\; dt,qquad\; Re(z)\; ge\; 0.$

The function Ei is related with E₁ as follows:

:$lim\_\{delta\; o0pm\}\{\; m\; Ei\}(-x+\{\; m\; i\}delta)\; =\; -\; \{\; m\; E\}\_1(x)\; mp\; \{\; m\; i\}\; pi,qquad\; x>0,$

:$-\{\; m\; Ei\}(x)\; =\; frac\{1\}\{2\}lim\_\{delta\; o0+\}\{\; m\; E\}\_1(-x+\{\; m\; i\}delta)\; +\; frac\{1\}\{2\}lim\_\{delta\; o0-\}\{\; m\; E\}\_1(-x+\{\; m\; i\}delta),qquad\; x>0.$

The extension of Ei to the complex plane may have cut at the negative values of argument. Then, area of analyticity of function Ei is complementary to that of E₁.

Properties

Several properties of the exponential integral below, in certain cases, allow to avoid its explicit evaluation through the definition above.

Convergent series

E₁ has logarithmic peculiarity at zero. The extraction allows to write the exponential integral in terms of convergent series: :$mbox\{Ei\}(x)\; =\; gamma+ln\; x+\; sum\_\{k=1\}^\{infty\}\; frac\{x^k\}\{k;\; k!\}\; ,,\; ~~~~~x>0$

:$E\_1(z)\; =-gamma-ln\; z+\; sum\_\{k=1\}^\{infty\}\; frac\{(-1)^\{k+1\}\; z^k\}\{k;\; k!\}\; ,,~~~~~~~~\; \{\; m\; Re\}(z)>0$

where $~gammaapprox\; 0.5772156649015328606...~$ is the Euler gamma constant.The series converges at any complex value of the argument, but definition of Ei requires that$~x!>!0~$.

Asymptotic (divergent) series

At large values of the argument, evaluation of exponential integral with convergent series above is difficult, if at all possible. For this case, there exist so-called divergent, or asymptotic series::$E\_1(z)=frac\{exp(-z)\}\{z\}\; left(sum\_\{n=0\}^\{N-1\}\; frac\{n!\}\{(-z)^n\}\; +\; \{mathcal\{Oleft(\; frac\{N!\}\{z^N\}\; ight)\; ight)$The truncated sum can be used for evaluation at $~\{\; m\; Re\; \}(z)!gg!\; 1~$.The more terms are taken into account in the sum, the larger should be the real part of the argument in order to make the truncated sum useful for the evaluation.

The relative error of the approximation above is plotted on the figure.With truncated series, the function$E\_1(z,N)=frac\{exp(-z)\}\{z\}\; sum\_\{n=0\}^\{N-1\}\; frac\{n!\}\{(-z)^n\}$approximates $E\_1(z)$ at $~\{\; m\; Re\}(z)\; gg\; 1$.The relative error$R(x)=frac\{\; E\_1(x,N)-E\_1(x)\; \}\{E\_1(x)\}$ is plotted versus $~x~$ for
$~N=1~$ (red),
$~N=2~$ (green),
$~N=3~$ (yellow),
$~N=4~$ (blue),
$~N=5~$ (pink).
The larger is $~x~$, the more terms in this expansion can be taken into account for approximation of the function (and the better is the resulting approximation). At $~x!>!40~$, the approximation$E\_1(z,40)$ is preciser, that the direct evaluation of the initial integral with standard (double) variables in most of programming languages. The modulus of the last term taken into account characterizes the error of such approximation.

Exponential and logarithmic behavior: bracketing

From the two series suggested in previous subsecitons, it follows, that $~\; E\_1~$ behaves similar to an exponential at large values of the argument and as logarithm at small values. In the range of positive values of argument, $~E\_1~$ can be bracketed with elementary functions as follows::The left-hand side of this inequality is shown in the Figure with red curve.The central part $~\{\; m\; E\}\_1(x)~$ is shown with the black curve.The right-hand side is shown with blue curve.

Relation with other functions

The exponential integral is closely related to the logarithmic integral function li("x"),

:li("x") = Ei (ln ("x")) for all positive real "x" ≠ 1.

Also closely related is a function which integrates over a different range:

:$\{\; m\; E\}\_1(x)\; =\; int\_1^infty\; frac\{e^\{-tx\{t\},\; dt\; =\; int\_x^infty\; frac\{e^\{-t\{t\}\; dt.$

This function may be regarded as extending the exponential integral to thenegative reals by

: $\{\; m\; Ei\}(-x)\; =\; -\; \{\; m\; E\}\_1(x).,$

We can express both of them in terms of an entire function,

:$\{\; m\; Ein\}(x)\; =\; int\_0^x\; (1-e^\{-t\})frac\{dt\}\{t\}=\; sum\_\{k=1\}^infty\; frac\{(-1)^\{k+1\}x^k\}\{k;\; k!\}.$

Using this function, we then may define, using the logarithm,

:

and

:$\{\; m\; Ei\}(x)\; ,=,\; gamma+ln\; x\; -\; \{\; m\; Ein\}(-x),~~~~~~x>0.$

The exponential integral may also be generalized to

:$\{\; m\; E\}\_n(x)\; =\; int\_1^infty\; frac\{e^\{-xt\{t^n\},\; dt,$

which can be written as a special case of the incomplete gamma function:

: $\{\; m\; E\}\_n(x)\; =x^\{n-1\}Gamma(1-n,x).,$

The generalized form is sometimes called the Misra function $varphi\_m(x)$, defined as

:$varphi\_m(x)=\{\; m\; E\}\_\{-m\}(x).,$

Derivatives

Functions $~\{\; m\; E\}\_n~$ are simply related with derivatives of $~\{\; m\; E\}\_1~$::$$ m E}_n} '(z) = -{ m E}_{n-1}(z),~~~~~~~~(|{ m Arg}(z)|0) However, $~n~$ is supposed to be integer; the generalization for complex $~n~$ is not yet reported in the literature, although the definition of $~\{\; m\; E\}\_n~$ through the integral could allow such a generalization.

Exponential integral of imaginary argument

---

right|200px|thumb|">$\{\; m\; E\}\_1(\; \{\; m\; i\}!~\; x)$ versus $~x~$,real part(black) and imaginary part (red).From the representation

:$\{\; m\; E\}\_1(z)\; =\; int\_1^inftyfrac\{exp(-zt)\}\{t\},\{\; m\; d\}\; t,~~~~~~(\{\; m\; Re\}(z)\; ge\; 0)$

the relation of the exponential integral with integral sinus (Si) and integral cosinus (Ci) is seen:

:$\{\; m\; E\}\_1(\; \{\; m\; i\}!~\; x)=-frac\{pi\}\{2\}+\{\; m\; Si\}(x)-\{\; m\; i\}cdot\; \{\; m\; Ci\}(x),~~~~~~~~~(x>0)$

Real and imaginary parts of $~\{\; m\; E\}\_1(x)~$ are plotted in Figure with black and red curves. The real part has logarithmic peculiarity at zero (As the Exponential integral of the real argument).

Integrals with exponential integral

Applications

* Time-dependent heat transfer
* Nonequilibrium groundwater flow in the Theis solution (called a "well function")
* Radiative transfer in stellar atmospheres
* Radial Diffusivity Equation for transient or unsteady state flow with line sources and sinks

Notes

References

*Press, William H. et al. Numerical Recipes (FORTRAN). Cambridge University Press, New York: 1989.
* Milton Abramowitz and Irene A. Stegun, eds. "Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables." New York: Dover, 1972. " [http://www.math.sfu.ca/~cbm/aands/page_227.htm (See Chapter 5)] "

* R. D. Misra, Proc. Cambridge Phil. Soc. 36, 173 (1940)
*S. Chandrasekhar, Radiative transfer, reprinted 1960, Dover

External links

*MathWorld|urlname=ExponentialIntegral|title=Exponential Integral
*MathWorld|urlname=En-Function|title="En"-Function
* [http://functions.wolfram.com/GammaBetaErf/ExpIntegralEi/ Formulas and identities for Ei]