cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 18 results. Next

A079006 Expansion of q^(-1/4) * (eta(q) * eta(q^4)^2 / eta(q^2)^3)^2 in powers of q.

Original entry on oeis.org

1, -2, 5, -10, 18, -32, 55, -90, 144, -226, 346, -522, 777, -1138, 1648, -2362, 3348, -4704, 6554, -9056, 12425, -16932, 22922, -30848, 41282, -54946, 72768, -95914, 125842, -164402, 213901, -277204, 357904, -460448, 590330, -754368, 960948, -1220370
Offset: 0

Views

Author

Michael Somos, Dec 22 2002

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
The Lagrange series reversion of Sum_{n >= 1} a(n-1)*x^n is Sum_{n >= 1} A002103(n-1)*x^n. See the example in A002103. - Wolfdieter Lang, Jul 09 2016

Examples

			G.f. A(x) = 1 - 2*x + 5*x^2 - 10*x^3 + 18*x^4 - 32*x^5 + 55*x^6 - 90*x^7 + 144*x^8 + ...
G.f. B(q) = q * A(q^4) = q - 2*q^5 + 5*q^9 - 10*q^13 + 18*q^17 - 32*q^21 + 55*q^25 - 90*q^29 + ...
		

References

  • A. Cayley, A memoir on the transformation of elliptic functions, Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 9, p. 128.
  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; Eq. (34.3).

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ Product[(1 + x^(k + 1)) / (1 + x^k), {k, 1, n, 2}]^2, {x, 0, n}]; (* Michael Somos, Jul 08 2011 *)
    a[ n_] := With[ {m = InverseEllipticNomeQ[ q]}, SeriesCoefficient[ (m / 16 / q)^(1/4), {q, 0, n}]]; (* Michael Somos, Jul 08 2011 *)
    QP = QPochhammer; s = (QP[q]*(QP[q^4]^2/QP[q^2]^3))^2 + O[q]^40; CoefficientList[s, q] (* Jean-François Alcover, Nov 23 2015 *)
    nmax = 50; CoefficientList[Series[Product[(1+x^(2*k))^4 / (1+x^k)^2, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jul 04 2016 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ x^4]^2 / QPochhammer[ -x]^2, {x, 0, n}]; (* Michael Somos, Apr 19 2017 *)
  • PARI
    {a(n) = my(N, A); if( n<0, 0, N = (sqrtint(16*n + 1) + 1)\2; A = contfracpnqn( matrix(2, N, i, j, if( i==1, if( j<2, 1 + O(x^(N^2 + N)), (x^(j-1) + x^(3*j - 3))^2), 1 - x^(4*j - 2)))); polcoeff( A[2,1] / A[1,1], 4*n))}; /* Michael Somos, Sep 01 2005 */
    
  • PARI
    {a(n) = my(A, m); if( n<0, 0, A = 1 + O(x); m = 1; while( m<=n, m*=2; A = subst(A, x, x^2); A = sqrt(A / (1 + 4 * x*A^2))); polcoeff(A, n))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x + A) * eta(x^4 + A)^2 / eta(x^2 + A)^3)^2, n))};

Formula

a(n) = (2/n)*Sum_{k=1..n} (-1)^k*A046897(k)*a(n-k). - Vladeta Jovovic, Dec 24 2002
Expansion of q^(-1/4) * (1/2) * k^(1/2) in powers of q, where k^2 is the parameter and q the Jacobi nome of elliptic functions.
Expansion of (1/(2*q)) * (1 - sqrt(k')) / (1 + sqrt(k')) in powers of q^4, where k'^2 is the complementary parameter and q the Jacobi nome of elliptic functions. See the Fricke reference.
Expansion of psi(x^2) / phi(x) = psi(x)^2 / phi(x)^2 = psi(x^2)^2 / psi(x)^2 = psi(-x)^2 / phi(-x^2)^2 = chi(-x)^2 / chi(-x^2)^4 = 1 / (chi(x)^2 * chi(-x^2)^2) = 1 / (chi(x)^4 * chi(-x)^2) = f(-x^4)^2 / f(x)^2 in powers of x where phi(), psi(), chi(), f() are Ramanujan theta functions.
Euler transform of period 4 sequence [-2, 4, -2, 0, ...].
G.f. A(x) satisfies A(x)^2 = A(x^2) / (1 + 4 * x * A(x^2)^2). - Michael Somos, Mar 19 2004
Given g.f. A(x), then B(q) = q * A(q^4) satisfies 0 = f(B(q), B(q^2)) where f(u, v) = u^2 * (1 + 4 * v^2) - v. - Michael Somos, Jul 09 2005
Given g.f. A(x), then B(q) = q * A(q^4) satisfies 0 = f(B(q), B(q^2), B(q^3), B(q^6)) where f(u1, u2, u3, u6) = u1*u3 * (u6 + u2)^2 - u2*u6. - Michael Somos, Jul 09 2005
G.f.: (Product_{k>0} (1 + x^(2*k)) / (1 + x^(2*k-1)))^2 = (Product_{k>0} (1 - x^(4*k)) / (1 - (-x)^k))^2.
Expansion of continued fraction 1 / (1 - x^2 + (x^1 + x^3)^2 / (1 - x^6 + (x^2 + x^6)^2 / (1 - x^10 + (x^3 + x^9)^2 / ...))) in powers of x^4. - Michael Somos, Sep 01 2005
Given g.f. A(x), then B(q) = 2 * q * A(q^4) satisfies 0 = f(B(q), B(q^3)) where f(u, v) = (1 - u^4) * (1 - v^4) - (1 - u*v)^4 . - Michael Somos, Jan 01 2006
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = (1/2) g(t) where q = exp(2 Pi i t) and g() is g.f. for A189925.
Convolution inverse of A029839. Convolution square of A083365. a(n) = (-1)^n * A001936(n).
G.f.: 1/Q(0), where Q(k)= 1 - x^(k+1/2) + (x^((k+1)/4) + x^((3*k+3)/4))^2/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 02 2013
a(n) ~ (-1)^n * exp(Pi*sqrt(n)) / (2^(7/2)*n^(3/4)). - Vaclav Kotesovec, Jul 04 2016
Given g.f. A(x), and B(x) is the g.f. for A008441, then A(x) = B(x^2) / B(x) and A(x) * A(x^2) * A(x^4) * ... = 1 / B(x). - Michael Somos, Apr 20 2017
Expansion of continued fraction 1 / (1 - x^1 + x^1*(1 + x^1)^2 / (1 - x^3 + x^2*(1 + x^2)^2 / (1 - x^5 + x^3*(1 + x^3)^2 / ...))) in powers of x^2. - Michael Somos, Apr 20 2017
a(n) = A208933(4*n+1) - A215348(4*n+1) (conjectured). - Thomas Baruchel, May 14 2018
A(x^4) = (1/(m*x)) * ( chi(x)^m - chi(-x)^m ) / ( chi(x)^m + chi(-x)^m ) at m = 2, where chi(x) = Product_{i >= 0} (1 + x^(2*i+1)) is the g.f. of A000700. The formula gives generating functions related to A092869 when m = 1 and A001938 (also A093160) when m = 4. - Peter Bala, Sep 23 2023

A001938 Expansion of k/(4*q^(1/2)) in powers of q, where k defined by sqrt(k) = theta_2(0, q)/theta_3(0, q).

Original entry on oeis.org

1, -4, 14, -40, 101, -236, 518, -1080, 2162, -4180, 7840, -14328, 25591, -44776, 76918, -129952, 216240, -354864, 574958, -920600, 1457946, -2285452, 3548550, -5460592, 8332425, -12614088, 18953310, -28276968, 41904208, -61702876, 90304598, -131399624
Offset: 0

Views

Author

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
k^2 is the parameter and q the Jacobi nome of elliptic functions. See, e.g., Fricke, p. 11, eq. (8) with p. 10. eq. (1). - Wolfdieter Lang, Jul 04 2016

Examples

			G.f. = 1 - 4*x + 14*x^2 - 40*x^3 + 101*x^4 - 236*x^5 + 518*x^6 - 1080*x^7 + ...
G.f. of B(q) = q * A(q^2): q - 4*q^3 + 14*q^5 - 40*q^7 + 101*q^9 - 236*q^11 + 518*q^13 - 1080*q^15 + ...
		

References

  • A. Cayley, A memoir on the transformation of elliptic functions, Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 9, p. 128.
  • E. T. Copson, An Introduction to the Theory of Functions of a Complex Variable, 1935, Oxford University Press, p. 385.
  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; Eq. (34.3).
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ 1 / (QPochhammer[ -x, x^2] QPochhammer[ x^2, x^4])^4, {x, 0, n}]; (* Michael Somos, Sep 24 2011 *)
    a[ n_] := SeriesCoefficient[ (QPochhammer[ x^4] / QPochhammer[ -x])^4, {x, 0, n}]; (* Michael Somos, Sep 24 2011 *)
    a[ n_] := SeriesCoefficient[ (Product[ 1 - x^k, {k, 4, n, 4}] / Product[ 1 - (-x)^k, {k, n}])^4, {x, 0, n}]; (* Michael Somos, Sep 24 2011 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 2, 0, q^(1/2)] / (2 EllipticTheta[ 3, 0, q]))^4, {q, 0, n + 1/2}]; (* Michael Somos, Sep 24 2011 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 2, 0, q] / EllipticTheta[ 2, 0, q^(1/2)])^4, {q, 0, n + 1/2}]; (* Michael Somos, Sep 24 2011 *)
  • PARI
    {a(n) = my(A, A2, m); if( n<0, 0, n = 2*n + 1; A = x + O(x^3); m=2; while( mMichael Somos, Mar 26 2004 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x + A) * eta(x^4 + A)^2 / eta(x^2 + A)^3)^4, n))}; /* Michael Somos, Mar 26 2004 */

Formula

Expansion of (psi(x^2) / phi(x))^2 = (psi(x) / phi(x))^4 = (psi(x^2) / psi(x))^4 = (psi(-x) / psi(-x^2))^4 = (chi(-x) / chi(-x^2)^2)^4 = (chi(x)^2 * chi(-x))^-4 = (chi(x) * chi(-x^2))^-4 = (f(-x^4) / f(x))^4 in powers of x where phi(), psi(), chi(), f() are Ramanujan theta functions. - Michael Somos, Feb 26 2012
G.f. A(x) satisfies 1 = (1 - 16 * x * A(x)^2) * (1 + 16 * x * A(-x)^2). - Michael Somos, Mar 26 2004
Expansion of q^(-1/2) * (eta(q) * eta(q^4)^2 / eta(q^2)^3)^4 in powers of q.
Euler transform of period 4 sequence [ -4, 8, -4, 0, ...].
Given g.f. A(x), then B(q) = q * A(q^2) satisfies 0 = f(B(q), B(q^2)) where f(u, v) = v - (u * (1 + 4*v))^2. - Michael Somos, Mar 26 2004
G.f. A(q) satisfies A(q) = sqrt(A(q^2)) / (1 + 4*q*A(q^2)); together with limit_{n->infinity} A(x^n) = 1 this gives a fast algorithm to compute the series. - Joerg Arndt, Aug 06 2011
G.f.: (Product_{k>0} (1 + x^(2*k)) / (1 + x^(2*k - 1)))^4.
a(n) = (-1)^n * A093160(n). Convolution square of A079006.
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = 1/4 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A139820. - Michael Somos, Jun 04 2015
G.f.: ((Sum_{n >= 0} x^(n*(n+1))) / (1 + Sum_{n >= 1} x^(n^2)))^4 (from the sum representation of the Jacobi theta functions evaluated at vanishing argument). - Wolfdieter Lang, Jul 04 2016
a(n) ~ (-1)^n * exp(sqrt(2*n)*Pi) / (32 * 2^(1/4) * n^(3/4)). - Vaclav Kotesovec, Nov 15 2017

Extensions

Edited by N. J. A. Sloane, Mar 31 2007

A092877 Expansion of (eta(q^4) / eta(q))^8 in powers of q.

Original entry on oeis.org

1, 8, 44, 192, 718, 2400, 7352, 20992, 56549, 145008, 356388, 844032, 1934534, 4306368, 9337704, 19771392, 40965362, 83207976, 165944732, 325393024, 628092832, 1194744096, 2241688744, 4152367104, 7599231223, 13749863984
Offset: 1

Views

Author

Michael Somos, Mar 19 2004

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = q + 8*q^2 + 44*q^3 + 192*q^4 + 718*q^5 + 2400*q^6 + 7352*q^7 + 20992*q^8 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 0, 0, SeriesCoefficient[ -InverseEllipticNomeQ[ -x] / 16, {x, 0, n}]]; (* Michael Somos, Jun 13 2011 *)
    a[ n_] := If[ n < 0, 0, SeriesCoefficient[ With[ {lambda = ModularLambda[ Log[x] / ( Pi I)]}, lambda / (16 * (1 - lambda))], {x, 0, n}]]; (* Michael Somos, Jun 13 2011 *)
    a[ n_] := SeriesCoefficient[ q (QPochhammer[ q^4] / QPochhammer[ q])^8, {q, 0, n}]; (* Michael Somos, Aug 09 2015 *)
    a[1] = 1; a[n_] := a[n] = (8/(n-1))*Sum[DivisorSum[k, Identity, Mod[#, 4] != 0&]*a[n-k], {k, 1, n-1}]; Array[a, 26] (* Jean-François Alcover, Mar 01 2018, after Seiichi Manyama *)
    eta[q_]:= q^(1/6) QPochhammer[q]; a[n_]:=SeriesCoefficient[(eta[q^4] / eta[q])^8, {q, 0, n}]; Table[a[n], {n, 4, 35}] (* Vincenzo Librandi, Oct 18 2018 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( x * prod(k=1, (n+1)\2, (1 + x^(2*k)) / (1 - x^(2*k-1)), 1 + x * O(x^n))^8, n))};
    
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( (eta(x^4 + A) / eta(x + A))^8, n))};

Formula

Expansion of q * (psi(-q) / phi(-q))^8 = q * (psi(q^2) / psi(-q))^8 = q * (psi(q) / phi(-q^2))^8 = q * (psi(q^2) / phi(-q))^4 = q * (chi(q) / chi(-q^2)^2)^8 = q / (chi(-q) * chi(-q^2))^8 = q / (chi(q) * chi(-q)^2)^8 = q * (f(-q^4) / f(-q))^8 in powers of q where phi(), psi(), chi(), f() are Ramanujan theta functions. - Michael Somos, Jun 13 2011
Euler transform of period 4 sequence [ 8, 8, 8, 0, ...].
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v) = u^2 - v - 16*u*v - 16*v^2 - 256*u*v^2.
G.f.: x * (Product_{k>0} (1 + x^(2*k)) / (1 - x^(2*k - 1)))^8.
G.f.: theta_2^4 / (16*theta_4^4) = lambda / (16 * (1 - lambda)).
G.f.: exp( Integral theta_3(x)^4/x dx ). - Paul D. Hanna, May 03 2010
a(n) = (-1)^n * A005798(n).
a(2*n) = 8 * A014103(n). - Michael Somos, Aug 09 2015
Convolution inverse of A124972, 8th power of A001935, 4th power of A001936, square of A093160. - Michael Somos, Aug 09 2015
a(n) ~ exp(2*Pi*sqrt(n))/(512*n^(3/4)). - Vaclav Kotesovec, Sep 07 2015
a(1) = 1, a(n) = (8/(n-1))*Sum_{k=1..n-1} A046897(k)*a(n-k) for n > 1. - Seiichi Manyama, Apr 01 2017

A296044 a(n) = [x^n] Product_{k>=1} ((1 + x^(2*k))/(1 - x^(2*k-1)))^n.

Original entry on oeis.org

1, 1, 5, 22, 101, 481, 2330, 11425, 56549, 281911, 1413465, 7120136, 36006362, 182681916, 929461993, 4740491107, 24229115109, 124069449335, 636376573943, 3268955179686, 16814509004601, 86593280920756, 446437797872016, 2303948443259841, 11900990745759578, 61526182236027756
Offset: 0

Views

Author

Ilya Gutkovskiy, Dec 03 2017

Keywords

Crossrefs

Programs

  • Mathematica
    Table[SeriesCoefficient[Product[((1 + x^(2 k))/(1 - x^(2 k - 1)))^n, {k, 1, n}], {x, 0, n}], {n, 0, 25}]
    Table[SeriesCoefficient[Product[((1 - x^(4 k))/(1 - x^k))^n, {k, 1, n}], {x, 0, n}], {n, 0, 25}]
    Table[SeriesCoefficient[(EllipticTheta[2, 0, x]/EllipticTheta[2, Pi/4, x^(1/2)]/(16 x)^(1/8))^n, {x, 0, n}], {n, 0, 25}]
    (* Calculation of constant d: *) With[{k = 4}, 1/r /. FindRoot[{s == QPochhammer[(r*s)^k] / QPochhammer[r*s], k*(-(s*QPochhammer[r*s]*(Log[1 - (r*s)^k] + QPolyGamma[0, 1, (r*s)^k]) / Log[(r*s)^k]) + (r*s)^k * Derivative[0, 1][QPochhammer][(r*s)^k, (r*s)^k]) == s*QPochhammer[r*s] + s^2*(-(QPochhammer[r*s]*(Log[1 - r*s] + QPolyGamma[0, 1, r*s]) / (s*Log[r*s])) + r*Derivative[0, 1][QPochhammer][r*s, r*s])}, {r, 1/5}, {s, 1}, WorkingPrecision -> 70]] (* Vaclav Kotesovec, Jan 17 2024 *)

Formula

a(n) = [x^n] Product_{k>=1} ((1 - x^(4*k))/(1 - x^k))^n.
a(n) ~ c * d^n / sqrt(n), where d = 5.2749356339591798618290252741994029798069148326559... and c = 0.2726256757090475625917361066565981461455343437... - Vaclav Kotesovec, Dec 05 2017

A001937 Expansion of (psi(x^2) / psi(-x))^3 in powers of x where psi() is a Ramanujan theta function.

Original entry on oeis.org

1, 3, 9, 22, 48, 99, 194, 363, 657, 1155, 1977, 3312, 5443, 8787, 13968, 21894, 33873, 51795, 78345, 117312, 174033, 255945, 373353, 540486, 776848, 1109040, 1573209, 2218198, 3109713, 4335840, 6014123, 8300811, 11402928, 15593702, 21232521, 28790667, 38884082
Offset: 0

Views

Author

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
The Cayley reference is actually to A187053. - Michael Somos, Jul 26 2012

Examples

			1 + 3*x + 9*x^2 + 22*x^3 + 48*x^4 + 99*x^5 + 194*x^6 + 363*x^7 + 657*x^8 + ...
q^3 + 3*q^11 + 9*q^19 + 22*q^27 + 48*q^35 + 99*q^43 + 194*q^51 + 363*q^59 + ...
		

References

  • A. Cayley, A memoir on the transformation of elliptic functions, Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 9, p. 128.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Maple
    g100:= mul((1+x^(2*k))/(1-x^(2*k-1)),k=1..50)^3:
    S:= series(g100,x,101):
    seq(coeff(S,x,j),j=0..100); # Robert Israel, Nov 30 2015
  • Mathematica
    CoefficientList[ Series[Product[(1 - x^k)^(-3*Boole[Mod[k, 4] != 0]), {k, 1, 101}], {x, 0, 100}], x] (* Olivier GERARD, May 06 2009 *)
    QP = QPochhammer; s = (QP[q^4]/QP[q])^3 + O[q]^40; CoefficientList[s, q] (* Jean-François Alcover, Nov 30 2015, adapted from PARI *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^4 + A) / eta(x + A))^3, n))} /* Michael Somos, Mar 06 2011 */

Formula

Expansion of q^(-3/8) * (eta(q^4) / eta(q))^3 in powers of q. - Michael Somos, Jul 26 2012
Euler transform of period 4 sequence [ 3, 3, 3, 0, ...]. - Michael Somos, Mar 06 2011
Convolution cube of A001935. A187053(n) = (-1)^n * a(n). - Michael Somos, Mar 06 2011
G.f.: (Product_{k>0} (1 + x^(2*k)) / (1 - x^(2*k-1)))^3.
a(n) ~ 3^(1/4) * exp(sqrt(3*n/2)*Pi) / (16*2^(3/4)*n^(3/4)). - Vaclav Kotesovec, Nov 15 2017

Extensions

Corrected and extended by Simon Plouffe
Checked and more terms from Olivier GERARD, May 06 2009

A001940 Absolute value of coefficients of an elliptic function.

Original entry on oeis.org

1, 6, 27, 98, 309, 882, 2330, 5784, 13644, 30826, 67107, 141444, 289746, 578646, 1129527, 2159774, 4052721, 7474806, 13569463, 24274716, 42838245, 74644794, 128533884, 218881098, 368859591, 615513678, 1017596115, 1667593666, 2710062756, 4369417452
Offset: 0

Views

Author

Keywords

References

  • A. Cayley, A memoir on the transformation of elliptic functions, Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 9, p. 128.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    nn = 4*10; b = Flatten[Table[{6, 6, 6, 0}, {nn/4}]]; CoefficientList[x*Series[Product[1/(1 - x^m)^b[[m]], {m, nn}], {x, 0, nn}], x] (* T. D. Noe, Aug 17 2012 *)
    nmax = 40; CoefficientList[Series[Product[((1 - x^(4*k)) / (1 - x^k))^6, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Nov 15 2017 *)

Formula

G.f.: Product ( 1 - x^k )^(-c(k)), c(k) = 6, 6, 6, 0, 6, 6, 6, 0, ....
a(n) ~ 3^(1/4) * exp(sqrt(3*n)*Pi) / (128*sqrt(2)*n^(3/4)). - Vaclav Kotesovec, Nov 15 2017
G.f.: Product_{k>=1} ((1 + x^(2*k))/(1 - x^(2*k-1)))^6. - Ilya Gutkovskiy, Dec 04 2017

Extensions

Extended and corrected by Simon Plouffe

A001941 Absolute values of coefficients of an elliptic function.

Original entry on oeis.org

1, 7, 35, 140, 483, 1498, 4277, 11425, 28889, 69734, 161735, 362271, 786877, 1662927, 3428770, 6913760, 13660346, 26492361, 50504755, 94766875, 175221109, 319564227, 575387295, 1023624280, 1800577849, 3133695747, 5399228149, 9214458260, 15584195428
Offset: 0

Views

Author

Keywords

References

  • A. Cayley, A memoir on the transformation of elliptic functions, Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 9, p. 128.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    nn = 4*10; b = Flatten[Table[{7, 7, 7, 0}, {nn/4}]]; CoefficientList[x*Series[Product[1/(1 - x^m)^b[[m]], {m, nn}], {x, 0, nn}], x] (* T. D. Noe, Aug 17 2012 *)
    nmax = 40; CoefficientList[Series[Product[((1 - x^(4*k)) / (1 - x^k))^7, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Nov 15 2017 *)

Formula

G.f.: Product ( 1 - x^k )^-c(k), c(k) = 7, 7, 7, 0, 7, 7, 7, 0, ....
a(n) ~ 7^(1/4) * exp(sqrt(7*n/2)*Pi) / (256*2^(3/4)*n^(3/4)). - Vaclav Kotesovec, Nov 15 2017
G.f.: Product_{k>=1} ((1 + x^(2*k))/(1 - x^(2*k-1)))^7. - Ilya Gutkovskiy, Dec 04 2017

A123655 Expansion of q * psi(q^8) / phi(-q) in powers of q where psi(), phi() are Ramanujan theta functions.

Original entry on oeis.org

1, 2, 4, 8, 14, 24, 40, 64, 101, 156, 236, 352, 518, 752, 1080, 1536, 2162, 3018, 4180, 5744, 7840, 10632, 14328, 19200, 25591, 33932, 44776, 58816, 76918, 100176, 129952, 167936, 216240, 277476, 354864, 452392, 574958, 728568, 920600, 1160064
Offset: 1

Views

Author

Michael Somos, Oct 04 2006

Keywords

Comments

Ramanujan theta functions: phi(q) (A000122), psi(q) (A010054).
Number 12 of the 14 eta-quotients listed in Table 2 of Moy 2013. - Michael Somos, Sep 19 2013

Examples

			G.f. = q + 2*q^2 + 4*q^3 + 8*q^4 + 14*q^5 + 24*q^6 + 40*q^7 + 64*q^8 + 101*q^9 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 2, 0, q^4] / EllipticTheta[ 4, 0, q] / 2, {q, 0, n}]; (* Michael Somos, Sep 19 2013 *)
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^16 + A)^2 / (eta(x + A)^2 * eta(x^8 + A)), n))};

Formula

Expansion of eta(q^2) * eta(q^16)^2 / (eta(q)^2 * eta(q^8)) in powers of q.
Euler transform of period 16 sequence [ 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 0, ...].
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v) = u^2 - v * (1 + 4*u) * (1 + 2*v).
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = 1/8 * g(t) where q = exp(2 Pi i t) and g() is g.f. for A185338.
a(n) is odd iff n is an odd square. If n>2 is a power of 2 then the highest power of 2 dividing a(n) is (n/2)^3. - Michael Somos, Feb 18 2007
4 * a(n) = A007096(n) unless n=0. -(-1)^n * a(n) = A208605(n). Convolution inverse of A185338.
G.f.: x * Product_{k>0} (1 + x^k)^2 * (1 + x^(2*k)) * (1 + x^(4*k)) * (1 + x^(8*k))^2. Michael Somos, Sep 19 2013
a(2*n) = 2 * A107035(n). a(2*n + 1) = A093160(n). - Michael Somos, Sep 19 2013
a(n) ~ exp(sqrt(n)*Pi) / (2^(9/2) * n^(3/4)). - Vaclav Kotesovec, Nov 15 2017

A285928 Expansion of (Product_{k>0} (1 - x^(5*k)) / (1 - x^k))^5 in powers of x.

Original entry on oeis.org

1, 5, 20, 65, 190, 501, 1240, 2890, 6440, 13775, 28502, 57205, 111880, 213670, 399620, 733128, 1321850, 2345340, 4100700, 7072520, 12045005, 20272465, 33746060, 55595635, 90706390, 146638756, 235016940, 373580735, 589238640, 922537655, 1434232510, 2214817165
Offset: 0

Views

Author

Seiichi Manyama, Apr 28 2017

Keywords

Comments

In general, if m > 1 and g.f. = Product_{k>=1} ((1 - x^(m*k)) / (1 - x^k))^m, then a(n, m) ~ exp(Pi*sqrt(2*(m-1)*n/3)) * (m-1)^(1/4) / (2^(5/4) * 3^(1/4) * m^(m/2) * n^(3/4)). - Vaclav Kotesovec, Apr 30 2017

Crossrefs

(Product_{k>0} (1 - x^(m*k)) / (1 - x^k))^m: A022567 (m=2), A285927 (m=3), A093160 (m=4), this sequence (m=5).

Programs

  • Mathematica
    nmax = 40; CoefficientList[Series[Product[((1 - x^(5*k)) / (1 - x^k))^5, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Apr 30 2017 *)

Formula

a(0) = 1, a(n) = (5/n)*Sum_{k=1..n} A116073(k)*a(n-k) for n > 0.
a(n) ~ exp(2*Pi*sqrt(2*n/3)) / (2^(3/4) * 3^(1/4) * 5^(5/2) * n^(3/4)). - Vaclav Kotesovec, Apr 30 2017

A112128 Expansion of phi(q^4) / phi(q) in powers of q where phi() is a Ramanujan theta function.

Original entry on oeis.org

1, -2, 4, -8, 16, -28, 48, -80, 128, -202, 312, -472, 704, -1036, 1504, -2160, 3072, -4324, 6036, -8360, 11488, -15680, 21264, -28656, 38400, -51182, 67864, -89552, 117632, -153836, 200352, -259904, 335872, -432480, 554952, -709728, 904784, -1149916, 1457136
Offset: 0

Views

Author

Michael Somos, Aug 27 2005

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 - 2*q + 4*q^2 - 8*q^3 + 16*q^4 - 28*q^5 + 48*q^6 - 80*q^7 + 128*q^8 + ...
		

Crossrefs

Programs

  • Mathematica
    QP = QPochhammer; s = QP[q]^2*(QP[q^8]^5/QP[q^2]^5/QP[q^16]^2) + O[q]^40; CoefficientList[s, q] (* Jean-François Alcover, Nov 30 2015, adapted from PARI *)
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q^4] / EllipticTheta[ 3, 0, q], {q, 0, n}]; (* Michael Somos, Dec 11 2016 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^8 + A)^5 / (eta(x^2 + A)^5 * eta(x^16 + A)^2), n))};

Formula

Expansion of (eta(q) / eta(q^16))^2 * (eta(q^8) / eta(q^2))^5 in powers of q.
Euler transform of period 16 sequence [ -2, 3, -2, 3, -2, 3, -2, -2, -2, 3, -2, 3, -2, 3, -2, 0, ...].
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v) = u^2 - (1 - 2*u + 2*u^2) * (1 - 2*v + 2*v^2).
G.f.: (Sum_{k in Z} x^(4*k^2)) / (Sum_{k in Z} x^(k^2)) = theta_3(0, x^4) / theta_3(0, x).
G.f.: Product_{k>0} ((1 + x^(2*k)) * (1 + x^(4*k)))^3 / ((1 + x^k) * (1 + x^(8*k)))^2.
Expansion of continued fraction 1 / (1 + 2*x / (1 - x^2 + (x^1 + x^3)^2 / (1 - x^6 + (x^2 + x^6)^2 / (1 - x^10 + (x^3 + x^9)^2 / ...)))).
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = 1/2 * g(t) where q = exp(2 Pi i t) and g() is the g.f. for A208724.
(-1)^n * a(n) = A208933(n). a(2*n) = A131126(n). a(2*n + 2) = -2 * A093160(n). - Michael Somos, Dec 11 2016
Convolution inverse of A208274. - Michael Somos, Dec 11 2016
a(n) ~ (-1)^n * exp(sqrt(n)*Pi) / (2^(7/2) * n^(3/4)). - Vaclav Kotesovec, Nov 15 2017
Showing 1-10 of 18 results. Next