In case the input is an h-vector the script calls BEdegrees to compute the degree pattern of the (unique) minimal size skew-symmetric matrix whose pfaffians define a Gorenstein ideal with that h-vector. We then choose a random skew-symmetric matrix of forms of the given BEdegrees.
i1 : S=ZZ/10009[x_0..x_2]
o1 = S
o1 : PolynomialRing
|
i2 : hV1=hVector 9; hV=hV1|reverse hV1
o3 = {1, 3, 5, 5, 3, 1}
o3 : List
|
i4 : (J,M)=pickGorensteinIdeal(S,hV);
|
i5 : betti res J
0 1 2 3
o5 = total: 1 3 3 1
0: 1 . . .
1: . 1 . .
2: . 2 . .
3: . . 2 .
4: . . 1 .
5: . . . 1
o5 : BettiTally
|
i6 : hVector J
o6 = {1, 3, 5, 5, 3, 1}
o6 : List
|