I'm trying to make sense of Hankel's expansion for the first Bessel function $J_0$. According to NIST, it is 10.17(i) here. Using their notation, we have $\nu=0.$ Hence the linked formula becomes:
Let
$$a_k = \frac{(-1^k)1^2*3^2* . . . *(2k-1)^2}{8^kk!}$$
and
$$\omega = z - \pi/4.$$
Then
$$J_0(z) \approx \sqrt{\frac{2}{\pi z}}\left(\cos(\omega)\sum_{k=0}^{\infty}{(-1)^k \frac{a_{2k}}{z^{2k}} - \sin(\omega)\sum_{k=0}^{\infty}{(-1)^k \frac{a_{2k+1}}{z^{2k+1}}}}\right)$$
But that makes no sense, as the series on the right both diverge! The problem is that the absolute value of the numerator of $a_k$ is greater than $(k!)^2$. Therefore, the absolute value of $a_k$ is greater than $\frac{k!}{8^k}$; hence for large $k$, it dominates the $z^{(2k+1)}$ in the denominator.
Can anyone rectify the apparent contradiction? (NIST goofed, maybe?)