plot (x-3)^2 + (y-1)^2 = 1