Make the (x+1)^2 expansion work (fixes #38)

We just needed to take into account the overall coefficient in Mul and Add.
A test was added.