w2:=sqrt(2); w5:=sqrt(5); w10:=w2*w5; r1:=sqrt(w5+1); I:=sqrt(-1); phi:=[ 8, ((w10+w2-2*w5-2)*r1-4)*I+(w10+w2)*r1+4*w2-4, (8*w2-8)*I, ((w10+w2-2*w5-2)*r1+4)*I+(w10+w2)*r1-4*w2+4 ];