x,y = var('x y')
f=1/sqrt(x^2+y^2)
grad= f.gradient([x,y])
plot_vector_field(grad, (x,-1,1), (y,-1,1)).show()
gives me this
As you can see the errors very close to the pole "cross out" the whole pic. Obviously they are very large vectors, but they do harm the overall image. i have other plots where its worse. How can i work around that? can i blank out the direkt proximity of the pole, somehow? can i limit the vektors' length? How?