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]); [y^2+x, x*y-1, x^2+y]