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!)
A152600 a(n)=0^n+sum{k=0..n-1, C(n+k-1,2k)*A000108(k)*3^k*2^(n-k)} 2
1, 2, 10, 80, 790, 8720, 103060, 1275680, 16326190, 214280720, 2868504460, 39014154080, 537592643740, 7488960021920, 105295566289960, 1492291482505280, 21296015905884190, 305755507155234320 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

Hankel transform is 2^n*3^C(n+1,2)*5^C(n,2). A152601(n)=a(n+1)/2.

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 0..200

FORMULA

a(n) = 2^n * (4*(n+1)*LegendreP(n+1,4) - (31*n+16)*LegendreP(n,4))/(3*n*(n-1)) for n>1. - Mark van Hoeij, May 27 2010

Recurrence: n*a(n) = 8*(2*n-3)*a(n-1) - 4*(n-3)*a(n-2). - Vaclav Kotesovec, Oct 20 2012

a(n) ~ sqrt(8*sqrt(15)-30)*(8+2*sqrt(15))^n/(6*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Oct 20 2012

MATHEMATICA

Flatten[{1, 2, Table[2^n*(4*(n+1)*LegendreP[n+1, 4]-(31*n+16)*LegendreP[n, 4])/(3*n*(n-1)), {n, 2, 20}]}] (* Vaclav Kotesovec, Oct 20 2012 *)

PROG

(PARI) a(n)=if(n>1, (4*(n+1)*pollegendre(n+1, 4) - (31*n+16)*pollegendre(n, 4))/(3*n*(n-1))<<n, n+1) \\ Charles R Greathouse IV, Mar 19 2017

CROSSREFS

Sequence in context: A152168 A003578 A274276 * A220112 A048286 A227463

Adjacent sequences:  A152597 A152598 A152599 * A152601 A152602 A152603

KEYWORD

easy,nonn

AUTHOR

Paul Barry, Dec 09 2008

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 August 9 01:18 EDT 2022. Contains 356016 sequences. (Running on oeis4.)