[en] We have obtained a simple representation to the observed invariant cross-section for the production of neutral pions in proton-proton collisions. Making use of this representation, we have calculated the differential and integral production spectra of gamma rays in the Galaxy from interactions of cosmic ray nuclei with interstellar gas. It is shown that the uncertainties in deducing interstellar proton spectrum by demodulating the observed spectrum do not affect very much the gamma ray spectrum. We have also determined the gamma ray production spectrum through bremsstrahlung process for a typical interstellar electron spectrum derived from the radio spectrum in the Galaxy. From these, the total gamma ray production spectrum resulting from the interaction of cosmic rays with interstellar matter is compared with the observed gamma ray spectrum in the Galaxy and some inferences have been obtained. We also point out the possible uncertainty in the present calculation and suggest the improvements needed. (orig.)