t=2^64-1;
for i=0:t
 a=m64(i);
 display(i)
 for j=i:t
  if mod(j,16384) == 0
   display(j)
  end 
  b=m64(j);
  c=a'*b+b'*a;
  d=a'*b-b'*a;
  if norm(c)==0 
   fprintf(' %20d, %20d:\n',i,j);
   display(a)
   display(b)
   display(d)
   fprintf('\n');
  end
 end
end
