[solution]
**(a) First iterated integral.**
**Step 1 (Inner integral in $y$).** Observe that $\frac{x^2 - y^2}{(x^2 + y^2)^2} = -\frac{\partial}{\partial y}\left(\frac{y}{x^2 + y^2}\right)$. Verification: $\frac{\partial}{\partial y}\frac{y}{x^2 + y^2} = \frac{(x^2 + y^2) - y \cdot 2y}{(x^2 + y^2)^2} = \frac{x^2 - y^2}{(x^2 + y^2)^2}$. Therefore:
\begin{align*}
\int_0^1 \frac{x^2 - y^2}{(x^2 + y^2)^2} \, dy = \left[\frac{y}{x^2 + y^2}\right]_{y=0}^{y=1} \cdot (-1) \quad \text{...}
\end{align*}
Wait — we have $f = \frac{\partial}{\partial y}\frac{y}{x^2 + y^2}$, not $-\frac{\partial}{\partial y}$. So:
\begin{align*}
\int_0^1 f(x,y) \, dy = \left[\frac{y}{x^2 + y^2}\right]_0^1 = \frac{1}{x^2 + 1} - 0 = \frac{1}{x^2 + 1}.
\end{align*}
**Step 2 (Outer integral in $x$).**
\begin{align*}
\int_0^1 \frac{1}{x^2 + 1} \, dx = [\arctan(x)]_0^1 = \frac{\pi}{4}.
\end{align*}
**(b) Second iterated integral.**
**Step 3 (By antisymmetry).** Observe $f(x, y) = \frac{x^2 - y^2}{(x^2 + y^2)^2} = -f(y, x)$. Therefore:
\begin{align*}
\int_0^1\left(\int_0^1 f(x,y) \, dx\right) dy = \int_0^1\left(-\int_0^1 f(y, x) \, dx\right) dy = -\int_0^1 \frac{1}{y^2 + 1} \, dy = -\frac{\pi}{4}.
\end{align*}
The two iterated integrals are $\pi/4$ and $-\pi/4$ — they differ.
**(c) Why no contradiction with Fubini.**
**Step 4 (Check integrability of $|f|$).** By Tonelli's theorem, if $\int \int |f| < \infty$, then Fubini applies and the iterated integrals must agree. Since they differ, $\int_0^1 \int_0^1 |f(x,y)| \, dx \, dy = \infty$. Indeed, near $(0, 0)$: $|f(x, y)| \approx |x^2 - y^2|/(x^2 + y^2)^2$. In polar coordinates $x = r\cos\theta$, $y = r\sin\theta$: $|f| \approx |\cos(2\theta)|/r^2$, and $\int_0^\varepsilon r^{-2} \cdot r \, dr = \int_0^\varepsilon r^{-1} \, dr = \infty$.
Fubini's theorem requires $f \in L^1$. Since $|f| \notin L^1((0,1)^2)$, the hypotheses fail and unequal iterated integrals are permitted. $\blacksquare$
[/solution]