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