3

I tried to use Buchberger's Algorithm (see this or this) to compute a Groebner basis for an ideal in rational field. Here is my GAP script:

F := Rationals;
R := PolynomialRing( F, [ "x", "y", "z" ]);
x := IndeterminatesOfPolynomialRing(R)[1];
y := IndeterminatesOfPolynomialRing(R)[2];
z := IndeterminatesOfPolynomialRing(R)[3];
I := Ideal (R, [x^2+2*x*y^2, x*y + 2*y^3 - 1]);
ord := MonomialLexOrdering(x,y,z);

GroebnerBasis( I, ord );

but the result is always this:

[ 2*x*y^2+x^2, 2*y^3+x*y-1, -x, -4*y^4+2*y, 2*y^3-1 ]

Obviously, the fourth can be completely divided by the last basis, the first and second can be completely divided by the third basis. The expected result should be like this:

 [ -x, 2*y^3-1 ]

So my question is how to get the simplified Groebner basis in GAP?

Mike Pierce
  • 1,390
  • 1
  • 12
  • 35
jscoot
  • 2,019
  • 3
  • 24
  • 27

1 Answers1

4

Try the ReducedGroebnerBasis command:

gap> ReducedGroebnerBasis(I, ord);
[ y^3-1/2, x ]
Mike Pierce
  • 1,390
  • 1
  • 12
  • 35
rcs
  • 67,191
  • 22
  • 172
  • 153