The groebnerbasis command computes the reduced Groebner basis for a set of polynomials A in the variables X. The command uses the F4 algorithm in graded reverse lexicographical order. An optional argument P is the modulus, which should be a machine size prime. Progress information is available by setting info().
A = [x^2+y, x*y-1]; groebnerbasis(A,[x,y]); [x+y^2, -1+x*y, y+x^2]