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

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