Abstract
The energy levels and wave functions of the Schrödinger equation involving the potential x2+λx2/(1+gx2) are calculated by the variational method, for any range of λ and g, without having to resort to numerical quadrature. Using properly scaled (in λ and g) harmonic oscillator functions as a basis set, an easy to compute analytical expression of the current Hamiltonian matrix element is derived. Perturbative results are also given.