A method proposed earlier by the author for Fermi systems is applied to systems of interacting Bose particles. The energy of the ground state is represented in the form of expansions for weakly nonideal systems and in the form of approximating expressions for strongly nonideal systems. The results are illustrated by a one-dimensional model with delta-functional repulsive potential between the particles, for which the energy and binary distribution function are calculated.