Power Series Expansion for Exponential Integral Function/Formulation 1

Theorem
Let $\Ei: \R_{>0} \to \R$ denote the exponential integral function:


 * $\map \Ei x = \ds \int_{t \mathop = x}^{t \mathop \to +\infty} \frac {e^{-t} } t \rd t$

Then:

where $\gamma$ denotes the Euler-Mascheroni constant.