/*
www-ATLAS of Group Representations.
2.L2(31) represented as 16 x 16 matrices over GF(5).
*/

F:=GF(5);

x:=CambridgeMatrix(1,F,16,[
"0100000000000000",
"4000000000000000",
"0001000000000000",
"0040000000000000",
"0000001000000000",
"0000000100000000",
"0000400000000000",
"0000040000000000",
"0000000000100000",
"0000000000001000",
"0000000040000000",
"0000000000000001",
"0000000004000000",
"0024020430130201",
"3414422241133333",
"0000000000040000"]);

y:=CambridgeMatrix(1,F,16,[
"0010000000000000",
"4110000000000000",
"0000100000000000",
"0000010000000000",
"1000000000000000",
"0444440000000000",
"0000000010000000",
"0000000001000000",
"0000000000010000",
"0000000000000100",
"0000000000000010",
"0000001000000000",
"1402103344220021",
"0000000100000000",
"1024010100033342",
"4212344241014442"]);

G<x,y>:=MatrixGroup<16,F|x,y>;
print "Group G is 2.L2(31) < GL(16,GF(5))";