The linear second order ordinary differential equation of type

\[{{x^2}y^{\prime\prime} + xy’ }+{ \left( {{x^2} – {v^2}} \right)y }={ 0}\]

is called the Bessel equation. The number \(v\) is called the order of the Bessel equation.

The given differential equation is named after the German mathematician and astronomer Friedrich Wilhelm Bessel who studied this equation in detail and showed (in \(1824\)) that its solutions are expressed in terms of a special class of functions called cylinder functions or Bessel functions.

Concrete representation of the general solution depends on the number \(v.\) Further we consider separately two cases:

- The order \(v\) is non-integer;
- The order \(v\) is an integer.

### Case \(1.\) The Order \(v\) is Non-Integer

Assuming that the number \(v\) is non-integer and positive, the general solution of the Bessel equation can be written as

\[{y\left( x \right) }={ {C_1}{J_v}\left( x \right) + {C_2}{J_{ – v}}\left( x \right),}\]

where \({C_1},\) \({C_2}\) are arbitrary constants and \({J_v}\left( x \right),\) \({J_{ – v}}\left( x \right)\) are Bessel functions of the first kind.

The Bessel function can be represented by a series, the terms of which are expressed using the so-called Gamma function:

\[{{J_v}\left( x \right) \text{=}}\kern0pt{ \sum\limits_{p = 0}^\infty {\frac{{{{\left( { – 1} \right)}^p}}}{{\Gamma \left( {p + 1} \right)\Gamma \left( {p + v + 1} \right)}} \cdot}\kern0pt{ {{\left( {\frac{x}{2}} \right)}^{2p + v}}} }\]

The Gamma function is the generalization of the factorial function from integers to all real numbers. It has, in particular, the following properties:

\[

{\Gamma \left( {p + 1} \right) = p!,\;\;}\kern-0.3pt

{\Gamma \left( {p + v + 1} \right) }={ \left( {v + 1} \right)\left( {v + 2} \right) \cdots }\kern0pt{ \left( {v + p} \right)\Gamma \left( {v + 1} \right).}

\]

The Bessel functions of the negative order \(\left(-v\right)\) (assuming that \(v \gt 0\)) are written in similar way:

\[{{J_{ – v}}\left( x \right) \text{=}}\kern0pt{ \sum\limits_{p = 0}^\infty {\frac{{{{\left( { – 1} \right)}^p}}}{{\Gamma \left( {p + 1} \right)\Gamma \left( {p – v + 1} \right)}} \cdot}\kern0pt{ {{\left( {\frac{x}{2}} \right)}^{2p – v}}} }\]

The Bessel functions can be calculated in most mathematical software packages as well as in MS Excel. For example, the Bessel functions of the \(1\)st kind of orders \(v = 0\) to \(v = 4\) are shown in Figure \(1.\)

### Case \(2.\) The Order \(v\) is an Integer

If the order \(v\) of the Bessel differential equation is an integer, the Bessel functions \({J_v}\left( x \right)\) and \({J_{ – v}}\left( x \right)\) can become dependent from each other. In this case the general solution is described by another formula:

\[{y\left( x \right) }={ {C_1}{J_v}\left( x \right) + {C_2}{Y_v}\left( x \right),}\]

where \({Y_v}\left( x \right)\) is the Bessel function of the second kind. Sometimes this family of functions is also called Neumann functions or Weber functions.

The Bessel function of the second kind \({Y_v}\left( x \right)\) can be expressed in terms of the Bessel functions of the first kind \({J_v}\left( x \right)\) and \({J_{ – v}}\left( x \right):\)

\[{{Y_v}\left( x \right) }={ \frac{{{J_v}\left( x \right)\cos \pi v – {J_{ – v}}\left( x \right)}}{{\sin \pi v}}.}\]

The graphs of the functions \({Y_v}\left( x \right)\) for several first orders \(v\) are shown in Figure \(2.\)

#### Note:

Actually the general solution of the differential equation expressed in terms of Bessel functions of the first and second kind is valid for non-integer orders as well.

### Some Differential Equations Reducible to Bessel’s Equation

- One of the well-known equations tied with the Bessel’s differential equation is the modified Bessel’s equation that is obtained by replacing \(x\) with \(ix.\) This equation has the form:
\[{{x^2}y^{\prime\prime} + xy’ }-{ \left( {{x^2} + {v^2}} \right)y }={ 0.}\]The general solution of this equation can be expressed in terms of the so-called modified Bessel functions of the first and second kind:\[ {y\left( x \right) }={ {C_1}{J_v}\left( { – ix} \right) + {C_2}{Y_v}\left( { – ix} \right) } = {{C_3}{I_v}\left( x \right) + {C_4}{K_v}\left( x \right),} \]where \({I_v}\left( x \right)\) and \({K_v}\left( x \right)\) are modified Bessel functions of the \(1\)st and \(2\)nd kind, respectively, and \({C_1},\) \({C_2},\) \({C_3},\) \({C_4}\) are arbitrary constants.
- The Airy differential equation known in astronomy and physics has the form:
\[y^{\prime\prime} – xy = 0.\]It can also be reduced to the Bessel equation. Its solution is given by the Bessel functions of the fractional order \( \pm {\large\frac{1}{3}\normalsize}:\)\[ {y\left( x \right) } = {{C_1}\sqrt x {J_{\large\frac{1}{3}\normalsize}}\left( {\frac{2}{3}i{x^{\large\frac{3}{2}\normalsize}}} \right) } + {{C_2}\sqrt x {J_{ – \large\frac{1}{3}\normalsize}}\left( {\frac{2}{3}i{x^{\large\frac{3}{2}\normalsize}}} \right).} \]
- The differential equation of type
\[{{x^2}y^{\prime\prime} + xy’ }+{ \left( {{a^2}{x^2} – {v^2}} \right)y }={ 0}\]differs from the Bessel equation only by a factor \({a^2}\) before \({x^2}\) and has the general solution in the form:\[{y\left( x \right) }={ {C_1}{J_v}\left( {ax} \right) + {C_2}{Y_v}\left( {ax} \right).}\]
- The similar differential equation
\[{{x^2}y^{\prime\prime} + axy’ }+{ \left( {{x^2} – {v^2}} \right)y }={ 0}\]is reduced to the Bessel equation\[{{x^2}z^{\prime\prime} + xz’ }+{ \left( {{x^2} – {n^2}} \right)z }={ 0}\]by using the substitution\[y\left( x \right) = {x^{\large\frac{{1 – a}}{2}\normalsize}}z\left( x \right).\]Here the parameter \({n^2}\) denotes:\[{n^2} = {v^2} + \frac{1}{4}{\left( {a – 1} \right)^2}.\]As a result, the general solution of the differential equation is given by\[{y\left( x \right) }={ {x^{\large\frac{{1 – a}}{2}\normalsize}}\left[ {{C_1}{J_n}\left( x \right) }\right.}+{\left.{ {C_2}{Y_n}\left( x \right)} \right].}\]

The special Bessel functions are widely used in solving problems of theoretical physics, for example in investigating

- wave propagation;
- heat conduction;
- vibrations of membranes

in the systems with cylindrical or spherical symmetry.

## Solved Problems

Click a problem to see the solution.