syms a b g
Rz1 = [cos(a) -sin(a) 0; sin(a) cos(a) 0; 0 0 1];
Ry1 = [cos(b) 0 sin(b); 0 1 0; -sin(b) 0 cos(b)];
Rz2 = [cos(g) -sin(g) 0; sin(g) cos(g) 0; 0 0 1];
Rzyz1 = Rz1*Ry1*Rz2

%expand(Rzyz1(1,1)^2)
%expand(Rzyz1(1,2)^2)
%expand(Rzyz1(2,1)^2)
%expand(Rzyz1(2,2)^2)
simplify(Rzyz1(1,1)^2+Rzyz1(1,2)^2)

for i = 3:1
    y=0
    for j = 3:1
        x=0
        for n = 3:1
            y=y+1
            for m = 3:1
                A = simplify(Rzyz1(i,j)^2+Rzyz1(n,m)^2);
                x=x+1
            end
        end
    end
end
A