SIMBAD references

The exact solution for any eccentricity of the Peters-Mathews model of binary stars emitting gravitational waves leads to hypergeometric functions of Appell and Gauss type. Efficient (fast and accurate) numerical computations of such functions can be obtained by using the Padé approximation. Padé rational approximants, for equal numbers of coefficients, are more accurate than Taylor series representation. Unfortunately, the algorithm computing Padé coefficients numerically from Taylor series is known to be unstable. This numerical instability lowers the accuracy of the Padé approximant. In this paper, double precision accuracy Padé approximant coefficients are computed symbolically. Advantages of the symbolic computation with respect to the numerical treatment are shown, providing accurate representations of the Appell and Gauss hypergeometric functions. The symbolic approach can be extended to higher post-Newtonian models.