The basic goal of this study (as the first step) is to collect the appropriate set of the fundamental astronomic-geodetics parameters for their further use to obtain the components of the density distributions for the terrestrial and outer planets of the Solar system (in the time interval of more than 10 years). The initial data were adopted from several steps of the general way of the exploration of the Solar system by iterations through different spacecraft. The mechanical and geometrical parameters of the planets allow finding the solution of the inverse gravitational problem (as the second stage) in the case of the continued Gaussian density distribution for the Moon, terrestrial planets (Mercury, Venus, Earth, Mars) and outer planets (Jupiter, Saturn, Uranus, Neptune). This law of Gaussian density distribution or normal density was chosen as a partial solution of the Adams-Williamson equation and the best approximation of the piecewise radial profile of the Earth, including the PREM model based on independent seismic velocities. Such conclusion already obtained for the Earth’s was used as hypothetic in view of the approximation problem for other planets of the Solar system where we believing to get the density from the inverse gravitational problem in the case of the Gaussian density distribution for other planets because seismic information, in that case, is almost absent. Therefore, if we can find a stable solution for the inverse gravitational problem and corresponding continue Gaussian density distribution approximated with good quality of planet’s density distribution we come in this way to a stable determination of the gravitational potential energy of the terrestrial and giant planets. Moreover to the planet’s normal low, the gravitational potential energy, Dirichlet’s integral, and other planets’ parameters were derived. It should be noted that this study is considered time-independent to avoid possible time changes in the gravitational fields of the planets.