/*
www-ATLAS of Group Representations.
4a.L3(4):2a represented as 16 x 16 matrices over GF(7).
*/

F:=GF(7);

x:=CambridgeMatrix(1,F,16,[
"0100000000000000",
"1000000000000000",
"0001000000000000",
"0010000000000000",
"0000001000000000",
"0000000100000000",
"0000100000000000",
"0000010000000000",
"0000000000010000",
"0000000000000100",
"0101166642560230",
"0000000010000000",
"6154121615651566",
"0000000001000000",
"2323126233650320",
"5305004065220516"]);

y:=CambridgeMatrix(1,F,16,[
"0010000000000000",
"5660000000000000",
"0000100000000000",
"0000010000000000",
"1060100000000000",
"0000000010000000",
"0000000001000000",
"0000000000100000",
"0000000000001000",
"0000000000000010",
"6002200466606140",
"0000000000000001",
"0001000000000000",
"2011034020544216",
"4044143125632142",
"4042063032263152"]);

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