function findsperm(k)
pp=factorial(k);
nn=2^k;
e=eye(k);
for p1=1:pp
 for n1=1:nn
  A1=spermat(k,p1,n1);
  for p2=p1:pp
   for n2=1:nn
    A2=spermat(k,p2,n2);
    if norm(A1*A2'-A2*A1')==0
     fprintf('Amicable: (%d,%d), (%d,%d)\n',p1,n1,p2,n2);
    end
    if norm(A1*A2'+A2*A1')==0
     fprintf('Fraxious: (%d,%d), (%d,%d)\n',p1,n1,p2,n2);
    end
   end
  end
 end
end
  