A1=randperm(4);
B=eye(4);
C=[B(:,3:4),B(:,1:2)];
D=[B(:,2:4),B(:,1)];
A2=A1*C;
A3=A2*D;
A4=A3*C;
A=[A1;A2;A3;A4]

    Source: geocities.com/hk/tamwingyin0415

               ( geocities.com/hk)