As is well known, the algorithm of the kinetic force method, based on the kinetic equation for the two-particle distribution function, guarantees conservation of energy and momentum during a "quasicollision" of each separate pair of quasiparticles. But in the modeling of relaxation of the entire gas an energy deficit nevertheless appears at the step of reconstructing the one-particle distribution function. A new variant of the algorithm of the kinetic force method is proposed in which the energy of the system of molecules is conserved at any step of the modeling of relaxation without additional corrections. By virtue of this, the advantages of modeling on the basis of the two-particle distribution function of pairs of quasiparticles are used more efficiently.