dsolve (y'(x) + x y(x) = x^2 (x+1), y, x)