  The Chisholm rational approximant is a natural generalization to two
variables of the well-known single variable Pad\'e approximant, and has the
advantage of reducing to the latter when one of the variables is set equals to
0. We present, to our knowledge, the first automated Mathematica package to
evaluate diagonal Chisholm approximants of two variable series. For the moment,
the package can only be used to evaluate diagonal approximants i.e. the maximum
powers of both the variables, in both the numerator and the denominator, is
equal to some integer $M$. We further modify the original method so as to allow
us to evaluate the approximants around some general point $(x,y)$ not
necessarily $(0,0)$. Using the approximants around general point $(x,y)$,
allows us to get a better estimate of the result when the point of evaluation
is far from $(0,0)$. Several examples of the elementary functions have been
studied which shows that the approximants can be useful for analytic
continuation and convergence acceleration purposes. We continue our study using
various examples of two variable hypergeometric series,
$\mathrm{Li}_{2,2}(x,y)$ etc that arise in particle physics and in the study of
critical phenomena in condensed matter physics. The demonstration of the
package is discussed in detail and the Mathematica package is provided as an
ancillary file.
