There are essentially two ways to approach this. The first is
The other is to use partial fractions:
Now, tanh^(-1) and ln have an overlapping domain and in this region everything is just fine and C = 0. The problem is that the ln solution has a much larger domain than the tanh one.
How can two solutions of an integral have different domains?