OPERATOR R; MATRIX RR(3,3); FOR I:=1:3 DO FOR J:=1:3 DO RR(I,J):=R(I,J); V:=MAT((V1),(V2),(V3)); RV:=RR*V; Q:=MAT( (V3,V1+i*V2),(V1-i*V2,-V3)); % Now the unitary matrix. a and c are complex numbers, % ab and cb are their complex conjugates U:=MAT( (a,-cb),(c,ab)); % Ub is the Heritian conjugate of U Ub:= MAT( (ab, cb),(-c,a)); Qt:=U*Q*Ub; Vt3:=Qt(1,1); Vt1:=(Qt(1,2)+Qt(2,1))/2; Vt2:=(Qt(1,2)-Qt(2,1))/(2*i); COEFFN(RV(1,1)-Vt1,v1,1); COEFFN(RV(1,1)-Vt1,v2,1); COEFFN(RV(1,1)-Vt1,v3,1); COEFFN(RV(2,1)-Vt2,v1,1); COEFFN(RV(2,1)-Vt2,v2,1); COEFFN(RV(2,1)-Vt2,v3,1); COEFFN(RV(3,1)-Vt3,v1,1); COEFFN(RV(3,1)-Vt3,v2,1); COEFFN(RV(3,1)-Vt3,v3,1); END; In "H:/reduce myfiles/un1.red";