a[ n_] := SeriesCoefficient[ InverseEllipticNomeQ @ x, {x, 0, n}];
a[ n_] := If[ n < 0, 0, SeriesCoefficient[ ModularLambda[ Log[q] / (Pi I)], {q, 0, n}]];
a[ n_] := SeriesCoefficient[ (EllipticTheta[ 2, 0, q] / EllipticTheta[ 3, 0, q])^4, {q, 0, n}];
a[ n_] := SeriesCoefficient[ 1/16 (EllipticTheta[ 2, 0, q] / EllipticTheta[ 3, 0, q^2])^8, {q, 0, n}]; (*
Michael Somos, May 26 2016 *)