function E=Efield(x,y,z)
E=[-x,-y,-z];