To see that it is the same, let's compute the integral with hyperbolic substitution:
[tex]I=\int{C}osh^{2}(v)dv=Sinh(v)Cosh(v)-\int{S}inh^{2}(v)dv=Sinh(v)Cosh(v)-\int(Cosh^{2}(v)-1)dv=[/tex]
[tex]Sinh(v)Cosh(v)+v-\int{C}osh^{2}(v)dv=Sinh(v)Cosh(v)+v-I[/tex]
That is, we have the expression:
[tex]I=Sinh(v)Cosh(v)+v-I[/tex]
yielding:
[tex]I=\frac{1}{2}Sinh(v)Cosh(v)+\frac{1}{2}v[/tex]
Now, we have [tex]v=Sinh^{-1}(u)[/tex]
Using the identity:
[tex]Cosh^{2}(y)-Sinh^{2}(y)=1[/tex]
we get [tex]Cosh(Sinh^{-1}(u))=\sqrt{1+u^{2}}[/tex]
since we have, of course [tex]Sinh(Sinh^{-1}(u))=u[/tex]
Thus, our first term may be rewritten as:
[tex]\frac{1}{2}Sinh(v)Cosh(v)=\frac{u}{2}\sqrt{1+u^{2}}[/tex]
Using the exponential representation of the the hyperbolic sine, you should have little trouble representing the inverse function of hypsine in terms of the natural logarithm.
That constitutes the second term in the formula you found.