The method of the geopotential parameters determination with the use of the gradiometry data is considered. The second derivative of the gravitational potential in the correction equation on the rectangular coordinates x, y, z is used as a measured variable. For the calculated value of the measured quantity required for the formation of a free member of the correction equation, the the Cunningham polynomials were used. We give algorithms for computing the second derivatives of the Cunningham polynomials on rectangular coordinates x, y, z, which allow to calculate the second derivatives of the geopotential at the rectangular coordinates x, y, z.Then we convert derivatives obtained from the Cartesian coordinate system in the coordinate system of the gradiometer, which allow to calculate the free term of the correction equation. Afterwards the correction equation coefficients are calculated by differentiating the formula for calculating the second derivative of the gravitational potential on the rectangular coordinates x, y, z. The result is a coefficient matrix of the correction equations and corrections vector of the free members of equations for each component of the tensor of the geopotential. As the number of conditional equations is much more than the number of the specified parameters, we go to the drawing up of the system of normal equations, from which solutions we determine the required corrections to the harmonic coefficients.
A nonlinear second derivative method with a variable step-size based on continued fractions for singular initial value problems
