Procedure from library
rinvar.lib (see rinvar_lib).
- HilbertSeries(I, w); ideal I, intvec wt
- compute the polynomial p of the Hilbert Series, represented by p/q, of
the ring K[t_1,...,t_m,y_1,...,y_r]/I1 where 'w' are the weights of
the variables, computed, e.g., by 'HilbertWeights', 'I1' is of the
form I - y_1,...,I[r] - y_r and is quasihomogenous w.r.t. 'w'
- the leading 0 of the result does not belong to p, but is needed in
the Hilbert driven 'std'.