The eliminationbasis command computes the reduced Groebner basis for a set of polynomials A. The command uses the F4 algorithm in an elimination ordering with X > Y. 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]; eliminationbasis(A,[x],[y]); [1+y^3, x+y^2]