w2:=sqrt(2); r1:=sqrt(2*w2+1); I:=sqrt(-1); phi:=[ w2*r1*I+-w2-2, 2, 2, -r1*I+w2-1, 2, -r1*I+w2-1, -r1*I+w2-1 ];