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

F:=GF(2);

x:=CambridgeMatrix(1,F,16,[
"0100000000000000",
"1000000000000000",
"0001000000000000",
"0010000000000000",
"1111100000000000",
"0000000100000000",
"0011011100000000",
"0000010000000000",
"0000000010000000",
"0101110100010000",
"1010100001110000",
"0101110101000000",
"0000000000001000",
"0000000000000100",
"0110110101010010",
"1001110101010001"]);

y:=CambridgeMatrix(1,F,16,[
"0010000000000000",
"0110000000000000",
"0000100000000000",
"0000010000000000",
"0000001000000000",
"0000000010000000",
"0000000001000000",
"0000000000100000",
"0000000000010000",
"0000000000001000",
"0000000000000100",
"0000000000000010",
"0000000000000001",
"1111011000100101",
"1111010100001000",
"0001100100101110"]);

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