In pre-equilibrium nuclear reactions, the geometry-dependent hybrid model is applied with the use of the neutron and proton densities to investigate the effect of initial exciton numbers on the nucleon emission spectra. The initial exciton numbers calculated with the theoretical neutron and proton densities have been obtained within the Skryme-Hartree-Fock method with SKM* and SLy4 forces on target nuclei in the Fe-54,Fe-56(p, xp) reaction at 61.5-MeV incident proton energy by using a new calculationmethod of Tel et al. Also, the differences between the initial exciton numbers for protons and neutrons as a function of nuclear radius, focusing on systematic discrepancies correlated to differences in the proton and neutron densities have been investigated.