The bound state solutions to the radial Schrödinger equation are obtained in three-dimensional space using the series expansion method within the framework of a general interaction potential. The energy eigenvaluesof the pseudoharmonic and Kratzer potentials are given as special cases. The obtained analytical results are applied to several diatomic molecules, i.e. $\rm{N}_{2}$,CO,NO and CH. In order to check the accuracy of the present method, a comparison is made with similar results obtained in the literature by using other techniques.
