The project of European Space Agency and a recent achievement in satellite geodesy, the GOCE satellite mission (Gravity of Field and Steady-State of Ocean Circulation Explorer) exploits a method of satellite gradientometry. Gravitational field of the Earth is usually represented as a finite series of spherical harmonic functions, the model containing a finite number of coefficients, Cnm, Snm. The coefficients Cnm, Snm are derived in our work, based on the second method of Neumann and the Gauss-Legandre quadrature decomposition.

- Sneeuw N. Global spherical harmonic analysis by least-squares and numerical quadrature methods in historical perspective, Delft University of Technology, Faculty of Geodetic Engineering, Geophys. J. Int.,1994. 118, – P.707-716.
- GOCE. Level 2. Product Data Handbook
- Lundberg J.B., Schutz B.E. (1988) Recursion formulas for Legendre functions for use with nonsingular geopotential models // J. Guidance Control Dyn. 11(1), – P. 31-38.
- Gofman–Vellengof B., Moric G. Fizicheskaja geodezija. – M.:, MiiGAiK, – 2007. – S. 285-286