THERMAL: A routine designed to calculate neutron thermal scattering

THERMAL is designed to calculate neutron thermal scattering that is isotropic in the center of mass system. At low energy thermal motion will be included. At high energies the target nuclei are assumed to be stationary. The point of transition between low and high energies has been defined to insure a smooth transition. It is assumed that at low energy the elastic cross section is constant in the center of mass system. At high energy the cross section can be of any form. You can use this routine for all energies where the elastic scattering is isotropic in the center of mass system. In most materials this will be a fairly high energy