(setf *grobner-debug* t)

(setf vars '(e1 e2 e3 e4 h1 h2 h3 h4 x1 x2 x3 x4 y1 y2 y3 y4))

(setf ord #'lex>)

(string-ideal-polysaturation
 "[x1^2+e1+e2*x1+e3*y1+e4*x1*y1,
 y1^2+h1+h2*x1+h3*y1+h4*x1*y1,
 x2^2+e1+e2*x2+e3*y2+e4*x2*y2,
 y2^2+h1+h2*x2+h3*y2+h4*x2*y2,
 x3^2+e1+e2*x3+e3*y3+e4*x3*y3,
 y3^2+h1+h2*x3+h3*y3+h4*x3*y3,
 x4^2+e1+e2*x4+e3*y4+e4*x4*y4,
 y4^2+h1+h2*x4+h3*y4+h4*x4*y4]"
 (list "[x1-x2,y1-y2]"
       "[x1-x3,y1-y3]"
       "[x1-x4,y1-y4]"
       "[x2-x3,y2-y3]"
       "[x2-x4,y2-y4]"
       "[x3-x4,y3-y4]")
 vars :order ord)