a[ n_] := With[ {A = q (QPochhammer[ q^7] / QPochhammer[ q])^4}, SeriesCoefficient[ (1 + 7 A)^2 / A, {q, 0, n}]]; (* Michael Somos, Mar 30 2015 *)
PROG
(PARI) {a(n) = my(A); if( n<-1, 0, n++; A = x * O(x^n); A = (eta(x^7 + A) / eta(x + A))^4; polcoeff( (1 + 7 * x * A)^2 / A, n))}; /* Michael Somos, Feb 02 2012 */