load(atrig1)$load(ntrig)$

f(x) := sqrt(x^2+x+1)/(x+1);

integrate(f(x),x);

trigsimp(%);