Closed Newton–Cotes trigonometrically-fitted formulae of high order for the numerical integration of the Schrödinger equation