用matlab数值求解电偶极子的等电势图和电场线图

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个电