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