Evaluar: $$I=\iint_{[0,1]^2}\frac{dxdy}{(1+x^2+y^2)^{3/2}}.$$
El intento. Trabajando en la sustitución de $x=\sqrt{1+y^2}\,\sinh u$ tenemos $$I=\int\limits_{0}^{1}\frac{x}{(1+y^2)\sqrt{x^2+y^2+1}}\Bigg|_{0}^{1}\,dy=\int\limits_{0}^{1}\frac{dy}{(1+y^2)\sqrt{2+y^2}}$$ and then working on $y=\sqrt{2}\,\bronceado u$ obtenemos: $$I=\arctan\left(\frac{y}{\sqrt{y^2+2}}\right)\bigg|_{0}^{1}=\arctan\left(\frac{1}{\sqrt{3}}\right)=\frac{\pi}{6}.$$
Por otro lado, las coordenadas polares:
$$I=\int\limits_{0}^{\pi/4}\int\limits_{0}^{\frac{1}{\cos \phi}}\frac{r}{(1+r^2)^{3/2}}dr d\phi+\int\limits_{\pi/4}^{\pi/2}\int\limits_{0}^{\frac{1}{\sin \phi}}\frac{r}{(1+r^2)^{3/2}}dr d\phi$$ lo que nos lleva a: $$I=\int\limits_{0}^{\pi/4}\left(1-\frac{\cos\phi}{\sqrt{1+\cos^2\phi}}\right)d\phi+\int\limits_{\pi/4}^{\pi/2}\left(1-\frac{\sin\phi}{\sqrt{1+\sin^2\phi}}\right)d\phi$$ y los cálculos que parecen difíciles de manejar.
Hay una manera más fácil abordar el cálculo de esta integral?
Gracias de antemano.