Example (2) Find the GF for p   (n).
==========                   3,3

==============================================================================
> x:=mul((1-q^(4*k))/(1-q^k),k=1..3); 
                                   4        8        12
                             (1 - q ) (1 - q ) (1 - q  )
                        x := ---------------------------
                                            2        3
                              (1 - q) (1 - q ) (1 - q )

> normal(x);
  8    7    6    2            7    6    5    4    3    2
(q  - q  + q  + q  - q + 1) (q  + q  + q  + q  + q  + q  + q + 1)

      3    2
    (q  + q  + q + 1)

> expand(%);
           4      8      2      12      3      10      6      7      5      15
1 + q + 3 q  + 5 q  + 2 q  + 5 q   + 3 q  + 5 q   + 5 q  + 5 q  + 4 q  + 3 q

          14      13      11      9    18    17      16
     + 3 q   + 4 q   + 5 q   + 6 q  + q   + q   + 2 q

> sort(%);
 18    17      16      15      14      13      12      11      10      9      8
q   + q   + 2 q   + 3 q   + 3 q   + 4 q   + 5 q   + 5 q   + 5 q   + 6 q  + 5 q

          7      6      5      4      3      2
     + 5 q  + 5 q  + 4 q  + 3 q  + 3 q  + 2 q  + q + 1
==============================================================================

The last output is the required GF.