gens:=[
    Matrix(GF(3),4,4,[1, 0, 0, 0, 0, 1, 0, 0, 2, 2, 2, 0, 0, 0, 0, 2 ]),

    Matrix(GF(3),4,4,[1, 1, 1, 1, 2, 0, 1, 1, 2, 2, 2, 1, 2, 2, 1, 0 ]),

    Matrix(GF(3),4,4,[0, 0, 1, 1, 2, 1, 1, 1, 2, 2, 0, 1, 2, 0, 2, 2 ])
];



NH:=PermutationGroup<22|(1, 3)(4, 16)(6, 13)(7, 8)(9, 14)(12, 22)(15, 19)(18, 20),(2, 9, 12, 5)(3, 6, 17, 20)(4, 13, 18, 16)(7, 10)(8, 15, 21, 19)(14, 22), (2, 8, 12, 21)(3, 4, 17, 18)(5, 15, 9, 19)(6, 16, 20, 13)(7, 22)(10, 14)>;


V49d:=GModule(NH,gens);


print("Green correspondent V49d of the simple M22-module D49d=D49^* 
in NH (normalizer of vertex P in M22) over GF(3)");