I have to solve the following double integral
$$\iint_{A}\frac{dx\,dy}{(1+x^2)(1+x^2 y^2)}$$
with $A= \left[0,+\infty\right[ \times [0,1].$ So far I've tried to solve it integrating w.r.t. $y$ first.
$$\iint_0^1\frac{dy\,dx}{(1+x^2)(1+x^2 y^2)} = \int_0^\infty\frac{1}{1+x^2}\int_0^1\frac{dy}{1+x^2 y^2} \, dx. $$
I've solved the internal integral by substitution, remembering that $\int\frac{du}{1+u^2}=\arctan u$ Substitution:
$$x^2 y^2= u^2 \to y=\frac{1}{x}u \to dy=\frac{1}{x}du.$$ $$y=0 \to u=0, \qquad y=1 →u=x$$
So:
\begin{align} & \int_0^∞\frac{1}{1+x^2} \left( \int_0^x\frac{1}{x}\frac{1}{1+u^2}\,du \right) \, dx \\[8pt] = {} & \int_0^\infty\frac{1}{x}\frac{1}{1+x^2}[\arctan u] \, dx \\[8pt] = {} & \int_0^\infty\frac{\arctan x}{x(1+x^2)} \, dx. \end{align}
Now I have to solve this last integral with the Fubini's theorem but I don't know how to do it.
Solution 1. Let $I$ denote the double integral. Then
\begin{align*} I &= \int_0^1 \int_0^\infty \frac{1}{(1+x^2)(1+x^2 y^2)} \, \mathrm{d}x\,\mathrm{d}y \\ &= \int_0^1 \int_0^\infty \frac{1}{1-y^2} \left( \frac{1}{1+x^2} - \frac{y^2}{1+x^2y^2} \right) \, \mathrm{d}x\,\mathrm{d}y \\ &= \int_{0}^{1} \frac{1}{1-y^2} \left( \frac{\pi}{2} - \frac{\pi y}{2} \right) \, \mathrm{d}y \\ &= \frac{\pi}{2} \int_0^1 \frac{1}{1+y} \, \mathrm{d}y \\ &= \frac{\pi}{2} \log 2. \end{align*}
Solution 2.
\begin{align*} I = \int_0^\infty \int_0^1 \frac{1}{(1+x^2)(1+x^2 y^2)} \, \mathrm{d}y \, \mathrm{d}x = \int_0^\infty \frac{\arctan x}{(1+x^2)x} \, \mathrm{d}x. \end{align*}
Substituting $x = \tan\theta$, then
\begin{align*} \require{cancel} I &= \int_0^{\frac{\pi}{2}} \frac{\theta}{\tan\theta} \, \mathrm{d}\theta = \cancel{\left[ \theta \log \sin\theta \right]_0^{\frac{\pi}{2}}} - \int_0^{\frac{\pi}{2}} \log \sin\theta \, \mathrm{d}\theta. \end{align*}
Regarding the last integral, we observe that the following holds:
$$ I = -\int_0^{\frac{\pi}{2}} \log \sin\theta \, \mathrm{d}\theta = -\int_0^{\frac{\pi}{2}} \log \cos\theta \, \mathrm{d}\theta = -\int_0^{\frac{\pi}{2}} \log \sin(2\theta) \, \mathrm{d}\theta. $$
Using this, we get
$$ 2I = -\int_0^{\frac{\pi}{2}} \log (\sin\theta\cos\theta) \, \mathrm{d}\theta = -\int_0^{\frac{\pi}{2}} \log \left(\frac{\sin2\theta}{2}\right) \, \mathrm{d}\theta = I + \frac{\pi}{2}\log 2. $$
Therefore
$$ I = \frac{\pi}{2}\log 2. $$