A variational calculation of the ground state and of the first excited state of a bound polaron is given for all values of the electron-phonon coupling and of the Coulomb potential. The applied unitary transformation formalism is adapted in order to take account of correlation between the longitudinal optical phonons.