matlab 矩形线圈磁场仿真

clear
clc
Nt=510;i0=1;u0=4*pi*1e-7;  %常数参数的设定
n=11;
a=[-0.139/20.139/2];
b=[-0.139/20.139/2];
c1=0.105;	
c2=-0.105;
k=Nt*i0*u0/（4*pi）;%常数参数的设定


xmax=0.0125;
ymax=0.0125;% 工作范围的限定
ngrid1=25;
x=linspace（-xmaxymaxngrid1）; %计算点分布
y=linspace（-xmaxymaxngrid1）; %计算点分布
Bx1（1:ngrid11:ngrid1）=zeros;By1（1:ngrid11:ngrid1）=zeros; 

for h=1:ngrid1
    for m=1:ngrid1
       
        for i=1:2
             for j=1:2
                 R=（-1）.^（i+j+1）*sqrt（（x（h）-c1）.^2+（y（m）-a（j））.^2+（0-b（i））.^2）;