I am interested in determining a closed expression for the n-th derivative of the Bessel function of the first kind $J_0(x)$, centered in $x=0$: \begin{equation} \left.\left(\frac{\mathrm d}{\mathrm d x} \right)^n J_0(x)\right|_{x=0} \end{equation} Can I compute it? If yes, how?
Thanks in advance!
The defining ODE of the zeroth order Bessel function is
$$x^2J_0''(x) + xJ_0'(x) + x^2J_0(x) = 0.$$
Solve this using the power-series method (Frobenius method): take the ansatz $J_0(x) = \sum_{n=0}^\infty a_n x^n$ and insert it into the ODE to get a recurrence relation for the $a_n$'s and solve this. With this solution in hand note that $$\left.\frac{d^n}{dx^n}J_0(x)\right|_{x=0} = n! a_n$$ which gives you all the numbers you seek.