  Previously, an explicit solution for the time evolution of the Wigner
function was presented in terms of auxiliary phase space coordinates which obey
simple equations that are analogous with, but not identical to, the classical
equations of motion. They can be solved easily and their solutions can be
utilized to construct the time evolution of the Wigner function. In this paper,
the usefulness of this explicit solution is demonstrated by solving a numerical
example in which the Wigner function has strong spatial and temporal variations
as well as regions with negative values. It is found that the explicit solution
gives a correct description of the time evolution of the Wigner function. We
examine next the pseudoparticle approximation which uses classical trajectories
to evolve the Wigner function. We find that the pseudoparticle approximation
reproduces the general features of the time evolution, but there are
deviations. We show how these deviations can be systematically reduced by
including higher-order correction terms in powers of $\hbar^2$.
