Read("types.gap");
Read("perms2matrices.gap");

# n = number of points on which gg acts
# types = collection of conjugacy classes  

ntypes:=Size(types);;
gg:=Representative(types[ntypes]);;


Print("\nPsi ",psi,"\n");

for symnumber in [1..ntypes] do
  Print("\nSymmetry type S",symnumber,"\n");
  S:=types[symnumber];;              
  for s in S do
    Print(s,"\n");
    int:=IdentityMat(n);
    for g in s do
      mg:=perm2matrix(g,n);
      nul:=NullspaceMat(psi*mg-psi);
      int:=SumIntersectionMat(int,nul)[2];
    od;
    Print(" Basis ",int,"\n");
  od;
od;