Evolution of Observables in a Numerical Kinetic Model
  The numerical solutions of the nonrelativistic and relativistic Boltzmann equations have been studied at various initial conditions. Particularly, the known analytical solution of the nonrelativistic Boltzmann equation at spherically symmetric Gaussian initial conditions has been recovered and the corresponding conservation of particle spectra and correlation functions has been confirmed. Similar conservation properties have been found also for a more general class of solutions describing the nonrelativistic and relativistic systems at anisotropic initial conditions. 
N.S. Amelin, R. Lednicky, L. V. Malinina, T. A. Pocheptsov, Yu. M.   Sinyukov