  We present the Mathematica package $\texttt{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 arXiv:2208.01000v2. 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 https://github.com/souvik5151/MultiHypExp.
