GV

GeezerVison

The World Wide Web, restored to 800 pixels of glory.

Home|About|Gaming cabinet|Archives|How it works|Starter sites

On this site

  • Algebra
  • Applied Mathematics
  • Calculus and Analysis
  • Discrete Mathematics
  • Foundations of Mathematics
  • Geometry
  • History and Terminology
  • Number Theory
  • Probability and Statistics
  • Recreational Mathematics
  • Topology
  • Alphabetical Index
  • Random Entry
  • New in MathWorld

This view

  • Original page
  • New address

GeezerVison is showing mathworld.wolfram.com — Bessel Function of the First Kind -- from Wolfram MathWorld

Search

Bessel Function of the First Kind

DOWNLOAD Mathematica NotebookDownload Wolfram Notebook
BesselJ

The Bessel functions of the first kind J_m(x) are conventionally normalized solutions to the Bessel differential equation

x^2(d^2y)/(dx^2)+x(dy)/(dx)+(x^2-m^2)y=0.
(1)

For nonnegative integer order m, they are nonsingular at the origin. They are sometimes also called cylinder functions or cylindrical harmonics. The above plot shows J_n(x) for n=0, 1, 2, ..., 5. The notation J_(z,n) was first used by Hansen (1843) and subsequently by Schlömilch (1857) to denote what is now written J_n(2z) (Watson 1966, p. 14). However, Hansen's definition of the function itself in terms of the generating function

e^(z(t-1/t)/2)=sum_(n=-infty)^inftyt^nJ_n(z).
(2)

is the same as the modern one (Watson 1966, p. 14). Bessel used the notation I_k^h to denote what is now called the Bessel function of the first kind (Cajori 1993, vol. 2, p. 279).

For integer n, the Bessel function J_n(z) can also be defined by the contour integral

J_n(z)=1/(2pii)∮e^((z/2)(t-1/t))t^(-n-1)dt,
(3)

where the contour encloses the origin and is traversed in a counterclockwise direction (Arfken 1985, p. 416).

The Bessel function of the first kind is implemented in the Wolfram Language as BesselJ[nu, z].

To solve the Bessel differential equation for x>0, keep the order m fixed and apply Frobenius method using a series solution of the form

y=x^ksum_(n=0)^inftya_nx^n=sum_(n=0)^inftya_nx^(n+k).
(4)

Substituting into the Bessel differential equation yields

x^2sum_(n=0)^infty(k+n)(k+n-1)a_nx^(k+n-2)+xsum_(n=0)^infty(k+n)a_nx^(k+n-1)+x^2sum_(n=0)^inftya_nx^(k+n)-m^2sum_(n=0)^inftya_nx^(n+k)=0
(5)
sum_(n=0)^infty(k+n)(k+n-1)a_nx^(k+n)+sum_(n=0)^infty(k+n)a_nx^(k+n) +sum_(n=2)^inftya_(n-2)x^(k+n)-m^2sum_(n=0)^inftya_nx^(n+k)=0.
(6)

The indicial equation, obtained by equating the coefficient of x^k to zero, is

a_0[k(k-1)+k-m^2]=a_0(k^2-m^2)=0.
(7)

Since a_0 is defined as the first nonzero coefficient, k^2-m^2=0, so k=+/-m. Now, if k=m,

sum_(n=0)^infty[(m+n)(m+n-1)+(m+n)-m^2]a_nx^(m+n)+sum_(n=2)^inftya_(n-2)x^(m+n)=0
(8)
sum_(n=0)^infty[(m+n)^2-m^2]a_nx^(m+n)+sum_(n=2)^inftya_(n-2)x^(m+n)=0
(9)
sum_(n=0)^inftyn(2m+n)a_nx^(m+n)+sum_(n=2)^inftya_(n-2)x^(m+n)=0
(10)
a_1(2m+1)x^(m+1)+sum_(n=2)^infty[a_nn(2m+n)+a_(n-2)]x^(m+n)=0.
(11)

First, look at the special case m=-1/2, then (11) becomes

sum_(n=2)^infty[a_nn(n-1)+a_(n-2)]x^(m+n)=0,
(12)

so

a_n=-1/(n(n-1))a_(n-2).
(13)

Now let n=2l, where l=1, 2, ....

a_(2l)=-1/(2l(2l-1))a_(2l-2)
(14)
=((-1)^l)/([2l(2l-1)][2(l-1)(2l-3)]...[2·1·1])a_0
(15)
=((-1)^l)/(2^ll!(2l-1)!!)a_0,
(16)

which, using the identity 2^ll!(2l-1)!!=(2l)!, gives

a_(2l)=((-1)^l)/((2l)!)a_0.
(17)

Similarly, letting n=2l+1,

a_(2l+1)=-1/((2l+1)(2l))a_(2l-1)=((-1)^l)/([2l(2l+1)][2(l-1)(2l-1)]...[2·1·3][1])a_1,
(18)

which, using the identity 2^ll!(2l+1)!!=(2l+1)!, gives

a_(2l+1)=((-1)^l)/(2^ll!(2l+1)!!)a_1=((-1)^l)/((2l+1)!)a_1.
(19)

Plugging back into (4) with k=m=-1/2 gives

y=x^(-1/2)sum_(n=0)^(infty)a_nx^n
(20)
=x^(-1/2)[sum_(n=1,3,5,...)^(infty)a_nx^n+sum_(n=0,2,4,...)^(infty)a_nx^n]
(21)
=x^(-1/2)[sum_(l=0)^(infty)a_(2l)x^(2l)+sum_(l=0)^(infty)a_(2l+1)x^(2l+1)]
(22)
=x^(-1/2)[a_0sum_(l=0)^(infty)((-1)^l)/((2l)!)x^(2l)+a_1sum_(l=0)^(infty)((-1)^l)/((2l+1)!)x^(2l+1)]
(23)
=x^(-1/2)(a_0cosx+a_1sinx).
(24)

The Bessel functions of order +/-1/2 are therefore defined as

J_(-1/2)(x)=sqrt(2/(pix))cosx
(25)
J_(1/2)(x)=sqrt(2/(pix))sinx,
(26)

so the general solution for m=+/-1/2 is

y=a_0^'J_(-1/2)(x)+a_1^'J_(1/2)(x).
(27)

Now take m>=0. Equation (11) requires

a_1(2m+1)=0
(28)
[a_nn(2m+n)+a_(n-2)]x^(m+n)=0
(29)

for n=2, 3, ..., so

a_1=0
(30)
a_n=-1/(n(2m+n))a_(n-2)
(31)

for n=2, 3, .... Let n=2l+1, where l=1, 2, ..., then

a_(2l+1)=-1/((2l+1)[2(m+l)+1])a_(2l-1)
(32)
=...=f(l,m)a_1=0,
(33)

where f(l,m) is the function of l and m obtained by iterating the recurrence relation down to a_1. Now let n=2l, where l=1, 2, ..., so

a_(2l)=-1/(2l(2m+2l))a_(2l-2)
(34)
=-1/(4l(m+l))a_(2l-2)
(35)
=((-1)^l)/([4l(m+l)][4(l-1)(m+l-1)]...[4·(m+1)])a_0.
(36)

Plugging back into (4),

y=sum_(n=0)^(infty)a_nx^(n+m)=sum_(n=1,3,5,...)^(infty)a_nx^(n+m)+sum_(n=0,2,4,...)^(infty)a_nx^(n+m)
(37)
=sum_(l=0)^(infty)a_(2l+1)x^(2l+m+1)+sum_(l=0)^(infty)a_(2l)x^(2l+m)
(38)
=a_0sum_(l=0)^(infty)((-1)^l)/([4l(m+l)][4(l-1)(m+l-1)]...[4(m+1)])x^(2l+m)
(39)
=a_0sum_(l=0)^(infty)((-1)^lGamma(m+1))/(2^(2l)l!Gamma(m+l+1))x^(2l+m).
(40)

Now define

J_m(x)=sum_(l=0)^infty((-1)^l)/(2^(2l+m)l!Gamma(m+l+1))x^(2l+m),
(41)

where Gamma(z) is the gamma function. This series also defines J_m(x) for any complex number m when x>0. For the solution derived above, the normalization gives

y=a_02^mGamma(m+1)J_m(x)=a_0^'J_m(x).
(42)

Changing the sign of m leaves the Bessel differential equation unchanged, so J_m and J_(-m) solve the same equation. For a nonnegative integer m, we can relate them by writing

J_(-m)(x)=sum_(l=0)^infty((-1)^l)/(2^(2l-m)l!(l-m)!)x^(2l-m).
(43)

Now let l=l^'+m. Then

J_(-m)(x)=sum_(l^'+m=0)^(infty)((-1)^(l^'+m))/(2^(2l^'+m)(l^'+m)!l^'!)x^(2l^'+m)
(44)
=sum_(l^'=-m)^(-1)((-1)^(l^'+m))/(2^(2l^'+m)l^'!(l^'+m)!)x^(2l^'+m)+sum_(l^'=0)^(infty)((-1)^(l^'+m))/(2^(2l^'+m)l^'!(l^'+m)!)x^(2l^'+m).
(45)

Here reciprocal factorials are interpreted using the gamma function, so 1/Gamma(l^'+1)=0 for l^'=-m, ..., -1. The finite sum therefore vanishes, giving

J_(-m)(x)=sum_(l=0)^(infty)((-1)^(l+m))/(2^(2l+m)l!(l+m)!)x^(2l+m)
(46)
=(-1)^mJ_m(x).
(47)

Note that the Bessel differential equation is second-order, so there must be two linearly independent solutions. We have found both only for |m|=1/2. For a general nonintegral order, the linearly independent solutions are J_m and J_(-m). When m is an integer, the real general solution is of the form

Z_m=C_1J_m(x)+C_2Y_m(x),
(48)

where J_m is a Bessel function of the first kind, Y_m (a.k.a. N_m) is the Bessel function of the second kind (a.k.a. Neumann function or Weber function), and C_1 and C_2 are constants. Complex solutions are given by the Hankel functions (a.k.a. Bessel functions of the third kind).

For a>0 and real nu>=0, the Bessel functions are orthogonal in [0,a] according to

int_0^aJ_nu(alpha_(num)rho/a)J_nu(alpha_(nun)rho/a)rhodrho=1/2a^2[J_(nu+1)(alpha_(num))]^2delta_(mn),
(49)

where alpha_(num) is the mth positive zero of J_nu and delta_(mn) is the Kronecker delta (Arfken 1985, p. 592).

Except when 2m is a negative integer,

J_m(z)=(z^(-1/2))/(2^(2m+1/2)i^(m+1/2)Gamma(m+1))M_(0,m)(2iz),
(50)

where Gamma(x) is the gamma function and M_(0,m) is a Whittaker function.

In terms of a confluent hypergeometric limit function, the Bessel function is written

J_nu(z)=((1/2z)^nu)/(Gamma(nu+1))_0F_1(nu+1;-1/4z^2).
(51)

A derivative identity for expressing higher order Bessel functions in terms of J_0(z) is

J_n(z)=i^nT_n(id/(dz))J_0(z),
(52)

where T_n(z) is a Chebyshev polynomial of the first kind. For fixed m>=0 and positive real z, asymptotic forms for the Bessel functions are

J_m(z) approx 1/(Gamma(m+1))(z/2)^m
(53)

for z<<1 and

J_m(z) approx sqrt(2/(piz))cos(z-(mpi)/2-pi/4)
(54)

for z>>|m^2-1/4| (correcting the condition of Abramowitz and Stegun 1972, p. 364).

A derivative identity is

d/(dx)[x^mJ_m(x)]=x^mJ_(m-1)(x).
(55)

An integral identity is

int_0^uu^'J_0(u^')du^'=uJ_1(u).
(56)

Some sum identities are

sum_(k=-infty)^inftyJ_k(x)=1
(57)

(which follows from the generating function (2) with t=1),

1=[J_0(x)]^2+2sum_(k=1)^infty[J_k(x)]^2
(58)

(Abramowitz and Stegun 1972, p. 363),

1=J_0(x)+2sum_(k=1)^inftyJ_(2k)(x)
(59)

(Abramowitz and Stegun 1972, p. 361),

0=sum_(k=0)^(2n)(-1)^kJ_k(z)J_(2n-k)(z)+2sum_(k=1)^inftyJ_k(z)J_(2n+k)(z)
(60)

for n>=1 (Abramowitz and Stegun 1972, p. 361),

J_n(2z)=sum_(k=0)^nJ_k(z)J_(n-k)(z)+2sum_(k=1)^infty(-1)^kJ_k(z)J_(n+k)(z)
(61)

(Abramowitz and Stegun 1972, p. 361), and the Jacobi-Anger expansion

e^(izcostheta)=sum_(n=-infty)^inftyi^nJ_n(z)e^(intheta),
(62)

which can also be written

e^(izcostheta)=J_0(z)+2sum_(n=1)^inftyi^nJ_n(z)cos(ntheta).
(63)

The Bessel function addition theorem states

J_n(y+z)=sum_(m=-infty)^inftyJ_m(y)J_(n-m)(z).
(64)

For positive integers m, the sums of k^m[J_k(x)]^2 over k>=1 give Bessel moment polynomials in x^2/4 for even m and finite expressions in J_0(x) and J_1(x) for odd m.

Various integrals can be expressed in terms of Bessel functions

J_n(z)=1/piint_0^picos(zsintheta-ntheta)dtheta,
(65)

which is Bessel's first integral,

J_n(z)=(i^(-n))/piint_0^pie^(izcostheta)cos(ntheta)dtheta
(66)
J_n(z)=1/(2pii^n)int_0^(2pi)e^(izcosphi)e^(inphi)dphi
(67)

for n=1, 2, ...,

J_n(z)=2/pi(z^n)/((2n-1)!!)int_0^(pi/2)sin^(2n)ucos(zcosu)du
(68)

for n=1, 2, ...,

J_n(x)=1/(2pii)int_gammae^((x/2)(z-1/z))z^(-n-1)dz
(69)

for n>-1/2. The Bessel functions are normalized so that

int_0^inftyJ_n(x)dx=1
(70)

for positive integral (and real) n. Integrals involving J_1(x) include

int_0^infty[(J_1(x))/x]^2dx=4/(3pi)
(71)
int_0^infty[(J_1(x))/x]^2xdx=1/2.
(72)

Ratios of Bessel functions of the first kind have continued fraction

(J_(n-1)(z))/(J_n(z))=(2n)/z-(z/(2(n+1)))/(1-(((z/2)^2)/((n+1)(n+2)))/((1-((z/2)^2)/((n+2)(n+3)))/(1-...)))
(73)

(Wall 1948, p. 349).

BesselJ0ReIm
BesselJ0Contours

The special case of n=0 gives J_0(z) as the series

J_0(z)=sum_(k=0)^infty(-1)^k((1/4z^2)^k)/((k!)^2)
(74)

(Abramowitz and Stegun 1972, p. 360), or the integral

J_0(z)=1/piint_0^pie^(izcostheta)dtheta.
(75)

See also

Bessel Function of the Second Kind, Bessel Function Zeros, Debye's Asymptotic Representation, Dixon-Ferrar Formula, Hansen-Bessel Formula, Kapteyn Series, Kneser-Sommerfeld Formula, Mehler's Bessel Function Formula, Modified Bessel Function of the First Kind, Modified Bessel Function of the Second Kind, Nicholson's Formula, Poisson's Bessel Function Formula, Rayleigh Function, Schläfli's Formula, Schlömilch's Series, Sommerfeld's Formula, Sonine-Schafheitlin Formula, Watson-Nicholson Formula, Watson's Formula, Weber's Discontinuous Integrals, Weber's Formula, Weber-Schafheitlin Integral, Weber-Sonine Formula, Weyrich's Formula Explore this topic in the MathWorld classroom

Related Wolfram sites

https://functions.wolfram.com/Bessel-TypeFunctions/BesselJ/

Explore with Wolfram|Alpha

WolframAlpha

More things to try:

  • Anger function
  • Bessel functions
  • Bessel function of the first kind of order zero

References

Abramowitz, M. and Stegun, I. A. (Eds.). "Bessel Functions J and Y." §9.1 in Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th printing. New York: Dover, pp. 358-364, 1972.Arfken, G. "Bessel Functions of the First Kind, J_nu(x)" and "Orthogonality." §11.1 and 11.2 in Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, pp. 573-591 and 591-596, 1985.Cajori, F. A History of Mathematical Notations, Vols. 1-2. New York: Dover, 1993.Hansen, P. A. "Ermittelung der absoluten Störungen in Ellipsen von beliebiger Excentricität und Neigung, I." Schriften der Sternwarte Seeberg. Gotha, Germany, 1843.Le Lionnais, F. Les nombres remarquables. Paris, France: Hermann, 1983.Lehmer, D. H. "Arithmetical Periodicities of Bessel Functions." Ann. Math. 33, 143-150, 1932.Morse, P. M. and Feshbach, H. Methods of Theoretical Physics, Part I. New York: McGraw-Hill, pp. 619-622, 1953.Schlömilch, O. X. "Ueber die Bessel'schen Function." Z. für Math. u. Phys. 2, 137-165, 1857.Spanier, J. and Oldham, K. B. "The Bessel Coefficients J_0(x) and J_1(x)" and "The Bessel Function J_nu(x)." Chs. 52-53 in An Atlas of Functions. Washington, DC: Hemisphere, pp. 509-520 and 521-532, 1987.Wall, H. S. Analytic Theory of Continued Fractions. New York: Chelsea, 1948.Watson, G. N. A Treatise on the Theory of Bessel Functions, 2nd ed. Cambridge, England: Cambridge University Press, 1966.

Referenced on Wolfram|Alpha

Bessel Function of the First Kind

Cite this as:

Weisstein, Eric W. "Bessel Function of the First Kind." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/BesselFunctionoftheFirstKind.html

Subject classifications

You are visitor 01402503 — GeezerVison does not store the pages you look up.
Best viewed with Netscape Navigator 4.7 or Internet Explorer 6 at 800×600. No cookies required.