We have used 2 libraries here:
# [] in the below function is x^(n)'s coeffs for n=0 ... N
p = Polynomial( [1,3,5] )
#
q = fromroots( [1,-4,4] )
p(x) # evaluate for x
NOTE: Basic maths operation work without “.”:
p + q # works
2*p - q # workds
derivative(p)
integrate(p)
root(q) # find root of poly_type var