f := randpoly([x,y,z],terms=30);
S := {seq(i=i*h, i=indets(f))};
g := subs(S, f);
g := coeff(g, h, 3);

[coeffs(g, indets(g))];