In this work, we investigate the role of interactions in the process of thermalization of a weakly interacting Bose gas. The system of kinetic equations based on the 'Fermi's golden rule' is solved numerically using special transformation for calculation efficiency. We study the distribution function for particles in various conditions, including interaction with phonon subsystem, i.e. energy exchange with thermal bath. The possibility to achieve the state of Bose-Einstein condensation with specific values of parameters, is also discussed.