login
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
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
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 * A371460 A220112 A367851
KEYWORD
easy,nonn
AUTHOR
Paul Barry, Dec 09 2008
STATUS
approved