x = var('x')
y = var('y')

T = arctan(y/x)
contour_plot(T, (x, -500, 500), (y, -500,500))