option(redSB);
ring r=2,(x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16),lp;
poly f1=x2*(1+x16)-x1; 
poly f2=(x5+1)*(x3+1)*(1+x16)-x2;
poly f3= (x2*x3+x2+x3)*(1+x16)-x3;
poly f4= x15*(1+x16)-x4 ;
poly f5= x4*(1+x16)-x5;
poly f6= (x7+1)*(1+x16)-x6;
poly f7= x15*(1+x16)-x7;
poly f8= (x15+1)*x6*(1+x16)-x8;
poly f9= (x6*x8*x11+x6*x8+x6*x11+x6+x8)*(1+x16)-x9;
poly f10= (x9*x12*x13+x9*x12+x12*x13+x9+x12)*(1+x16)-x10;
poly f11= (x9+1)*(1+x16)-x11 ;
poly f12= (x14+1)*(1+x16)-x12;
poly f13= (x12+1)*(1+x16)-x13 ;
poly f14= (x9+1)*(1+x16)-x14;
poly f15= (x8+1)*(1+x16)-x15;
poly f16= (x10*x16+x10+x16)-x16;
poly p1=x1^2-x1;
poly p2=x2^2-x2;
poly p3=x3^2-x3;
poly p4=x4^2-x4;
poly p5=x5^2-x5;
poly p6=x6^2-x6;
poly p7=x7^2-x7;
poly p8=x8^2-x8;
poly p9=x9^2-x9;
poly p10=x10^2-x10;
poly p11=x11^2-x11;
poly p12=x12^2-x12;
poly p13=x13^2-x13;
poly p14=x14^2-x14;
poly p15=x15^2-x15;
poly p16=x16^2-x16;
ideal i=f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,f13,f14,f15,f16,p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14,p15,p16;
ideal j=groebner(i);
j;