A 50 molecule sample in contact with a heat bath of similar molecules was used as a model to calculate the rotational collision numbers, Zrot, for HCI and DCI at 300 and SOOoK, by the method of molecular dynamics. The mean rotational energy of the sample exhibits an exponential approach to a steady state value; however, there was considerable scatter about the exponential curve. The results for Zrot generally agree with previous theoretical calculations. They are intermediate to experimental acoustic data and the earlier perturbation results. second, and fourth moments of angular momentum and Boltzmann's H function, the second moment being essentially the kinetic energy. They experienced difficulties in extracting relaxation times because of large fluctuations. An interesting aspect of this work is the demonstration of the existence of a set of initial 4523