反函数纽比·波利瓦尔()

Question

我想知道是否有一个方便的反函数np.polyval公司（），其中我给出y值，它解出x？在

我知道我可以这样做的一个方法是：

import numpy as np

# Set up the question
p = np.array([1, 1, -10])
y = 100

# Solve
p_temp = p
p_temp[-1] -= y
x = np.roots(p_temp)

不过，我想大多数人都会同意这段代码的可读性很低。有什么建议吗？在

Answer 1

这样的怎么样？在

In [19]: p = np.poly1d([1, 1, -10])  # Use a poly1d to represent the polynomial.

In [20]: y = 100

In [21]: (p - y).roots
Out[21]: array([-11.,  10.])

poly1d对象实现算术运算以返回新的poly1d对象，因此p - y是一个新的poly1d：

^{pr2}$

poly1d的roots属性返回预期值。在