w2:=sqrt(2);
w5:=sqrt(5);
I:=sqrt(-1);

s1:=sqrt(1+w5);

w10:=w2*w5;
zeta_8:=(1+I)/w2; 

fid4a:=[
8*w2-8,
((w10+w2)*s1+4*w2-4)-((-w10-w2+2*w5+2)*s1+4)*I,
8*I,
(-(w10+w2)*s1+4*w2-4)-((w10+w2-2*w5-2)*s1+4)*I
];