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

F<w>:=GF(9);

x:=CambridgeMatrix(1,F,16,[
"0100000000000000",
"1000000000000000",
"0001000000000000",
"0010000000000000",
"0000001000000000",
"0000000100000000",
"0000100000000000",
"0000010000000000",
"0000000000010000",
"0000000000000100",
"0000000000000001",
"0000000010000000",
"6357834613021600",
"0000000001000000",
"7142203627811428",
"0000000000100000"]);

y:=CambridgeMatrix(1,F,16,[
"1200000000000000",
"0010000000000000",
"0000100000000000",
"0000010000000000",
"0220200000000000",
"0000000010000000",
"0000000001000000",
"0000000000100000",
"0000000000001000",
"0000000000000010",
"8121336277064716",
"4438381227772276",
"0001000000000000",
"1873148470374525",
"0001000000002001",
"0001021000000000"]);

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