¡@

Bessel Functions

Download Maple Worksheet

Bessel function of order m is a solution to the differential equation (Bessel equation)

x^2 d*y^2/(d^2*x)+x d*y/(d*x)+(x^2-m^2)*y = 0

The Bessel functions are named after the German astronomer Friedrich Bessel (1784-1846). These functions first arose when Bessel solved Kepler's equation for describing planetary motion. These functions have been applied in many different physical situation, including the temperature distribution in a circular plate and the shape of a vibrating drumhead.

The Bessel function of order 0 is defined by

J[0](x) = sum((-1)^n*x^(2*n)/(2^(2*n)*n!^2),n = 0 ....

Let  a[n] = (-1)^n*x^(2*n)/(2^(2*n)*n!^2) .

Then

abs(a[n+1]/a[n]) = abs((-1)^(n+1)*x^(2*(n+1))/(2^(2...

 

 = x^(2*n+2)*2^(2*n)*n!^2/(2^(2*n+2)*(n+1)^2*n!^2*x^(2...

 

             = x^2/(4*(n+1)^2) approaches to 0 < 1 for all x

 

Thus , by Ratio test , the series converges for all values of x .   In other words, the domain of the Bessel function J[0] is (-infinity, infinity) = R .

 

Here are the first four partial sums and their graphs :

1-1/4*x^2

1-1/4*x^2+1/64*x^4

1-1/4*x^2+1/64*x^4-1/2304*x^6

1-1/4*x^2+1/64*x^4-1/2304*x^6+1/147456*x^8


[Maple Plot]

 

[Maple Plot]

We see that J[0],  the Bessel function of order 0, is an even function with J[0](0) = 1.   Moreover, by the Alternating Series Estimation Theorem, we have

 

abs(J[0](x)-sum((-1)^k*x^(2*k)/(2^(2*k)*k!^2),k = 0...

 

 

The animation below shows how the first 25 partial sums approach the Bessel function of order 0.

[Maple Plot]

 

The Bessel function of order 1 is defined by

J[1](x) = sum((-1)^n*x^(2*n+1)/(2^(2*n+1)*n!*(n+1)!...

Let b[n] = (-1)^n*x^(2*n+1)/(2^(2*n+1)*n!*(n+1)!) .

Then

abs(b[n+1]/b[n]) = abs((-1)^(n+1)*x^(2*(n+1)+1)/(2^...

 

= abs(x^(2*n+3)*2^(2*n+1)*n!*(n+1)!/(2^(2*n+3)*(n+1)!...                 

 

= x^2/(4*(n+1)*(n+2)) approaches to 0 < 1 for all x          

 

Thus, by Ratio test, the series converges for all values of x . In other words, the domain of the Bessel function J[1]is also (-infinity, infinity) = R .

 

Here is the graph of J[1](x) where abs(x) <= 20 .

[Maple Plot]

 

We see that J[1] , the Bessel function of order 1, is an odd function with J[1](0) = 0 .

The animation below shows how the first 25 partial sums approach J[1] .

[Maple Plot]
 

Here we graph J[0](x) and J[1](x) on a common screen, can you see the relation between them ?

[Maple Plot]

It seems that d/(d*x) J[0](x) = -J[1](x),  let's check.

d/(d*x)J[0](x) = sum((-1)^n/(2^(2*n)*n!^2),n = 0 .. infini...d/(d*x)x^(2*n)

                     =sum((-1)^n*2*n/(2^(2*n)*n!^2),n = 1 .. infinity)x^(2*n-1)

                                      =sum((-1)^n/(2^(2*n-1)*(n-1)!*n!),n = 1 .. infinity)...x^(2*n-1)

                                                        =sum((-1)^(n+1)/(2^(2*n+1)*n!*(n+1)!),n = 0 .. infin...x^(2*n+1) = -J[1](x)

 

To show that J[0](x) satisfies the equation

x^2 d*y^2/(d^2*x)+x d*y/(d*x)+x^2*y = 0

 

Let y = J[0](x) , we have

x^2*y=x^2sum((-1)^n*x^(2*n)/(2^(2*n)*n!^2),n = 0 .. infinity...
             =
sum((-1)^n*x^(2*n+2)/(2^(2*n)*n!^2),n = 0 .. infini...

¡@

xd*y/(d*x)= xsum((-1)^(n+1)/(2^(2*n+1)*n!*(n+1)!),n = 0 .. infin... x^(2*n+1)
            =
sum((-1)^(n+1)/(2^(2*n+1)*n!*(n+1)!),n = 0 .. infin... x^(2*n+2)

 

x^2 d*y^2/(d^2*x) = x^2 sum((-1)^(n+1)*(2*n+1)/(2^(2*n+1)*n!*(n+1)!),n = 0 ... x^(2*n)
                     =
sum((-1)^(n+1)*(2*n+1)/(2^(2*n+1)*n!*(n+1)!),n = 0 ... x^(2*n+2)

Thus,

x^2d*y^2/(d^2*x)+xd*y/(d*x)+x^2*y= sum([(-1)^(n+1)*(2*n+1)/(2^(2*n+1)*n!*(n+1)!)+(-1)^...x^(2*n+2)

and

(-1)^(n+1)*(2*n+1)/(2^(2*n+1)*n!*(n+1)!)+(-1)^(n+1)...

= ((-1)^(n+1)*(2*n+1)*n!+(-1)^n*n!+(-1)^n*2*(n+1)!)/(... = 0

 

Therefore, x^2d*y^2/(d^2*x)+xd*y/(d*x)+x^2*y = 0.

 

Similarly, we can show that J[1](x)satisfies the equation

x^2 d*y^2/(d^2*x)+x d*y/(d*x)+(x^2-1)*y = 0

 

 

For more information about the Bessel functions and the Bessel equations, please refer to

¡@

http://mathworld.wolfram.com/BesselFunctionoftheFirstKind.html

¡@


Download Maple Worksheet