The evolution equations for Green functions of the Reggeon Field Theory (RFT) are equivalent to those of the inclusive distributions for the reaction–diffusion system of classical particles. We use this equivalence to obtain numerically Green functions and amplitudes of the RFT with all loop contributions included. The numerical realization of the approach is described and some important applications including total and elastic proton–proton cross sections are studied. It is shown that the loop diagram contribution is essential but can be imitated in the eikonal cross section description by changing the Pomeron intercept. A role of the quartic Pomeron coupling which is an inherent part of the stochastic model is shown to be negligible for available energies.