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