Python调用SymPy求解方程输出结果与预期不符应如何调整？

问题原因 #

• 参数p取值错误：你需要代入的是概率值0.6而非填写的60，该公式的p为0~1区间的概率变量，代入60会导致所有计算项数值异常放大，得到完全不符合预期的负数结果。
• 代码冗余：当前方程中x直接等于右侧表达式，不需要调用sympy的solve函数解方程，直接计算右侧表达式即可得到结果，多余的求解操作没有实际意义。

修正后代码 #

p = 0.6
x = (p**4)*((15-(4*p))-((10*(p**2))/(1-((2*p)*(1-p)))))
print(x)

运行结果 #

运行上述代码输出为0.7357292307692305，和你预期的0.735729231完全匹配。

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