local z, r, result;

result := rec();

result.comment := "L2(67) as 68 x 68 monomial matrices over Z(z33)\n";

# Change the value of r to any number between 1 and 16
# to get the complete set of inequivalent faithful irreducible 68-dimensional
# representations of L2(67)
r := 1;
z := E(33)^r;

result.symmetricforms := [ ]; 
result.antisymmetricforms := [ ]; 
result.hermitianforms := [ IdentityMat(68) ];
result.centralizeralgebra := [ IdentityMat(68) ];

result.generators := [
DiagonalMat([z^17,z^28,z^10,z^24,z^21,z^30,z^26,z^32,z^26,z^6,z^14,
  z^29,z^6,z,z^21,z^8,1,z^30,z^13,z^18,z^30,z^17,z^11,z^27,z^29,z^3,
  z^15,z^27,z^31,z^20,z^4,z^28,z^14,z^19,z^12,z^15,z^5,z^7,z^25,z^9,
  z^20,z^13,z^3,z^10,z^31,z^32,z^9,z^22,z^23,z^23,z^8,z^4,z^12,z^16,
  z^11,z^2,z^5,z^19,z^18,z^16,z^2,z^24,z^22,1,z^25,z,z^7,z^3]) *
PermutationMat(
( 1,54)( 2,37)( 3,49)( 4,40)( 5,35)( 6,26)( 7,38)( 8,66)( 9,67)(10,24)(11,58)
(12,52)(13,28)(14,46)(15,53)(16,39)(17,64)(18,43)(19,30)(20,27)(21,68)(22,60)
(23,63)(25,31)(29,61)(32,57)(33,34)(36,59)(41,42)(44,50)(45,56)(47,62)(48,55)
(51,65), 68)
,
DiagonalMat([z^30,z^27,z^9,z^19,z^14,z,z^21,z^13,z^10,z^26,z^31,z^26,
  z^30,z^24,z^7,z^9,z^16,z^28,z^7,z^11,1,z^17,z^8,z^27,z^32,z^25,z^28,
  z^18,z^20,z^5,z^22,1,z^26,z^16,z^23,z^11,z^12,z^22,z^2,z^3,z^3,z^21,
  z^17,z^25,z^29,z^24,z^15,z^18,z^32,z^31,z^29,z^15,z^4,z^2,z^13,z^5,
  z^14,z^20,z^8,z^23,z^4,z,z^10,z^6,z^12,z^6,z^19,z^7]) *
PermutationMat(
( 1,55,35)( 2,53,54)( 3,40,42)( 4,26,38)( 5,18,46)( 6,24,30)( 7,45,34)
( 8,39,48)( 9,17,15)(10,44,47)(11,51,66)(12,68,32)(13,14,65)(16,56,67)
(19,28,59)(20,23,57)(21,49,62)(22,60,33)(25,27,64)(29,41,63)(37,43,61)
(50,58,52), 68)];

return result;