[x,y]=meshgrid(-2:0.1:2,-2:0.1:2);
%以0.1為步長建立平面數據網格
z=1./sqrt((x-1).^2+y.^2+0.01)... %寫出電勢表達式
-1./sqrt((x+1).^2+y.^2+0.01);
[px,py]=gradient(z);
%求電勢在x,y方向的梯度即電場強度
contour(x,y,z,[-12,-8,-5,-3,-1,... %畫出等勢線
-0.5,-0.1,0.1,0.5,1,3,5,8,12])
hold on %作圖控制
quiver(x,y,px,py,'k') %畫出各點上電場的大小和方向