All cross sections of neutron-induced reactions, angular distributions, double-differential cross sections, and energy spectra for neutron, proton, deuteron, triton, helium, and alpha-particle emissions are consistently calculated and analyzed for 59Co with incident neutron energies below 20 MeV by using nuclear theoretical models, which are the optical model, preequilibrium and equilibrium reaction theories, and the distorted wave Born approximation theory. Theoretical calculated results are compared with existing experimental data and with evaluated data from ENDF/B-VII and JENDL-3.3.