added a check that expand(x^((a+b)^2-a^2-b^2-2*a*b)) -> 1