|
|
A065919
|
|
Bessel polynomial y_n(4).
|
|
6
|
|
|
1, 5, 61, 1225, 34361, 1238221, 54516085, 2836074641, 170218994545, 11577727703701, 880077524475821, 73938089783672665, 6803184337622361001, 680392371852019772765, 73489179344355757819621, 8525425196317119926848801, 1057226213522667226687070945
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
|
|
REFERENCES
|
J. Riordan, Combinatorial Identities, Wiley, 1968, p. 77.
|
|
LINKS
|
|
|
FORMULA
|
y_n(x) = Sum_{k=0..n} (n+k)!*(x/2)^k/((n-k)!*k!).
Recurrence relation: a(0) = 1, a(1) = 5, a(n) = 4*(2*n-1)*a(n-1) + a(n-2) for n >= 2. Sequence A143412(n) satisfies the same recurrence relation.
1/sqrt(e) = 1 - 2*Sum_{n = 0..inf} (-1)^n/(a(n)*a(n+1)) = 1 - 2*( 1/(1*5) - 1/(5*61) + 1/(61*1225) - ... ). (End)
G.f.: 1/Q(0), where Q(k)= 1 - x - 4*x*(k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 17 2013
a(n) = exp(1/4)/sqrt(2*Pi)*BesselK(n+1/2,1/4). - Gerry Martens, Jul 22 2015
a(n) = 1/n!*Integral_{x = 0..inf} x^n*(1 + 2*x)^n dx.
E.g.f.: d/dx( exp(x*c(2*x)) ) = 1 + 5*x + 61*x^2/2! + 1225*x^3/3! + ..., where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108. (End)
G.f.: (1/(1-x))*hypergeometric2f0(1,1/2; - ; 8*x/(1-x)^2). - G. C. Greubel, Aug 16 2017
|
|
MAPLE
|
seq(simplify(2^n*KummerU(-n, -2*n, 1/2)), n=0..16); # Peter Luschny, May 10 2022
|
|
MATHEMATICA
|
Table[Sum[(n+k)!*2^k/((n-k)!*k!), {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Jul 22 2015 *)
|
|
PROG
|
(PARI) for (n=0, 100, if (n>1, a=4*(2*n - 1)*a1 + a2; a2=a1; a1=a, if (n, a=a1=5, a=a2=1)); write("b065919.txt", n, " ", a) ) \\ Harry J. Smith, Nov 04 2009
(PARI) a(n) = sum(k=0, n, (n+k)!*2^k/((n-k)!*k!) ); \\ Joerg Arndt, May 17 2013
(Magma)
A065919:= func< n | (&+[Binomial(n, k)*Factorial(n+k)*2^k/Factorial(n): k in [0..n]]) >;
(SageMath)
def A065919(n): return sum(binomial(n, k)*factorial(n+k)*2^k/factorial(n) for k in range(n+1))
|
|
CROSSREFS
|
Polynomial coefficients are in A001498.
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|