A matrix method for finding the energy eigenvalues of one-dimensional systems is given. The method utilizes the correspondence principle to find the separation of every two successive eigenvalues and the uncertainty principle to find the lowest energy. The harmonic oscillator and the infinite square well are solved exactly, and the WKB approximation for bound state energies is deduced.