How to compute the $$\int_{0}^{1} \int_{0}^{1} \frac{(\log(1+x^2)-\log(1+y^2))^2 }{|x-y|^{2}}dx dy.$$ Is it possible to compute the integral analytically upto some terms. I believe it should involve hypergeometric series. Any ideas are welcome.