Data on the numerical solution of a system of kinetic Boltzmann equations for a homogeneous multicomponent mixture of reacting gases with molecules of different “colors” that change in the reactions are given. The solution is obtained using a well-known version of the direct statistical simulation (Monte-Carlo) method, namely, the majorant frequency method, under conditions when the molecules belonging to the high-velocity “tails” of the corresponding distribution functions enter into the color change reaction. The properties of the numerical solution are compared with solutions obtained within the framework of the usual perturbation methods. It is shown that to obtain correct solutions over the range of threshold molecular velocities it is necessary to modify substantially the procedure of the perturbation method, while the traditional approach can be used only on the range of thermal particle velocities. Earlier, this was definitely established for distributions of the reacting molecules over their internal degrees of freedom and for the distributions of reactant-molecules participating in a high-threshold reaction.