The univariatepoly command computes the univariate polynomial in Y 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]); univariatepoly(y,G,[x,y]); 1+y^3