contourplot arg(-4 (tan(x) + i tan(y))^3 + 1/2) from x = -pi/2 to pi/2 and y = -pi/2 to pi/2