In the report, the algorithm of the Kinetic Force Method founded on a kinetic equation for auxiliary two-particle distribution function of quasiparticle pairs is considered. It is known that in the interaction of quasiparticle pairs, energy and momentum are conserved for each individual pair. However, in the course of the numerical simulation of relaxation of the entire ensemble, the energy deficit still appears. In the report, the new algorithm of the Kinetic Force Method is proposed: in this algorithm the energy of a system of molecules is conserved without any additional corrections. Thus, the advantages of modeling on the base of two-particle distribution function of quasiparticle pairs are used more effectively.