A MATLAB routine to compute numerically the solution to the 1D stationary Schrödinger equation in the highly oscillatory regime.