Derivation of an Integral Expression of the Euler-Mascheroni Constant

\begin{equation}
\int\limits_{0}^{\infty} \mathrm{e}^{-x} \mathrm{ln}(x) \mathrm{d} x = -\gamma
\label{eq:1608111}
\tag{1}
\end{equation} There are of course an enormous number of integral expressions of the Euler-Mascheroni Constant. This particular expression appeared on the cover of Irresistible Integrals: Symbolics, Analysis and Experiments in the Evaluation of Integrals by George Boros and Victor Moll and was derived in section 9.3 of the book. Here, I will show this derivation and fill in the steps as well as parts that were not shown in the book.

Let us begin by providing a basic definition of the Euler-Mascheroni Constant
\begin{equation}
\gamma := \lim_{n \to \infty} H_{n} – \mathrm{ln}(n)
\label{eq:1608112}
\tag{2}
\end{equation}
where $$H_{n}$$ are the harmonic numbers defined as
\begin{equation}
H_{n} = \displaystyle\sum_{k=1}^{n} \frac{1}{k}
\label{eq:1608113}
\tag{3}
\end{equation}

We will break up the derivation into parts to make it easier to follow.

1. Exercise 9.3.1 – Prove
\begin{equation}
\int\limits_{0}^{1} \frac{1-(1-x)^{n}}{x} \mathrm{d} x = H_{n}
\label{eq:1608114}
\tag{4}
\end{equation}

\begin{align}
\label{eq:1608115}
\tag{5}
\int\limits_{0}^{1} \frac{1-(1-x)^{n}}{x} \mathrm{d} x & = \int\limits_{0}^{1} \frac{1-y^{n}}{1-y} \mathrm{d} y \\
& = \int\limits_{0}^{1} \frac{1}{1-y} \mathrm{d} y \, – \int\limits_{0}^{1} \frac{y^{n}}{1-y} \mathrm{d} y \\
& = \int\limits_{0}^{1} \displaystyle\sum_{k=0}^{\infty} y^{k} \mathrm{d} y \,\, – \int\limits_{0}^{1} y^{n} \displaystyle\sum_{k=0}^{\infty} y^{k} \mathrm{d} y \\
& = \displaystyle\sum_{k=0}^{\infty} \frac{y^{k+1}}{k+1} \big|_{0}^{1} \,\, – \displaystyle\sum_{k=0}^{\infty} \frac{y^{k+n+1}}{k+n+1} \big|_{0}^{1} \\
& = (1 + \frac{1}{2} + \frac{1}{3} + \dots) \, – (\frac{1}{n+1} + \frac{1}{n+2} +\frac{1}{n+3} + \dots) \\
& = 1 + \frac{1}{2} + \frac{1}{3} + \dots + \frac{1}{n} = H_{n}
\end{align}

Notes:

• $$|x| < 1$$
• Let $$y=1-x$$
• Expand $$1/(1-y)$$ as a geometric series \begin{equation} \frac{1}{1-y} = \displaystyle\sum_{k=0}^{\infty} y^{k} \end{equation}
• The integrals are improper, we have dropped the limits for clarity.
• We assumed that the integrals converged uniformly so that we could integrate term by term. We acknowledge that this should be proven, but the assumption is warranted based on the form of the integrand and experience.

2. Proposition 9.3.1 – The Euler constant is given by
\begin{equation}
\gamma = \int\limits_{0}^{1} \frac{1-\mathrm{e}^{-x}}{x} \mathrm{d} x \,\, – \int\limits_{1}^{\infty} \frac{\mathrm{e}^{-x}}{x} \mathrm{d} x
\label{eq:1608116}
\tag{6}
\end{equation}
The proof in the book is fairly straightforward but there are some missing steps.

Make the substitution $$x = \frac{y}{n}$$ in \eqref{eq:1608114} to obtain
\begin{align}
H_{n} & = \int\limits_{0}^{n} \frac{1-(1-\frac{y}{n})^{n}}{y} \mathrm{d} y \\
& = \int\limits_{0}^{1} \frac{1-(1-\frac{y}{n})^{n}}{y} \mathrm{d} y \,\, + \int\limits_{1}^{n} \frac{1-(1-\frac{y}{n})^{n}}{y} \mathrm{d} y \\
& = \int\limits_{0}^{1} \frac{1-(1-\frac{y}{n})^{n}}{y} \mathrm{d} y \,\, + \int\limits_{1}^{n} \frac{1}{y} \mathrm{d} y \,\, – \int\limits_{1}^{n} \frac{(1-\frac{y}{n})^{n}}{y} \mathrm{d} y \\
& = \int\limits_{0}^{1} \frac{1-(1-\frac{y}{n})^{n}}{y} \mathrm{d} y \,\, + \mathrm{ln}(n) \, – \int\limits_{1}^{n} \frac{(1-\frac{y}{n})^{n}}{y} \mathrm{d} y
\label{eq:1608117}
\tag{7}
\end{align}

Applying $$\lim n \to \infty$$ to equation \eqref{eq:1608117}, noting the limit definition of
\begin{equation}
\mathrm{e}^{\pm x} = \lim_{n \to \infty} \left(1 \pm \frac{x}{n} \right)^{n}
\end{equation}
and invoking equation \eqref{eq:1608112} yields our result.

3. Proposition 9.3.2 – Euler’s constant is given by equation \eqref{eq:1608111}.

Here we require one obvious alteration of the integral in light of \eqref{eq:1608117} and a non obvious trick.

\begin{equation}
\int\limits_{0}^{\infty} \mathrm{e}^{-x} \mathrm{ln}(x) \mathrm{d} x = \int\limits_{0}^{1} \mathrm{e}^{-x} \mathrm{ln}(x) \mathrm{d} x + \int\limits_{1}^{\infty} \mathrm{e}^{-x} \mathrm{ln}(x) \mathrm{d} x = \mathrm{I}_{1} + \mathrm{I}_{2}
\label{eq:1608118}
\tag{8}
\end{equation}

Integrating $$\mathrm{I}_{2}$$ by parts and applying the appropriate limit yields
\begin{equation}
\mathrm{I}_{2} = \int\limits_{1}^{\infty} \frac{\mathrm{e}^{-x}}{x} \mathrm{d} x
\label{eq:1608119}
\tag{9}
\end{equation}

For $$\mathrm{I}_{1}$$ use the substitution
\begin{equation}
\mathrm{e}^{-x} = -\frac{\mathrm d}{\mathrm d x} (\mathrm{e}^{-x} – 1),
\end{equation}
integrate by parts and apply the appropriate limit to obtain
\begin{equation}
\mathrm{I}_{1} = -\int\limits_{0}^{1} \frac{\mathrm{e}^{-x} – 1}{x} \mathrm{d} x
\label{eq:16081110}
\tag{10}
\end{equation}
Substituting equations \eqref{eq:16081110} and \eqref{eq:1608119} into \eqref{eq:1608118} and using \eqref{eq:1608116} yields our final result.