restart;
loadPackage ("FourTiTwo");

-- load matrix, set up the ring
M = getMatrix("h0.zinhom");  -- h1 or h2
nvar = numColumns(M)//2;
R = QQ[p_0 .. p_(nvar-1)];

-- turn a vector into a monomial
toMonomial = V -> product for i from 0 to (nvar-1) list (R_i)^(V#i);
-- turn a matrix into a monomial ideal (each row is a monomial)
toMonomialIdeal = M -> monomialIdeal for i from 0 to numRows(M)-1 list toMonomial flatten entries (M^{i});

I = toMonomialIdeal M;
SPI = standardPairs I;

-- output for h1: {{1, {p_11, p_8, p_5, p_13, p_14, p_3, p_0, p_6}}}
-- output for h2: {{1, {p_10, p_9, p_4, p_12, p_15, p_1, p_2, p_7}}}