Cross sections for production of gamma rays due to neutron interactions with iron have been measured as a function of both neutron and gamma-ray energy. Two experimental configurations were used to obtain the data: a Nal-spectrometer system using the Oak Ridge Linear Accelerator as the neutron source and a Ge(Li)-spectrometer system using a pulsed Van de Graaff and the D( d, n) reaction as the neutron source. The Nal-spectrometer system, described completely in this report, was used to acquire data for 0.8 ≤ En ≤ 20 MeV and θγ = 125 deg, which were unfolded to obtain d2σ/dωdE values for gamma-ray energies between 0.7 and 10 MeV. The Ge(Li) system was used to obtain high resolution information on the production of discrete-line dσ/dω values for 4.85 ≤ En ≤ 9.0 MeV and θγ = 55, 75, and 90 deg. Our data are compared with previously reported experimental data and with the current ENDF/B evaluation. Although there is generally reasonable (20%) agreement, important differences among these data are discussed.