A new approach to the epsilon expansion of generalized hypergeometric functions

Assuming that the parameters of a generalized hypergeometric function depend linearly on a small variable ε$\epsilon$, the successive derivatives of the function with respect to that small variable are evaluated at ε=0$\epsilon =0$ to obtain the coefficients of the ε$\epsilon$-expansion of the function. The procedure, which is quite naive, benefits from simple explicit expressions of the derivatives, to any order, of the Pochhammer and reciprocal Pochhammer symbols with respect to their argument. The algorithm may be used algebraically, irrespective of the values of the parameters. It reproduces the exact results obtained by other authors in cases of especially simple parameters. Implemented numerically, the procedure improves considerably, for higher orders in ε$\epsilon$, the numerical expansions given by other methods.