Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #14 Jul 31 2023 07:23:36
%S 0,1,-3,-27,147,-183,-54,-423,1434,1323,-4141,-2580,-972,20394,-1929,
%T -20169,-15597,-7746,58212,-27964,76398,-19683,-117585,-94050,-60696,
%U 474840,-103989,-89910,-72487,148359,100440,-1107363,43554,629289,1257516,517584,-745848,-1430837,-748068
%N Fourier coefficients of the modular form (1/t_{6a}^3) * (1-12*sqrt(-3) / t_{6a})^(1/6) * F_{6a}^10.
%C Here, F_{6a} is the hypergeometric function F(1/3, 1/2; 1; 12*sqrt(-3)/t_{6a}). The definition given on page 23 in the linked manuscript has a minor typo where "t_{3A}" should be "t_{6a}". - _Robin Visser_, Jul 31 2023
%H Masao Koike, <a href="/A004016/a004016.pdf">Modular forms on non-compact arithmetic triangle groups</a>, Unpublished manuscript [Extensively annotated with OEIS A-numbers by N. J. A. Sloane, Feb 14 2021. Sloane wrote 2005 on the first page but the internal evidence suggests 1997.] See page 30.
%o (Sage)
%o def a(n):
%o if n==0: return 0
%o theta2 = sum([1]+[2*x^(k^2/2) for k in range(1, n+1)])
%o theta3 = sum([2*x^((k^2 + k + 1/4)/2) for k in range(n)])
%o phi0 = theta2(x=x^2)*theta2(x=x^6) + theta3(x=x^2)*theta3(x=x^6)
%o phi1 = theta2(x=x^2)*theta3(x=x^6) + theta3(x=x^2)*theta2(x=x^6)
%o phi02, phi12 = phi0(x=x^2), phi1(x=x^2)
%o f = phi0*(phi12^3*(phi02^2 - phi12^2)^3)/8
%o return f.taylor(x, 0, n+1).coefficient(x^(n+1/2)) # _Robin Visser_, Jul 31 2023
%K sign
%O 0,3
%A _Robert C. Lyons_, Feb 15 2021
%E More terms from _Robin Visser_, Jul 31 2023