用数值方法求解非线性不可分的ODEs二阶系统？

Question

我正在尝试对非常特殊的泛函进行数值求解Euler-Lagrange方程，但是不幸的是，所得微分方程的形式似乎无法操纵为可输入ODE求解器的形式。在以下等式中，A（t）是n维向量，Q是nxn维矩阵。

有人知道一个能够处理这种计算的ODE求解器吗？

谢谢！