Various cross sections of the (p + 206,207,208Pb), (p + natPb), (p + 209Bi), and (p + 197Au) reactions have been calculated at energies between 20 and 1600 MeV, using the Cascade-Exciton model, Geometry-Dependent Hybrid model, and Hybrid model. The neutron multiplicity has been predicted. The calculated results have been compared with the experimental data taken from the literature. The accordance of the calculated results with the experimental data has been examined.