A High-Order Two-Step Phase-Fitted Method for the Numerical Solution of the Schrödinger Equation