The Exponential Integral
Not to be confused with the Exponential Anti-derivative (ie., the integral of exp(z)), the exponential integral, Ei(z), cannot be evaluated in closed form (with finite terms of elementary functions) and, as such, it is a special function. It is surprising when considering the apparent simplicity of the integrands, but the integration does not lead to a trivial solution. Like other special functions, the exponential integral can be evaluated by approximation to an arbitrary accuracy using infinite series. For a more detailed description see Wolfram functions.
Definition
The Exponential Integral is defined by the integral,
E1(x)=-∫exp(-x)xdx
Series Expansion
The function can be approximation using the series expansion, which is useful for evaluation in software,
E1(z)=-γ-lnz-∞∑n=1(-1)nznn!n
Hypergeometric Form
It is possible to express the Exponential integral as a hypergeometric function,
E1(z)=-γ-lnz+ 2F2(1,1;2,2;-z)z
Graph
The Polynomial Denominator
Another useful variant is the integral with a square denominator
∫exp(-x)x2dx
We can go about evaluating this by finding the exponential integral by integration by parts. Expressing the exponential integral (with x1 in the denominator) in the form,
∫exp(-x)xdx=∫u′vdx=uv-∫v′udx
We have, u′=exp(-x) and v=1x, and so u=-exp(x) and v′=-1x2. Hence,
∫exp(-x)xdx=-exp(-x)x-∫exp(-x)x2
Now we recognise the integral on the left as being the definition -E1(x), and hence we can solve it for the integral on the right. ie.,
∫exp(-x)x2dx=E1(x)-exp(-x)x
The Double Integral
Another property which can be useful is the integral of the exponential integral,
∫∫exp(-x)xdxdx=∫E1(x)dx=xE1(x)-exp(-x)