圆柱体磁异常正演

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