  We develop a variational method to obtain accurate bounds for the
eigenenergies of H = -Delta + V in arbitrary dimensions N>1, where V(r) is the
nonpolynomial oscillator potential V(r) = r^2 + lambda r^2/(1+gr^2), lambda in
(-infinity,\infinity), g>0. The variational bounds are compared with results
previously obtained in the literature. An infinite set of exact solutions is
also obtained and used as a source of comparison eigenvalues.
