The maximum entropy approach is used to solve the classical moment problem of stellar kinematics.
If an extended set of moments is available, the current method provides a
linear estimation algorithm, which is given by a Gramian system of
equations, that leads to a fast and suitable estimation of the velocity
distribution. In particular, it can be used as an alternative approach for modelling
multimodal distributions that can not be described through gaussian mixtures.