We propose a new method for the energy eigenvalue problem for boson systems with either weak or strong two body interactions. It is also shown that this method is particularly simple and effective for obtaining the explicit analytical expressions of energy eigenvalues and eigenstates for interacting harmonic oscillators especially whenthe number of oscillators is huge.