/*
www-ATLAS of Group Representations.
2.O9(3) represented as 16 x 16 matrices over GF(3).
*/

F:=GF(3);

x:=CambridgeMatrix(1,F,16,[
"1000000000000000",
"0001000000000000",
"0000100000000000",
"0100000000000000",
"0010000000000000",
"0000000001000000",
"0000002000000000",
"0000000200000000",
"0000000000000100",
"0000010000000000",
"1000000000200000",
"1020120101210020",
"1201001020201120",
"0000000010000000",
"1000000000000020",
"0102002000100021"]);

y:=CambridgeMatrix(1,F,16,[
"0100000000000000",
"0010000000000000",
"1000000000000000",
"0000010000000000",
"0000000100000000",
"0000000000100000",
"0001000000000000",
"0000000000001000",
"0000100000000000",
"1220102002021100",
"0000000000000001",
"0000001000000000",
"0010110000020000",
"0000000002000000",
"0000000010000000",
"2201021120221022"]);

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