Exponentially fitted open Newton-Cotes differential methods as multilayer symplectic integrators.

ABSTRACT

Classical open and closed Newton-Cotes differential methods possessing the characteristics of multilayer symplectic structures have been constructed in the past. In this paper, we study the exponentially fitted open Newton-Cotes differential methods of order two, four, and six. It is shown that these integrators, just as their classical counterparts, preserve the volume in the phase space of a Hamiltonian system. They can be converted into a multilayer symplectic structure so that volume-preserving integrators of a Hamiltonian system are obtained. A numerical example has been carried out to show the effectiveness of the present differential method.