MultiHypExp: A Mathematica package for expanding multivariate hypergeometric functions in terms of multiple polylogarithms

We present the Mathematica package MultiHypExp that allows for the expansion of multivariate hypergeometric functions (MHFs), especially those likely to appear as solutions of multi-loop, multi-scale Feynman integrals, in the dimensional regularization parameter. The series expansion of MHFs can be carried out around integer values of parameters to express the series coefficients in terms of multiple polylogarithms. The package uses a modified version of the algorithm prescribed in [1]. In the present work, we relate a given MHF to a Taylor series expandable MHF by a differential operator. The Taylor expansion of the latter MHF is found by first finding the associated partial differential equations (PDEs) from its series representation. We then bring the PDEs to the Pfaffian system and further to the canonical form, and solve them order by order in the expansion parameter using appropriate boundary conditions. The Taylor expansion so obtained and the differential operators are used to find the series expansion of the given MHF. We provide examples to demonstrate the algorithm and to describe the usage of the package, which can be found in this url https://github.com/souvik5151/MultiHypExp.