Supongamos que la integral se $I=\int_0^{\infty} e^{-y^2-\frac{x^2}{y^2}}dy$. A continuación, tomamos nota de que $y^2+\frac{x^2}{y^2}=\left(y-\frac{|x|}{y}\right)^2+2|x|$.
Por lo tanto, hemos
$$I=e^{-2|x|}\int_0^{\infty}e^{-\left(y-\frac{|x|}{y}\right)^2}dy \tag 1$$
Ahora, suplente $y\to |x|/y$, de modo que $dy\to -|x|dy/y^2$. A continuación,
$$I=e^{-2|x|}\int_0^{\infty}\frac{|x|}{y^2}e^{-\left(y-\frac{|x|}{y}\right)^2}dy \tag 2$$
Si añadimos $(1)$ y (2), nos encontramos con
$$\begin{align}
I&=\frac12\,e^{-2|x|}\int_0^{\infty}\left(1+\frac{|x|}{y^2}\right)e^{-\left(y-\frac{|x|}{y}\right)^2}dy \\\\
&=\frac12\,e^{-2|x|}\int_{-\infty}^{\infty}e^{-y^2}dy\\\\
&=e^{-2|x|}\frac{\sqrt{\pi}}{2}
\end{align}$$
Así que, aunque no bastante "Feynmann" truco, es una manera eficaz de evaluación.