Approximate numerical integration suggests a possible closed form:
I was not able to rigorously establish the equality, but the value is correct up to at least 900 decimal digits.
Is it the correct exact value of the integral I?
I will attempt to simplify the integral. It might take me a while and I might edit this answer a few times, so this is not its final form yet.
Let D be the differential operator.
We know ∫dx√x2−1=ln(2(x+√x2−1)).
So by using integration by parts we reduce the problem to solving A=∫∞1D[arccot(1+2πarcothx−arccscx)]ln(2(x+√x2−1))dx.
We know that D[arccot(x)]=−11+x2 hence by the chain rule for derivatives we get
Because we know that D[arccot(x)]=−11+x2 and also that ln(2z)=ln(2)+ln(z) we can use integration by parts again and arrive at
Now we can turn this into an infinite sum because we set U=2πarcoth(x)−arccsc(x) and use the Taylor expansion for z1+(1+z)2.
In other words we substitute z=U.
( Remember that ∫Σ=Σ∫ )
Finally the core problem is reduced mainly to solving:
by induction we get the need to solve for C1 and then are able to get the others.
At this point, I must admit that I have ignored convergeance issues.
Those issues can be solved by taking limits.
For instance C1 does not actually converge by itself.
For all clarity the problem is not resolved.
In fact it might require a new bounty. Still thinking about it …