The quasiexactly solvable potential method is used to determine the energies and the corresponding exact eigenfunctions for a system of N particles with equal mass interacting via an anharmonic potential. For systems with five and seven particles, we compute the ground state and the first excited state only, and compare the spectrums with the results obtained by Ritz approximation method.
