# Computing Stieltjes constants using complex integration

1 LFANT - Lithe and fast algorithmic number theory
IMB - Institut de Mathématiques de Bordeaux, Inria Bordeaux - Sud-Ouest
Abstract : The Stieltjes constants $\gamma_n$ are the coefficients appearing in the Laurent series of the Riemann zeta function at $s = 1$. We give a simple and efficient method to compute a $p$-bit approximation of $\gamma_n$ with rigorous error bounds. Starting from an integral representation due to Blagouchine, we shift the contour to eliminate cancellation. The integral is then evaluated numerically in ball arithmetic using the Petras algorithm, with the use of a Taylor expansion for bounds near the saddle point. This appears to be the first algorithm for Stieltjes constants with uniformly low complexity with respect to both $n$ and $p$. An implementation is provided in the Arb library. We can, for example, compute $\gamma_n$ to 1000 digits in a minute for any $n \le 10^{100}$.
Document type :
Preprints, Working Papers, ...
Domain :

Cited literature [19 references]

https://hal.inria.fr/hal-01758620
Contributor : Fredrik Johansson <>
Submitted on : Wednesday, April 4, 2018 - 4:06:46 PM
Last modification on : Tuesday, June 5, 2018 - 1:20:12 AM

### Files

stieltjes.pdf
Files produced by the author(s)

### Identifiers

• HAL Id : hal-01758620, version 1
• ARXIV : 1804.01679

### Citation

Fredrik Johansson. Computing Stieltjes constants using complex integration. 2018. ⟨hal-01758620v1⟩

Record views