  In this paper, we introduce a novel analytical solution to
Tolman-Oppenheimer-Volkoff (TOV) equation, which is ultimately a hydrostatic
equilibrium equation derived from the general relativity in the framework of
relativistic isothermal spheres. Application of the traditional power series
expansions on solving TOV equation results in a limited physical range to the
convergent power series solution. To improve the convergence radii of the
obtained series solutions, a combination of the two techniques of Euler-Abel
transformation and Pade approximation has done. The solutions are given in
\exi-\theta and \exi-\mu phase planes taking into account the general
relativistic effects \sigma= 0.1, 0.2 and 0.3. An Application to a neutron star
has done. A Comparison between the results obtained by the suggested approach
in the present paper and the numerical one indicates a good agreement with a
maximum relative error of order 10-3, which establishes the validity and
accuracy of the method. The procedure we have applied accelerated the power
series solution with about ten times than of traditional one.
