Explain and express general solution of Bessel's equation



Bessel's equation is a second-order linear differential equation that appears in various physical contexts, such as heat conduction, wave propagation, and electromagnetism. The general form of Bessel's equation is:

x2d2ydx2+xdydx+(x2n2)y=0x^2 \frac{d^2 y}{dx^2} + x \frac{dy}{dx} + (x^2 - n^2) y = 0

where nn is a constant, which can be any real or complex number. The solutions to this equation are known as Bessel functions.

To derive the general solution of Bessel's equation, we start with the standard form of Bessel's equation of order 

nn:

x2d2ydx2+xdydx+(x2n2)y=0x^2 \frac{d^2 y}{dx^2} + x \frac{dy}{dx} + (x^2 - n^2) y = 0

We'll use the method of Frobenius to find the solution, which involves looking for a solution in the form of a power series:

y(x)=k=0akxk+ry(x) = \sum_{k=0}^{\infty} a_k x^{k+r}

where aka_k are the coefficients to be determined, and rr is the exponent to be determined from the indicial equation.

Step 1: Substitute the Power Series into the Differential Equation

First, we need the derivatives of y(x)y(x):

dydx=k=0ak(k+r)xk+r1\frac{dy}{dx} = \sum_{k=0}^{\infty} a_k (k+r) x^{k+r-1} 


d2ydx2=k=0ak(k+r)(k+r1)xk+r2\frac{d^2 y}{dx^2} = \sum_{k=0}^{\infty} a_k (k+r)(k+r-1) x^{k+r-2}

Substituting these into the Bessel's equation gives:

x2k=0ak(k+r)(k+r1)xk+r2+xk=0ak(k+r)xk+r1+(x2n2)k=0akxk+r=0x^2 \sum_{k=0}^{\infty} a_k (k+r)(k+r-1) x^{k+r-2} + x \sum_{k=0}^{\infty} a_k (k+r) x^{k+r-1} + (x^2 - n^2) \sum_{k=0}^{\infty} a_k x^{k+r} = 0

x2k=0ak(k+r)(k+r1)xk+r2+xk=0ak(k+r)xk+r1+(x2n2)k=0akxk+r=0

Simplify each term:

k=0ak(k+r)(k+r1)xk+r+k=0ak(k+r)xk+r+k=0akxk+r+2n2k=0akxk+r=0\sum_{k=0}^{\infty} a_k (k+r)(k+r-1) x^{k+r} + \sum_{k=0}^{\infty} a_k (k+r) x^{k+r} + \sum_{k=0}^{\infty} a_k x^{k+r+2} - n^2 \sum_{k=0}^{\infty} a_k x^{k+r} = 0

k=0ak(k+r)(k+r1)xk+r+k=0ak(k+r)xk+r+k=0akxk+r+2n2k=0akxk+r=0

Step 2: Combine Like Terms

Combine the series:

k=0ak(k+r)(k+r1)xk+r+k=0ak(k+r)xk+r+k=0akxk+r+2n2k=0akxk+r=0\sum_{k=0}^{\infty} a_k (k+r)(k+r-1) x^{k+r} + \sum_{k=0}^{\infty} a_k (k+r) x^{k+r} + \sum_{k=0}^{\infty} a_k x^{k+r+2} - n^2 \sum_{k=0}^{\infty} a_k x^{k+r} = 0

k=0ak(k+r)(k+r1)xk+r+k=0ak(k+r)xk+r+k=0akxk+r+2n2k=0akxk+r=0

Rewrite the third sum to match the powers of xx:

k=0ak(k+r)(k+r1)xk+r+k=0ak(k+r)xk+r+k=2ak2xk+rn2k=0akxk+r=0\sum_{k=0}^{\infty} a_k (k+r)(k+r-1) x^{k+r} + \sum_{k=0}^{\infty} a_k (k+r) x^{k+r} + \sum_{k=2}^{\infty} a_{k-2} x^{k+r} - n^2 \sum_{k=0}^{\infty} a_k x^{k+r} = 0

k=0ak(k+r)(k+r1)xk+r+k=0ak(k+r)xk+r+k=2ak2xk+rn2k=0akxk+r=0

Now we have:

k=0[ak(k+r)(k+r1)+ak(k+r)+ak2n2ak]xk+r=0\sum_{k=0}^{\infty} \left[ a_k (k+r)(k+r-1) + a_k (k+r) + a_{k-2} - n^2 a_k \right] x^{k+r} = 0

Step 3: Indicial Equation and Recurrence Relation

The indicial equation is obtained by setting k=0k = 0:

a0r(r1)+a0rn2a0=0a_0 r(r-1) + a_0 r - n^2 a_0 = 0

r(r1)+rn2=0r(r-1) + r - n^2 = 0

r2n2=0r^2 - n^2 = 0

r=±nr = \pm n

So, the roots are r1=nr_1 = n and r2=nr_2 = -n.

Next, we find the recurrence relation by setting the coefficients of xk+rx^{k+r} to zero for k2k \ge 2:

ak[(k+r)(k+r1)+(k+r)n2]+ak2=0a_k \left[ (k+r)(k+r-1) + (k+r) - n^2 \right] + a_{k-2} = 0

ak[k2+2kr+r2n2]+ak2=0a_k \left[ k^2 + 2kr + r^2 - n^2 \right] + a_{k-2} = 0

ak[(k+r)2n2]+ak2=0a_k \left[ (k+r)^2 - n^2 \right] + a_{k-2} = 0

Step 4: Solving the Recurrence Relation

For r=nr = n:

ak[(k+n)2n2]+ak2=0a_k \left[ (k+n)^2 - n^2 \right] + a_{k-2} = 0

ak[k2+2kn]+ak2=0a_k \left[ k^2 + 2kn \right] + a_{k-2} = 0

akk(k+2n)+ak2=0a_k k(k+2n) + a_{k-2} = 0

ak=ak2k(k+2n)a_k = -\frac{a_{k-2}}{k(k+2n)}

For r=nr = -n:

ak[(kn)2n2]+ak2=0a_k \left[ (k-n)^2 - n^2 \right] + a_{k-2} = 0

ak[k22kn]+ak2=0a_k \left[ k^2 - 2kn \right] + a_{k-2} = 0

akk(k2n)+ak2=0a_k k(k-2n) + a_{k-2} = 0

ak=ak2k(k2n)a_k = -\frac{a_{k-2}}{k(k-2n)}

Step 5: General Solution

The general solution to Bessel's equation is the linear combination of the two independent solutions y1(x)y_1(x) and y2(x)y_2(x) corresponding to r=nr = n and r=nr = -n:

y(x)=AJn(x)+BJn(x)y(x) = A J_n(x) + B J_{-n}(x)

where Jn(x)J_n(x) and Jn(x)J_{-n}(x) are the Bessel functions of the first kind of order nn and n-n, respectively. For non-integer nn, Jn(x)J_{-n}(x) is independent of Jn(x)J_n(x). For integer nn, Jn(x)J_{-n}(x) is not independent, and the second linearly independent solution is given by the Bessel function of the second kind, Yn(x)Y_n(x):

y(x)=AJn(x)+BYn(x)y(x) = A J_n(x) + B Y_n(x)

Thus, we have derived the general solution of Bessel's equation.