使用2020，IGRF-13球谐系数，计算点位磁正常场。igrf.m，输入点位三维坐标(Use 2020, igrf-13 spherical harmonic coefficient to calculate point magnetic normal field,IGRF. M, input three-dimensional coordinates of point position)