fprintf('graph G {\n');
n=3;
N=2^(n*n);
e=n*eye(n);
for i=0:(N-1)
 fprintf(' %2d [fixedsize=true,height=0.1,width=0.1,label=""]\n',i);
 a=mat(n,i);
 for j=(i+1):(N-1)
  b=mat(n,j);
  if norm(a*a'+b*b'-2*e)==0
   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
end
fprintf('}\n');

