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）;
????????????