k=1;q=1;d=0.01;h=0.8;
[xy]=meshgrid（[-h:0.1:h]）;%建立二維網格
r2=sqrt（（x-d/2）.^2+y.^2+0.001）;r1=sqrt（（x+d/2）.^2+y.^2+0.001）;
phi=k*q*（1./r2-1./r1）;%電勢計算公式
[ExEy]=gradient（-phi）;%求電勢phi對x和y的梯度
AE=sqrt（Ex.^2+Ey.^2）;Ex=Ex./AE;Ey=Ey./AE;%場強歸一化，使箭頭等長
CV=linspace（min（min（phi））max（max（phi））50）;%產生49個電