A Rational Approximant for the Digamma Function

Ernst Joachim Weniger

DOI: https://doi.org/10.1023/A%3A1025517617217

2003-06-20

Abstract:Power series representations for special functions are computationally satisfactory only in the vicinity of the expansion point. Thus, it is an obvious idea to use instead Padé approximants or other rational functions constructed from sequence transformations. However, neither Padé approximants nor sequence transformation utilize the information which is avaliable in the case of a special function -- all power series coefficients as well as the truncation errors are explicitly known -- in an optimal way. Thus, alternative rational approximants, which can profit from additional information of that kind, would be desirable. It is shown that in this way a rational approximant for the digamma function can be constructed which possesses a transformation error given by an explicitly known series expansion.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to construct a more effective rational approximation method to calculate the digamma function (also known as the psi function). Specifically, the author points out that the power series representation converges slowly or diverges far from the expansion point, resulting in its inefficiency in numerical calculations. Therefore, it is necessary to find alternative rational approximation methods.

### Problem Background
1. **Limitations of Power Series**
- The power series representation is very important in the calculation of special functions, but it converges slowly or diverges far from the expansion point.
- This makes it difficult to calculate special functions efficiently and reliably relying only on power series.

2. **Deficiencies of Padé Approximation and Other Sequence Transformations**
- Although Padé approximation is commonly used, it does not fully utilize the known power series coefficients and truncation error information.
- Other sequence transformation methods (such as Levin transformation) are more effective than Padé approximation, but there is still room for improvement.

### Goals of the Paper
- **Optimize Rational Approximation**: By using all known power series coefficients and truncation error information, construct a new rational approximation method to improve the efficiency and accuracy of calculating special functions (especially the digamma function).
- **Explicit Expressions and Error Estimates**: Ensure that the new method can not only provide better approximation results but also give an explicit truncation error expression, thereby further accelerating convergence.

### Solution
The paper proposes a new rational approximation method based on sequence transformation, especially for the calculation of the digamma function \( \psi(z) \). This is achieved through the following steps:

1. **Power Series Decomposition**: Decompose the power series representation of the digamma function into a partial sum and a remainder term.
\[
Z(z) = Z_n(z) + R_n(z)
\]
where
\[
Z_n(z) = \sum_{\nu = 0}^n \zeta(\nu + 2)(-z)^\nu, \quad R_n(z) = (-z)^{n + 1} \sum_{\nu = 0}^\infty \zeta(n + \nu + 3)(-z)^\nu
\]

2. **Construct Rational Approximation**: Use the known power series coefficients and truncation error information to construct an explicit rational approximation formula.
\[
T(k)_n(Z_n(z), \ldots, Z_{n + k}(z))=\frac{\sum_{\kappa = 0}^k e(k)_{k - \kappa}(q_1, \ldots, q_k)E^\kappa Z_n(z)}{\sum_{\kappa = 0}^k e(k)_\kappa(q_1, \ldots, q_k)}
\]
where \( E \) is the shift operator and \( e(k)_\kappa \) is the elementary symmetric polynomial.

3. **Error Analysis**: Give an explicit truncation error expression and verify the effectiveness of the new method through numerical experiments.

### Conclusion
The paper successfully constructs a new rational approximation method for calculating the digamma function. This method not only improves the computational efficiency but also provides an explicit truncation error expression, which is helpful for further optimization and accelerating convergence.