%% 此函数为计算 点偶极子 周围的磁场的函数　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*3*3的矩阵，为每个计算点处的张量矩阵
% I 为Num_obs*2的矩阵，第一列为每个计算点的 张量不变量I1，第二列为每个计算点的张量不变量I2
%%
function [B,G]=cioujizitensor(M,r0,r)