clc
clear
//input
r=5500//radius
g1=6.7*10^-11
g=7//acceleration due to gravity
//calculation of mean density
p=3*g/(4*%pi*r*10^3*g1)//mean density
//output
printf("the mean density of planet is %3.3f kgm^-3",p)