e=2*eye(2);
for i=0:15
 a=m4(i);
 display(i)
 for j=i:15
  b=m4(j);
  c=a*b'+b*a';
  f=a*a'+b*b'-2*e;
  if norm(c)==0 && norm(f)==0
   d=a*b'-b*a';
   g=a*a'-b*b';
   fprintf(' %2d, %2d:\n',i,j);
   display(a)
   display(b)
   display(d)
   display(g)
   fprintf('\n');
  end
 end
end
