Плохо не знать небесной механики.....
Существует прекрасная формула Мак-Калуфа для определения потенциала гравитационного поля
обьекта произвольной формы на больших расстояниях от его центра (см. например В.А. Жаров "Сферическая астрономия"
стр.136-138, формулы (3.16), (3.17)
http://iaaras.ru/media/library/zharov_sf.pdf.
Формула основана на том же самом мультипольном разложении и из нее вроде как следуют результаты, рассмотренные выше.
Вопрос только в том, нельзя ли из этой формулы "вытянуть" что-то еще более изощренное чем две точечных массы?
Может нужно учитывать возможные ненулевые центробежные моменты инерции тела, т.е. "уйти от вращательной симметрии тела"....