 Solution

For $y^{\prime \prime }+P(x)y^{\prime }+Q(x)y=0$, a singular point $x_0$ is a regular singular point if $(x-x_0)P(x)$ and $(x-x_0{)}^{2}Q(x)$ are analytic there. The Bessel differential equation has $xP=1$ and $x^{2}Q=x^2-{\nu }^2$, so zero is regular singular.

Use the Frobenius method with $y=x^s{\sum }_{n\ge 0}{a}_n{x}^n$, $a_0\ne 0$. Matching coefficients gives $(s^2-{\nu }^2)a_0=0$ and $[(s+n{)}^2-{\nu }^2]a_n+a_{n-2}=0$, with missing negative-index coefficients zero. The regular root is $s=\nu$. Then $a_1=0$, all odd coefficients vanish, and the Frobenius series for Bessel functions has recurrenceFor noninteger $\nu$, the root $s=-\nu$ supplies an independent even-power series $y_2=x^{-\nu }{\sum }_{m\ge 0}{b}_{2m}{x}^{2m}$ with $b_{2m}=-b_{2m-2}/[4m(m-\nu )]$ and $b_0\ne 0$. None of these denominators vanishes. With standard normalization this is $J_{-\nu }(x)$, independent of the Bessel function of the first kind $J_{\nu }(x)$. The real powers here are understood for $x>0$.

For integer $\nu \ge 0$ and the stated normalization $a_0=1/(2^{\nu }\nu !)$, iteration givesThe ratio of consecutive term magnitudes is $|x{|}^2/[4(m+1)(m+\nu +1)]\to 0$. Thus the series converges absolutely for every finite complex $x$, defining an entire function for integer $\nu$. At integer order, $J_{-\nu }=(-1{)}^{\nu }{J}_{\nu }$; an independent solution is instead the Bessel function of the second kind.