  We assume that a finite set of moments of a random vector is given. Its
underlying density is unknown. An algorithm is proposed for efficiently
calculating Dirac mixture densities maintaining these moments while providing a
homogeneous coverage of the state space.
