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

F<w>:=GF(16);

x:=CambridgeMatrix(3,F,9,\[
0,1,0,0,0,0,0,0,0,
1,0,0,0,0,0,0,0,0,
0,0,0,1,0,0,0,0,0,
0,0,1,0,0,0,0,0,0,
0,0,0,0,0,1,0,0,0,
0,0,0,0,1,0,0,0,0,
8,8,0,0,10,10,1,0,0,
0,0,0,0,0,0,0,0,1,
0,0,0,0,0,0,0,1,0]);

y:=CambridgeMatrix(3,F,9,\[
0,0,1,0,0,0,0,0,0,
0,1,1,0,0,0,0,0,0,
1,0,1,0,0,0,0,0,0,
0,0,0,0,1,0,0,0,0,
0,0,0,0,0,0,1,0,0,
0,0,0,0,0,0,0,1,0,
0,0,0,1,0,0,0,0,0,
11,8,1,1,13,6,12,1,1,
10,0,3,13,12,7,1,7,1]);

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