求面密度为 长短边为a*b的均匀矩形带电面在空间产生的电场强度分布

clc;
format long
rou = 2;
A = 3;
B = 1;
e0 = 8.854187817e-12;
d = 0.1;
N = A./d;
M = B./d;
z=0;
ex=0;
ey=0;
ez=0;

for i=0:N
    for j=0:M
        ex=ex+(rou/(4*pi*e0)).*(((x-i.*d).*d.^2)./((x-i.*d).^2+(y-j.*d).^2+z.^2).^1.5);
        ey=ey+(rou/(4*pi*e0)).*(((y-y.*d).*d.^2)./((x-i.*d).^2+(y-j.*d).^2+z.^2).^1.5);
        ez=ez+(rou/(4*pi*e0)).*((z.*d.^2)./((x-i.*d).^2+(y-j.*d).^2+z.^2).^1.5);
    end
end
E=(ex.^2+ey.^2+ez.^2).^0.5


%surf(x,y,E)

当Z轴不取0时只需改变z的取值即可。