scipy处理多项式方法

# -*- coding:utf-8 -*-
from scipy import poly1d
import numpy as np

p = poly1d([1, 2, 3])
print (p)
print (p*p)

print (p.integ(k=1)) #求p(x)的不定积分，指定常数项为1
print (p.deriv())    #求p(x)的一阶导数
p([4, 5])    #计算每个值代入p(x)的结果

def addsubtract(a, b):    # 按照原始定义，仅接受可比较的数字作为参数
    if a > b:
        return a - b
    else:
        return b - a
vec_addsubtract = np.vectorize(addsubtract) #将标运算转换为向量运算。
print (vec_addsubtract([0, 3, 6, 9], [1, 3, 5, 7]))