plot3d Re[cot(1/(1.5x+i y))] from x = -1/pi to 1/pi, y = -1/pi to 1/pi