clear?all;
clc;
u0=4*pi*1E-7;?%絕對磁導率
K=1E-3;
xc=50;
zc=30;?%柱體中心深度
z=0;
r=8;?%柱體半徑
B0=50000.0;?%nT
B0=B0*1E-9;?%T
G=6.67*1E-11;
I0=45;
D0=5;
A=12;
Ir=60;
Dr=30;
Mr=30;?%nT
DEG2ARC=pi/180.0;
ARC2DEG=180.0/pi;
I0=I0*DEG2ARC;
D0=D0*DEG2ARC;
A=A*DEG2ARC;
Ir=Ir*DEG2ARC;
Dr=Dr*DEG2ARC;
x=0:100;
iNumx=size（x2）;

B0x=B0*cos（I0）*cos（D0-A）;
Mix=K*B0x/u0
B0z=B0*sin（I0）;
Miz=K*B0z/u0
Mrx=Mr*cos（Ir）*cos（Dr-A）
Mrz=Mr*sin（Ir）

Bx=zeros（1iNumx）;
Bz=zeros（1iNumx）;
dB=zeros（1iNumx）;
%dB=zeros（iNumxiNumy）;
for?i=1:iNumx
????Bx（i）=1E9*u0*r*r*（（Mix+Mrx）*（（xc-x（i