FeynGKZ: A Mathematica package for solving Feynman integrals using GKZ hypergeometric systems

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öbner 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 GitHub.