The interior potential of a spherical planetary body (and the corresponding density distribution) is expressed as the sum of two parts, the first given uniquely by the external gravity field and the surface density and the second depending on an arbitrary function. The first part of the potential (density distribution) is shown to be a 3-harmonic (biharmonic) function, while the second part can be expressed as a series whose i-th term (for i ≥ 0) is an (i + 4)-harmonic ((i + 3)-harmonic) function. From this general solution a single solution is then chosen: this is done by imposing certain natural conditions on it, among others that this particular solution is an n-harmonic function for n as small as possible.
Read full abstract