用matlab画出电偶极子的电场线

q=1;d=2;e0=8.854187817*10.^-12;%真空介电常数x=-3:0.1:3;y=-3:0.1:3;[x,y]=meshgrid(x,y);z=q.*(1./sqrt((y-1).^2+x.^2+0.01)-1./sqrt((y+1).^2+x.^2+0.01))./(4*pi*e0);%求电势，运用公式E=q4πε0(1𝑟1−1𝑟2)contour(x,y,z);[px,py]=gradient(z);%求电势在x，y方向的梯度（即电场强度）holdonstreamslice(x,y,px,py,'k')%得到结果图