  An alternative Monte Carlo estimator for the one-body density rho(r) is
presented. This estimator has a simple form and can be readily used in any type
of Monte Carlo simulation. Comparisons with the usual regularization of the
delta-function on a grid show that the statistical errors are greatly reduced.
Furthermore, our expression allows accurate calculations of the density at any
point in space, even in the regions never visited during the Monte Carlo
simulation. The method is illustrated with the computation of accurate
Variational Monte Carlo electronic densities for the Helium atom (1D curve) and
for the water dimer (3D grid containing up to 51x51x51=132651 points).
