Uncertainty quantification of multidimensional dynamical systems based on adaptive numerical solutions of the Liouville equation