3305.sagews

epsilon=.1
q=1
k=(1/(4*pi*epsilon))
E(x,y)=[k*x/((x^2+y^2)*sqrt(x^2+y^2)),k*y/((x^2+y^2)*sqrt(x^2+y^2))]