Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A026620
a(n) = A026615(2*n-1, n-2).
16
1, 9, 37, 147, 576, 2244, 8723, 33891, 131716, 512278, 1994202, 7770734, 30310320, 118343970, 462501135, 1809134115, 7082699580, 27750808470, 108812919270, 426966196410, 1676471166240, 6586744582080, 25894139638302, 101852815940622, 400840469986376, 1578280410414204
OFFSET
2,2
LINKS
FORMULA
From G. C. Greubel, Jun 13 2024: (Start)
a(n) = (7*n^2 - 11*n + 6)*binomial(2*n, n)/(4*(n+1)*(2*n-1)) - [n=2].
G.f.: ( (2 - 7*x + 3*x^2) - (2 - 3*x + x^2 + 2*x^3)*sqrt(1-4*x) )/(2*x*sqrt(1-4*x)).
E.g.f.: (1/2)*exp(2*x)*( 3*(-1 + x)*BesselI(0, 2*x) + (4 - 3*x)*BesselI(1, 2*x) ) + (1/2)*(3 - x - x^2). (End)
MATHEMATICA
Table[(7*n^2-11*n+6)*Binomial[2*n, n]/(4*(n+1)*(2*n-1))-Boole[n==2], {n, 2, 40}] (* G. C. Greubel, Jun 13 2024 *)
PROG
(Magma) [n eq 2 select 1 else (7*n^2-11*n+6)*Catalan(n)/(4*(2*n-1)): n in [2..40]]; // G. C. Greubel, Jun 13 2024
(SageMath) [(7*n^2-11*n+6)*binomial(2*n, n)/(4*(n+1)*(2*n-1)) -int(n==2) for n in range(2, 41)] # G. C. Greubel, Jun 13 2024
KEYWORD
nonn
STATUS
approved