We compute the probability of producing n particles from few colliding particles in the unbroken (3 + 1)-dimensional λϕ 4 theory. To this end we numerically implement semiclassical method of singular solutions which works at n ≫ 1 in the weakly coupled regime λ ≪ 1. For the first time, we obtain reliable results in the region of exceptionally large final-state multiplicities n ≫ λ −1 where the probability decreases exponentially with n, $P\left(\mathrm{few}\to n\right)~exp\left(f_{\infty }\left(\epsilon \right)n\right)$ $$ \mathcal{P}\left(\textrm{few}\to n\right)\sim \exp \left\{{f}_{\infty}\left(\varepsilon \right)n\right\} $$ , and its slope f ∞ < 0 depends on the mean kinetic energy ε of produced particles. In the opposite case n ≪ λ −1 our data match well-known tree-level result, and they interpolate between the two limits at n ~ λ −1. Overall, this proves exponential suppression of the multiparticle production probability at n ≫ 1 and arbitrary ε in the unbroken theory. Using numerical solutions, we critically analyze the mechanism for multiple Higgs boson production suggested in the literature. Application of our technique to the scalar theory with spontaneously broken symmetry can eradicate (or confirm) it in the nearest future.