>>> np.poly([1,-2,3])
array([ 1, -2, -5, 6])
>>> np.roots([ 1, -2, -5, 6])
array([-2., 3., 1.])
f(x) = x3 – 2x2 -5x + 6