SymPy化简问题：代入c2=0后结果含无穷大，如何进一步化简？

问题原因 #

你用SymPy求解三次方程得到的解析解基于卡丹公式，当直接代入c2=0时，表达式中会出现0/0型不定式，SymPy直接计算会返回无穷大，但实际上这个不定式的极限就是原方程退化为线性方程时的解x1=F_C/c1。

解决方案 #

以下是两种可靠的处理方式：

方法1：计算c2趋近于0的极限

通过极限运算处理不定式，得到正确的化简结果：

import sympy as sp
c1, c2, x1, F_C = sp.symbols("c1, c2, x1, F_C", real=True)
F_C_ex = c1*x1 + c2*x1**3
sol = sp.solve(F_C_ex - F_C, x1)[0]

# 计算c2→0时的极限
correct_sol = sp.limit(sol, c2, 0)
print(correct_sol)  # 输出 F_C/c1

方法2：先代入c2=0再求解

如果只需要c2=0时的解，可以先将参数代入原方程，再求解线性方程，避免处理三次方程的复杂解析解：

import sympy as sp
c1, c2, x1, F_C = sp.symbols("c1, c2, x1, F_C", real=True)
F_C_ex = c1*x1 + c2*x1**3

# 先替换c2=0，再解方程
sol_c2_zero = sp.solve(F_C_ex.subs(c2, 0) - F_C, x1)[0]
print(sol_c2_zero)  # 输出 F_C/c1

补充说明 #

直接subs(c2, 0)失效是因为三次方程的解析解结构中，c2同时出现在分子和分母的项中，形成0/0型的不定式。SymPy的subs只会做直接替换，不会自动处理不定式，而limit函数可以正确计算这种情况下的极限值。

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