๐Ÿงฎ Brain Teaser

The Fresnel-Flavored Integral

Compute the definite integral

I=โˆซ0โˆžsinโกxxโ€‰dx.I = \int_0^\infty \frac{\sin x}{\sqrt{x}}\, dx.

(You may use the known Gaussian integral โˆซ0โˆžeโˆ’t2โ€‰dt=ฯ€2\int_0^\infty e^{-t^2}\,dt = \dfrac{\sqrt{\pi}}{2}.)

Gaussian integralLaplace transformdouble integralFubiniFresnel

Answer: Fresnel-Flavored Integral

Key Idea / Intuition

The trick is to get rid of the x\sqrt{x} in the denominator by writing 1x\frac{1}{\sqrt{x}} as a Gaussian integral: 1x=2ฯ€โˆซ0โˆžeโˆ’t2xโ€‰dt\frac{1}{\sqrt{x}} = \frac{2}{\sqrt{\pi}}\int_0^\infty e^{-t^2 x}\,dt. This converts the problem into a double integral where the xx-integral becomes a standard Laplace-transform type integral of sinโกx\sin x, which is elementary. The two integrals then decouple beautifully.


Formal Proof / Solution

Step 1: Represent 1/x1/\sqrt{x} as a Gaussian.

From the Gaussian integral, substituting u=txu = t\sqrt{x}: โˆซ0โˆžeโˆ’t2xโ€‰dt=ฯ€2x,\int_0^\infty e^{-t^2 x}\,dt = \frac{\sqrt{\pi}}{2\sqrt{x}}, so 1x=2ฯ€โˆซ0โˆžeโˆ’t2xโ€‰dt.\frac{1}{\sqrt{x}} = \frac{2}{\sqrt{\pi}}\int_0^\infty e^{-t^2 x}\,dt.

Step 2: Substitute into II.

I=โˆซ0โˆžsinโกxโ‹…2ฯ€โˆซ0โˆžeโˆ’t2xโ€‰dtโ€…โ€Šdx=2ฯ€โˆซ0โˆžโˆซ0โˆžeโˆ’t2xsinโกxโ€…โ€Šdxโ€…โ€Šdt.I = \int_0^\infty \sin x \cdot \frac{2}{\sqrt{\pi}}\int_0^\infty e^{-t^2 x}\,dt\;dx = \frac{2}{\sqrt{\pi}}\int_0^\infty \int_0^\infty e^{-t^2 x}\sin x\;dx\;dt.

(Fubini is justified since the double integral converges absolutely after a small regularization argument.)

Step 3: Evaluate the inner xx-integral.

For fixed t>0t > 0, use the standard Laplace transform: โˆซ0โˆžeโˆ’axsinโกxโ€…โ€Šdx=11+a2,a>0.\int_0^\infty e^{-ax}\sin x\;dx = \frac{1}{1+a^2}, \quad a > 0.

With a=t2a = t^2: โˆซ0โˆžeโˆ’t2xsinโกxโ€…โ€Šdx=11+t4.\int_0^\infty e^{-t^2 x}\sin x\;dx = \frac{1}{1+t^4}.

Step 4: Evaluate the outer tt-integral.

I=2ฯ€โˆซ0โˆždt1+t4.I = \frac{2}{\sqrt{\pi}}\int_0^\infty \frac{dt}{1+t^4}.

Now compute J=โˆซ0โˆždt1+t4\displaystyle J = \int_0^\infty \frac{dt}{1+t^4}. Factor 1+t4=(t2+2โ€‰t+1)(t2โˆ’2โ€‰t+1)1+t^4 = (t^2+\sqrt{2}\,t+1)(t^2-\sqrt{2}\,t+1) and use partial fractions, or use the residue theorem. The standard result is: J=ฯ€22.J = \frac{\pi}{2\sqrt{2}}.

(Quick derivation: by the substitution tโ†ฆ1/tt \mapsto 1/t and symmetry, one can show J=ฯ€22J = \frac{\pi}{2\sqrt 2} via residues at eiฯ€/4e^{i\pi/4} and e3iฯ€/4e^{3i\pi/4}.)

Step 5: Assemble.

I=2ฯ€โ‹…ฯ€22=ฯ€2ฯ€=ฯ€2.I = \frac{2}{\sqrt{\pi}}\cdot\frac{\pi}{2\sqrt{2}} = \frac{\pi}{\sqrt{2\pi}} = \sqrt{\frac{\pi}{2}}.

I=ฯ€2.\boxed{I = \sqrt{\dfrac{\pi}{2}}.}

Why this is beautiful: Three seemingly unrelated objects โ€” the Gaussian, Laplace transforms of sinโก\sin, and a rational integral โ€” combine to give the clean answer ฯ€/2\sqrt{\pi/2}, the same constant that appears in the normal distribution.

Type: IntegrationEdit on GitHub โ†—