๐Ÿงฎ Brain Teaser

Lobachevsky's Folding Trick

Evaluate

I=โˆซ0โˆžsinโกxxโ‹…11+cosโก2xโ€‰dx.I=\int_0^\infty \frac{\sin x}{x}\cdot\frac{1}{1+\cos^2 x}\,dx .

Bonus. Show more generally that if ff is continuous, ฯ€\pi-periodic, and satisfies f(ฯ€โˆ’x)=f(x)f(\pi-x)=f(x), then

โˆซ0โˆžsinโกxxโ€‰f(x)โ€‰dx=โˆซ0ฯ€/2f(x)โ€‰dx.\int_0^\infty \frac{\sin x}{x}\,f(x)\,dx=\int_0^{\pi/2} f(x)\,dx .

As a check, fโ‰ก1f\equiv 1 gives the Dirichlet integral ฯ€/2\pi/2.

Hint. You may use the partial fraction expansion of the cosecant:

1sinโกu=โˆ‘mโˆˆZ(โˆ’1)mu+mฯ€,\frac{1}{\sin u}=\sum_{m\in\mathbb Z}\frac{(-1)^m}{u+m\pi},

where the sum is taken symmetrically.

integration beeLobachevsky integral formulaperiodic foldingpartial fractions of cosecantMittag-Lefflerconditionally convergent integral

Answer: Lobachevsky's Folding Trick

Key Idea / Intuition

The factor f(x)f(x) repeats every ฯ€\pi. The only non-periodic part is sinโกxx\dfrac{\sin x}{x}, and sinโกx\sin x itself just flips sign on each period.

So cut [0,โˆž)[0,\infty) into intervals of length ฯ€\pi and stack them all on top of [0,ฯ€][0,\pi]. Then use the mirror symmetry f(ฯ€โˆ’x)=f(x)f(\pi-x)=f(x) to fold [0,ฯ€][0,\pi] onto [0,ฯ€/2][0,\pi/2].

After stacking and folding, the factor 1x\dfrac1x becomes an alternating sum over all its shifted copies 1u+mฯ€\dfrac{1}{u+m\pi}. By Euler's partial fractions, that sum is exactly 1sinโกu\dfrac{1}{\sin u}. It cancels the sinโกu\sin u, and what remains is simply โˆซ0ฯ€/2f\int_0^{\pi/2} f.

The infinite oscillatory integral collapses to a finite, elementary one.


Formal Proof / Solution

Step 1 (Convergence). The function g(x)=sinโกxโ€‰f(x)g(x)=\sin x\, f(x) satisfies g(x+ฯ€)=โˆ’g(x)g(x+\pi)=-g(x). So gg has period 2ฯ€2\pi and zero mean over a period, which means its antiderivative is bounded. Since 1/xโ†“01/x\downarrow 0, Dirichlet's test shows that โˆซ1โˆžg(x)/xโ€‰dx\int_1^\infty g(x)/x\,dx converges. Near 00 the integrand is bounded. Hence

I=limโกNโ†’โˆžโˆซ0Nฯ€sinโกxxf(x)โ€‰dx.I=\lim_{N\to\infty}\int_0^{N\pi}\frac{\sin x}{x}f(x)\,dx .

Step 2 (Stack the periods). On [kฯ€,(k+1)ฯ€][k\pi,(k+1)\pi], substitute x=t+kฯ€x=t+k\pi. This uses sinโก(t+kฯ€)=(โˆ’1)ksinโกt\sin(t+k\pi)=(-1)^k\sin t and f(t+kฯ€)=f(t)f(t+k\pi)=f(t):

โˆซkฯ€(k+1)ฯ€sinโกxxf(x)โ€‰dx=(โˆ’1)kโˆซ0ฯ€sinโกtโ€‰f(t)t+kฯ€โ€‰dt.\int_{k\pi}^{(k+1)\pi}\frac{\sin x}{x}f(x)\,dx=(-1)^k\int_0^{\pi}\frac{\sin t\, f(t)}{t+k\pi}\,dt .

Step 3 (Fold the period in half). Split [0,ฯ€][0,\pi] into [0,ฯ€/2]โˆช[ฯ€/2,ฯ€][0,\pi/2]\cup[\pi/2,\pi]. On the second half put t=ฯ€โˆ’ut=\pi-u, which uses sinโก(ฯ€โˆ’u)=sinโกu\sin(\pi-u)=\sin u and f(ฯ€โˆ’u)=f(u)f(\pi-u)=f(u). The kk-th piece becomes

(โˆ’1)kโˆซ0ฯ€/2sinโกuโ€‰f(u)[1u+kฯ€+1(k+1)ฯ€โˆ’u]du.(-1)^k\int_0^{\pi/2}\sin u\, f(u)\left[\frac{1}{u+k\pi}+\frac{1}{(k+1)\pi-u}\right]du .

Step 4 (Sum over kk). Set ak(u)=1u+kฯ€+1(k+1)ฯ€โˆ’ua_k(u)=\dfrac{1}{u+k\pi}+\dfrac{1}{(k+1)\pi-u}. For uโˆˆ(0,ฯ€/2]u\in(0,\pi/2] we have ak(u)>0a_k(u)>0, ak(u)a_k(u) decreasing in kk, and ak(u)โ†’0a_k(u)\to0.

So โˆ‘k(โˆ’1)kak(u)\sum_k(-1)^k a_k(u) is a Leibniz series. Its partial sums lie between 00 and a0(u)=1u+1ฯ€โˆ’ua_0(u)=\frac1u+\frac1{\pi-u}. Therefore

โˆฃsinโกuโ€‰f(u)โˆ‘k=0Nโˆ’1(โˆ’1)kak(u)โˆฃโ‰คโˆฅfโˆฅโˆžโ€‰sinโกu(1u+1ฯ€โˆ’u)โ‰คC.\Big|\sin u\, f(u)\sum_{k=0}^{N-1}(-1)^k a_k(u)\Big|\le \|f\|_\infty\,\sin u\left(\frac1u+\frac{1}{\pi-u}\right)\le C .

This bound is integrable on [0,ฯ€/2][0,\pi/2], so dominated convergence lets us swap the sum and the integral:

I=โˆซ0ฯ€/2sinโกuโ€‰f(u)โ€‰S(u)โ€‰du,S(u)=โˆ‘kโ‰ฅ0(โˆ’1)k[1u+kฯ€+1(k+1)ฯ€โˆ’u].I=\int_0^{\pi/2}\sin u\, f(u)\,S(u)\,du,\qquad S(u)=\sum_{k\ge0}(-1)^k\left[\frac{1}{u+k\pi}+\frac{1}{(k+1)\pi-u}\right].

Step 5 (Recognize the cosecant). In the second sum, reindex with m=โˆ’(k+1)m=-(k+1), so mโ‰คโˆ’1m\le -1. Then

(โˆ’1)k(k+1)ฯ€โˆ’u=โˆ’(โˆ’1)mโˆ’(u+mฯ€)=(โˆ’1)mu+mฯ€.\frac{(-1)^k}{(k+1)\pi-u}=\frac{-(-1)^m}{-(u+m\pi)}=\frac{(-1)^m}{u+m\pi}.

The terms with kโ‰ฅ0k\ge0 (that is, mโ‰ฅ0m\ge0) and with mโ‰คโˆ’1m\le-1 combine into the symmetric full sum. By Euler's expansion ฯ€sinโกฯ€z=โˆ‘m(โˆ’1)mzโˆ’m\dfrac{\pi}{\sin \pi z}=\sum_m\dfrac{(-1)^m}{z-m} with z=u/ฯ€z=u/\pi,

S(u)=โˆ‘mโˆˆZ(โˆ’1)mu+mฯ€=1sinโกu.S(u)=\sum_{m\in\mathbb Z}\frac{(-1)^m}{u+m\pi}=\frac{1}{\sin u}.

So the sinโกu\sin u cancels, and we obtain Lobachevsky's formula:

โˆซ0โˆžsinโกxxf(x)โ€‰dx=โˆซ0ฯ€/2f(u)โ€‰du.\int_0^\infty\frac{\sin x}{x}f(x)\,dx=\int_0^{\pi/2}f(u)\,du .

Step 6 (Apply to our ff). The function f(x)=11+cosโก2xf(x)=\dfrac{1}{1+\cos^2x} is continuous and ฯ€\pi-periodic, and cosโก2(ฯ€โˆ’x)=cosโก2x\cos^2(\pi-x)=\cos^2x, so the formula applies. With v=tanโกuv=\tan u, dv=secโก2uโ€‰dudv=\sec^2u\,du, and 1+cosโก2u=secโก2u+1secโก2u=v2+2secโก2u1+\cos^2u=\dfrac{\sec^2u+1}{\sec^2 u}=\dfrac{v^2+2}{\sec^2 u}:

โˆซ0ฯ€/2du1+cosโก2u=โˆซ0โˆždvv2+2=12โ‹…ฯ€2.\int_0^{\pi/2}\frac{du}{1+\cos^2u}=\int_0^{\infty}\frac{dv}{v^2+2}=\frac{1}{\sqrt2}\cdot\frac{\pi}{2}. โ€‰I=ฯ€22โ€‰\boxed{\,I=\dfrac{\pi}{2\sqrt2}\,}

Remarks.

  • With fโ‰ก1f\equiv1 this recovers โˆซ0โˆžsinโกxxโ€‰dx=ฯ€2\int_0^\infty \frac{\sin x}{x}\,dx=\frac{\pi}{2}.
  • With f(x)=sinโก2xx2โ‹…x2sinโก2xf(x)=\dfrac{\sin^2x}{x^2}\cdot\dfrac{x^2}{\sin^2 x} the formula does not apply directly. However, a variant of the same folding gives โˆซ0โˆžsinโก2xx2f(x)โ€‰dx=โˆซ0ฯ€/2f\int_0^\infty\frac{\sin^2x}{x^2}f(x)\,dx=\int_0^{\pi/2}f as well. This variant uses โˆ‘m1(u+mฯ€)2=cscโก2u\sum_m\frac{1}{(u+m\pi)^2}=\csc^2u.
  • The structural lesson: periodic times 1/x1/x can be summed by folding, and the fold turns 1/x1/x into a Mittag-Leffler series.
Type: IntegrationEdit on GitHub โ†—