Show that the function $$f(x,y) = \frac{x^2-y^2}{(x^2+y^2)^2}$$ is not Lebesgue integrable on the unit square $[0,1]^2$.
It's been asked here before (Is the function $f(x,y) = \frac{x^2 - y^2}{(x^2+y^2)^2}$ lebesgue integrable in $[0,1]\times[0,1]$?). However I don't really understand the answers.
I'm not very comfortable with product measures, but it seems Fubini's theorem might be useful.
Edit: I was thinking about something along the same lines as in this example from Wikipedia. Is that correct?
A different thing you can try for this particular integral is the following. First, note that since $x^2=(-x)^2$, we have $x^2-y^2$ is invariant under flipping across the $x$ or $y$ axes. Thus its Lebesgue integrable iff the integral over $[-1,1]^2$ is finite. And for this it suffices to show that the (unsigned) integral over a small ball around 0, $B(0,r)$ is infinite.
For integrals on balls around singularities, you can count the degree of the polynomial as a heuristic. On top, you have a degree 2 polynomial, and at the bottom you have a degree 4 polynomial. Therefore, at best, you have a pole of order 2, and in $\mathbb R^d$, only poles of order $<d$ are integrable. This is something you can verify for the simpler integral $\int_{B(0,r)} \frac{d\vec x}{|\vec x|^{d}}$ using radial coordinates, and you can prove each case via calculations similar to the other answer.
Regarding your edit, yes, you can absolutely do it that way.
Fubini's Theorem says that if a function $f(x,y)$ is Lebesgue integrable on the product space, then the two iterated integrals must be the same. Wikipedia uses the contrapositive- they show that the two iterated integrals are different, which implies that $f(x,y)= -\partial_x\partial_y \arctan(y/x)$ is not Lebesgue integrable.
Well, they didn't exactly show it. But its not that hard, if you knew that $f$ miraculously had this formula! My experience is that this is quite rare, and the other methods work more generally.
But lets work through it anyway. You can proceed by first computing derivatives and working backwards: since $$\int_0^1 -\partial_x\partial_y \arctan(y/x) dx = -\partial_y \arctan(y/x)|_{x=0}^1 = \left.\frac{-1}x\frac1{y^2/x^2 + 1}\right|_{x=0}^1 = \left.\frac{-x}{y^2+x^2}\right|_{x=0}^1 = \frac{-1}{y^2+1}$$ we have $$\int_0^1 \int_0^1 -\partial_x\partial_y \arctan(y/x) dxdy =\int_0^1\int_0^1 \frac{x^2 - y^2}{(x^2+y^2)^2} dxdy = \int_0^1 \frac{-1}{y^2 + 1} dy < 0$$ On the other hand - $$\int_0^1 -\partial_x\partial_y \arctan(y/x) dy = -\partial_x \arctan(y/x)|_{y=0}^1 = -1\times \left.\frac{-y}{x^2}\frac1{y^2/x^2 + 1}\right|_{y=0}^1 = \left.\frac{y}{y^2+x^2}\right|_{y=0}^1 = \frac{1}{x^2+1}$$ which implies $$\int_0^1 \int_0^1 -\partial_x\partial_y \arctan(y/x) dydx = \int_0^1\int_0^1 \frac{x^2 - y^2}{(x^2+y^2)^2} dydx = \int_0^1 \frac{+1}{x^2 + 1} dy > 0$$ A related observation- note that $f(x,y) = -f(y,x)$. Therefore
$$ \int_0^1 \int_0^1 f(x,y) dxdy = \int_0^1 \int_0^1 -f(y,x) dxdy$$ But in the right hand side, $y$ and $x$ are dummy variables that are integrated over the same region. So we can switch all $y$s with all $x$s to get $$ \int_0^1 \int_0^1 f(x,y) dxdy = - \int_0^1 \int_0^1 f(x,y) dy dx$$ This means that if $f$ is integrable, its integral must be $0$. So for example, $\sin\frac{x^2 - y^2}{(x^2+y^2)^2}$ is Lebesgue integrable on $[0,1]^2$ with integral 0.