Paso 1. Deje $s > 0$$\beta \in \mathbb{C}$. A continuación, utilizando la integral de gauss $\int_{-\infty}^{\infty} e^{-sx^2} \, dx = \sqrt{\pi/s}$, tenemos
\begin{align*}
\int_{-\infty}^{\infty} e^{-s(x-\beta)^2} \, dx
&= \int_{-\infty}^{\infty} \bigg( e^{-sx^2} + \int_{0}^{1} \overbrace{ 2s\beta(x-\beta t) e^{-s(x-\beta t)^2} }^{= \frac{\partial}{\partial t} e^{-s(x-\beta t)^2}} \, dt \bigg) \, dx \\
&= \sqrt{\frac{\pi}{s}} -\beta \int_{0}^{1} \int_{-\infty}^{\infty} (-2s)(x-\beta t) e^{-s(x-\beta t)^2} \, dx dt \quad {\small(\because\text{Fubini})}\\
&= \sqrt{\frac{\pi}{s}} -\beta \int_{0}^{1} \left[ e^{-s(x-\beta t)^2} \right]_{x=-\infty}^{x=\infty} \, dt
= \bbox[border:1px dashed green,6px]{ \sqrt{\frac{\pi}{s}} }.
\end{align*}
Luego de conectar $\beta=\pm\mathrm{i}/2s$, podemos comprobar que
$$ \int_{-\infty}^{\infty} e^{-sx^2}\cos(x) \, dx = \sqrt{\frac{\pi}{s}} e^{-1/4s}. $$
Paso 2. Mediante el paso anterior,
\begin{align*}
I := \int_{-\infty}^{\infty} \frac{\cos x}{x^2+1} \, dx
&= \int_{-\infty}^{\infty} \cos x \left( \int_{0}^{\infty} e^{-(x^2+1)s} \, ds \right) \, dx \\
&= \int_{0}^{\infty} \left( \int_{-\infty}^{\infty} e^{-sx^2}\cos(x) \, dx \right) e^{-s} \, ds \quad {\small(\because\text{Fubini})}\\
&= \int_{0}^{\infty} \sqrt{\frac{\pi}{s}} e^{-\left( s + \frac{1}{4s}\right)} \, ds \\
&= \int_{0}^{\infty} \sqrt{2\pi} e^{-\frac{1}{2}\left( t^2 + \frac{1}{t^2}\right)} \, dt. \quad {\small(2s=t^2)}
\end{align*}
Finalmente aquí está un muy hábil manera de calcular la última integral. La aplicación de la sustitución de $t\mapsto 1/t$ muestra que
$$ \int_{0}^{\infty} e^{-\frac{1}{2}\left( t^2 + \frac{1}{t^2}\right)} \, dt = \int_{0}^{\infty} \frac{1}{t^2} e^{-\frac{1}{2}\left( t^2 + \frac{1}{t^2}\right)} \, dt. $$
Así que el promedio,
$$
I = \sqrt{\frac{\pi}{2}} \int_{0}^{\infty} \left(1 + \frac{1}{t^2}\right) e^{-\frac{1}{2}\left( t - \frac{1}{t}\right)^2 - 1} \, dt.
$$
Por último, la aplicación de la sustitución de $u = t - \frac{1}{t}$ demuestra
$$
Yo
= \sqrt{\frac{\pi}{2}} \int_{-\infty}^{\infty} e^{-\frac{u^2}{2} - 1} \, du
= \frac{\pi}{e}. $$