The quality and reliability of the computational simulation of a macroscopic nuclear device are directly related to the quality of the underlying basic nuclear data. To meet these needs, according to advanced nuclear models that account for details of nuclear structure and the quantum nature of nuclear reaction and the experimental data of total, nonelastic, and elastic scattering cross sections, and elastic scattering angular distributions of Pb and its isotopes, all cross sections of neutron-induced reaction, angular distributions, energy spectra, especially the double-differential cross sections for neutron, proton, deuteron, triton, helium, and alpha emissions are calculated and analyzed for n + 204,206,207,208,natPb at incident neutron energies below 20 MeV by using the UNF nuclear model code. At neutron incident energies 20 < En 250 MeV, MEND codes are used. Theoretical calculations are compared with existing experimental data and other evaluated data from ENDF/B-VI and JENDL-3.