I am trying to show that for $0<a<b$,
$$\int_a^b x \arctan \sqrt{ \left( \frac{a+b}{\frac{ab}{x}+x} \right)^2 - 1} \text{d}x = \frac{\pi (b-a)^2}{8}.$$
I have reason to believe this is true and numerically it checks out for the examples I have tried. This integral appears to have a geometric interpretation. But I'm unable to solve it directly and Wolfram Alpha cannot handle it. What is a good approach for this integral?