import sympy as sym

# Defining symbols:

x, y, a, b = sym.symbols('x, y, a, b')

# Factorization:

E_fact = sym.factor(x**2 + 2*x, x)
print('E_fact =', E_fact)

# Expansion:

E_exp = sym.expand((x + a)*(x + b), x)
print('E_exp =', E_exp)

# Simplification:

E_simp = sym.simplify(y*(sym.sin(x)**2 + sym.cos(x)**2)) 
print('E_simp =', E_simp)

# Partial decomposition:

E_part = sym.apart(2*x/((x + 1)*(x + 2)))
print('E_part =', E_part)