Automatic integration of Euler-Lagrange equations with constraints

摘要

Numerical difficulties in the integration of Euler-Lagrange and similar equations are discussed. A technique for reducing their index from three to two is introduced and it is shown that variable-order, variable-step BDF methods converge for these index two problems. The practical application of this reduction in a numerical setting is examined.