Runge Kutta DAE 求解器：微分代数系统可以使用经过调整的 Runge Kutta 方法求解。-matlab开发

这是一种适用的 Runge Kutta 方法（理论背景，例如在 Hairer、Lubich、Roche“通过runge-kutta 方法对微分代数系统的数值解”）来解决 DAE。 已经实现的是第 1、2 和 3 阶段的 Radau II A 方法，但基本上每个 Butcher 画面都可以实现（详见代码）。 提供了一个简单的例子（数学钟摆）来说明用法。 非线性系统求解器是牛顿法，但也可以互换。