In this paper, we propose a discrete model for the quantum harmonic oscillator. The eigenfunctions and eigenvalues for the corresponding Schrödinger equation are obtained through the factorization method. It is shown that this problem is also connected with the equation for Meixner polynomials.
