integrate(dx/(x + sqrt(x^2 + x + 1)))