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

T =y/(x^2 + y^2)
contour_plot(T, (x, -500, 500), (y, -500,500))