On the Jacobi formula for Bivariate Pade Approximants of Rectangular Type

In this paper a recursive algorithm is presented for evaluating multivariate Padé approximants (of the rectangular type described in the work of Lutterodt) which is analogous to the Jacobi formula for univariate Padé approximants. This algorithm is then applied to a (singular) Riccati differential equation to generate fast and accurate approximate solutions.