with(qseries);
  [aqprod, changes, etamake, etaq, findcong, findhom, findhomcombo,
        findhomcombomodp, findhommodp, findlincombo, findlincombomodp,
        findnonhom, findnonhomcombo, findpoly, jac2prod, jac2series,
        jacprod, jacprodmake, mprodmake, packageversion, prodmake,
        qbin, qdegree, qetamake, qfactor, qs2jaccombo, quinprod, sift,
        theta, theta2, theta3, theta4, tripleprod, winquist, zqfactor
        ]

> x:=add(q^(n*(n+1)/2),n=0..20);

                3    6    10    15    21    28    36    45    55
  x := 1 + q + q  + q  + q   + q   + q   + q   + q   + q   + q

            66    78    91    105    120    136    153    171    190
         + q   + q   + q   + q    + q    + q    + q    + q    + q

            210
         + q

> etamake(x,q,20);

                                       2
                             eta(2 tau)
                           ---------------
                            (1/8)
                           q      eta(tau)

> qetamake(x,q,20);

                                   2 2
                               _E(q )
                               -------
                                _E(q)

> P:=series(1/etaq(q,1,1001),q,1001):
> findcong(P,1000);

                              [4, 5, 5]


                              [5, 7, 7]


                             [6, 11, 11]


                             [24, 25, 25]


          {[6, 11, 11], [24, 25, 25], [4, 5, 5], [5, 7, 7]}

> qdegree(P);
                                     1000