A Variational Symplectic Scheme Based on Lobatto’s Quadrature
摘要
We present a variational integrator based on the Lobatto quadrature for the time integration of dynamical systems issued from the least action principle. This numerical method uses a cubic interpolation of the states and the action is approximated at each time step by Lobatto’s formula. Numerical analysis is performed on both a harmonic oscillator and a nonlinear pendulum. The geometric scheme is conditionally stable, sixth-order accurate, and symplectic. It preserves an approximate energy quantity. Simulation results illustrate the performance and the superconvergence of the proposed method. [GSI 2025, 26 June 2025.]