2.1.3.1 second order bessel ode

\[\begin {aligned} y^{\prime \prime }+\frac {y^{\prime }}{x}+y k^{2}&=0 \end {aligned}\]

Entering second order bessel ode solverWriting the ode as

\begin{align*} x^{2} y^{\prime \prime }+y^{\prime } x +k^{2} x^{2} y = 0\tag {1} \end{align*}

Bessel ode has the form

\begin{align*} x^{2} y^{\prime \prime }+y^{\prime } x +\left (-n^{2}+x^{2}\right ) y = 0\tag {2} \end{align*}

The generalized form of Bessel ode is given by Bowman (1958) as the following

\begin{align*} x^{2} y^{\prime \prime }+\left (1-2 \alpha \right ) x y^{\prime }+\left (\beta ^{2} \gamma ^{2} x^{2 \gamma }-n^{2} \gamma ^{2}+\alpha ^{2}\right ) y = 0\tag {3} \end{align*}

With the standard solution

\begin{align*} y&=x^{\alpha } \left (c_1 \operatorname {BesselJ}\left (n , \beta \,x^{\gamma }\right )+c_2 \operatorname {BesselY}\left (n , \beta \,x^{\gamma }\right )\right )\tag {4} \end{align*}

Comparing (3) to (1) and solving for \(\alpha ,\beta ,n,\gamma \) gives

\begin{align*} \alpha &= 0\\ \beta &= k\\ n &= 0\\ \gamma &= 1 \end{align*}

Substituting all the above into (4) gives the solution as

\begin{align*} y = c_1 \operatorname {BesselJ}\left (0, k x \right )+c_2 \operatorname {BesselY}\left (0, k x \right ) \end{align*}