Uncertainty Quantification for Multidimensional Dynamical Systems Based on Adaptive Numerical Solution of Liouville Equation