抄録
Precise and high-speed estimation of perturbation due to geopotential field is essentially important in orbital analyses of near Earth satellites, since perturbation due to geopotential field is much larger than any other perturbations and its estimation is one of the most time-consuming process in the analyses. However, prevalent algorithm has spurious singularity at poles. Here we report a new algorithm to estimate accerelation due to geopotential field at high speed, three times as fast as older ones, without any singularity.