print "4 by 4 matrices over GF(9).  U(4,3) with wrong indexing poly";
print "Order 1620 = 2^2 * 3^4 * 5";
print "Field: K; Primitive element: w; Group: G.";

K<w> := ExtensionField<GF(3), x | x^2 + x + 2>;
G := MatrixGroup<4, K |
    [ w^2, w^1, w^2, 0, 1, w^3, 1, 0, w^7, 0, w^3, 0, 0, 0, 0, 1 ],
    [ 1, w^2, w^6, 2, w^6, 2, 2, w^6, w^2, 2, 2, w^2, 2, w^2, w^6, 1 ]
>;
