My (old, PASCAL) programs for invariants are for linear groups over
finite fields, but the algorithms should work for linear groups over an
arbitrary field. To be more precise, given matrices generating a
subgroup G of GL(n,q), and a natural mnumber k, my programs calculate
a basis for the vector space of the n-variable degree k polynomials
invariant under G.

Regards, Leonard Soicher.