I want presentation of non-Abelian group $(C_3\times C_3)\rtimes C_3$ of order $27.$ Actually it is $$(C_3\times C_3)\rtimes C_3\cong \{\begin{pmatrix}1&a&b\\0&1&c\\0&0&1\end{pmatrix}\;|\; a,b,c\in \Bbb{Z}_3\}$$. But i want its presentation.
There is link Classify groups of order 27 which classifies all group of order $27.$ In this link presentation is given as $\langle a,b\mid a^3,b^3,[a,b]=[a,b]^a=[a,b]^b\rangle$ but i did't get this. Please help me so that i can write all $27$ elements of group. I also tried in GAP as below
gap> p:=3;
3
gap> N:=CyclicGroup(p);
gap> K:=DirectProduct(N,N);
gap> H:=CyclicGroup(p);
gap> hom:=AllHomomorphisms(H,AutomorphismGroup(K));
[ [ f1 ] -> [ IdentityMapping( ) ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1*f2^2, f2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1, f1^2*f2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1*f2, f2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1^2*f2, f1^2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1, f1*f2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f2^2, f1*f2^2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f2, f1^2*f2^2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1^2*f2^2, f1 ] ] ]
gap> hom:=AllHomomorphisms(H,AutomorphismGroup(K)); [ [ f1 ] -> [ IdentityMapping( ) ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1*f2^2, f2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1, f1^2*f2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1*f2, f2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1^2*f2, f1^2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1, f1*f2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f2^2, f1*f2^2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f2, f1^2*f2^2 ] ], [ f1 ] -> [ Pcgs([ f1, f2 ]) -> [ f1^2*f2^2, f1 ] ] ]
gap> g:=SemidirectProduct(H,hom[1],K);
gap> g:=SemidirectProduct(H,hom[1],K);
gap> StructureDescription(g);
"C3 x C3 x C3"
gap> g:=SemidirectProduct(H,hom[2],K);
gap> StructureDescription(g);
"(C3 x C3) : C3"
gap> IsAbelian(g);
false gap> Size(g);
27 gap> Elements(g);
[ of ..., f1, f2, f3, f1^2, f1*f2, f1*f3, f2^2, f2*f3, f3^2, f1^2*f2, f1^2*f3, f1*f2^2, f1*f2*f3, f1*f3^2, f2^2*f3, f2*f3^2, f1^2*f2^2, f1^2*f2*f3, f1^2*f3^2, f1*f2^2*f3, f1*f2*f3^2, f2^2*f3^2, f1^2*f2^2*f3, f1^2*f2*f3^2, f1*f2^2*f3^2, f1^2*f2^2*f3^2 ]
gap> t:=Elements(g); [ of ..., f1, f2, f3, f1^2, f1*f2, f1*f3, f2^2, f2*f3, f3^2, f1^2*f2, f1^2*f3, f1*f2^2, f1*f2*f3, f1*f3^2, f2^2*f3, f2*f3^2, f1^2*f2^2, f1^2*f2*f3, f1^2*f3^2, f1*f2^2*f3, f1*f2*f3^2, f2^2*f3^2, f1^2*f2^2*f3, f1^2*f2*f3^2, f1*f2^2*f3^2, f1^2*f2^2*f3^2 ]
So please help me to find a representation so that i can find all elements easily thanks.