A Runge–Kutta type four-step method with vanished phase-lag and its first and second derivatives for each level for the numerical integration of the Schrödinger equation