x = sympy.Symbol('x')
poly1 = (x + 1) * (x -2)
result1 = sympy.expand(poly1)
x = sympy.Symbol('x')
poly2 = x**2 - 5*x + 6
result2 = sympy.factor(poly2)
x = sympy.Symbol('x')
poly3 = x**2 + 2*x + 1
result3 = poly3.subs(x,5)
x = sympy.Symbol('x')
f = x**4 - 3*x**2 + 5
result4 = sympy.diff(f)
x = sympy.Symbol('x')
f = 2*x + 1
int_f = sympy.integrate(f)
area = sympy.integrate(f,(x,1,3))