sage: R.<a,b,c,d,e,A,B,C,D,E,F> = PolynomialRing(QQ, order='lex')
sage: I = ideal(-A+a^2+b^2+c^2+d^2+e^2,-B+a^2*b^2+a^2*c^2+a^2*d^2+a^2*e^2+b^2*c^2+b^2*d^2+b^2*e^2+c^2*d^
....: 2+c^2*e^2+d^2*e^2,-C+a^2*b^2*c^2+a^2*b^2*d^2+a^2*b^2*e^2+a^2*c^2*d^2+a^2*c^2*e^2+a^2*d^2*e^2+b^2*c
....: ^2*d^2+b^2*c^2*e^2+b^2*d^2*e^2+c^2*d^2*e^2,-D+a^2*b^2*c^2*d^2+a^2*b^2*c^2*e^2+a^2*c^2*d^2*e^2+a^2*
....: b^2*d^2*e^2+b^2*c^2*d^2*e^2,-E+a^2*b^2*c^2*d^2*e^2,-F+a+b+c+d+e)
....:
sage: G = I.groebner_basis()