Analytical computation of gravity effects for polyhedral bodies

On the basis of recent analytical results we derive new formulas for computing the gravity effects of polyhedral bodies which are expressed solely as function of the coordinates of the vertices of the relevant faces. We thus prove that such formulas exhibit no singularity whenever the position of the observation point is not aligned with an edge of a face. In the opposite case, the contribution of the edge to the potential to its first-order derivative and to the diagonal entries of the second-order derivative is deemed to be zero on the basis of some claims which still require a rigorous mathematical proof. In contrast with a common statement in the literature, it is proved that only the off-diagonal entries of the second-order derivative of the potential do exhibit a noneliminable singularity when the observation point is aligned with an edge of a face. The analytical provisions on the range of validity of the derived formulas have been fully confirmed by the Matlab® program which has been coded and thoroughly tested by computing the gravity effects induced by real asteroids at arbitrarily placed observation points.