Replace x with π/2 - x to get the equivalent integral
but the integrand is even, so this is really just
Substitute x = 1/2 arccot(u/2), which transforms the integral to
There are lots of ways to compute this. What I did was to consider the complex contour integral
where γ is a semicircle in the complex plane with its diameter joining (-R, 0) and (R, 0) on the real axis. A bound for the integral over the arc of the circle is estimated to be
which vanishes as R goes to ∞. Then by the residue theorem, we have in the limit
and it follows that
