We present the yield and spectra of stable secondary particles (γ, e±, νe, ν̄e, νμ, and v̄μ of p-p interaction in parameterized formulae to facilitate calculations involving them in astronomical environments. The formulae are derived from the up-to-date/p-p interaction model by Kamae et al. (2005), which incorporates the logarithmically rising inelastic cross section, the diffraction dissociation process, and the Feynman scaling violation. To improve fidelity to experimental data in lower energies, two baryon resonance contributions have been added: one representing Δ(1232) and the other representing multiple resonances around 1600 MeV/c-2. The parameterized formulae predict that all secondary particle spectra be harder by about 0.05 in power-law indices than that of the incident proton and their inclusive cross sections be larger than those predicted by p-p interaction models based on the Feynman scaling.