Cycle index polynomial of PGLk(q)

This routine computes the cycle index polynomial of the natural action of   PGLk(q) on   PGk-1(q).

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

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