Question:
I wish to find the inverse of the following function:
$f(x)=\frac{1}{2}\left(\arctan(x) + \ln\left(\sqrt{\frac{1+x}{1-x}}\right)\right)$
This is the equation for a radial null geodesic in a particular version of of the $AdS_2$-Schwarzschild metric. I want to invert the above function. I can plot as a parametric function $(f(x),x)$ (which I cannot post because of my non-existent reputation). It is monotonically increasing and bounded on $[0,\infty)$. I am sure that an inverse function exists but cannot be written in terms of elementary functions. What I want to obtain is at least the exact coefficients of the power series for the inverse function (an analytic expression for all coefficients would be nice but I am not that hopeful).
Attempt:
I tried to compute the coefficients of the power series of the inverse function using the Lagrange inversion theorem. I was able to find the first $100$ non-zero coefficients. However, the resulting power series was only convergent on $[0,1]$. I want a power series that converges over a larger interval. This should be possible since the inverse of $f(x)$ as plotted parametrically is bounded and without singularities on $[0,\infty)$. Is there a way to find such a power series?