deSolve R包中的ode()函数使用的默认数字方案是lsode方法,该方法实现BDF和隐式Adams线性多步方案来求解ODE系统。集成以可变步长完成,该步长通过估计每个步骤的本地截断误差(LTE)来控制。实际上,以使得LTE的估计小于某个预设值rtol和atol的方式来调整步长。 ode方法中的默认值为:
defaultValue
并且可以适应。
但是,实际评估的LTE总是高于这些值。有没有办法从求解器中提取此值?
下面是一个示例系统:
ode(rtol = 1e-6, atol = 1e-6)