e=4*eye(4);
t=2^16-1;
for i=0:t
 a=m16(i);
 %display(i)
 if norm(a*a'-e)==0
  for j=i:t
   b=m16(j);
   c=a*b'+b*a';
   %f=a*a'+b*b'-*e;
   f=b*b'-e;
   if norm(c)==0 && norm(f)==0
    d=a*b'-b*a';
    fprintf(' %6d, %6d:\n',i,j);
    display(a)
    display(b)
    display(d)
    fprintf('\n');
   end
  end
 end
end
