摘要

首先对带约束动力学中的辛算法作了改进,利用吴消元法求解多项式类型Euler-Lagrange方程.在辛算法的基础上,根据线性方程组理论和相容条件提出了一个求解约束的新算法.新算法的推导过程比辛算法严格,而且计算也比辛算法简单,并且多项式类型的Euler-Lagrange仍可以用吴消元法求解.另外,对于某些非多项式类型的Euler-Lagrange方程,可以先化为多项式类型,再用吴消元法求解.利用符号计算软件,上述算法都可以在计算机上实现.