This study is concerned with the calculations of double differential neutron cross-sections of the structural fusion materials of 56Fe and 90Zr isotopes that are bombarded with protons. Calculations were performed using the level density models of the TALYS 1.95 code and PHITS 3.22 Monte Carlo code. Constant Temperature Fermi Gas, Back Shifted Fermi Gas, and Generalized Super Fluid Models were employed for level density models. Calculations were performed at 22.2 MeV proton energies. Calculations were compared with the experimental data taken from Experimental Nuclear Reaction Data (EXFOR). In conclusion, the results showed that the level density model results of TALYS 1.95 codes for the double differential neutron cross-sections of 56Fe and 90Zr isotopes are consistent with experimental data. On the other hand, PHITS 3.22 results gave lower cross-section values than experimental data at 120 and 150°.