| i1 : R = QQ[x,y,z,w]
 o1 = R
 
 o1 : PolynomialRing
 | 
  
    | i2 : I = ideal(x*y-z,y^2-w-1,w^4-3)
 2           4
 o2 = ideal (x*y - z, y  - w - 1, w  - 3)
 
 o2 : Ideal of R
 | 
  
    | i3 : gb(I, PairLimit => 1)
 o3 = {0} | y2-w-1 |
 
 o3 : GroebnerBasis
 | 
  
    | i4 : gb(I, PairLimit => 2)
 o4 = {0} | y2-w-1 xy-z |
 
 o4 : GroebnerBasis
 | 
  
    | i5 : gb(I, PairLimit => 3)
 o5 = {0} | yz-xw-x y2-w-1 xy-z |
 
 o5 : GroebnerBasis
 |