solve y' = y + y^2 + y^3 + 1