python - 使用 Numpy 查找函数 a x^n + bx - c = 0 的根，其中 n 不是整数？

Question

我正在用 python 编写一个程序，我需要在其中找到一个函数的根:

a*x^n + b*x -c = 0 

其中 a 和 b 是程序中较早计算的常量，但它们有几千个。我需要对 a 和 b 的所有值重复此等式两次，一次是 n = 77/27 一次是 n = 3。

我如何在 python 中执行此操作？我检查了 numpy.roots(p)，我认为这适用于 n = 3 的情况。但是对于 n = 77/27 我怎么能做到呢？

Answer 1

我认为你的野兽选择是scipy.optimize.brentq() :

def f(x, n, a, b, c):
    return a * x**n + b * x - c

print scipy.optimize.brentq(
    f, 0.0, 100.0, args=(77.0/27.0, 1.0, 1.0, 10.0))

打印

2.0672035922580592