A new method for solving the Wigner-Liouville equation allowing a treatment of quantum dynamics has been developed. The method combines two widely known approaches: molecular dynamics and Monte Carlo. Numerical results have been obtained for one-, two- and three-dimensional potential wells of realistic type. The quantum dynamics method was tested by a comparison of the numerical results with analytical estimates for such values as the average position (x) over bar(t), average momentum (p) over bar(t), energy (E) over bar(t), position and momentum dispersions, and Heisenberg inequality.