  In the Lee-Pomeransky representation, Feynman integrals can be identified as
a subset of Euler-Mellin integrals, which are known to satisfy
Gel'fand-Kapranov-Zelevinsky (GKZ) system of partial differential equations.
Here we present an automated package to derive the associated GKZ system for a
given Feynman diagram and solve it in terms of hypergeometric functions using
two equivalent algorithms, namely the triangulation method and the Gr\"obner
deformation method. We present our code in the form of a Mathematica package
FeynGKZ.wl which requires the softwares polymake, Macaulay2 and TOPCOM, and the
packages AMBRE and Olsson.wl as dependencies. As applications of the package,
we find series solutions to the GKZ systems of several one-loop and two-loop
Feynman integrals. These are included in the file Examples.nb that can be
downloaded along with the package from https://github.com/anant-group/FeynGKZ.
