Здравствуйте, Петров Игорь Андреевич!
Поскольку так и не удалось выяснить, использование каких именно значений подразумевалось авторами задачи, рассмотрим несколько вариантов определения плотности.
Сперва варианты с использованием гравитационной постоянной
A. Ускорение свободного падения и диаметр
g=mG/R
2m=gR
2/G
[$961$]=m/((4/3)[$960$]R
3)=3g/(4G[$960$]R)=3g/(2G[$960$]D)
B. Ускорение свободного падения и масса
g=mG/R
2R
2=mG/g
[$961$]=m/((4/3)[$960$]R
3)=m/((4/3)[$960$](mG/g)
1,5)=3g
1,5/(4[$960$]G
1.5[$8730$](m))
Теперь варианты с использованием отношений и плотности Земли
C. Ускорение свободного падения и диаметр
[$961$]=[$961$]
З[$183$](g/g
З)/(D/D
З)
D. Ускорение свободного падения и масса
[$961$]=[$961$]
З[$183$](g/g
З)
1,5/[$8730$](m/m
З)
И в ещё пара вариантов непосредственно с использованием размеров планет
E. [$961$]=m/((4/3)[$960$]R
3)=6m/[$960$]D
3F. [$961$]=[$961$]
З[$183$](m/m
З)/(D/D
З)
3G. [$961$]=[$961$]
З[$183$](m/m
З)/(V/V
З)
Результаты расчёта (полная версия в приложении)
Как видите, в отдельных случаях точность исходных данных заметно влияет на результат. Также нужно учесть, что экваториальный диаметр может значительно отличаться от среднего диаметра из-за формы планет.