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!)
A352313 a(n) is -2^(2n) times the derivative of order 2n of the logarithm of J_0(x) (the Bessel function of the first kind of order zero) evaluated at zero. 1
0, 2, 6, 80, 2310, 114912, 8741040, 943113600, 136984998150, 25771171025600, 6096112602836256, 1770904261963549440, 619781244413394663600, 257206839532730597664000, 124885627021632853758096000, 70139282938382254349182510080, 45116676897361263845107791884550 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
The sequence results from some interesting combinatorics that come from the interaction of
1) the quotient rule,
2) the derivative recurrence relation for Bessel functions, and
3) the fact that J_n(0) = 0 unless n = 0 and J_0(0) = 1.
At each stage the numerator is the sum of products of Bessel functions. Each derivative introduces new terms which also (by the Bessel derivative recurrence relation) increment or decrement the Bessel function orders in the terms of the numerator. a(n) is the coefficient in front of the term (J_0(x))^(2^n) (if it exists). There are complicated combinatorics involved in what contributes to that coefficient.
LINKS
Eric Weisstein's World of Mathematics, Bessel Function of the First Kind
FORMULA
a(n) = -2^(2n)*(d^(2n)/dx^(2n))(log(J_0(x)))|_x=0.
An expression which does not require derivatives is:
a(n) = 2^(2n)*Sum_{k=1..2n} (-1)^k*(k-1)!*Y_{2n,k}(d_i),
where Y_{n,k} is the partial Bell polynomial with inputs given by
d_{2i} = (-1)^i*2^(-2i)*binomial(2i,i) and d_{2i+1} = 0.
EXAMPLE
For n = 0:
-1*(d^0/dx^0)log(J_0(x)) = -log(J_0(0)) = 0.
For n = 1:
-4*(d^2/dx^2)log(J_0(x)) = (-J_{-2}J_0 + 2J_0^2 - J_2J_0 + J_{-1}J_{-1} -
2J_{-1}J_1 + J_1 J_1)/(J_0)^4, which evaluates to 2 (because of 2J_0^2).
MAPLE
a := n -> -2^(2*n)*limit(diff(log(BesselJ(0, x)), [x$2*n]), x=0):
seq(a(n), n=0..5);
MATHEMATICA
Table[2^n*Sum[(-1)^k*Factorial[k - 1]*BellY[n, k, Flatten[Table[{0, (-1)^{(j + 1)/2}*2^{-(j + 1)}*Binomial[j + 1, (j + 1)/2]}, {j, 1, N - k + 1, 2}]]], {k, 1, n}], {n, 0, 32, 2}]
CROSSREFS
Cf. A352284 (signed version).
Sequence in context: A244084 A362581 A352284 * A357028 A245463 A325949
KEYWORD
nonn,easy
AUTHOR
Geoff Goehle, Mar 11 2022
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 1 10:18 EDT 2024. Contains 374816 sequences. (Running on oeis4.)