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

F:=GF(2);

x:=CambridgeMatrix(1,F,16,[
"1001011000000000",
"0101011000000000",
"0011011000000000",
"0001000000000000",
"0001111000000000",
"1010010000000000",
"1010001000000000",
"1011011100000000",
"0000000010010001",
"0000000010010100",
"0000000001101101",
"0000000001110010",
"0000000000011001",
"0000000010100011",
"0000000010011011",
"0000000001100011"]);

y:=CambridgeMatrix(1,F,16,[
"0000000010000011",
"0000000011011001",
"0000000010001101",
"0000000011110100",
"0000000001011011",
"0000000000011101",
"0000000000111010",
"0000000001100111",
"1001101000000000",
"0000011100000000",
"1000000100000000",
"1001011000000000",
"1010001100000000",
"1101100100000000",
"1111101000000000",
"1100000000000000"]);

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