e-expansion of multivariable hypergeometric functions appearing in Feynman integral calculus

We present a new methodology, suitable for implementation on computer, to perform the e-expansion of hypergeometric functions with linear e dependent Pochhammer parameters in any number of variables. Our approach allows one to perform Taylor as well as Laurent series expansion of multivariable hypergeo-metric functions. Each of the coefficients of e in the series expansion is expressed as a linear combination of multivariable hypergeometric functions with the same domain of convergence as that of the original hy-pergeometric function. We present illustrative examples of hypergeometric functions in one, two and three variables which are typical of Feynman integral calculus.(c) 2023 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons .org /licenses /by /4 .0/). Funded by SCOAP3.