球体磁梯度张量计算，根据文献中方法编写，可以快速计算各分量(Magnetic gradient tensor of sphere)

%% 此函数为计算 点偶极子 周围的磁场的函数 2013年5月6日 %% 计算坐标系及公式参考吴招才博士论文第4.2.1节（61页）的内容 % M为1*3的矩阵，为体磁矩，单位为A*m^2 % r0为1*3的矩阵，包含三个坐标值，r0(1)表示x0，r0(2)表示y0，r0(3)表示z0，为水平圆柱体的中心坐标为r0 % r为Num_obs*3的矩阵，其中存放的是要计算的点的坐标 %% % B 为Num_obs*3的矩阵，为每个计算点处的三分量磁场值 % G 为Num_obs*