用python求解负幂或非整数幂多项式方程

Question

我试着解这个方程

1600 = 0.41 + 6.31*d**-1.54 + 2.42*d**-3 

Mathematica可以在不到1秒内给我d=6.4673。但我无法用python符号计算得到答案。在

使用“解”的sympy采取永远。有没有办法用python符号计算来解这个方程？看来问题主要来自于非整数负幂。在

Answer 1

SymPy通过mpmath库提供数值计算；这包括通过nsolve进行数值寻根。在本例中，由于denomominator中有一个d，因此我们按照nsolve的docstring的建议，使用表达式的分子并给出一个初始猜测。很快就会找到已经被引用的同一个根：

>>> f
-6.31*d**(-1.54) + 1599.59 - 2.42/d**3
>>> nsolve(f.as_numer_denom()[0], 1)
0.119100503944930