login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A186100 Expansion of 2 * a(q^2)^2 - a(q)^2 in powers of q where a() is a cubic AGM theta function. 5
1, -12, -12, -12, -12, -72, -12, -96, -12, -12, -72, -144, -12, -168, -96, -72, -12, -216, -12, -240, -72, -96, -144, -288, -12, -372, -168, -12, -96, -360, -72, -384, -12, -144, -216, -576, -12, -456, -240, -168, -72, -504, -96, -528, -144, -72, -288 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Ramanujan's Eisenstein series: P(q) (see A006352), Q(q) (A004009), R(q) (A013973).
LINKS
J. M. Borwein and P. B. Borwein, A cubic counterpart of Jacobi's identity and the AGM, Trans. Amer. Math. Soc., 323 (1991), no. 2, 691-701.
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
FORMULA
Expansion of b(x) * b(x^2) - c(x) * c(x^2) in powers of x where b(), c() are cubic AGM functions.
Expansion of (phi(-x) * phi(-x^3))^2 - 8 * x * (psi(x) * psi(x^3))^2 in powers of x where phi(), psi() are Ramanujan theta functions.
Expansion of (P(q) - 2*P(q^2) - 3*P(q^3) + 6*P(q^6)) / 2 in powers of q where P() is a Ramanujan Eisenstein series. - Michael Somos, Jul 07 2015
a(n) = -12 * A186099(n) if n>0. a(2*n) = a(n). a(2*n + 1) = - A008653(2*n + 1). a(n) = 2 * A008653(n) - A008653(2*n) = A131946(n) - 8 * A111932(n) = A131943(n) - 9 * A121443(n).
a(3*n) = a(n). a(6*n + 5) = -72 * A098098(n).- Michael Somos, Jul 07 2015
EXAMPLE
G.f. = 1 - 12*q - 12*q^2 - 12*q^3 - 12*q^4 - 72*q^5 - 12*q^6 - 96*q^7 + ...
MATHEMATICA
a[ n_] := If[ n < 1, Boole[n == 0], -12 DivisorSum[ n, # Boole[ 1 == GCD[#, 6]] &]]; (* Michael Somos, Jul 07 2015 *)
a[ n_] := SeriesCoefficient[(EllipticTheta[ 4, 0, x] EllipticTheta[ 4, 0, x^3])^2 - 1/2 (EllipticTheta[ 2, 0, x^(1/2)] EllipticTheta[ 2, 0, x^(3/2)])^2, {x, 0, n}]; (* Michael Somos, Jul 07 2015 *)
PROG
(PARI) {a(n) = if( n<1, n==0, -12 * sumdiv( n, d, d * (1 == gcd( d, 6) ) ) )};
(PARI) {a(n) = if( n<1, n==0, -12 * direuler( p=2, n, 1 / (1 - X) / (1 - (p>3) * p * X)) [n])};
CROSSREFS
Sequence in context: A028996 A113595 A109053 * A334245 A134324 A268914
KEYWORD
sign
AUTHOR
Michael Somos, Feb 12 2011
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 04:13 EDT 2024. Contains 371235 sequences. (Running on oeis4.)