clear all; close all; % q1=2; % q2=-2; e0=8.85*1e-12; k=1/(4*pi*e0); q1s=[1,2,3]; q2s=[-3,-2,-1]; x=-5:5; y=-5:5; r1=[-2,0]'; r2=[2,0]' [X,Y]=meshgrid(x,y); a1=sqrt((X-r1(1)).^2+(Y-r1(2)).^2); a2=sqrt((X-r2(1)).^2+(Y-r2(2)).^2); for i1=1:l…

clear all; close all; q=2; x=-10:10; y=-10:10; [X,Y]=meshgrid(x,y); Z=sqrt(100-X.^2-Y.^2); Ex=100*X./(X.^2+Y.^2+Z.^2).^(3/2); Ey=100*Y./(X.^2+Y.^2+Z.^2).^(3/2); Ez=100*Z./(X.^2+Y.^2+Z.^2).^(3/2); Ex(isnan(Ex))=0; Ey(isnan(Ey))=0; Ez(isnan(…

clear all; close all; q=2; r=[[1,0]',[0,1]',[-1,0]',[0,-1]',[1.5,1.5]',[-1.5,1.5]',[-1.5,-1.5]',[1.5,-1.5]']; r2=arrayfun(@(x,y)x.^2+y.^2,r(1,:),r(2,:)) r1=sqrt(r2); e0=8.85*1e-12; E=1/(4*pi*e0)*q./r2./r1.*r; plot(0,0,'bo'); hold on; array…

clear all; close all; qA=-2; qB=1; rA=[2,1,0]'; rB=[0,1,0]'; r=[rA rB]; e0=8.85*1e-12; k = 1/4/pi/e0; F = k*qA*qB/sqrt(sum((rA-rB).^2)); Fdummy=1;%correction for view F=Fdummy;%The real F is too huge to be shown, so just leave it v=rA-rB; …