  We present an improved POTENT method for reconstructing the velocity and mass
density fields from radial peculiar velocities, test it with mock catalogs, and
apply it to the Mark III Catalog. Method improvments: (a) inhomogeneous
Malmquist bias is reduced by grouping and corrected in forward or inverse
analyses of inferred distances, (b) the smoothing into a radial velocity field
is optimized to reduce window and sampling biases, (c) the density is derived
from the velocity using an improved nonlinear approximation, and (d) the
computational errors are made negligible. The method is tested and optimized
using mock catalogs based on an N-body simulation that mimics our cosmological
neighborhood, and the remaining errors are evaluated quantitatively. The Mark
III catalog, with ~3300 grouped galaxies, allows a reliable reconstruction with
fixed Gaussian smoothing of 10-12 Mpc/h out to ~60 Mpc/h. We present maps of
the 3D velocity and mass-density fields and the corresponding errors. The
typical systematic and random errors in the density fluctuations inside 40
Mpc/h are \pm 0.13 and \pm 0.18. The recovered mass distribution resembles in
its gross features the galaxy distribution in redshift surveys and the mass
distribution in a similar POTENT analysis of a complementary velocity catalog
(SFI), including the Great Attractor, Perseus-Pisces, and the void in between.
The reconstruction inside ~40 Mpc/h is not affected much by a revised
calibration of the distance indicators (VM2, tailored to match the velocities
from the IRAS 1.2Jy redshift survey). The bulk velocity within the sphere of
radius 50 Mpc/h about the Local Group is V_50=370 \pm 110 km/s (including
systematic errors), and is shown to be mostly generated by external mass
fluctuations. With the VM2 calibration, V_50 is reduced to 305 \pm 110 km/s.
