pi=4*a(1)
w1=2*c(2*pi/7)
w2=2*c(4*pi/7)
w3=2*c(8*pi/7)

w1^3 + w1^2 - 2*w1 - 1
w2^3 + w2^2 - 2*w2 - 1 
w3^3 + w3^2 - 2*w3 - 1 

w2-(w1^2-2)
w3-(w2^2-2)
w1-(w3^2-2)

w3-(-(w1^2)-w1+1)
w1-(-(w2^2)-w2+1)
w2-(-(w3^2)-w3+1)

s1= sqrt(7)
s2= -sqrt(7)

limit=5;

for(e00=-limit;e00<=limit;e00++)
  for(e01=-limit;e01<=limit;e01++)
    for(e10=-limit;e10<=limit;e10++)
      for(e11=-limit;e11<=limit;e11++)
        for(e20=-limit;e20<=limit;e20++)
          for(e21=-limit;e21<=limit;e21++) {
            x11=(e00+e01*s1)+(e10+e11*s1)*w1+e20*w1^2+e21*(4+3*w1+w1^2)*s1/7
            x12=(e00+e01*s2)+(e10+e11*s2)*w1+e20*w1^2+e21*(4+3*w1+w1^2)*s2/7
            x21=(e00+e01*s1)+(e10+e11*s1)*w2+e20*w2^2+e21*(4+3*w2+w2^2)*s1/7
            x22=(e00+e01*s2)+(e10+e11*s2)*w2+e20*w2^2+e21*(4+3*w2+w2^2)*s2/7
            x31=(e00+e01*s1)+(e10+e11*s1)*w3+e20*w3^2+e21*(4+3*w3+w3^2)*s1/7
            x32=(e00+e01*s2)+(e10+e11*s2)*w3+e20*w3^2+e21*(4+3*w3+w3^2)*s2/7
	    if(x12<0 || x22<0 || x32<0)continue
	    if(x11>0 && x21>0 && x31>0)continue
	    if(x11<0 && x21<0 && x31<0)continue
            print "[", \
              e00,",",e01,",",e10,",",e11,",",e20,",",e21, \
              "],\n";
          }