In this paper, we develop a consistent extension of RPA for bosonic systems. In order to illustrate the method, we consider the case of the anharmonic oscillator. We compare our results with those obtained in mean-field and standard RPA approaches, with the exact ones and show that they are very close to the exact ones.