plot z = 1/(2 pi) arctan(sqrt((1 - x^2)/(1 - y^2))) for x=-1..1, y=-1..1