Integral forms of the Boltzmann-like transport equations are derived to treat neutron emissions from heavy ion collisions. The derivations are based on the probability of contribution (score) from nucleon occupied phase space bins to the double differential emission spectra. Two-body interactions among the nucleons in hot spot and cold spot regions of the composite system resulting in creation, annihilation, or energy redistribution of the colliding particles are considered. The resulting equations are solved numerically in energy-angle bins to study thermal and kinematic equilibrium of the target 1 projectile (165Ho+20Ne) excitedsystem at different projectile energies (10 MeV/A to 30 MeV/A). ©2000 The American Physical Society.