Abstract
We present a novel numerical method for solving ODEs while preserving polynomial first integrals. The method is based on introducing multiple quadratic auxiliary variables to reformulate the ODE as an equivalent but higher-dimensional ODE with only quadratic integrals to which the midpoint rule is applied. The quadratic auxiliary variables can subsequently be eliminated yielding a midpoint-like method on the original phase space. The resulting method is shown to be a novel discrete gradient method. Furthermore, the averaged vector field method can be obtained as a special case of the proposed method. The method can be extended to higher-order through composition and is illustrated through a number of numerical examples.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
