High order closed Newton-Cotes exponentially and trigonometrically fitted formulae as multilayer symplectic integrators and their application to the radial Schrödinger equation