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-6 of 6 results.

A051577 a(n) = (2*n + 3)!!/3 = A001147(n+2)/3.

Original entry on oeis.org

1, 5, 35, 315, 3465, 45045, 675675, 11486475, 218243025, 4583103525, 105411381075, 2635284526875, 71152682225625, 2063427784543125, 63966261320836875, 2110886623587616875, 73881031825566590625, 2733598177545963853125, 106610328924292590271875
Offset: 0

Views

Author

Keywords

Comments

Row m = 3 of the array A(3; m,n) := (2*n+m)!!/m!!, m >= 0, n >= 0.

Crossrefs

Cf. A000165, A001147, A002866(n+1) (m=0..2 rows of A(3; m,n)).

Programs

  • GAP
    F:=Factorial;; List([0..25], n-> F(2*n+4)/(12*2^n*F(n+2)) ); # G. C. Greubel, Nov 12 2019
  • Magma
    F:=Factorial; [F(2*n+4)/(12*2^n*F(n+2)): n in [0..25]]; // G. C. Greubel, Nov 12 2019
    
  • Maple
    seq( doublefactorial(2*n+3)/3,n=0..10) ; # R. J. Mathar, Sep 29 2013
  • Mathematica
    Table[(2*n + 3)!!/3!!, {n, 0, 25}] (* G. C. Greubel, Jan 22 2017 *)
    a[n_] := Sum[(-1)^k*Binomial[2*n + 1, n + k]*StirlingS1[n + k + 1 ,k], {k , 1, n + 1}]; Flatten[Table[a[n], {n, 0, 18}]] (* Detlef Meya, Jan 17 2024 *)
  • PARI
    vector(26, n, (2*n+2)!/(6*2^n*(n+1)!) ) \\ G. C. Greubel, Nov 12 2019
    
  • Sage
    f=factorial; [f(2*n+4)/(12*2^n*f(n+2)) for n in (0..25)] # G. C. Greubel, Nov 12 2019
    

Formula

a(n) = (2*n + 3)!!/3!!.
E.g.f.: 1/(1 - 2*x)^(5/2).
a(n) ~ (4/3) * sqrt(2) * n^2 * 2^n * e^(-n) * n^n *{1 + (47/24)*n^(-1) + ...}. - Joe Keane (jgk(AT)jgk.org), Nov 23 2001
Ramanujan polynomials -psi_n(n, x) evaluated at 0. - Ralf Stephan, Apr 16 2004
a(n) = 2^(2 + n) * Gamma(n + 5/2)/(3 * sqrt(Pi)). - Gerson Washiski Barbosa, May 05 2010
From Peter Bala, May 26 2017: (Start)
D-finite with recurrence a(n+1) = (2*n + 5)*a(n) with a(0) = 1.
O.g.f. A(x) satisfies the Riccati differential equation 2*x^2*A(x)' = (1 - 5*x)*A(x) - 1 with A(0) = 1.
G.f. as an S-fraction: A(x) = 1/(1 - 5*x/(1 - 2*x/(1 - 7*x/(1 - 4*x/(1 - 9*x/(1 - 6*x/(1 - ... - (2*n+3)*x/(1 - 2*n*x/(1 - ...))))))))) (by Stokes, 1982).
Reciprocal as an S-fraction: 1/A(x) = 1/(1 + 5*x/(1 - 7*x/(1 - 2*x/(1 - 9*x/(1 - 4*x/(1 - 11*x/(1 - 6*x/(1 - ... - (2*n + 5)*x/(1 - 2*n*x/(1 - ...)))))))))). (End)
From Amiram Eldar, Dec 11 2022: (Start)
Sum_{n>=0} 1/a(n) = 3*(sqrt(e*Pi/2) * erf(1/sqrt(2)) - 1), where erf is the error function.
Sum_{n>=0} (-1)^n/a(n) = 3*(1 - sqrt(Pi/(2*e)) * erfi(1/sqrt(2))), where erfi is the imaginary error function. (End)
a(n) = Sum_{k=1..n+1} (-1)^k*binomial(2*n + 1, n + k)*Stirling1(n + k + 1, k). - Detlef Meya, Jan 17 2024

A051578 a(n) = (2*n+4)!!/4!!, related to A000165 (even double factorials).

Original entry on oeis.org

1, 6, 48, 480, 5760, 80640, 1290240, 23224320, 464486400, 10218700800, 245248819200, 6376469299200, 178541140377600, 5356234211328000, 171399494762496000, 5827582821924864000, 209792981589295104000, 7972133300393213952000, 318885332015728558080000
Offset: 0

Views

Author

Keywords

Comments

Row m=4 of the array A(3; m,n) := (2*n+m)!!/m!!, m >= 0, n >= 0.

Crossrefs

Cf. A000165, A001147(n+1), A002866(n+1), A051577 (rows m=0..3), A051579, A051580, A051581, A051582, A051583.
Cf. A052587 (essentially the same).

Programs

  • GAP
    List([0..20], n-> 2^(n-1)*Factorial(n+2) ); # G. C. Greubel, Nov 11 2019
  • Magma
    [2^(n-1)*Factorial(n+2): n in [0..20]]; // G. C. Greubel, Nov 11 2019
    
  • Maple
    a:= proc(n) option remember; `if`(n=0, 1, 2*(n+2)*a(n-1)) end:
    seq(a(n), n=0..20);  # Alois P. Heinz, Apr 29 2019
    seq(2^(n-1)*(n+2)!, n=0..20); # G. C. Greubel, Nov 11 2019
  • Mathematica
    Table[2^(n-1)(n+2)!, {n,0,20}] (* Jean-François Alcover, Oct 05 2019 *)
    Table[(2n+4)!!/8,{n,0,20}] (* Harvey P. Dale, Apr 06 2023 *)
  • PARI
    vector(21, n, 2^(n-2)*(n+1)! ) \\ G. C. Greubel, Nov 11 2019
    
  • PARI
    apply( {A051578(n)=(n+2)!<<(n-1)}, [0..18]) \\ M. F. Hasler, Nov 10 2024
    
  • Sage
    [2^(n-1)*factorial(n+2) for n in (0..20)] # G. C. Greubel, Nov 11 2019
    

Formula

a(n) = (2*n+4)!!/4!!.
E.g.f.: 1/(1-2*x)^3.
a(n) ~ 2^(-1/2)*Pi^(1/2)*n^(5/2)*2^n*e^-n*n^n*{1 + 37/12*n^-1 + ...}. - Joe Keane (jgk(AT)jgk.org), Nov 23 2001
a(n) = (n+2)!*2^(n-1). - Zerinvary Lajos, Sep 23 2006. [corrected by Gary Detlefs, Apr 29 2019]
a(n) = 2^n*A001710(n+2). - R. J. Mathar, Feb 22 2008
From Peter Bala, May 26 2017: (Start)
a(n+1) = (2*n + 6)*a(n) with a(0) = 1.
O.g.f. satisfies the Riccati differential equation 2*x^2*A(x)' = (1 - 6*x)*A(x) - 1 with A(0) = 1.
G.f. as an S-fraction: A(x) = 1/(1 - 6*x/(1 - 2*x/(1 - 8*x/(1 - 4*x/(1 - 10*x/(1 - 6*x/(1 - ... - (2*n + 4)*x/(1 - 2*n*x/(1 - ...))))))))) (by Stokes 1982).
Reciprocal as an S-fraction: 1/A(x) = 1/(1 + 6*x/(1 - 8*x/(1 - 2*x/(1 - 10*x/(1 - 4*x/(1 - 12*x/(1 - 6*x/(1 - ... - (2*n + 6)*x/(1 - 2*n*x/(1 - ...)))))))))). (End)
From Amiram Eldar, Dec 11 2022: (Start)
Sum_{n>=0} 1/a(n) = 8*sqrt(e) - 12.
Sum_{n>=0} (-1)^n/a(n) = 8/sqrt(e) - 4. (End)
a(n) = A052587(n+2) for n > 0. - M. F. Hasler, Nov 10 2024

A051580 a(n) = (2*n+6)!!/6!!, related to A000165 (even double factorials).

Original entry on oeis.org

1, 8, 80, 960, 13440, 215040, 3870720, 77414400, 1703116800, 40874803200, 1062744883200, 29756856729600, 892705701888000, 28566582460416000, 971263803654144000, 34965496931549184000, 1328688883398868992000
Offset: 0

Views

Author

Keywords

Comments

Row m=6 of the array A(3; m,n) := (2*n+m)!!/m!!, m >= 0, n >= 0.

Crossrefs

Cf. A000165, A001147(n+1), A002866(n+1), A051577, A051578, A051579 (rows m=0..5), A051581, A051582, A051583.

Programs

  • GAP
    List([0..20], n-> Product([1..n], j-> 2*j+6) ); # G. C. Greubel, Nov 11 2019
  • Magma
    [1] cat [(&*[2*j+6: j in [1..n]]): n in [1..20]]; // G. C. Greubel, Nov 11 2019
    
  • Maple
    seq( mul(2*j+6, j=1..n), n=0..20); # G. C. Greubel, Nov 11 2019
  • Mathematica
    Table[2^n*Pochhammer[4, n], {n,0,20}] (* G. C. Greubel, Nov 11 2019 *)
    Table[(2n+6)!!/6!!,{n,0,20}] (* Harvey P. Dale, Mar 03 2022 *)
  • PARI
    vector(20, n, prod(j=1,n-1, 2*j+6) ) \\ G. C. Greubel, Nov 11 2019
    
  • Sage
    [product( (2*j+6) for j in (1..n)) for n in (0..20)] # G. C. Greubel, Nov 11 2019
    

Formula

a(n) = (2*n+6)!!/6!!.
E.g.f.: 1/(1-2*x)^4.
a(n) = n!*2^(n-4)/3, n>=3. - Zerinvary Lajos, Sep 23 2006
G.f.: G(0)/2, where G(k)= 1 + 1/(1 - x/(x + 1/(2*k+8)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 02 2013
From Peter Bala, May 26 2017: (Start)
a(n+1) = (2*n + 8)*a(n) with a(0) = 1.
O.g.f. satisfies the Riccati differential equation 2*x^2*A(x)' = (1 - 8*x)*A(x) - 1 with A(0) = 1.
G.f. as an S-fraction: A(x) = 1/(1 - 8*x/(1 - 2*x/(1 - 10*x/(1 - 4*x/(1 - 12*x/(1 - 6*x/(1 - ... - (2*n + 6)*x/(1 - 2*n*x/(1 - ...))))))))) (by Stokes 1982).
Reciprocal as an S-fraction: 1/A(x) = 1/(1 + 8*x/(1 - 10*x/(1 - 2*x/(1 - 12*x/(1 - 4*x/(1 - 14*x/(1 - 6*x/(1 - ... - (2*n + 8)*x/(1 - 2*n*x/(1 - ...)))))))))). (End)
From Amiram Eldar, Dec 11 2022: (Start)
Sum_{n>=0} 1/a(n) = 48*sqrt(e) - 78.
Sum_{n>=0} (-1)^n/a(n) = 30 - 48/sqrt(e). (End)

A051579 a(n) = (2*n+5)!!/5!!, related to A001147 (odd double factorials).

Original entry on oeis.org

1, 7, 63, 693, 9009, 135135, 2297295, 43648605, 916620705, 21082276215, 527056905375, 14230536445125, 412685556908625, 12793252264167375, 422177324717523375, 14776206365113318125, 546719635509192770625
Offset: 0

Views

Author

Keywords

Comments

Row m=5 of the array A(3; m,n) := (2*n+m)!!/m!!, m >= 0, n >= 0.

Crossrefs

Cf. A000165, A001147(n+1), A002866(n+1), A051577, A051578 (rows m=0..4).

Programs

  • GAP
    List([0..20], n-> Product([0..n-1], j-> 2*j+7) ); # G. C. Greubel, Nov 12 2019
  • Magma
    [1] cat [(&*[2*j+7: j in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Nov 12 2019
    
  • Maple
    df:=doublefactorial; seq(df(2*n+5)/df(5), n = 0..20); # G. C. Greubel, Nov 12 2019
  • Mathematica
    Table[2^n*Pochhammer[7/2, n], {n,0,20}] (* G. C. Greubel, Nov 12 2019 *)
  • PARI
    vector(20, n, prod(j=1,n-1, 2*j+5) ) \\ G. C. Greubel, Nov 12 2019
    
  • Sage
    [product( (2*j+7) for j in (0..n-1)) for n in (0..20)] # G. C. Greubel, Nov 12 2019
    

Formula

a(n) = (2*n+5)!!/4!!.
E.g.f.: 1/(1-2*x)^(7/2).
a(n) ~ 8/15*sqrt(2)*n^3*2^n*e^-n*n^n*(1 + 107/24*n^-1 + ...). - Joe Keane (jgk(AT)jgk.org), Nov 23 2001
G.f.: G(0)/(10*x) -1/(5*x), where G(k)= 1 + 1/(1 - x*(2*k+5)/(x*(2*k+5) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 13 2013
From Peter Bala, May 26 2017: (Start)
a(n+1) = (2*n + 7)*a(n) with a(0) = 1.
O.g.f. satisfies the Riccati differential equation 2*x^2*A(x)' = (1 - 7*x)*A(x) - 1 with A(0) = 1.
G.f. as an S-fraction: A(x) = 1/(1 - 7*x/(1 - 2*x/(1 - 9*x/(1 - 4*x/(1 - 11*x/(1 - 6*x/(1 - ... - (2*n + 5)*x/(1 - 2*n*x/(1 - ...))))))))) (by Stokes 1982).
Reciprocal as an S-fraction: 1/A(x) = 1/(1 + 7*x/(1 - 9*x/(1 - 2*x/(1 - 11*x/(1 - 4*x/(1 - 13*x/(1 - 6*x/(1 - ... - (2*n + 7)*x/(1 - 2*n*x/(1 - ...)))))))))). (End)
From Amiram Eldar, Dec 11 2022: (Start)
Sum_{n>=0} 1/a(n) = 15 * sqrt(e*Pi/2) * erf(1/sqrt(2)) - 20, where erf is the error function.
Sum_{n>=0} (-1)^n/a(n) = 15 * sqrt(Pi/(2*e)) * erfi(1/sqrt(2)) - 10, where erfi is the imaginary error function. (End)

A051582 a(n) = (2*n+8)!!/8!!, related to A000165 (even double factorials).

Original entry on oeis.org

1, 10, 120, 1680, 26880, 483840, 9676800, 212889600, 5109350400, 132843110400, 3719607091200, 111588212736000, 3570822807552000, 121407975456768000, 4370687116443648000, 166086110424858624000, 6643444416994344960000
Offset: 0

Views

Author

Keywords

Comments

Row m=8 of the array A(3; m,n) := (2*n+m)!!/m!!, m >= 0, n >= 0.

Crossrefs

Cf. A000165, A001147(n+1), A002866(n+1).
Cf. A051577, A051578, A051579, A051580, A051581 (rows m=0..7), A051583.

Programs

  • GAP
    F:=Factorial;; List([0..20], n-> 2^n*F(n+4)/F(4) ); # G. C. Greubel, Nov 12 2019
  • Magma
    F:=Factorial; [2^n*F(n+4)/F(4): n in [0..20]]; // G. C. Greubel, Nov 12 2019
    
  • Maple
    seq(2^n*pochhammer(5, n), n=0..20); # G. C. Greubel, Nov 12 2019
  • Mathematica
    (2Range[0,20]+8)!!/8!! (* Harvey P. Dale, Feb 03 2013 *)
    Table[2^n*Pochhammer[5, n], {n,0,20}] (* G. C. Greubel, Nov 12 2019 *)
  • PARI
    vector(20, n, n--; (n+4)!*2^(n-1)/12) \\ Michel Marcus, Feb 09 2015
    
  • Sage
    f=factorial; [2^n*f(n+4)/f(4) for n in (0..20)] # G. C. Greubel, Nov 12 2019
    

Formula

a(n) = (2*n+8)!!/8!!.
E.g.f.: 1/(1-2*x)^5.
a(n) = (n+4)!*2^(n-1)/12. - Zerinvary Lajos, Sep 23 2006
From Peter Bala, May 26 2017: (Start)
a(n+1) = (2*n + 10)*a(n) with a(0) = 1.
O.g.f. satisfies the Riccati differential equation 2*x^2*A(x)' = (1 - 10*x)*A(x) - 1 with A(0) = 1.
G.f. as an S-fraction: A(x) = 1/(1 - 10*x/(1 - 2*x/(1 - 12*x/(1 - 4*x/(1 - 14*x/(1 - 6*x/(1 - ... - (2*n + 8)*x/(1 - 2*n*x/(1 - ...))))))))) (by Stokes 1982).
Reciprocal as an S-fraction: 1/A(x) = 1/(1 + 10*x/(1 - 12*x/(1 - 2*x/(1 - 14*x/(1 - 4*x/(1 - 16*x/(1 - 6*x/(1 - ... - (2*n + 10)*x/(1 - 2*n*x/(1 - ...)))))))))). (End)
From Amiram Eldar, Dec 11 2022: (Start)
Sum_{n>=0} 1/a(n) = 384*sqrt(e) - 632.
Sum_{n>=0} (-1)^n/a(n) = 384/sqrt(e) - 232. (End)

A051583 a(n) = (2*n+9)!!/9!!, related to A001147 (odd double factorials).

Original entry on oeis.org

1, 11, 143, 2145, 36465, 692835, 14549535, 334639305, 8365982625, 225881530875, 6550564395375, 203067496256625, 6701227376468625, 234542958176401875, 8678089452526869375, 338445488648547905625, 13876265034590464130625
Offset: 0

Views

Author

Keywords

Comments

Row m=9 of the array A(3; m,n) := (2*n+m)!!/m!!, m >= 0, n >= 0.

Crossrefs

Cf. A000165, A001147(n+1), A002866(n+1), A178647.
Cf. A051577, A051578, A051579, A051580, A051581, A051582 (rows m=0..8).

Programs

  • GAP
    List([0..20], n-> Product([0..n-1], j-> 2*j+11) ); # G. C. Greubel, Nov 12 2019
  • Magma
    [1] cat [(&*[2*j+11: j in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Nov 12 2019
    
  • Maple
    seq(2^n*pochhammer(11/2,n), n = 0..20); # G. C. Greubel, Nov 12 2019
  • Mathematica
    (2*Range[0,20]+9)!!/945 (* Harvey P. Dale, Apr 10 2019 *)
    Table[2^n*Pochhammer[11/2, n], {n,0,20}] (* G. C. Greubel, Nov 12 2019 *)
  • PARI
    vector(20, n, prod(j=0,n-2, 2*j+11) ) \\ G. C. Greubel, Nov 12 2019
    
  • Sage
    [product( (2*j+11) for j in (0..n-1)) for n in (0..20)] # G. C. Greubel, Nov 12 2019
    

Formula

a(n) = (2*n+9)!!/9!!.
E.g.f.: 1/(1-2*x)^(11/2).
From Peter Bala, May 26 2017: (Start)
a(n+1) = (2*n + 11)*a(n) with a(0) = 1.
O.g.f. satisfies the Riccati differential equation 2*x^2*A(x)' = (1 - 11*x)*A(x) - 1 with A(0) = 1.
G.f. as an S-fraction: A(x) = 1/(1 - 11*x/(1 - 2*x/(1 - 13*x/(1 - 4*x/(1 - 15*x/(1 - 6*x/(1 - ... - (2*n + 9)*x/(1 - 2*n*x/(1 - ...))))))))) (by Stokes 1982).
Reciprocal as an S-fraction: 1/A(x) = 1/(1 + 11*x/(1 - 13*x/(1 - 2*x/(1 - 15*x/(1 - 4*x/(1 - 17*x/(1 - 6*x/(1 - ... - (2*n + 11)*x/(1 - 2*n*x/(1 - ...)))))))))). (End)
From Amiram Eldar, Dec 11 2022: (Start)
Sum_{n>=0} 1/a(n) = 945 * sqrt(e*Pi/2) * erf(1/sqrt(2)) - 1332, where erf is the error function.
Sum_{n>=0} (-1)^n/a(n) = 945 * sqrt(Pi/(2*e)) * erfi(1/sqrt(2)) - 684, where erfi is the imaginary error function. (End)
Showing 1-6 of 6 results.