The "trick" with this problem is that is that the transformation x = (sin u ) / ( cos v ) and y = (sin v ) / ( cos u )
transforms the square 0 ≤ x ≤ 1 , 0 ≤ y ≤ 1 into the triangle 0 ≤ u ≤ π/2 , 0 ≤ v ≤ π/2− u. This can be done
by transforming each side of the square separately.
Notice that (sin u ) / ( cos v ) =1 suggests sin u = cos v or equivalently sin u = sin ( π /2 − v )
thus u + v = π/2
Then as Dayv O. correctly observed the Jacobian is the reciprocal of the integrant and therefore
the ∫10 ∫10 1 / ( 1- x2 y2 ) d x d y = ∫0π/2 ∫0π/2 - u dv du = ∫0π/2 ( π/2 −u )du =[ ( π/2) u − u2/2 ] |0π/2 =
π2/4 − π2/8 = π2/8

Dayv O.
I see better now. If u is allowed to be zero, which is required for x=0 bound, then v must be allowed to range for its bound to up to π/2. If you graph f(x,y) as descibed you can see The answer must be above 1 cubic unit, and your answer does that.07/31/21