Sympy的2nd_power_series_ordinary ODE求解器阶项错误问题咨询

问题1：级数解阶数不一致是否为Bug #

这是Sympy幂级数求解模块的可复现Bug，核心原因如下：

• 按dsolve的幂级数接口定义，n=3时结果应该包含所有阶数低于3的项（即0次、1次、2次项），最终返回带O(s**3)的截断结果。你遇到的n=3结果缺漏s²项，是因为求解器在处理带符号参数的二阶常微分方程时，系数递推逻辑的提前终止判定有缺陷，没有完成2次项系数的计算就直接做了截断。
• 该问题在含符号参数的幂级数求解场景下出现概率较高，你可以将复现代码提交到Sympy官方Issue队列跟进修复进度。

临时规避方案：
使用比你需要的阶数高2的n值调用dsolve，拿到结果后手动截断到目标阶数即可。

问题2：不指定hint触发递归深度超限 #

这是Sympy 1.9版本迭代ODE分类器时引入的回归问题（Regression）：

• 1.8版本的分类器逻辑没有递归判定分支，不会触发该问题；1.9之后的分类器在匹配你构造的这类分子分母包含相同高阶导项的方程时，进入了无终止条件的递归判定逻辑，最终触发最大递归深度限制报错。

临时规避方案：

• 求解该类方程时始终显式指定hint="2nd_power_series_ordinary"参数
• 若需要使用自动分类功能，可临时降级到Sympy 1.8版本使用

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