Cycle index polynomial of GLk(q)

This routine computes the cycle index polynomial of the natural action of  GLk(q) on GF(q)k.

please enter the dimension
k=
please enter the size of the finite field
q=

Please use small parameters, for example: k=3,q=5.