import sympy as sp

x = sp.symbols("x")
f = x**3 - 3*x + 2
print("Polynomial:", f)
print("Derivative:", sp.diff(f, x))
print("Integral on [0, 1]:", sp.integrate(f, (x, 0, 1)))
print("Factorization:", sp.factor(f))
assert sp.expand(sp.factor(f)) == f