This is a Mathematica package to compute the Gross-Keating invariant and the Siegel series of a half-integral matrix over Zp.
Load the package into Mathematica as follows (assuming that computeGK.m is in an appropriate location)
<< computeGK.m;
Some examples can be found in the file computeGK_samples.nb.