solve_ivp中LSODA算法生成dense_output采用的是何种算法？

• BDF求解器：采用与当前求解阶数匹配的可变阶赫米特多项式插值，插值精度和求解步的计算精度完全对齐，不会额外引入精度损耗。该逻辑完全沿用了ODEPACK库中原始BDF求解器的配套插值实现。
• LSODA求解器：作为自适应切换求解格式的刚性/非刚性兼容求解器，插值方法随当前激活的求解格式动态调整：
  • 处理非刚性区间、激活Adams线性多步法时，采用对应阶数的Adams插值多项式
  • 处理刚性区间、激活BDF格式时，采用和独立BDF求解器完全一致的可变阶赫米特多项式插值

scipy官方文档没有标注这两个求解器的插值方法，是因为二者都是从经典Fortran ODEPACK库移植的成熟实现，文档仅专门标注了后续自行实现的RK45、Radau等求解器的特殊插值逻辑。

内容的提问来源于stack exchange，提问作者noir1993