
print "The unitary group U( 3, 3 ) represented as a permutation";
print "group of degree 28. ";
print "Order: 6,048 = 2^5 * 3^3 * 7;  Base: 1, 2, 3. ";
print "Group: G";

G := PermutationGroup<28 |
    \[ 1,  4,  3,  6,  8,  9, 12, 11,  5, 16, 19, 21, 18, 23, 10, 14,
 7, 17,  2, 22, 20, 13, 26, 15, 24, 25, 28, 27],
    \[ 2,  5,  3,  7,  1, 10, 13, 12, 15, 17, 20, 22,  4, 23, 24, 25,
 6, 21, 14, 26, 27,  8, 19,  9, 28, 11, 18, 16],
    \[ 3,  5,  1,  4,  2, 11, 14,  8,  9, 18,  6, 15, 23,  7, 12, 20,
25, 10, 19, 16, 26, 24, 13, 22, 17, 21, 28, 27]>;
AssertAttribute(G, "Order", 2^5 * 3^3 * 7);
