I will discuss a method of solving the full 6-dimensional Vlasov-Poisson equation describing the Liouville dynamics of the system of gravitating particles via a system of auxiliary Schroedinger-like equations and demonstrate the validity of the technique on several simple test cases.