w5:=sqrt(5); s1:=sqrt(2*w5+1); I:=sqrt(-1); phi:=[ 2*(w5+1), -s1*I+1, s1*I+1, s1*I+1, -s1*I+1, -s1*I+1, -s1*I+1, -s1*I+1, s1*I+1, -s1*I+1, s1*I+1, -s1*I+1, s1*I+1, s1*I+1, s1*I+1, s1*I+1, -s1*I+1, -s1*I+1, s1*I+1 ];