#We find all primitive subgroups of $GL(9,7)$ with parameters $e=9$, $a=1$
#In the notations of Theorem 5.1, if $H$ is a primitive solvable subgroup, then it possesses a series

#$$1<U<F<A<H,$$

#where $U$ is cyclic of order $6$, $F=U\circ E$ and $E$ is extraspecial of order $3^{1+4}$, in this case $H=A$. 

#First we make general variables:


gap> z:=Z(7);;
gap> v:=[z,z^2,z^3,0*z,z^5,z,z^4,z^3,z^5];;
gap> result:=[];;


#Now we run the following loop, where gens is the list of generators from generators-GL(9,7).txt. Each entry of the list result has three items: the first is the
#number of the group in the list gens, the second is the order of the group, and the third is the order of the stabilizer of v in H. 
#If the order of the stabilizer is 1, then |v^H=|H| and Corollary 2.5 can be applied.

gap> for j in [1..Size(gens)] do
> H:=Subgroup(GL(9,7),gens[j]);;
> n:=Size(H);;
> m:=Size(v^H);
> result[j]:=[j,n,n/m];
> Display(result[j]);
> od;
[ 1, 2430, 1 ]
[ 2, 4860, 1 ]
[ 3, 3888, 1 ]
[ 4, 3888, 1 ]
[ 5, 9720, 1 ]
[ 6, 7776, 1 ]
[ 7, 7776, 1 ]
[ 8, 7776, 1 ]
[ 9, 7776, 1 ]
[ 10, 7776, 1 ]
[ 11, 7776, 1 ]
[ 12, 7776, 1 ]
[ 13, 11664, 1 ]
[ 14, 19440, 1 ]
[ 15, 15552, 1 ]
[ 16, 15552, 1 ]
[ 17, 15552, 1 ]
[ 18, 15552, 1 ]
[ 19, 15552, 1 ]
[ 20, 15552, 1 ]
[ 21, 23328, 1 ]
[ 22, 23328, 1 ]
[ 23, 23328, 1 ]
[ 24, 23328, 1 ]
[ 25, 23328, 1 ]
[ 26, 31104, 1 ]
[ 27, 31104, 1 ]
[ 28, 31104, 1 ]
[ 29, 46656, 1 ]
[ 30, 46656, 1 ]
[ 31, 46656, 1 ]
[ 32, 46656, 1 ]
[ 33, 46656, 1 ]
[ 34, 77760, 1 ]
[ 35, 62208, 1 ]
[ 36, 93312, 1 ]
[ 37, 93312, 1 ]
[ 38, 93312, 1 ]
[ 39, 93312, 1 ]
[ 40, 155520, 1 ]
[ 41, 186624, 1 ]
[ 42, 186624, 1 ]
[ 43, 279936, 1 ]
[ 44, 559872, 1 ]