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

Abstract

We present a new methodology, suitable for implementation on computer, to perform the $\epsilon$-expansion of hypergeometric functions with linear $\epsilon$ dependent Pochhammer parameters in any number of variables. Our approach allows one to perform Taylor as well as Laurent series expansion of multivariable hypergeometric functions. Each of the coefficients of $\epsilon$ 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 hypergeometric function. We present illustrative examples of hypergeometric functions in one, two and three variables which are typical of Feynman integral calculus.