A discretization method referred to as the Quadrature Discretization Method (QDM) is introduced for the solution of the Schrodinger equation. The method has been used previously for the solution of Fokker-Planck equations. The Fokker-Planck equation can be transformed to a Schrodinger equation with a potential of the form that occurs in supersymmetric quantum mechanics, For this class of potentials, the groundstate wave function is known. The QDM is based on the discretization of the wave function on a grid of points that coincide with the points of a quadrature. The quadrature is based on a set of nonclassical polynomials orthogonal with respect to a weight function determined by the potential function in the Schrodinger equation. For the Fokker-Planck operator, the weight function that provides rapid convergence of the eigenvalues are the steady distributions at infinite time, that is, the ground state wave functions. In the present paper, the weight functions used in an analogous solution of the Schrodinger equation are related to the ground state wave functions if known, or some approximate form. Calculations are carried out for a model systems, the Morse potential, and for the vibrational levels of O-2 and Ar-Xe with realistic pair potentials. For O-2, the wave functions are used to calculate the vibrationally inelastic transition amplitudes for a Morse potential and compared with exact analytic results. The eigenvalues of a two-dimensional Schrodinger equation with the Henon-Heiles potential are also calculated. The rate of convergence of the eigenvalues and the eigenfunctions of the Schrodinger equation is very rapid with this approach.