The fglm command computes the reduced lexicographic Groebner basis from a grevlex basis A in the variables X. The system must have a finite number of solutions, i.e. its hilbert dimension should be zero. An optional argument P is the modulus, which should be a machine size prime. Progress information is available by setting info().
sys = [x^2+y, x*y-1]; G = groebnerbasis(sys,[x,y]); fglm(G,[x,y],[x,y]); [1+y^3, x+y^2]