Numerical implementation (using mpmath) of the series explained in Some geometric series for Euler’s constant (2026-03-31).
This serves as a brief mathematical introduction to Kempner-Irwin harmonic series, and to the formulas by the author regarding them, which allow computational algorithms distinct from the Baillie and Schmelzer-Baillie ones (links to the latter two are given on the page).
It also serves as front page to the repository where one finds some SageMath scripts implementing these formulas.
And it provides links to the further papers by the author (limited to year 2024) on this and related topics.
Links directly (no front page) to the GitLab repository which provides Python code using mpmath (for arbitrary bases and sets of excluded digits) in addition to legacy Maple code (in its 2024 version with English comments) as used by the author some years ago.
The Maple code is superseded by the April 2026 Python+mpmath code, and was already so earlier by the May 2025 SageMath scripts at the Irwin repository (only for the case of a single excluded digit).