  The direct gravimetry problem is solved using the subdivision of each body of
a deposit into a set of vertical adjoining bars, and in the inverse problem
each body of a deposit is modeled by a uniform ellipsoid of revolution
(spheroid). Well-known formulas for z-component of gravitational intensity of a
spheroid are transformed to a convenient form. Parameters of a spheroid are
determined by minimizing the Tikhonov smoothing functional using constraints on
the parameters. This makes the ill-posed inverse problem by unique and stable.
The Bulakh algorithm for initial estimating the depth and mass of a deposit is
modified. The technique is illustrated by numerical model examples of deposits
in the form of two and five bodies. The inverse gravimetry problem is
interpreted as a gravitational tomography problem or the intravision of the
Earth's crust and mantle.
