Nuclear reactions govern major aspects of the chemical evolution of galaxies and stars. Analytic study of the reaction rates and reaction probability integrals is attempted here. Exact expressions for the reaction rates and reaction probability integrals for nuclear reactions in the cases of nonresonant, modified nonresonant, screened nonresonant and resonant cases are given. These are expressed in terms of H-functions, G-functions and in computable series forms. Computational aspects are also discussed.