p = polynom([1, 1]);
q = polynom([1, 2]);
r = (p+1)*(q+2)
[a, b] = r/[1, 1]