The methods used for computing the Coulomb functions described below are similar to those in §13.29.
The power-series expansions of §§33.6 and 33.19 converge for all finite values of the radii and , respectively, and may be used to compute the regular and irregular solutions. Cancellation errors increase with increases in and , and may be estimated by comparing the final sum of the series with the largest partial sum. Use of extended-precision arithmetic increases the radial range that yields accurate results, but eventually other methods must be employed, for example, the asymptotic expansions of §§33.11 and 33.21.
When numerical values of the Coulomb functions are available for some radii, their values for other radii may be obtained by direct numerical integration of equations (33.2.1) or (33.14.1), provided that the integration is carried out in a stable direction (§3.7). Thus the regular solutions can be computed from the power-series expansions (§§33.6, 33.19) for small values of the radii and then integrated in the direction of increasing values of the radii. On the other hand, the irregular solutions of §§33.2(iii) and 33.14(iii) need to be integrated in the direction of decreasing radii beginning, for example, with values obtained from asymptotic expansions (§§33.11 and 33.21).
In a similar manner to §33.23(iii) the recurrence relations of §§33.4 or 33.17 can be used for a range of values of the integer , provided that the recurrence is carried out in a stable direction (§3.6). This implies decreasing for the regular solutions and increasing for the irregular solutions of §§33.2(iii) and 33.14(iii).
Bardin et al. (1972) describes ten different methods for the calculation of and , valid in different regions of the ()-plane.
Thompson and Barnett (1985, 1986) and Thompson (2004) use combinations of series, continued fractions, and Padé-accelerated asymptotic expansions (§3.11(iv)) for the analytic continuations of Coulomb functions.
WKBJ approximations (§2.7(iii)) for are presented in Hull and Breit (1959) and Seaton and Peach (1962: in Eq. (12) should be ). A set of consistent second-order WKBJ formulas is given by Burgess (1963: in Eq. (16) should be ). Seaton (1984) estimates the accuracies of these approximations.