# Evaluate $\int_0^\infty \frac{(\mathrm{e}^{-ax}-\mathrm{e}^{-bx})^2}{x^2}dx=\ln\frac{(2a)^{2a}(2b)^{2b}}{(a+b)^{2(a+b)}}$

(1) Prove that$$\int_0^\infty \frac{(\mathrm{e}^{-ax}-\mathrm{e}^{-bx})^2}{x^2}\mathrm{d}x=\ln\frac{(2a)^{2a}(2b)^{2b}}{(a+b)^{2(a+b)}}$$ (2) Prove that$$\int_0^\infty \frac{\mathrm{e}^{-ax^2}-\mathrm{e}^{-bx^2}}{x^2}\mathrm{d}x=\frac{1}{2}\ln\frac{b}{a}$$

I just learned improper integral with parameters and worked out a basic example below $$\int_0^\infty \frac{\mathrm{e}^{-ax}-\mathrm{e}^{-bx}}{x}\mathrm{d}x=\ln\frac{b}{a}$$ But I do not know how to compute the two integrals above. I think the basic example may be of some help and tried integration by parts, inspired by $$\int_0^\infty\frac{\sin x}{x}\mathrm{d}x$$ and $$\int_0^\infty\frac{\sin^2 x}{x^2}\mathrm{d}x$$.

Any hint is appreciated.

• Differentiating under the integral sign would be definitely helpful here. Assume the integral to be $f(a,b)$ and then partially differentiate wrt one variable, then again integrate it. It's the easiest and fastest method!
– V.G
Dec 11, 2020 at 18:15

Write $$\frac{1}{x^2}=\int_0^\infty ye^{-xy}dy$$ so the LHS of (1) is\begin{align}&\int_0^\infty dy\:y\int_0^\infty dx(e^{-(2a+y)x}-2e^{-(a+b+y)x}+e^{-(2b+y)x})\\&=\int_0^\infty dy\:y\left(\frac{1}{2a+y}-\frac{2}{a+b+y}+\frac{1}{2b+y}\right).\end{align}You can handle (2) the same way, thought it helps to first replace $$\int_0^\infty$$ with $$\tfrac12\int_{-\infty}^\infty$$, so you just have to evaluate Gaussians on $$\Bbb R$$.

• Very nice! idea! Dec 11, 2020 at 12:19

Integrate (1) by parts with $$\frac{dx}{x^2}=-d(\frac1x)$$, \begin{align} &\int_0^\infty \frac{(\mathrm{e}^{-ax}-\mathrm{e}^{-bx})^2}{x^2}dx\\ = &2a \int_0^\infty \frac{{e}^{-(a+b)x}-{e}^{-2ax}}{x}dx + 2b \int_0^\infty \frac{{e}^{-(a+b)x}-{e}^{-2bx}}{x}dx\\ = &2a \ln \frac{2a}{a+b} + 2b \ln \frac{2b}{a+b} =\ln\frac{(2a)^{2a}(2b)^{2b}}{(a+b)^{2(a+b)}} \end{align} where $$\int_0^\infty \frac{\mathrm{e}^{-px }-\mathrm{e}^{-qx }}{x}\mathrm{d}x=\ln\frac{q}{p}$$ is used.

Integrate (2) by parts $$\int_0^\infty \frac{{e}^{-ax^2}-{e}^{-bx^2}}{x^2}{d}x= 2\int_0^\infty (b{e}^{-bx^2}-a{e}^{-ax^2})dx= \sqrt{\pi}(\sqrt b-\sqrt a)$$ where $$\int_0^\infty e^{-px^2} =\frac{\sqrt{\pi}}{2\sqrt p}$$ is used

As fjaclot says, $$(2)$$ is wrong. Making the substitution $$u = x\sqrt{a}$$ yields $$\sqrt{a}\int_{0}^{\infty}\frac{e^{-u^{2}}-e^{-tu^{2}}}{u^{2}}du$$ where $$t = \frac{b}{a}$$. Let $$f(t)$$ be the integral. Then $$f'(t) = \int_0^{\infty} e^{-tu^2} du = \frac{\sqrt{\pi}}{2\sqrt{t}}$$

Since $$f(1) = 0$$, $$(2)$$ is $$\sqrt{a}\left( \sqrt{\pi t}-\sqrt{\pi} \right) = \sqrt{\pi}\left(\sqrt{b}-\sqrt{a}\right)$$

Similarly for $$(1)$$, make the substitution $$u=ax$$ and let $$f(t) = \int_{0}^{\infty}\frac{\left(e^{-u}-e^{-tu}\right)^{2}}{u^{2}}du$$. Then $$f'(t) = 2\int_{0}^{\infty}\frac{e^{-\left(1+t\right)u}-e^{-2tu}}{u}du$$

Then you can apply that formula you found so that $$f'(t) = 2 \ln\left( \frac{2t}{1+t} \right) = 2\ln\left(2t\right)-2\ln\left(1+t\right)$$

The integral is then $$f(t) = 2t\ln\left( \frac{2t}{t+1} \right)-2\ln(t+1)+2\ln(2)$$. Therefore the answer for $$(1)$$ is $$af\left(\frac{b}{a}\right)$$ which simplifies to $$\ln\left(\frac{(2a)^{2a}(2b)^{2b}}{(a+b)^{2(a+b)}}\right)$$

The answer for integral (2) is not correct. It should be :

     (Pi)^(1/2) [ (b)^(1/2) - (a)^(1/2) ]


fjaclot

$$I=\int_0^\infty\frac{(e^{-ax}-e^{-bx})^2}{x^2}dx=\int_0^\infty(e^{-ax}-e^{-bx})^2d\frac{1}{x}=\frac{e^{-ax}-e^{-bx}}{x}\Big|_0^\infty-\int_0^\infty\frac{2(e^{-ax}-e^{-bx})(-ae^{-ax}+be^{-bx})}{x}dx$$

$$=-\int_0^\infty\frac{2(e^{-ax}-e^{-bx})(-ae^{-ax}+be^{-bx})}{x}dx =-2\int_0^\infty\frac{1}{x}(-ae^{-2ax}+ae^{-(a+b)x}+be^{-(a+b)x}-be^{-2bx})dx$$

$$=2(a\int_0^\infty\frac{1}{x}(e^{-(a+b)x}-e^{-2ax})dx+b\int_0^\infty\frac{1}{x}(e^{-2bx}-e^{-(a+b)x})dx$$

by Frullani integral

$$\int_0^\infty\frac{1}{x}(e^{-(a+b)x}-e^{-2ax})dx=\ln(\frac{2a}{a+b})$$

$$\int_0^\infty\frac{1}{x}(e^{-2bx}-e^{-(a+b)x})=\ln(\frac{a+b}{2b})$$

so $$I=2(a\ln(\frac{2a}{a+b})-b\ln(\frac{a+b}{2b}))=\ln\frac{(2a)^{2a}(2b)^{2b}}{(a+b)^{a+b}}$$

• How is this different from Quanto's answer below? (You have some typos as well.)
– Gary
Jan 6 at 9:34