Position-dependent mass Schrödinger equations allowing harmonic oscillator (HO) eigenvalues