In general, the deep understanding of proton-induced reactions is a crucial step for the further development of nuclear reactions theory. However there has been an interesting focus in nuclear physics. Some applications require accurate nuclear reaction data of common cross sections and especially need the data of neutron and proton induced energy-angle correlated spectra of secondary particles, as well as double differential cross sections. Double-differential nucleon-production cross-sections of Fe-56, Cu-63 and Zr-90 targets, bombarded with protons are calculated based on the nuclear theoretical models. Monte Carlo calculations with the TALYS 1.6 nuclear reaction simulation code are performed. Theoretical calculated results are compared with existing experimental data in EXFOR library.