我正在编写一个用于计算牛顿多项式的代码片段。在循环期间for，我需要得到 3 个术语（每个术语都必须在括号中）：
(x - 30) (x - 20) (x - 10) + (x - 20)*(x - 10) + (x - 10)
当我浏览代码时，我得到了这个：
x + (x - 30) (x - 20) (x - 10) + (x - 20)*(x - 10) - 10
拥有 3 个术语对我来说很重要，因为将来在另一个周期中，它们中的每一个都将乘以系数。
如何实现？
编码：
from sympy import symbols
x = symbols('x')
dig_list = [10, 20, 30]
poly_0 = x-dig_list[0]
summa = poly_0
for i in range(1, len(dig_list)):
poly_n = poly_0 * (x-dig_list[i])
summa += poly_n
poly_0 = poly_n
print(summa)
使用 sympy，我得到了一个函数的导数：
>>> import sympy as sp
>>> x = sp.Symbol('x')
>>> f = -0.01 * sp.sin(0.05 * x - 0.7) + 0.01 * sp.sin(6.2 * x + 4.5)
>>> f.diff(x)
−0.0005cos(0.05𝑥−0.7)+0.062cos(6.2𝑥+4.5)
我构建了一个函数及其导数的图表：
import matplotlib.pyplot as plt
import numpy as np
x = np.arange(100)
plt.plot(x, -0.01 * np.sin(0.05 * x - 0.7) + 0.01 * np.sin(6.2 * x + 4.5), 'b-')
plt.twinx()
plt.plot(x, -0.0005 * np.cos(0.05 * x - 0.7) + 0.062 * np.cos(6.2 * x + 4.5), 'g-')
plt.plot(plt.xlim(), [0, 0], 'k--')
plt.show()
导数为零的地方似乎不是极值。我不明白为什么。