In the first step we try to guess the result basing on our personal experience
and check if the derivative gives back the integrand. If we miss, we can correct our guess and see what happens. In case we are unable to make a good guess, we try to simplify the integrand. In our case
$$x^2+x=\left (x+{1\over 2}\right)^2-{1\over 4}={1\over 4}[(2x+1)^2-1]$$
It seems convenient to perform the substitution $u=2x+1.$ In this way we get
$${1\over 2}\int {\sqrt{u^2-1}\over u-1}\,du$$ The integrated function is defined on two disjoint intervals $u<-1$ and $u>1.$ For $u>1$ it is tempting to make the substitution $u=\cosh t={1\over 2}(e^t+e^{-t}),$ for $t>0,$ having in mind that $\cosh ^2t-\sinh^2 t=1$ (where $\sinh t={1\over 2}[e^t-e^{-t}]).$
Then $du=\sinh t\,dt$ and we end up with the integral
$${1\over 2}\int {\sinh ^2t\over \cosh t-1}\,dt={1\over 2}\int {\cosh^2t-1\over \cosh t-1}\,dt ={1\over 2}\int [\cosh t+1]\,dt={1\over 2}\sinh t +{1\over 2}t$$ Next we have to get back to the original variable $x.$ Solving $$u=\cosh t ={1\over 2}[e^t+e^{-t}],\qquad t>0$$ leads to the quadratic equation for $e^t$
$${1\over 2}\,(e^{t})^2-u\,e^t+{1\over 2}=0$$ from which we get $$e^t=u+\sqrt{u^2-1},\quad \sinh t=\sqrt{u^2-1},\quad t=\log (u+\sqrt{u^2-1}) $$
Remark There are two roots of the equation and their product is equal $1.$ We choose the greater greater root, as $e^t>1$ for $t>0.$
Eventually, as $u=2x+1>1,
$ i.e. $x>0,$ we obtain $$\int {\sqrt{x^2+x}\over x}\,dx = {1\over 2}\log(2x+1+\sqrt{4x^2+4x})+{1\over 2}\sqrt{4x^2+4x}\\ ={1\over 2}\log (2x+1+2\sqrt{x^2+x})+\sqrt{x^2+x}\qquad \qquad (*)$$
It remains to consider the case left behind, when $u<-1,$ i.e. $x<-1.$ But since formula $(*)$ holds for $x>1$ and $2x+1+2\sqrt{x^2+x}<0,$ for $x<-1,$
we get the final result
$$\int {\sqrt{x^2+x}\over x}\,dx ={1\over 2}\log \mid 2x+1+2\sqrt{x^2+x}\mid+\sqrt{x^2+x}\qquad x\notin [-1,0]$$
Remark Mind that if we are looking for a particular antiderivative $F(x),$ satisfying $F(1)=a $ and $F(-2)=b,$ we have to add to the solution above different constants, depending on the cases $x<-1$ or $x>0.$ Therefore I resist from writing $\int f(x)\,dx =F(x)+C,$ when $f(x)$ is defined on several disjoint open intervals. The formula $\int f(x)\,dx =F(x)$ is just equivalent to $F'(x)=f(x).$ For example $\int 2x\,dx =x^2$ and also $\int 2x\,dx =x^2+1,$ but it does not imply $0=1.$