|
Search: id:A144301
|
|
|
| A144301 |
|
a(0) = a(1) = 1; thereafter a(n) = (2*n-3)*a(n-1) + a(n-2). |
|
+0 6
|
|
| 1, 1, 2, 7, 37, 266, 2431, 27007, 353522, 5329837, 90960751, 1733584106, 36496226977, 841146804577, 21065166341402, 569600638022431, 16539483668991901, 513293594376771362, 16955228098102446847
(list; graph; listen)
|
|
|
OFFSET
|
0,3
|
|
|
COMMENT
|
A variant of A001515, which is the main entry.
|
|
REFERENCES
|
E. Grosswald, Bessel Polynomials, Lecture Notes Math., Vol. 698, 1978.
|
|
FORMULA
|
a(n) = A001515(n-1) for n>= 1.
E.g.f.: A(x) = exp(1-sqrt(1-2*x)).
A'(x) = A(x)/(1-sqrt(1-2*x)).
Hence a(n+1) = Sum_{ k = 0..n } a(n-k) binomial(n,k)*(2*k)!/(k!*2^k).
A''(x) = (A'(x)/(1-2*x))*(1 + 1/sqrt(1-2*x)).
A''(x) = 2*x*A''(x) + A'(x) + A(x), which is equivalent to the recurrence in the definition.
a(n) = Sum_{ k = 0..n-1 } binomial(n-1-k,2*k)*(2*k)!/(k!*2^k). [See Grosswald, p. 6, Eq. (8).]
a(n) ~ exp(1)*(2n-1)!/(n!*2^n) as n -> oo. [See Grosswald, p. 124]
|
|
CROSSREFS
|
See A001515 for much more about this sequence.
See A144498 for first differences.
Sequence in context: A125515 A135920 A001515 this_sequence A083659 A036247 A107877
Adjacent sequences: A144298 A144299 A144300 this_sequence A144302 A144303 A144304
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
David Applegate and N. J. A. Sloane (njas(AT)research.att.com), Dec 07 2008
|
|
|
Search completed in 0.002 seconds
|