OPERATOR R; MATRIX RR(3,3); FOR I:=1:3 DO FOR J:=1:3 DO RR(I,J):=R(I,J); R(1,1) := (a^2 + ab^2 - c^2 - cb^2)/2$ R(1,2) := i*(a^2 - ab^2 - c^2 + cb^2)/2$ R(1,3) := a*cb + ab*c$ R(2,1) := i*(ab^2 - a^2 + cb^2 - c^2)/2$ R(2,2) := (a^2 + ab^2 + c^2 + cb^2)/2$ R(2,3) := i*(ab*c - a*cb)$ R(3,1) := -a*c - ab*cb$ R(3,2) := i*(ab*cb - a*c)$ R(3,3) := a*ab - c*cb$ Let c*cb=1-a*ab; RR*TP(RR); Det(RR); END; In "H:/reduce myfiles/un2.red";