High-order integration formulas for the numerical integration of the one-dimensional Schrodinger equation have been generated by means of the symbolic computational language MAPLE. The formulas have been used to evaluate energies of the bound states for a series of model systems. The calculations have revealed that they possess high stability and accuracy. When performing calculations within the framework of MAPLE, practically arbitrary accuracy is attainable with a modest amount of computing time.