This work is primarily concerned with the calculation of neutron capture cross sections and capture gamma-ray spectra, in the framework of the Hauser-Feshbach statistical model and for neutrons from the resonance region up to several MeV. An argument is made that, for applied purposes such as constructing evaluated cross-section libraries, nonstatistical capture mechanisms may be completely neglected at low energies and adequately approximated at high energies in a simple way. The use of gamma-ray strength functions to obtain radiation widths is emphasized. Using the reaction /sup 89/Y + n as an example, the problems encountered in trying to construct a case that could be run equivalently on two different nuclear reaction codes are illustrated, and the effects produced by certain parameter variations are discussed.