Remainder Pade Approximants for the Exponential Function

Following our earlier research, we propose a new method for obtaining the complete Pade table of the exponential function. It is based on an explicit construction of certain Pade approximants, not for the usual power series for exp at 0 but for a formal power series related in a simple way to the remainder term of the power series for exp. This surprising and nontrivial coincidence is proved more generally for type II simultaneous Pade approximants for a family \((\exp(a_jz))_{j=1,\ldots, r}\) with distinct complex a's and we recover Hermite's classical formulas. The proof uses certain discrete multiple orthogonal polynomials recently introduced by Arvesu, Coussement, and van Assche, which generalize the classical Charlier orthogonal polynomials.