title('fontname{ 隸書 }fontsize 偶極子的場 '),hold on
clear;clf;q=2e-6;k=9e9;a=1.5;b=-1.5;x=-6:0.6:6;y=x;
[X,Y]=meshgrid(x,y); % 設置坐標網點
rp=sqrt((X-a).^2+(Y-b).^2);rm=sqrt((X+a).^2+(Y+b).^2);
V=q*k*(1./rp-1./rm); % 計算電勢
[Ex,Ey]=gradient(-V); % 計算場強
AE=sqrt(Ex.^2+Ey.^2);Ex=Ex./AE;Ey=Ey./AE;% 場強歸一化，使箭頭等長
cv=linspace(min(min(V)),max(max(V)),49);% 產生 49 個電位值
contourf(X,Y,V,cv,'k-') % 用黑實線畫填色等位線圖
%axis('square') % 在 Notebook 中，此指令不用
title('fontname{ 隸書 }fontsize 偶極子的場 '),hold on
quiver(X,Y,Ex,Ey,0.7) % 第五輸入宗量 0.7 使場強箭頭長短適中。
plot(a,b,'wo',a,b,'w+') % 用白線畫正電荷位置
plot(-a,-b,'wo',-a,-b,'w-') % 用白線畫負電荷位置
xlabel('x');ylabel('y'),hold off