Simulation of Proton Double Differential Cross 113, 256, 597 Induced Neutron-production Sections on Fe, Pb, Al at and 800 MeV