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

Abstract

In this paper we study the connection between: (i) closed Newton-Cotes formulae of high order, (ii) trigonometrically-fitted and exponentially-fitted differential methods, (iii) symplectic integrators. Several one step symplectic integrators have been produced based on symplectic geometry during the last decades (see the relevant literature and the references here). However, the study of multistep symplectic integrators is very poor. In this paper we investigate the High Order Closed Newton-Cotes Formulae and we write them as symplectic multilayer structures. We develop trigonometrically-fitted and exponentially-fitted symplectic methods which are based on the closed Newton-Cotes formulae. We apply the symplectic schemes in order to solve the resonance problem of the radial Schrodinger equation. Based on the theoretical and numerical results, conclusions on the efficiency of the new obtained methods are given.