 Solution

The inverse change of variables is $x=v/(1+u)$, $y=uv/(1+u)$. Thus the triangular region becomesThe upper limit depends on $u$: the new region is not a rectangle. The absolute inverse Jacobian determinant is $x^2/(x+y)$, so the change of variables formula cancels the rational factor in the integrand:Although the original integrand is undefined at the vertex $x=y=0$, its nonnegative improper integral exists. One can apply the change of variables formula on $x\ge \epsilon$ and let $\epsilon \downarrow 0$ by the monotone convergence theorem; the transformed integrand has the finite integral just computed.