上传者: 38576045
|
上传时间: 2022-03-17 16:20:42
|
文件大小: 20KB
|
文件类型: -
这是一种适用的 Runge Kutta 方法(理论背景,例如在 Hairer、Lubich、Roche“通过runge-kutta 方法对微分代数系统的数值解”)来解决 DAE。 已经实现的是第 1、2 和 3 阶段的 Radau II A 方法,但基本上每个 Butcher 画面都可以实现(详见代码)。 提供了一个简单的例子(数学钟摆)来说明用法。 非线性系统求解器是牛顿法,但也可以互换。