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

A162443 Numerators of the BG1[ -5,n] coefficients of the BG1 matrix.

Original entry on oeis.org

5, 66, 680, 2576, 33408, 14080, 545792, 481280, 29523968, 73465856, 27525120, 856162304, 1153433600, 18798870528, 86603988992, 2080374784, 2385854332928, 3216930504704, 71829033058304, 7593502179328, 281749854617600
Offset: 1

Views

Author

Johannes W. Meijer, Jul 06 2009

Keywords

Comments

The BG1 matrix coefficients are defined by BG1[2m-1,1] = 2*beta(2m) and the recurrence relation BG1[2m-1,n] = BG1[2m-1,n-1] - BG1[2m-3,n-1]/(2*n-3)^2 with m = .. , -2, -1, 0, 1, 2, .. and n = 1, 2, 3, .. . As usual beta(m) = sum((-1)^k/(1+2*k)^m, k=0..infinity). For the BG2 matrix, the even counterpart of the BG1 matrix, see A008956.
We discovered that the n-th term of the row coefficients can be generated with BG1[1-2*m,n] = RBS1(1-2*m,n)* 4^(n-1)*((n-1)!)^2/ (2*n-2)! for m >= 1. For the BS1(1-2*m,n) polynomials see A160485.
The coefficients in the columns of the BG1 matrix, for m >= 1 and n >= 2, can be generated with GFB(z;n) = ((-1)^(n+1)*CFN2(z;n)*GFB(z;n=1) + BETA(z;n))/((2*n-3)!!)^2 for n >= 2. For the CFN2(z;n) and the Beta polynomials see A160480.
The BG1[ -5,n] sequence can be generated with the first Maple program and the BG1[2*m-1,n] matrix coefficients can be generated with the second Maple program.
The BG1 matrix is related to the BS1 matrix, see A160480 and the formulas below.

Examples

			The first few formulas for the BG1[1-2*m,n] matrix coefficients are:
BG1[ -1,n] = (1)*4^(n-1)*(n-1)!^2/(2*n-2)!
BG1[ -3,n] = (1-2*n)*4^(n-1)*(n-1)!^2/(2*n-2)!
BG1[ -5,n] = (1-8*n+12*n^2)*4^(n-1)*(n-1)!^2/(2*n-2)!
The first few generating functions GFB(z;n) are:
GFB(z;2) = ((-1)*(z^2-1)*GFB(z;1) + (-1))/1
GFB(z;3) = ((+1)*(z^4-10*z^2+9)*GFB(z;1) + (-11 + z^2))/9
GFB(z;4) = ((-1)*( z^6- 35*z^4+259*z^2-225)*GFB(z;1) + (-299 + 36*z^2 - z^4))/225
		

Crossrefs

A162444 are the denominators of the BG1[ -5, n] matrix coefficients.
The BG1[ -3, n] equal (-1)*A002595(n-1)/A055786(n-1) for n >= 1.
The BG1[ -1, n] equal A046161(n-1)/A001790(n-1) for n >= 1.
The cs(n) equal A046161(n-2)/A001803(n-2) for n >= 2.
The BETA(z, n) polynomials and the BS1 matrix lead to the Beta triangle A160480.
The CFN2(z, n), the t2(n, m) and the BG2 matrix lead to A008956.
Cf. A162443 (BG1 matrix), A162446 (ZG1 matrix) and A162448 (LG1 matrix).

Programs

  • Maple
    a := proc(n): numer((1-8*n+12*n^2)*4^(n-1)*(n-1)!^2/(2*n-2)!) end proc: seq(a(n), n=1..21);
    # End program 1
    nmax1 := 5; coln := 3; Digits := 20: mmax1 := nmax1: for n from 0 to nmax1 do t2(n, 0) := 1 od: for n from 0 to nmax1 do t2(n, n) := doublefactorial(2*n-1)^2 od: for n from 1 to nmax1 do for m from 1 to n-1 do t2(n, m) := (2*n-1)^2* t2(n-1, m-1) + t2(n-1, m) od: od: for m from 1 to mmax1 do BG1[1-2*m, 1] := euler(2*m-2) od: for m from 1 to mmax1 do BG1[2*m-1, 1] := Re(evalf(2*sum((-1)^k1/(1+2*k1)^(2*m), k1=0..infinity))) od: for m from -mmax1 +coln to mmax1 do BG1[2*m-1, coln] := (-1)^(coln+1)*sum((-1)^k1*t2(coln-1, k1)*BG1[2*m-(2*coln-1)+2*k1, 1], k1=0..coln-1)/doublefactorial(2*coln-3)^2 od;
    # End program 2
    # Maple programs edited by Johannes W. Meijer, Sep 25 2012

Formula

a(n) = numer(BG1[ -5,n]) and A162444(n) = denom(BG1[ -5,n]) with BG1[ -5,n] = (1-8*n+12*n^2)*4^(n-1)*(n-1)!^2/(2*n-2)!.
The generating functions GFB(z;n) of the coefficients in the matrix columns are defined by
GFB(z;n) = sum(BG1[2*m-1,n]*z^(2*m-2), m=1..infinity).
GFB(z;n) = (1-z^2/(2*n-3)^2)*GFB(n-1) - 4^(n-2)*(n-2)!^2/((2*n-4)!*(2*n-3)^2) for n => 2 with GFB(z;n=1) = 1/(z*cos(Pi*z/2))*int(sin(z*t)/sin(t),t=0..Pi/2).
The column sums cs(n) = sum(BG1[2*m-1,n]*z^(2*m-2), m=1..infinity) = 4^(n-1)/((2*n-2)*binomial(2*n-2,n-1)) for n >= 2.
BG1[2*m-1,n] = (n-1)!^2*4^(n-1)*BS1[2*m-1,n]/(2*n-2)!

A055786 Numerators of Taylor series expansion of arcsin(x). Also arises from arccos(x), arccsc(x), arcsec(x), arcsinh(x).

Original entry on oeis.org

1, 1, 3, 5, 35, 63, 231, 143, 6435, 12155, 46189, 88179, 676039, 1300075, 5014575, 9694845, 100180065, 116680311, 2268783825, 1472719325, 34461632205, 67282234305, 17534158031, 514589420475, 8061900920775, 5267108601573
Offset: 0

Views

Author

N. J. A. Sloane, Jul 13 2000

Keywords

Comments

Note that the sequence is not monotonic.

Examples

			arcsin(x) is usually written as x + x^3/(2*3) + 1*3*x^5/(2*4*5) + 1*3*5*x^7/(2*4*6*7) + ..., which is x + 1/6*x^3 + 3/40*x^5 + 5/112*x^7 + 35/1152*x^9 + 63/2816*x^11 + ... (A055786/A002595) when reduced to lowest terms.
arccos(x) = Pi/2 - (x + (1/6)*x^3 + (3/40)*x^5 + (5/112)*x^7 + (35/1152)*x^9 + (63/2816)*x^11 + ...) (A055786/A002595).
arccsc(x) = 1/x + 1/(6*x^3) + 3/(40*x^5) + 5/(112*x^7) + 35/(1152*x^9) + 63/(2816*x^11) + ... (A055786/A002595).
arcsec(x) = Pi/2 -(1/x + 1/(6*x^3) + 3/(40*x^5) + 5/(112*x^7) + 35/(1152*x^9) + 63/(2816*x^11) + ...) (A055786/A002595).
arcsinh(x) = x - (1/6)*x^3 + (3/40)*x^5 - (5/112)*x^7 + (35/1152)*x^9 - (63/2816)*x^11 + ... (A055786/A002595).
i*Pi/2 - arccosh(x) = i*x + (1/6)*i*x^3 + (3/40)*i*x^5 + (5/112)*i*x^7 + (35/1152)*i*x^9 + (63/2816)*i*x^11 + (231/13312)*i*x^13 + (143/10240)*i*x^15 + (6435/557056)*i*x^17 + ... (A055786/A002595).
0, 1, 0, 1/6, 0, 3/40, 0, 5/112, 0, 35/1152, 0, 63/2816, 0, 231/13312, 0, 143/10240, 0, 6435/557056, 0, 12155/1245184, 0, 46189/5505024, 0, ... = A055786/A002595.
a(4) = 35 = 3*5*7*9 / gcd( 3*5*7*9, (2*4*6*8) * (2*4+1))
		

References

  • Bronstein-Semendjajew, Taschenbuch der Mathematik, 7th German ed. 1965, ch. 4.2.6
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 88.
  • H. B. Dwight, Tables of Integrals and Other Mathematical Data, Macmillan, NY, 1968, Chap. 3.

Crossrefs

Cf. A002595.
a(n) / A002595(n) = A001147(n) / ( A000165(n) * (2*n+1))
Cf. A162443 where BG1[-3,n] = (-1)*A002595(n-1)/A055786(n-1) for n >= 1. - Johannes W. Meijer, Jul 06 2009

Programs

  • Magma
    [Numerator( (n+1)*Binomial(2*n+2,n+1)/(2^(2*n+1)*(2*n+1)^2) ): n in [0..25]]; // G. C. Greubel, Jan 25 2020
    
  • Maple
    seq( numer( (n+1)*binomial(2*n+2,n+1)/(2^(2*n+1)*(2*n+1)^2) ), n=0..25); # G. C. Greubel, Jan 25 2020
  • Mathematica
    Numerator/@Select[CoefficientList[Series[ArcSin[x],{x,0,60}],x], #!=0&]  (* Harvey P. Dale, Apr 29 2011 *)
  • PARI
    vector(25, n, numerator(2*n*binomial(2*n,n)/(4^n*(2*n-1)^2)) ) \\ G. C. Greubel, Jan 25 2020
    
  • Sage
    [numerator( (n+1)*binomial(2*n+2,n+1)/(2^(2*n+1)*(2*n+1)^2) ) for n in (0..25)] # G. C. Greubel, Jan 25 2020

Formula

a(n) / A052469(n) = A001147(n) / ( A000165(n) *2*n ). E.g., a(6) = 77 = 1*3*5*7*9*11 / gcd( 1*3*5*7*9*11, 2*4*6*8*10*12*12 ).
a(n) = numerator((2*n)!/(2^(2*n)*(n)!^2*(2*n+1))). - Johannes W. Meijer, Jul 06 2009

Extensions

Edited by Johannes W. Meijer, Jul 06 2009

A052469 Denominators in the Taylor series for arccosh(x) - log(2*x).

Original entry on oeis.org

4, 32, 96, 1024, 2560, 4096, 28672, 524288, 1179648, 5242880, 11534336, 100663296, 218103808, 939524096, 134217728, 68719476736, 146028888064, 206158430208, 1305670057984, 2199023255552, 7696581394432, 96757023244288, 202310139510784, 1125899906842624
Offset: 1

Views

Author

Keywords

Examples

			arccosh(x) = log(2x) - 1/(4*x^2) - 3/(32*x^4) - 5/(96*x^6) - ... for x>1.
		

References

  • Bronstein-Semendjajew, sprawotchnik po matematikje, 6th Russian ed. 1956, ch. 4.2.6.

Crossrefs

Cf. A002595.

Programs

  • GAP
    List([1..30], n-> DenominatorRat( Factorial(2*n-1)/(4^n*(Factorial(n))^2) )) # G. C. Greubel, May 18 2019
  • Magma
    [Denominator(Factorial(2*n-1)/( 2^(2*n)* Factorial(n)^2)): n in [1..30]]; // Vincenzo Librandi, Jul 10 2017
    
  • Mathematica
    a[n_] := Denominator[(2*n-1)!/(2^(2*n)*n!^2)]; Array[a, 21] (* Jean-François Alcover, May 17 2017 *)
  • PARI
    {a(n) = denominator((2*n-1)!/(4^n*(n!)^2))}; \\ G. C. Greubel, May 18 2019
    
  • Sage
    [denominator(factorial(2*n-1)/(4^n*(factorial(n))^2)) for n in (1..30)] # G. C. Greubel, May 18 2019
    

Formula

A052468(n) / a(n) = A001147(n) / ( A000165(n) *2*n )
From Johannes W. Meijer, Jul 06 2009: (Start)
a(n) = denom((2*n-1)!/( 4^n * (n!)^2)).
Equals 2*A162442(n+1) for n >= 1.
A052468(n)/a(n) = (1/(2*n))*A001790(n)/A046161(n) for n>=1.
(End)

Extensions

Updated by Frank Ellermann, May 22 2001

A162444 Denominators of the BG1[ -5,n] coefficients of the BG1 matrix.

Original entry on oeis.org

1, 1, 3, 5, 35, 9, 231, 143, 6435, 12155, 3553, 88179, 96577, 1300075, 5014575, 102051, 100180065, 116680311, 2268783825, 210388475, 6892326441, 67282234305, 17534158031, 39583801575, 8061900920775, 169906729083
Offset: 1

Views

Author

Johannes W. Meijer, Jul 06 2009

Keywords

Comments

For the numerators of the BG1[ -5,n] coefficients see A162443.
We observe that BG1[ -3,n] = (-1)*A002595(n-1)/A055786(n-1), i.e. they equal the inverted coefficients of the series expansion of arcsin(x), and that BG1[ -1,n] = A046161(n-1)/A001790(n-1), i.e. they equal the inverted coefficients of the series expansion of 1/sqrt(1-x).

Examples

			The first few formulas for the BG1[1-2*m,n] matrix coefficients are:
BG1[ -1,n] = (1)*4^(n-1)*(n-1)!^2/(2*n-2)!
BG1[ -3,n] = (1-2*n)*4^(n-1)*(n-1)!^2/(2*n-2)!
BG1[ -5,n] = (1-8*n+12*n^2)*4^(n-1)*(n-1)!^2/(2*n-2)!
BG1[ -7,n] = (1-2*n+60*n^2-120*n^3)*4^(n-1)*(n-1)!^2/(2*n-2)!
		

Crossrefs

A162443 are the numerators of the BG1[ -5, n] matrix coefficients.
The BG1[ -3, n] equal A002595(n-1)/A055786(n-1) for n =>1.
The BG1[ -1, n] equal A046161(n-1)/A001790(n-1) for n =>1.

Formula

a(n) = denom(BG1[ -5,n]) and A162443(n) = numer(BG1[ -5,n]) with BG1[ -5,n] = 4^(n-1)*(1-8*n+12*n^2)*(n-1)!^2/ (2*n-2)!.

A091154 Numerator of Maclaurin expansion of (t*sqrt(t^2+1) + arcsinh(t))/2, the arc length of Archimedes' spiral.

Original entry on oeis.org

1, 1, -1, 1, -5, 7, -21, 11, -429, 715, -2431, 4199, -29393, 52003, -185725, 334305, -3231615, 3535767, -64822395, 39803225, -883631595, 1641030105, -407771117, 11435320455, -171529806825, 107492012277, -1215486600363, 2295919134019
Offset: 1

Views

Author

Eric W. Weisstein, Dec 22 2003

Keywords

Comments

From Mikhail Gaichenkov, Feb 05 2013: (Start)
For Archimedean spiral (r=at) and the arc length s(t)= a(t*sqrt(t^2+1) + arcsinh(t))/2, the limit of s’’(t)=a, t- -> infinity. In other words, a point moves with uniform acceleration along the spiral while the spiral corresponds to the locations over time of a point moving away from a fixed point with a constant speed along a line that rotates with constant angular velocity.
The error of approximation for large t: |a-s’’(t)| ~ a/(2(1+t^2)) (Gaichenkov private research).
The arc of the Archimedean spiral is approximated by the differential equation in polar coordinates r’^2+r^2=(at)^2 (see A202407). (End)

Examples

			t + t^3/6 - t^5/40 + t^7/112 - (5*t^9)/1152 + (7*t^11)/2816 - ...
		

Crossrefs

Denominators are in A002595.
Showing 1-5 of 5 results.