up previous next
2.2.12 A Groebner Basis Example
A Groebner basis of an ideal I is calculated with the command GBasis(I), as illustrated in the following example.

Let r be a root of the equation x^7-x-1 over the rationals. The minimal polynomial of (4r-1)/r^3 can be found by computing the reduced Groebner basis of the ideal (x^7-x-1,x^3y-4x+1) with respect to the lexicographic term-ordering with x>y.

  Use R ::= QQ[x,y], Lex;
  Set Indentation;  -- to improve the appearance of the output
  G := GBasis(Ideal(x^7-x-1,x^3y-4x+1));

[ 1602818757152090759440/34524608236181199361x - 4457540/5875764481y^7
  - 47746460716124220/34524608236181199361y^6 +
  890175715271333840/34524608236181199361y^5 -
  3541992534667352220/34524608236181199361y^4 -
  55943894513139464160/34524608236181199361y^3 -
  56473654361333280980/34524608236181199361y^2 -
  27971979712025453040/34524608236181199361y -
  400704689288022689860/34524608236181199361, 1/16384y^7 - 5/16384y^6
  + 147/16384y^4 + 5/128y^3 - 31/16384y^2 + 17/128y - 20479/16384]
  F := 16384*G[2];  -- clear denominators
y^7 - 5y^6 + 147y^4 + 640y^3 - 31y^2 + 2176y - 20479
The Groebner basis is reported as a list with two elements. The second gives a univariate polynomial which is the minimal polynomial for r.

Note that the statement declaring the ring includes the modifier, Lex. Without this modifier, the default term-ordering, DegRevLex, is used. The command Set Indentation forces each polynomial of the Groebner basis to be printed on a new line.