10.07.2015 Views

University of Paderborn Department of Mathematics Diploma Thesis ...

University of Paderborn Department of Mathematics Diploma Thesis ...

University of Paderborn Department of Mathematics Diploma Thesis ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

4.1. COMPUTING HILBERT SERIES OF STABLE IDEALS 77compute_Hilbert_series:= proc(M, n)local di, li, monom, H, H_red, i, j;beginH:= poly(1,[t]);for i from 1 to nops(M) domonom:= M[i];for j from n+1 downto 1 do/* Determine index li <strong>of</strong> last variable in monom */if monom[j] 0 thenli:= j-1;break;end_if;end_for;/* Determine the degree <strong>of</strong> monom */di:= _plus(op(monom));H:= H - poly((1-t)^li * t^di,[t]);end_for;/* At this point H is the numerator <strong>of</strong> the non-reducedHilbert series. Now compute the numerator H_red <strong>of</strong>the reduced Hilbert series via gcd. */H_red:= H/gcd(H, poly((1-t)^(n+1), [t]));if H_red(0) = -1 thenH_red:= H_red * poly(-1, [t]);end_if;return(H, H_red);end_proc:To illustrate how to use the procedure, we consider a few examples where we compute theHilbert series <strong>of</strong> some stable ideals, which also occurred in some <strong>of</strong> the examples <strong>of</strong> thepreceding two chapters.Example 4.2. First, let R := K[x 0 , x 1 , x 2 ], i.e. n = 2. We will consider some “larger”examples below in (iii) and (iv).(i) Consider the stable ideal I 1 := (x 2 0, x 0 x 1 , x 2 1) ⊂ R. The set <strong>of</strong> generators is encoded bythe list <strong>of</strong> lists[[2, 0, 0], [1, 1, 0], [0, 2, 0]].The value <strong>of</strong> n is 2. Hence, we compute the Hilbert series <strong>of</strong> R/I 1 by the following simplecommand:

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!