L代码

LSODE（Livermore常微分方程求解器）求解形式为dy/dt=f（t，y）的刚性和非刚性系统。在刚性情况下，它将雅可比矩阵df/dy视为稠密（全）矩阵或带状矩阵，并将其视为用户提供的或由差商进行内部近似。在非刚性情况下采用Adams方法（预测-校正），在刚性情况下采用后向微分公式（BDF）方法（Gear方法）。线性系统的求解采用直接法（LU因子/求解）。LSODE取代了旧的GEAR和GEARB包，反映了用户界面和内部组织的完全重新设计，并对算法进行了一些改进。LSODE有单独的双精度和单精度版本，称为DLSODE和SLSODE。关于DLSODE/SLSODE用法的文档在源文件的注释行的初始块中提供，其中包括一个简单的示例。还提供演示程序（单独的双精度/单精度版本）。