This paper reports on the numerical implementation of a convergent scheme solving a general class of nonlinear Boltzmann-like equations for reacting fluids. The scheme is tested on a model of space-homogeneous multicomponent gas with binary chemical reactions. Although the model cannot be solved exactly, it provides analytical expressions for certain low-order moments of the one-particle distribution functions. Computing these moments by applying the numerical scheme, one obtains a good agreement with the values obtained analytically. An error analysis is performed and error control considerations complete the theoretical support of the scheme developed in previous works.