%I #6 Apr 05 2018 20:31:37
%S 1,276,11178,177400,1612875,10131156,48897678,193740408,658523925,
%T 1980143600,5386270686,13477895856,31425764410,68969957700,
%U 143635113000,285718115112,545796171084,1005775268868,1794713445350,3111031518000
%N Expansion of (psi(x)^24 + psi(-x)^24) / 2 in powers of x^2 where psi() is a Ramanujan theta function.
%H G. C. Greubel, <a href="/A169976/b169976.txt">Table of n, a(n) for n = 0..1000</a>
%F a(n) = (A013959(2*n + 3) - A000594(2*n + 3)) / 176896 = A014809(2*n).
%e 1 + 276*x + 11178*x^2 + 177400*x^3 + 1612875*x^4 + 10131156*x^5 + ...
%e q^3 + 276*q^5 + 11178*q^7 + 177400*q^9 + 1612875*q^11 + 10131156*q^13 + ...
%t QP:= Pochhammer; a[n_]:= SeriesCoefficient[(QP[q, q])^24*(QP[-q^(1/2), q^(1/2)]^24 + QP[q^(1/2), -q^(1/2)]^24)/2, {q, 0, n}]; Table[a[n], {n,0,50}] (* _G. C. Greubel_, Apr 04 2018 *)
%o (PARI) {a(n) = local(A); if( n<0, 0, n *= 2; A = x * O(x^n); polcoeff( (eta(x^2 + A)^2 / eta(x + A))^24, n))}
%o (PARI) {a(n) = if( n<0, 0, n = 2*n + 3; (sigma(n, 11) - polcoeff( x * eta(x + x * O(x^n))^24, n)) / 176896 )}
%K nonn
%O 0,2
%A _Michael Somos_, Aug 15 2010