A new method for computing the ellipsoidal correction for Stokes's formula

 This paper generalizes the Stokes formula from the spherical boundary surface to the ellipsoidal boundary surface. The resulting solution (ellipsoidal geoidal height), consisting of two parts, i.e. the spherical geoidal height N0 evaluated from Stokes's formula and the ellipsoidal correction N1, makes the relative geoidal height error decrease from O(e2) to O(e4), which can be neglected for most practical purposes. The ellipsoidal correction N1 is expressed as a sum of an integral about the spherical geoidal height N0 and a simple analytical function of N0 and the first three geopotential coefficients. The kernel function in the integral has the same degree of singularity at the origin as the original Stokes function. A brief comparison among this and other solutions shows that this solution is more effective than the solutions of Molodensky et al. and Moritz and, when the evaluation of the ellipsoidal correction N1 is done in an area where the spherical geoidal height N0 has already been evaluated, it is also more effective than the solution of Martinec and Grafarend.