function find2(p)
fprintf('graph G {\n');
n=2;
EP=2^(p*p);
I=eye(p);
for i=0:(EP-1)
 B1=mat(p,i);
 printf('%5d\n',i);
 fflush(1);
 for j=i:(EP-1)
  B2=mat(p,j);
  M=B1+B2;
  B1B2T=B1*B2';
  if norm(M*M'-n*p*I)==0
   if norm(B1B2T+B1B2T')==0
    fprintf(' %2d -- %2d [color=red,len=1];\n',i,j);
   end
  end
  N=B1-B2;
  if norm(M*N'+2*B2*B2'-n*p*I)==0
   if norm(B1B2T-B1B2T')==0
    fprintf(' %2d -- %2d [color=blue,len=1];\n',i,j);
   end
  end
 end
end
fprintf('}\n');

