Integrate[1/(1 - x^2*y^2), {x, 0, 1}, {y, 0, 1}]