In this article we study BoseEinstein condensation of the gas with interactions in a harmonic trap using a numerical method. The distribution of the energy states was given, by which the author calculated the shift of the critical temperature and made a comparison with the situation without interactions. We have obtained a positive shift of critical temperature for systems with repulsive interaction, and a negative one for systems with attractive interaction.