We recast the well-known Numerov method for solving Schr€odinger’s equation into a representation of the kinetic energy operator on a discrete lattice. With just a few lines of code in a high-level programming environment such as MATHEMATICA, it is simple to calculate and plot accurate eigenvalues and eigenvectors for a variety of potential problems. We illustrate the method by calculating high-accuracy solutions for the jxj potential.