In this work we consider exponentially fitted and trigonometrically fitted Runge–Kutta–Nyström methods. These methods integrate exactly differential systems whose solutions can be expressed as linear combinations of the set of functions exp( w x ) , exp(− w x ) , or sin( w x ) , cos( w x ) , w . We modify existing RKN methods of fifth and sixth order. We apply these methods to the computation of the eigenvalues of the Schrödinger equation with different potentials as the harmonic oscillator, the doubly anharmonic oscillator and the exponential potential.