Python数值计算（9）

继续接着前面多项式的操作来说。

4. 多项式的计算

除了进行值替换，多项式之间可以进行代数的基本运算，如加减乘除等，甚至还可以进行微分和积分等常见运算。

分别使用传统的运算法就可以了，唯一需要注意的是除法，应该使用//，而不是/，测试如下：

a=P([1,-2,0,3])
b=P([-2,1])
print(a-b) # 3.0 - 3.0·x + 0.0·x² + 3.0·x³
print(a+b) # -1.0 - 1.0·x + 0.0·x² + 3.0·x³
print(a*b) # -2.0 + 5.0·x - 2.0·x² - 6.0·x³ + 3.0·x⁴
print(a//b) # 10.0 + 6.0·x + 3.0·x²

通常两个多项式的除法不大可能刚好整除，因此，会有余项，计算余项很简单，使用%运算符，例如上例中：

print(a%b) # 21.0

由于b是一个一次式，因此余数为一个常数。

在需要同时获取商和余式的时，最好的做法是使用divmod函数：

quo, rem = divmod(a,b)

甚至还可以使用** 计算多项式的幂：

a=P([1,-2,0,3])
print(a*a)
print(a**2)
# both display: 1.0 - 4.0·x + 4.0·x² + 6.0·x³ - 12.0·x⁴ + 0.0·x⁵ + 9.0·x⁶

积分和微分，使用函数integ和deriv，积分和微分之后仍旧是一个多项式：

a=P([1,-2,0,3])
print(a.integ()) # 0.0 + 1.0·x - 1.0·x² + 0.0·x³ + 0.75·x⁴
print(a.deriv()) # -2.0 + 0.0·x + 9.0·x²