<<<>>>
p1:= poly(x^2 - x + 2*y^2, [x,y]):
p2:= poly(x + 2*y - 1, [x,y]):
groebner::gbasis([p1,p2], LexOrder);
<<<>>>
export(groebner):
gbasis([p1,p2], LexOrder, Reorder);
<<<>>>
export(groebner, normalf):
p:=  poly(x^2*y - 2*x*y + 1, [x,y]):
p1:= poly(x^2 - x + 2*y^2, [x,y]):
p2:= poly(x + 2*y - 1, [x,y]):
normalf(p, [p1,p2], LexOrder);
<<<>>>
export(groebner, spoly):
p1:= poly(x^2 - x + 2*y^2, [x,y]):
p2:= poly(x + 2*y - 1, [x,y]):
spoly(p1, p2, LexOrder);
<<<>>>
