Definir $I(s)$ por
$$ I(s) = \int_{0}^{\infty} \frac{e^{-su^2}}{1+u^2} \, du. $$
A continuación, $I(s)$ resuelve la siguiente ecuación:
$$ I(s) - I'(s)
= \int_{0}^{\infty} e^{-do^2} \, du
= \frac{1}{2}\sqrt{\frac{\pi}{s}}. $$
Este es un 1st-orden de la ecuación diferencial, que puede ser resuelto de manera sistemática por medio de la integración del factor. El resultado es que
$$ I(s) = e^s \left( \mathsf{C}-\int \frac{1}{2}\sqrt{\frac{\pi}{s}}e^{-s} \, ds \right) $$
para algunos la elección adecuada de los constantes $\mathsf{C}$. Junto con la condición de contorno $I(\infty) = 0$, resulta que
$$ I(s)
= e^s \int_{s}^{\infty} \frac{1}{2}\sqrt{\frac{\pi}{s}}e^{-s'} \, ds'
= \frac{\pi}{2}e^s \operatorname{erfc}(\sqrt{s}). $$
Conectar $s = 1$ da el valor
$$ I(1) = \frac{e\pi}{2}\operatorname{erfc}(1) \approx 0.67164671082336758522\cdots. $$