/*
www-ATLAS of Group Representations.
L2(16):4 represented as 17 x 17 matrices over GF(17).
*/

F:=GF(17);

x:=CambridgeMatrix(3,F,17,\[
0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,
0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,
0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,
12,7,14,2,10,11,15,11,8,6,8,16,6,11,16,0,0,
0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,
16,1,1,16,1,0,16,0,0,0,0,0,0,1,0,0,0,
1,16,16,1,16,0,1,0,0,0,0,0,0,0,1,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0]);

y:=CambridgeMatrix(3,F,17,\[
1,16,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,
0,16,16,0,16,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,
0,10,7,0,2,0,8,0,7,0,11,3,0,0,6,0,0,
0,16,2,2,15,15,2,0,1,1,14,14,1,16,3,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,
0,1,15,0,1,0,16,0,15,0,3,3,0,0,15,0,0,
0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,
1,0,1,2,0,0,1,0,0,1,0,0,1,0,1,0,16]);

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