# Prove the equation

Prove that $$\int_0^{\infty}\exp\left(-\left(x^2+\dfrac{a^2}{x^2}\right)\right)\text{d}x=\frac{e^{-2a}\sqrt{\pi}}{2}$$ Assume that the equation is true for $a=0.$

-

$$I(a):=\int_0^{\infty}e^{-\left(x^2+\dfrac{a^2}{x^2}\right)}\text{d}x$$ $$\frac{dI}{da}=\int_0^{\infty}e^{-\left(x^2+\dfrac{a^2}{x^2}\right)}\left(-\frac{2a}{x^2}\right)\text{d}x$$ Now substitute $y=\frac{a}{x}$, so $dy=-\frac{a}{y^2}$: $$\frac{dI}{da}=2\int_{\infty}^{0}e^{-\left(\dfrac{a^2}{y^2}+y^2\right)}\text{d}y=-2\int_{0}^{\infty}e^{-\left(\dfrac{a^2}{y^2}+y^2\right)}\text{d}y=-2I$$ To obtain $I$ you just have to solve the simple ODE: $$\frac{dI}{da}=-2I$$ with initial condition given by $$I(0)=\frac{\sqrt{\pi}}{2}\ .$$ This gives you $$I(a)=\frac{e^{-2a}\sqrt{\pi}}{2}\ .$$

-
+1. This is a nice answer. –  Tunk-Fey Jun 24 at 13:58
Thanks. It is very clever. But can you elaborate a bit more on how you came up with the solution translating the original problem into solving an ODE? –  lovelesswang Jun 24 at 14:06
Do you mean the explicit resolution of the ODE? –  Dario Jun 24 at 14:08
The idea behind the solution is the technique of differentiation under the integral sign(en.wikipedia.org/wiki/…): when you have an integral that depend on a parameter, you can derive the integral w.r.t. that parameter and sometimes you obtain a much simpler integral. In this particular case the trick is to notice that after the derivation, getting rid of the extra factor you obtain by a change of variable you get again the first integral. This allow you to find a relation between $I$ and its derivative, i.e. the ODE. –  Dario Jun 24 at 14:13
show 1 more comment

In general \begin{align} \int_{x=0}^\infty \exp\left(-ax^2-\frac{b}{x^2}\right)\,dx&=\int_{x=0}^\infty \exp\left(-a\left(x^2+\frac{b}{ax^2}\right)\right)\,dx\\ &=\int_{x=0}^\infty \exp\left(-a\left(x^2-2\sqrt{\frac{b}{a}}+\frac{b}{ax^2}+2\sqrt{\frac{b}{a}}\right)\right)\,dx\\ &=\int_{x=0}^\infty \exp\left(-a\left(x-\frac{1}{x}\sqrt{\frac{b}{a}}\right)^2-2\sqrt{ab}\right)\,dx\\ &=\exp(-2\sqrt{ab})\int_{x=0}^\infty \exp\left(-a\left(x-\frac{1}{x}\sqrt{\frac{b}{a}}\right)^2\right)\,dx\\ \end{align} The trick to solve the last integral is by setting $$I=\int_{x=0}^\infty \exp\left(-a\left(x-\frac{1}{x}\sqrt{\frac{b}{a}}\right)^2\right)\,dx.$$ Let $t=-\frac{1}{x}\sqrt{\frac{b}{a}}\;\rightarrow\;x=-\frac{1}{t}\sqrt{\frac{b}{a}}\;\rightarrow\;dx=\frac{1}{t^2}\sqrt{\frac{b}{a}}\,dt$, then $$I_t=\sqrt{\frac{b}{a}}\int_{t=0}^\infty \frac{\exp\left(-a\left(-\frac{1}{t}\sqrt{\frac{b}{a}}+t\right)^2\right)}{t^2}\,dt.$$ Let $t=x\;\rightarrow\;dt=dx$, then $$I_t=\int_{t=0}^\infty \exp\left(-a\left(t-\frac{1}{t}\sqrt{\frac{b}{a}}\right)^2\right)\,dt.$$ Adding the two $I_t$s yields $$2I=I_t+I_t=\int_{t=0}^\infty\left(1+\frac{1}{t^2}\sqrt{\frac{b}{a}}\right)\exp\left(-a\left(t-\frac{1}{t}\sqrt{\frac{b}{a}}\right)^2\right)\,dt.$$ Let $s=t-\frac{1}{t}\sqrt{\frac{b}{a}}\;\rightarrow\;ds=\left(1+\frac{1}{t^2}\sqrt{\frac{b}{a}}\right)dt$ and for $0<t<\infty$ is corresponding to $-\infty<s<\infty$, then $$I=\frac{1}{2}\int_{s=-\infty}^\infty e^{-as^2}\,ds=\frac{1}{2}\sqrt{\frac{\pi}{a}},$$ where $I$ is a Gaussian integral. Thus \begin{align} \exp(-2\sqrt{ab})\int_{x=0}^\infty \exp\left(-a\left(x-\frac{1}{x}\sqrt{\frac{b}{a}}\right)^2\right)\,dx &=\large\color{blue}{\frac12\sqrt{\frac{\pi}{a}}e^{-2\sqrt{ab}}}. \end{align} In our case, put $a=1$ and $b=a^2$.

-