/*
www-ATLAS of Group Representations.
U3(4) represented as 3 x 3 matrices over GF(16).
*/

F<w>:=GF(16);

x:=CambridgeMatrix(3,F,3,\[
1,0,0,
0,1,0,
11,13,1]);

y:=CambridgeMatrix(3,F,3,\[
0,1,0,
0,0,1,
1,0,0]);

G<x,y>:=MatrixGroup<3,F|x,y>;
print "Group G is U3(4) < GL(3,GF(16))";