/*
3"J3 presented in terms of its standard generators.
*/

G<x,y>:=Group<x,y|
x^2,
y^3,
(x,y)^9,
((x*y)^6*(x*y^-1)^5)^2,
((x*y*x*y*x*y^-1)^2*x*y*x*y^-1*x*y^-1*x*y*x*y^-1)^2,
x*y*x*y*(x*y*x*y^-1)^3*x*y*x*y*(x*y*x*y^-1)^4*x*y^-1*(x*y*x*y^-1)^3,
(x*y)^4*(x*y^-1)^2*x*y*(x*y*x*y*x*y^-1)^2*(x*y*x*y^-1*x*y)^2*x*y*(x*y^-1)^4*(x*y)^4*(x*y^-1)^3
// ((x*y)^5*x*y^-1*(x*y)^2*(x*y^-1)^5*x*y*(x*y^-1)^2)^2, // Redundant.
// ((x*y)^5*(x*y^-1*x*y*x*y^-1)^2)^3 // Redundant.
>;

A:=x;B:=y;


K1:=sub<G|x,x^y>;
// D18, of order 18 and index 8372160.
K2:=sub<G|x,y^(x*y)>;
// L2(16), of order 4080 and index 36936.


// `Maximal' subgroups [L2(19)s in same order as those that are produced from
// the word programs J3G1-max2/3W1].
M1:=sub<G|x,y*(x*y*x*y^-1)^2*(x*y^-1)^2*x*y*x*y>;
// L2(16):2, of order 8160 and index 18468.
M2:=sub<G|x,y*x*y*x*y*(x*y*x*y^-1)^2*(x,y*x*y^-1*x*y^-1)>;
// L2(19), of order 3420 and index 44064.
M3:=sub<G|x,(y^-1*x*y*x*y,x)*(y,x)^2*y*x*y*x*y>;
// L2(19), of order 3420 and index 44064.
M5:=sub<G|x,(x*y)^4*(x*y^-1)^3*(x*y)^4>;
// L2(17), of order 2448 and index 61560.
M6:=sub<G|x,(y*x)^4*(y^-1*x*y*x)^2*y*(x*y*x*y^-1)^2*x*y^-1>;
// (3 x A6):2b, of order 2160 and index 69768.
M8:=sub<G|x,(y*x*y^-1*x)^2*(y^-1*x*y^-1*x*y*x*y*x)^2*x^-1*y^-1>;
// 2^{1+4}:A5, of order 1920 and index 78489.


M4:=sub<G|x,y^-1*(x*y)^3*(x*y*x*y^-1)^2*x*y^-1*(x*y)^3*x*y^-1>;
// 2^4:(3 x A5), of order 2880 and index 52326.
M9:=sub<G|x,(x*y)^3*(x*y^-1)^2*x*y*(x*y^-1)^3*(x*y)^3>;
// 2^{2+4}:(3 x S3), of order 1152 and index 130815.


/*
To establish the order of G, it is probably necessary to do an enumeration
as hard as finding the 8372160 cosets of K1 = <x,x^y> = D18.
The following ought to work [WITHOUT the redundant relations].
It is a very hard enumeration though.
18*Index(G,K1:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=10*10^6);


Other enumerations should work as follows [without the redundant relations]:
Index(G,K2:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=100000);
Index(G,M1:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=10^6);
Index(G,M2:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=2*10^6);
Index(G,M3:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=2*10^6);
Index(G,M5:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=10^6);
Index(G,M6:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=2*10^6);
Index(G,M8:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=10^6);
Index(G,M4:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=2*10^6);
Index(G,M9:Print:=2,Hard:=true,Grain:=10^5,CosetLimit:=10^6);
*/