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!)
A089915 Special values of generalized Laguerre polynomials L_n^(alpha)(x): a(n) = 4^n *n! *L_n^(-1/4)(-1). 1
1, 7, 93, 1747, 41881, 1214079, 41146933, 1592909707, 69246230193, 3336315914359, 176318707191181, 10135074699232707, 629237102421392713, 41946693027934462447, 2987346130479191856741, 226298375000985686615419, 18164688323228843027295073 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
L_n^(alpha)(x) is orthogonal over the interval (0,infinity) with weight exp(-x)*x^alpha. L_1^(alpha)(x) = -x+1+alpha.
LINKS
FORMULA
From Vaclav Kotesovec, Jul 31 2013: (Start)
a(n) = (8*n-1)*a(n-1) - 4*(n-1)*(4*n-5)*a(n-2).
a(n) ~ 2^(2*n-1/2)*n^(n+1/8)*exp(2*sqrt(n)-n-1/2) * (1+97/(192*sqrt(n))).
(End)
MATHEMATICA
Table[4^n*n!*LaguerreL[n, -1/4, -1], {n, 0, 20}] (* Vaclav Kotesovec, Jul 31 2013 *)
RecurrenceTable[{a[n] == (8*n-1)*a[n-1] - 4*(n-1)*(4*n-5)*a[n-2], a[0] == 1, a[1] == 7}, a, {n, 0, 50}] (* G. C. Greubel, May 14 2018 *)
PROG
(Magma) I:=[7, 93]; [1] cat [n le 2 select I[n] else (8*n-1)*Self(n-1) - 4*(n-1)*(4*n-5)*Self(n-2): n in [1..30]]; // G. C. Greubel, May 14 2018
(PARI) m=30; v=concat([7, 93], vector(m-2)); for(n=3, m, v[n]=(8*n-1)*v[n-1]-4*(n-1)*(4*n-5)*v[n-2]); concat([1], v) \\ G. C. Greubel, May 14 2018
CROSSREFS
Sequence in context: A292227 A006178 A029808 * A370939 A318697 A337457
KEYWORD
nonn
AUTHOR
Karol A. Penson, Nov 14 2003
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 March 29 08:13 EDT 2024. Contains 371265 sequences. (Running on oeis4.)