PrecisionLauricella: Package for numerical computation of Lauricella functions depending on a parameter

We introduce the PrecisionLauricella package, a computational tool developed in Wolfram Mathematica for high-precision numerical evaluations of the Laurent expansion coefficients of Lauricella functions whose parameters depend linearly on a small regulator, ε. In practical multi-loop calculations, Lauricella functions are required only as series around ε = 0, and PrecisionLauricella is designed specifically to deliver such coefficients with arbitrary precision. The package leverages a method based on analytic continuation via Frobenius generalized power series, providing an efficient and accurate alternative to conventional approaches relying on multi-dimensional series expansions or Mellin–Barnes representations. This one-dimensional approach is particularly advantageous for high-precision calculations and facilitates further optimization through ε-dependent reconstruction from evaluations at specific numerical values, enabling efficient parallelization. The underlying mathematical framework for this method has been detailed in our previous work, while the current paper focuses on the design, implementation, and practical applications of the PrecisionLauricella package.