An important ingredient in the performance of accelerator driven systems for energy production, waste transmutation and other applications are the number of spallation neutrons produced per incident proton. The neutron multiplicities, angular and energy distributions are usually calculated using simulation codes. We have presented multiplicities of the neutrons ernitted in the interaction of a high-energy proton (1500 MeV) with a thin target Pb, Bi. In this study we have used the code ISABEL to calculate multiplicities of the neutron emitted. The results obtained have been compared with the available data.