Methods based on the density matrix for calculating the time-dependent probability density of a quantum system approaching equilibrium are presented. Explicit expressions are derived for the time-dependent probability density for a double-well potential. The effects of tunneling and transitions between energy levels on the probability density are discussed. For the case of closely spaced energy levels, a partial differential form of the density matrix equation is derived and used to calculate time-dependent probability densities.