-xy^2(x+y)^3+x(x+y)^2
=x(x+y)^2[(-y^2(x+y)+1)]
=x(x+y)^2(1-xy^2-y^3)