fprintf('graph G {\n');
n=4;
N=2^(n*n);
e=n*eye(n);
for i=0:(N-1)
 a=m16(i);
 for j=(i+1):(N-1)
  b=m16(j);
  % f=a*a'+b*b'-2*e;
  if norm(a*b'+b*a')==0
   fprintf(' %2d -- %2d [color=red];\n',i,j);
  end
  if norm(a*b'-b*a')==0
   fprintf(' %2d -- %2d [color=blue,len=2];\n',i,j);
  end
 end
end
fprintf('}\n');

