Evaluate $\int_{y=-\infty}^1\int_{x=0}^\infty 4xy\sqrt{x^2+y^2}\mathrm dx\mathrm dy$ Evaluate $$\int_{y=-\infty}^1\int_{x=0}^\infty 4xy\sqrt{x^2+y^2}\mathrm dx\mathrm dy$$
I tried to solve this for hours without any success. I tried substitution method of $u = x^2+y^2$ but doesn't seem to work.
Here is my attempt:

$$\int_{-\infty}^1\int_0^\infty (16x^2y^2)^{\frac{1}{2}}(x^2 + y^2)^\frac{1}{2} \; \mathrm{d}x \mathrm{d}{y} = \int_{-\infty}^1 \int_0^\infty (16x^3y^2 + 16x^2 y^3)^{\frac{1}{2}} \; \mathrm{d}x\mathrm{d}{y}.$$
Let
\begin{align}
u &= 16x^3 y^2 + 16x^2 y^3 \\
\mathrm{d}u &= (48x^2y^2 + 32xy^3) \; \mathrm{d}{x} \\
\mathrm{d}u &= 16xy^2(3x + 2y) \; \mathrm{d}{x} \\
\mathrm{d}u &= 4xy(12xy + 8y^2) \; \mathrm{d}{x} \\
4xy \; \mathrm{d}{x} &= \frac{\mathrm{d}{u}}{12xy + 8y^2}.
\end{align}
$$\int_{-\infty}^1 \int_0^\infty 4xy \sqrt{x^2 + y^2} \; \mathrm{d}x \mathrm{d}y = \int_{-\infty}^1 \int_0^\infty \frac{\mathrm{d}u}{(12xy + 8y^2)}.$$

Attaching image for reference


EDITS
Original Post asked integrate: 4xy sqrt(x^2+y^2) where 0<x, y<1, and added this image. .
Later edits misread 0<x,y<1 as if 0<x<1and 0<y<1

Using polar coordinates, the integrand is
$$4\rho^2\cos\theta\sin\theta\,\rho\,\rho\,d\rho\,d\theta.$$
By symmetry, we can integrate twice in the triangle between $\theta=0,\theta=\dfrac\pi4$ and $\rho\cos\theta=1$ (or $\rho=\dfrac1{\cos\theta}$).
After integration on $\rho$,
$$8\cos\theta\sin\theta\ \left.\frac{4\rho^5}{5}\right|_0^{1/\cos\theta}=\frac{32}5\frac{\sin\theta}{\cos^4\theta}.$$
Then
$$\int_0^{\pi/4}\frac{\sin\theta}{\cos^4\theta}\,d\theta=\int_0^{1/\sqrt2}\frac{dt}{(1-t^2)^2}\,d\theta.$$
This integral is soluble by decomposition in simple fractions.
