  We present the yield and spectra of stable secondary particles (gamma, e^+/-,
nu_e, nubar_e, nu_mu, and nubar_mu) of p-p interaction in parameterized
formulae to facilitate calculations involving them in astronomical
environments. The formulae are derived on the up-to-date p-p interaction model
by [Kamae05] 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 Delta(1232)
and the other multiple resonances around 1600MeV/c^2. The parametrized 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.
