intro bessel bk nov08

7/28/2019 Intro Bessel Bk Nov08

1/39

MATH306 SUPPLEMENTARY MATERIAL

A BRIEF INTRODUCTION TO BESSEL and

RELATED SPECIAL FUNCTIONS

Edmund Y. M. Chiang

c Draft date December 1, 2008


2/39

Contents

Contents i

1 Trigonometric and Gamma Functions 3

1.1 Trigonometric Functions . . . . . . . . . . . . . . . . . . . . . . . . . 3

1.2 Gamma Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

1.3 Chebyshev Polynomials . . . . . . . . . . . . . . . . . . . . . . . . . . 4

2 Bessel Functions 9

2.1 Power Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

2.2 Bessel Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

2.3 Basic Properties of Bessel Functions . . . . . . . . . . . . . . . . . . . 12

2.3.1 Zeros of Bessel Functions . . . . . . . . . . . . . . . . . . . . . 12

2.3.2 Recurrence Formulae for J . . . . . . . . . . . . . . . . . . . 12

2.3.3 Generating Function for Jn . . . . . . . . . . . . . . . . . . . . 13

2.3.4 Lommels Polynomials . . . . . . . . . . . . . . . . . . . . . . 14

2.3.5 Bessel Functions of half-integer Orders . . . . . . . . . . . . . 14

2.3.6 Formulae for Lommels polynomials . . . . . . . . . . . . . . . 14

2.3.7 Pythagoras Theorem for Bessel Function . . . . . . . . . . 15

2.3.8 Orthogonality of the Lommel Polynomials . . . . . . . . . . . 15

2.4 Integral formaulae . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

2.5 Asymptotic Behaviours of Bessel Functions . . . . . . . . . . . . . . . 17

2.5.1 Addition formulae . . . . . . . . . . . . . . . . . . . . . . . . . 18

2.6 Fourier-Bessel Series . . . . . . . . . . . . . . . . . . . . . . . . . . . 18

i


3/39

CONTENTS 1

3 An New Perspective 21

3.1 Hypergeometric Equations . . . . . . . . . . . . . . . . . . . . . . . . 21

3.2 Confluent Hypergeometric Equations . . . . . . . . . . . . . . . . . . 23

3.3 A Definition of Bessel Functions . . . . . . . . . . . . . . . . . . . . . 23

4 Physical Applications of Bessel Functions 27

5 Separation of Variables of Helmholtz Equations 31

5.1 What is not known? . . . . . . . . . . . . . . . . . . . . . . . . . . . 36


4/39

2 CONTENTS


5/39

Chapter 1

Trigonometric and GammaFunctions

1.1 Trigonometric Functions

Pythagoras Theorem:

(1.1) sin2 x + cos2 x = 1

holds for all real x.

Addition Theorem:(1.2) sin(A + B) = sin A cos B + sin B cos A

(1.3) cos(A + B) = cos A cos B sin A sin Bhold for all angles A, B.

1.2 Gamma Function

We recall that(1.1) k! = k (k 1) 3 2 1.Euler was able to give a correct definition to k! when k is not a positive integer. Heinvented the Euler-Gamma function in the year 1729. That is,

(1.2) (x) =

0

tx1et dt.

3


6/39

4 CHAPTER 1. TRIGONOMETRIC AND GAMMA FUNCTIONS

So the integral will converge for all positive real x. Since

(1.3) (1) =

0

et dt = 1,

so one can show

(1.4) (x + 1) = x(x),

and so for each positive integer n

(1.5) (n + 1) = n!.

In fact, the infinite integral will not only converge for all positive real x but canbe analytically continued into the whole complex plane C. We have

(1.6) (z) =

0

tz1et dt,

(z) > 0,

which is analytic in C except at the simple poles at the negative integers including0. So it is a meromorphic function.

Euler worked out that

(1.7)1

2

! =

1

2

.

So

(1.8)

51

2

! =

11

2

9

2,

7

2,

5

2,

3

2

12

!.

One can even compute negative factorial:(1.9) (1/2) = (1/2)(3/2)(5/2)(19/2)(11/2)(11/2).

1.3 Chebyshev Polynomials

Theorem 1.3.1 (De Moivres Theorem). Let z = r(cos + i sin ) be a complexnumber. Then for every non-negative integer n,

(1.1) zn = rn(cos n + i sin n).

Proof It is clear that that the formula holds when n = 1. We assume it holds at theinteger n. Then the addition formulae for the trigonometric theorem imply that

zn+1 = rn+1(cos + i sin )(cos n + i sin n)

= rn+1[cos(n + 1) + i sin(n + 1)].

Thus, the result follows by applying the principle of induction.


7/39

1.3. CHEBYSHEV POLYNOMIALS 5

Exercise 1.3.2. Justify the formula for negative integers n.

Remark 1.3.1. De Moivres formula actually works for all real and complex n. SeeMATH304.

Example 1.3.3 (Multiple angle formulae). Letn = 2 in De Movires formula yields

(1.2) z2 = cos2 + i sin2.

But

(1.3) z2 = (cos + i sin )2 = cos2 sin2 + i(2cos sin ).Comparing the real and imaginary parts of (1.2) and (1.3) yields

cos2 = cos2 sin2 = 2cos2 1= 1 2sin2

and

(1.4) sin 2 = 2cos sin .

Similarly, we have

(1.5) z3 = cos3 + i sin3.

But

z3 = (cos + i sin )3

= cos3 + 3 cos2 (i sin ) + 3 cos (i sin )2 + (i sin )3

= cos3 3cos sin2 + i(3cos2 sin sin3 ).(1.6)Comparing the real and imaginary parts of (1.5) and (1.6) yields

cos3 = cos3 3cos sin2 = 4cos3 3cos

and

sin3 = 3cos2 sin sin3 = 3sin 4sin3 .


8/39


More generally, we have

cos n + i sin n = zn = (cos + i sin )n

=n

k=0

n

k

cosnk (i sin )k

=n

k=0

n

k

ik cosnk sink .(1.7)

Equating the real parts of the last equation and noting that i2l = (1)l give

cos n = cosn

n

2

cosn2 sin2 +

n

4

cosn4 sin4 +

+ (1)[n/2] n2[n/2] cosn2[n/2] sin2[n/2] .(1.8)But each of the sine function above has even power. Thus the substitution ofsin2 = 1 cos2 in (1.8) yields

(1.9) cos n =

[n/2]l=0

n

2l

cosn2l

lk=0

(1)k

l

k

cos2k

.

Finally, it is possible to re-write the summation formula as

(1.10) cos n =

[n/2]k=0

(1)k

[n/2]j=k

n2jj

k

cosn

2k .

Writing = arccos x for 0 , and Tn(x) = cos(n arccos x). We deduce fromthe above formula that Tn(x) is a polynomial in x of degree n. The polynomial Tnis called the Chebyshev1 polynomial of the first kind of degree n. The firstfew polynomials are given by

T0(x) = 1;

T1(x) = x;

T2(x) = 2x2

1;

T3(x) = 4x3 3x;T4(x) = 8x

4 8x2 + 1;T5(x) = 16x

5 20x3 + 5x;

1P. L. Chebyshev (18211894) Russian mathematician. The collection of the polynomials{Tn(x)}0 is a complete orthogonal polynomial set on [1, 1]


9/39

1.3. CHEBYSHEV POLYNOMIALS 7

The set of Chebyshev polynomials {Tn(x)} can generate a L2[1, 1] space. That is,it is a Hilbert space.

Interested reader should consult the handbook [1] and the book [5] about theChebyshev polynomials by T. J. Rivlin. Our library has a copy of this book.


10/39



11/39

Chapter 2

Bessel Functions

2.1 Power Series

We define the Bessel function of first kind of order to be the complex functionrepresented by the power series

(2.1) J(z) =+k=0

(1)k(12z)+2k(+ k + 1) k!

= z+k=0

(1)k(12)+2k(+ k + 1) k!

z2k.

Here is an arbitrary complex constant and the notation () is the Euler Gamma

function defined by

(2.2) (z) =

0

tz1et dt, (z) > 0.

2.2 Bessel Equation

A common introduction to the series representation for Bessel functions of first kindof order is to consider Frobenius method: Suppose

y(z) =k=0

ckz+k,

where is a fixed parameter and ck are coefficients to be determined. Substitutingthe above series into the Bessel equation

(2.1) z2d2y

dz2+ x

dy

dz+ (z2 2)y = 0.

9


12/39

10 CHAPTER 2. BESSEL FUNCTIONS

and separating the coefficients yields:

0 = z2 d

2y

dz2 + x

dy

dz + (z2

2

)y

=k=0

ck( + k)( + k 1)z+k +k=0

ck( + k)z+k

+ (z2 2) k=0

ckz+k

=k=0

ck[( + k) 2]z+k +k=0

ckz+k+2

= c0(2

2)z +

k=0{

ck[( + k)2

2] + ck2

}z+k

The first term on the right side of the above expression is

c0(2 2)z

and the remaining are

c1[( + 1)2 2] + 0 = 0

c2[( + 2)2 2] + c0 = 0

c3[( + 3)2

2] + c1 = 0

ck[( + k)

2 2] + ck2 = 0(2.2) (2.3)

Hence a series solution could exist only if = . When k > 1, then werequire

ck[( + k)2 2] + ck2 = 0,

and this determines ck in terms of ck2 unless = 2 or + = 2 is aninteger. Suppose we discard these exceptional cases, then it follows from (2.2) that

c1 = c3 = c5 = = c2k+1 = = 0.

Thus we could express the coefficients c2k in terms of

c2k =(1)kc0

( + 2)( + 4) ( + 2k) ( + + 2)( + + 4) ( + + 2k) .


13/39

2.2. BESSEL EQUATION 11

If we now choose = , then we obtain

(2.4) c0z

1 +

k=0

(

1)k(1

2

z)2k

k!(+ 1)(+ 2) (+ k)

.

Alternatively, if we choose = , then we obtain

(2.5) c0z

1 +k=0

(1)k(12

z)2k

k!(+ 1)(+ 2) (+ k)

.

Since c0 and c

0 are arbitrary, so we choose them to be

(2.6) c0 =

1

2(+ 1) , c

0 =

1

2(+ 1)so that the two series (2.4) and (2.5) can be written in the forms:

(2.7) J(z) =+k=0

(1)k(12z)+2k(+ k + 1) k!

, J(z) =+k=0

(1)k(12z)+2k(+ k + 1) k!

and both are called the Bessel function of order if the first kind.

It can be shown that the Wronskian of J and J is given by (G. N. WatsonA Treatise On The Theory Of Bessel Functions, pp. 4243):

(2.8) W(J, J) = 2sin z

.

This shows that the J and J forms a fundamental set of solutions when isnot equal to an integer. In fact, when = n is an integer, then we can easily checkthat

(2.9) Jn(z) = (1)nJn(z).

Thus it requires an extra effort to find another linearly independent solution.

It turns out a second linearly independent solution is given by

(2.10) Y(z) =J(z)cos J(z)

sin

when is not an integer. The case when is an integer n is defined by

(2.11) Yn(z) = limn

J(z)cos J(z)sin

.


14/39


The Y so defined is linearly independent with J for all values of .

In particular, we obtain

Yn(z) =1

n1k=0

(n k 1)!k!

z2

2kn(2.12)

+1

k=0

(1)k(z/2)n+2kk!(n + k)!

2log

z

2 (k + 1) (k + n + 1)

(2.13)

for | arg z| < and n = 0, 1, 2, with the understanding that we set the sumto be 0 when n = 0. Here the (z) = (z)/(z). We note that the function isunbounded when z = 0.

2.3 Basic Properties of Bessel Functions

The general reference for Bessel functions is G. N. Watsons classic: A Treatise onthe Theory of Bessel Functions, published by Cambridge University Press in 1922[6].

2.3.1 Zeros of Bessel Functions

See A. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions withformulas, graphs and mathematical tables, 10th Edt., National Bureau of Stan-dards, 1964. [1]. [3] and [6]

2.3.2 Recurrence Formulae for J

We consider arbitrary complex .

d

dzzJ(z) =

d

dz

(1)kz2+2k2+2k k!(+ k + 1)

=d

dz

(1)kz21+2k21+2k k!(+ k)

= zJ1(z).

But the left side can be expanded and this yields

(2.1) zJ(z) + J(z) = zJ1(z).


15/39

2.3. BASIC PROPERTIES OF BESSEL FUNCTIONS 13

Similarly,

(2.2)

d

dzz

nu

J(z) = z

J+1(z).

and this yields

(2.3) zJ(z) J(z) = zJ+1(z).Substracting and adding the above recurrence formulae yield

J1(z) + J+1(z) =2

zJ(z)(2.4)

J1(z) J+1(z) = 2J(z).(2.5)

2.3.3 Generating Function for Jn

Jacobi in 1836 gave

(2.6) e1

2z(t1

t) =

+k=

tkJk(z).

Many of the forumulae derived above can be obtained from this expression.

(2.7) e1

2z(t1

t) =

+

k=

ck(z)tk

for 0 < |t| < . We multiply the power series

(2.8) ezt

2 = 1 +(z/2)

1!t +

(z/2)2

1!t2 +

and

(2.9) ezt

2 = 1 (z/2)1!

t +(z/2)2

1!t2

Multiplying the two series and comparing the coefficients of tk yield

cn(z) = Jn(z), n = 0, 1, (2.10)cn(z) = (1)nJn(z), n = 1, 2, .(2.11)

Thus

(2.12) e1

2z(t 1

t) = J0(z) +

+k=1

Jk[tk + (1)ktk].


16/39


2.3.4 Lommels Polynomials

Iterating the recurrence formula

(2.13) J+1(z) =2

zJ(z) J1

with respect to a number of times give

(2.14) J+k(z) = P(1/z)J(z)Q(1/z)J1.Lommel (1871) [6, pp. 294295] found that

(2.15) J+k(z) = Rk,(z)J(z)Rk1,+1J1.

2.3.5 Bessel Functions of half-integer Orders

One can check that

(2.16) J

1

2

(z) = 2

z

12

cos z, J12

(z) = 2

z

12

sin z.

Moreover,

(2.17) Jn+ 12

(z) = (1)n

2

zn+

1

2

dz dz

n sin zz

, n = 0, 1, 2, .

Thus applying a recurrence formula and using the Lommel polynomials yield

(2.18) Jn+ 12

(z) = Rn,(z)J12

(z)Rn1,+1J12

(z)

That is, we have

(2.19) Jn+ 12

(z) = Rn,(z) 2

z

12

sin zRn1,+1 2

z

12

sin z.

2.3.6 Formulae for Lommels polynomials

For each fixed , the Lommel polynomials are given by

(2.20) Rn(z) =

[n/2]k=0

(1)k(n k)!()nkk!(n 2k)!()k

2

z

n2k

where the [x] means the largest integer not exceeding x.

Lommel is a German mathematician who made a major contribution to Besselfunctions around 1870s.


17/39

2.3. BASIC PROPERTIES OF BESSEL FUNCTIONS 15

2.3.7 Pythagoras Theorem for Bessel Function

These Lommel polynomials have remarkable properties. Since

(2.21) J

1

2

(z) = 2

z

12

cos z, J12

(z) = 2

z

12

sin z

and sin2 x + cos2 x = 1; we now have

(2.22) J2n+ 12

(z) + J2n1

2

(z) = 2(1)nR2n, 1

2n(z)

z.

That is, we have

(2.23) J2n+ 1

2

(z) + J2n1

2

(z) =2

z

n

k=0

(2z)2n2k(2n k)!(2n 2k)![(n k)!]

2

k!

.

A few special cases are

1. J212

(z) + J2

1

2

(z) =2

z;

2. J232

(z) + J2

3

2

(z) =2

z

1 +

1

z2

;

3. J252 (z) + J

25

2 (z) =

2

z

1 +

3

z2 +

9

z4

;

4. J272

(z) + J2

7

2

(z) =2

z

1 +

6

z2+

45

z4+

225

z6

2.3.8 Orthogonality of the Lommel Polynomials

Let us set

(2.24) hn, (z) = Rn, (1/z) =

[n/2]k=0

(1)k

(n k)!()nkk!(n 2k)!()k

z2n2k

,

then the set {hn (z) is called the modified Lommel polynomials. Since the Besselfunctions J(z) satisfies a three-term recurrence relation, so the Lommel polynomialsinherit this property:

(2.25) 2z(n + )hn (z) = hn+1, (z) + hn1, (z)


18/39


with initial conditions

(2.26) h0, (z) = 1, h1, (z) = 2z.

If one start with a different set of initial conditions

(2.27) h0, (z) = 0, h

1, (z) = 2,

then the sequence {hn, (z)} generated by the (2.25) is called the associated Lom-mel polynomials. It is known that a three-term recurrence relation for polynomialswith the coefficients of as in (2.25) will generate a set of orthogonal polynomialson (, +). That is, there is a probability measure onb (, +) with+ d = 1

Theorem 2.3.1 (A. A. Markov, 1895). Suppose the set of {pn(z)}of orthogonalpolynomials with its measure supported on a bounded internal [a, b], then

(2.28) limn

pn (z)

pn (z)=

ba

d(t)

z tholds uniformly for z [a, b].

Since we know

Theorem 2.3.2 (Hurwitz). The limit

(2.29) limn

(z/2)+nRn, +1(z)

(n + + 1)= J(z),

holds uniformly on compact subsets of C.

So we have

Theorem 2.3.3. For > 0, the polynomials{hn (z) are orthogonal with respect toa discrete measure normalized to have

+ d(t) = 1, and

(2.30)

R

d(t)

z t = 2J(1/z)

J1(1/z).

One can use this formula to recover the measure (t). Thus the Lommel poly-nomials was found to have very distinct properties to be complete in L2(1, +1).


19/39

2.4. INTEGRAL FORMAULAE 17

2.4 Integral formaulae

(2.1)J(z) =

(z/2)

(12

)(+ 12

)

11

(1 t2)1/2 cos ztdt, () > 1/2, | arg z| < .

Or equivalently,(2.2)

J(z) =(z/2)

(12)(+12)

0

cos(zcos )sin2d, () > 1/2, | arg z| < ,

where t = cos .

Writing H(1) (z) = J(z) + iY(z). Then

Theorem 2.4.1 (Hankel 1869).

(2.3) H(1) (z) = 2

z

12 ei[z/2/4]

(+ 1/2)

exp(i)

0

euu1

2

1 +

iu

2z

12

du,

where || < /2.

2.5 Asymptotic Behaviours of Bessel Functions

Expanding the integrand of Henkels contour integral by binomial expansion andafter some careful analysis, we have

Theorem 2.5.1. For < arg z < 2,

(2.1) H(1) (x) =

2

x

12

ei(x1

21

4)

p1m=0

(12 )m(12 + )m

(2ix)mm!+ R(1)p (x)

,

and for2 < arg x < ,

(2.2) H(2) (x) = 2x

1

2

ei(x1

2 1

4)

p1

m=0

(12 )m(12 + )m

(2ix)mm!+ R(2)p (x),

where

(2.3) R(1)p (x) = O(xp) and R(2)p (x) = O(x

p),

as x +, uniformly in + < arg x < 2 and2 + < arg x < respectively.


20/39


The two expansions are valid simultaneously in < arg x < .We thus have J(z) =

12

(H(1) (z) + H

(2) (z)). So

Theorem 2.5.2. For | arg z| < , we have

J(z)

2

z

12

cos

z 12

14

k=0

(1)k(, 2k)(2z)2k

sin z 12

14

k=0

(1)k(, 2k + 1)(2z)2k+1

2.5.1 Addition formulae

Schlafli (1871) derived

(2.4) J(z+ t) =

k=

Jk(t)Jk(z).

Theorem 2.5.3 (Neumann (1867)). Let z, Z, R forms a triangle and let be theangle opposite to the side R, then

(2.5) J0

(Z2 + z2 2Zzcos ) = k=0

kJk(Z)Jk(z)cos k,

where 0 = 1, k = 2 for k 1.

We note that Z, z can assume complex values. There are generalizations to= 0.

2.6 Fourier-Bessel Series

Theorem 2.6.1. Suppose the real function f(r) is piecewise continuous in (0, a)and of bounded variation in every [r1, r2] (0, a). Then if

(2.1)

a0

r|f(r)| dr < ,


21/39

2.6. FOURIER-BESSEL SERIES 19

then the Fouier-Bessel series

(2.2) f(r) =

k=1

ckJ

xkr

a

converges to f(r) at every continuity point of f(r) and to

(2.3)1

2[f(r + 0) + f(r 0)]

at every discontinuity point of f(r). Here

(2.4) ck =2

J2(k)

10

xf(x)J(kx) dx

is the Fourier-Bessel coefficients of f.

The xk are the zero ofJ(x) k = 0, 1, . The discussion of the orthogonalityof the Bessel functions is omitted. You may consult G. P. Tolstovs Fourier Series,Dover 1962.

The situation here is very much like the classical Fourier series in terms of sineand cosine functions.


22/39



23/39

Chapter 3

An New Perspective

3.1 Hypergeometric Equations

We shall base our consideration of functions on the finite complex plane C in thischapter. Let p(z) and q(z) be meromorphic functions defined in C. A point z0belongs to C is called a regular singular point of the second order differential equation

(3.1)d2f(z)

dz2+ p(z)

df(z)

dz+ q(z)f(z) = 0

if

(3.2) limzz0

(z z0)p(z) = A, limzz0

(z z0)2q(z) = B

both exist and being finite. We start by recalling that any second order linearequation with three regular singular points in the C can always be transformed ([2,pp. 7375]) to the hypergeometric differential equation (due to Euler, Gaussand Riemann)

(3.3) z(1 z)d2y

dz2+ [c (a + b + 1)z]dy

dz ab y = 0

which has regular singular points at 0, 1, .

It is known that the hypergeometric equation has a power series solution

(3.4) 2F1

a, b

c

z

= 2F1 [a, b; c; z] =+k=0

(a)k(b)k(c)k k!

zk

21


24/39

22 CHAPTER 3. AN NEW PERSPECTIVE

which converges in |z| < 1. It is called the Gauss hypermetric series whichdepends on three parameters a, b, c. Here the notation is the standard Pochhammer

notation (1891)

(a)k = a (a + 1) (z+ 2) (a + k 1), ()0 = 1

for each integer k.

Theorem 3.1.1 (Euler 1769). Suppose(c) > (b) > 0, then we have

(3.5) 2F1

a, b

c

z

=

(c)

(b)(c b)10

tb1(1 t)cb1(1 xt)a dt,

defined in the cut-plane C\[1, +), arg t = arg(1 t) = 0, and that (1 xt)atakes values in the principal branch (, ).

Thus Eulers integal representation of 2F1 thus continues the function analyt-ically into the cut-plane.

Special cases of (3.4) include: ...

In 1812 Gauss gave a complete set of 15 contiguous relations ...

We note that one can write familiar functions in terms of pFq:

1. ex = 0F0

x

;

2. (1 x)a = 1F0

ax

;

3. log(1 + z) = x 2F1

1, 1

2

x

;

4. cos x = 0F1

1/2

x24

;

5. sin1 x = x 2F1

1/2, 1/2

3/2

x2

.


25/39

3.2. CONFLUENT HYPERGEOMETRIC EQUATIONS 23

3.2 Confluent Hypergeometric Equations

Writing the (3.1) in the variable x = z/b results in the equation

(3.1) x(1 x/b) d2y

dx2+ [c (a + b + 1)x/b] dy

dx a y = 0

which has regular singular points at 0, b, . Letting b tends to results in

(3.2) xd2y

dx2+ (c x) dy

dx a y = 0

which is called the confluent hypergeometric equation. But the equation hasa regular singular point at 0 and an irregular singular point at .

Substituting z/b in (3.4) and letting b to infinity results in a series

(3.3) (a, c; z) := 1F1 [a; c; z] =+k=0

(a)k(c)k k!

zk.

This series is called a Kummer series. We note that the Gauss hypergeometricseries has radius of convergence 1, and after the above substitution ofz/b results ina series that has radius of convergence equal to |b|. So the Kummer series has radiusof convergence after taking the limit b . We call all solutions to the equation(3.2) confluent hypergeometric functions because they are obtained from the

confluence of singularities of the hypergeometric equations. Further details can befound in [4], better known as the Bateman project volumn one.

3.3 A Definition of Bessel Functions

We define the Bessel function of first kind of order to be the complex functionrepresented by the power series

(3.1) J(z) =

+k=0

(1)k

(

1

2z)

+2k

(+ k + 1) k!= z

+k=0

(1)k

(

1

2)

+2k

(+ k + 1) k!z2k.

Here is an arbitrary complex constant and the notation () is the Euler Gammafunction defined by

(3.2) (z) =

0

tz1et dt, (z) > 0.


26/39


It can be analytically continued onto the whole complex plane C, and has simplepoles at the negative integers including 0. So it is a meromorphic function. In

particular, (k + 1) = k! = k(k 1) 2 1. So the Euler-Gamma function is ageneralization of the factorial notation k!.

Let us set a = + 12 , c = 2+ 1 and replace z by 2iz in the Kummer series.

That is,

(+ 1/2, 2+ 1; 2iz) =k=0

(+ 12)k

(2+ 1)k k!(2iz)k

= 1F2

+ 12

2+ 1

2iz

= eiz0F1

+ 1

z22

= eizk=0

(1)k(+ 1)k k!

z2k

22k

= eizk=0

(1)k(+ 1)(+ k + 1) k!

z2

2k

= eiz(+ 1)z

2

z

2

k=0

(1)k(+ k + 1) k!

z2

2k

= eiz(+ 1)z

2 J(z),

where we have applied Kummers second transformation

(3.3) 1F1

a

2a

4x

= e2x 0F1

a + 1

2

x2

in the third step above, and the identity ()k = (+ k)/() in the fourth step.

Since the confluent hypergeometric function ( + 1/2, 2 + 1; 2iz) satisfiesthe hypergeometric equation

(3.4) zd2ydz2

+ (c z) dydz a y = 0

with a = + 1/2 and c = 2+1. Substituting eiz(+ 1)z2

J(z) into the (3.4)

and simplfying lead to the equation

(3.5) z2d2y

dz2+ x

dy

dz+ (z2 2)y = 0.


27/39

3.3. A DEFINITION OF BESSEL FUNCTIONS 25

The above derivation shows that the Bessel function of the first kind with order is a special case of the confluent hypergeometric functions with specifically chosen

parameters.


28/39



29/39

Chapter 4

Physical Applications of BesselFunctions

We consider radial vibration of circular membrane. We assume that an elasticcircular membrane (radius ) can vibrate and that the material has a uniform density.Let u(x, y, t) denote the displacement of the membrane at time t from its equilibriumposition. We use polar coordinate in the xyplane by the change of variables:(4.1) x = r cos , y = r sin .

Then the corresponding equation

(4.2) 2u

t2= c2

2ux2

+ 2u

y2

can be transformed to the form:

(4.3)2u

t2= c2

2ur2

+1

r

2u

r+

1

r22u

2

.

Since the membrane has uniform density, so u = u(r, t) that is, it is independent ofthe . Thus we have

(4.4)2u

t2= c2

2ur2

+1

r

2u

r

.

The boundary condition is u(, t) = 0, and the initial conditions take the form

(4.5) u(r, t) = f(r),u(r, 0)

t= g(r).

Separation of variables method yields

(4.6) u(r, t) = R(r)T(t),

27


30/39

28 CHAPTER 4. PHYSICAL APPLICATIONS OF BESSEL FUNCTIONS

which satisfies the boundary condition u(, t) = 0. Thus

(4.7) RT

= c2

R

T +

1

r R

T

.

Hence

(4.8)R + (1/r)R

R+

T

c2T; = 2,

where the is a constant. Thus

R +1

rR + 2R = 0(4.9)

T + c22T = 0.(4.10)

We notice that the first equation is Bessel equation with = 0. Thus its generalsolution is given by

(4.11) R(r) = C1J0(r) + C2Y0(r).

. Since Y0 is unbounded when r = 0, so C2 = 0. Thus the boundary conditionimplies that

(4.12) J0() = 0,

implying that = is a zero of J0(). Setting C1 = 1, we obtain for each integern = 1, 2, 3, ,

(4.13) Rn(r) = J0(nr) = J0(n

r),

where n = n is the nth zero of J0(). Thus we have

(4.14) un(r, t) = (An cos cnt + Bn sin cnt) J0(nr),

for n = 1, 2, 3, . Thus the general solution is given by

(4.15)n=1

un(r, t) =n=1

(An cos cnt + Bn sin cnt) J0(nr),

and by

(4.16) f(r) = u(r, 0) =n=1

An J0(nr),


31/39

29

and

(4.17) g(r) =

u(r, 0)

tt=0 =

n=1

Bncn J0(nr).

Fourier-Bessel Series theory implies that

(4.18) An =2

2J21 (n)

0

rf(r)J0(nr) dr,

(4.19) Bn =2

cn2J21 (n)

0

rg(r)J0(nr) dr.


32/39

30 CHAPTER 4. PHYSICAL APPLICATIONS OF BESSEL FUNCTIONS


33/39

Chapter 5

Separation of Variables ofHelmholtz Equations

The Helmholtz equation named after the German physicst Hermann von Helmholtzrefers to second order (elliptic) partial differential equations of the form:

(5.1) (2 + k2) = 0,

where k is a constant. If k = 0, then it reduces to the Laplace equations.

In this discussion, we shall restrict ourselves in the Euclidean space R3. Oneof the most powerful theories developed in solving linear PDEs is the the method

of separation of variables. For example, the wave equation

(5.2)

2 1c2

2

t2

(r, t) = 0,

can be solved by assuming (r, t) = (r) T(t) where T(t) = eit. This yields

(5.3)

2 2

c2

(r) = 0,

which is a Helmholtz equation. The questions now is under what 3dimensionalcoordinate system (u1, u2, u3) do we have a solution that is in the separation ofvariables form

(5.4) (r) = 1(u1) 2(u2) 3(u3) ?

Eisenhart, L. P. (Separable Systems of Stackel. Ann. Math. 35, 284-305,1934) determines via a certain Stackel determinant is fulfilled (see e.g. Morse,P. M. and Feshbach, H. Methods of Theoretical Physics, Part I. New York:McGraw-Hill, pp. 125126, 271, and 509510, 1953).

31


34/39

32CHAPTER 5. SEPARATION OF VARIABLES OF HELMHOLTZ EQUATIONS

Theorem 5.0.1 (Eisenhart 1934). There are a total of eleven curvilinear coordi-nate systems in which the Helmholtz equation separates.

Each of the curvilinear coordinate is characterized by quadrics. That is,surfaces defined by

(5.5) Ax2 + By2 + Cz2 + Dxy + Exz + F yz+ Gx + Hy + lz+ J = 0.

One can visit http://en.wikipedia.org/wiki/Quadric for some of the quadricsurfaces. Curvilinear coordinate systems are formed by putting relevant orthogonalquadric surfaces. Wikipedia contains quite a few of these pictures. We list the elevencoordinate systems here:


35/39

33

Co-ordinate SystemsName of system Transformation formulae Degenerate Surfaces(1) Cartesian x = x, y = y, z = z

(2) Cylindricalx = cos , y = sin , z = z

0, <

(3) Spherical polarx = r sin cos , y = r sin sin ,z = r cos r 0, 0 ,

(4)Paraboliccylinder

x = u2 v2, y = 2uv,z = zu 0, < v < + Half-plane

(5)Ellipticcylinder

x = fcosh cos , y = fsinh sin ,z = z 0, < < +

Infinite strip;Plane withstraight aperture

(6)Rotationparaboloidal

x = 2uv cos , 2uv sin ,z = u2 v2u, v 0, < <

Half-line

(7)Prolatespheroidal

x = sinh u sin v cos ,y = sinh u sin v sin ,z = cosh u cos v,u 0, 0 v , <

Finite line ;segmentTwo half-lines

(8)Oblatespheroidal

x = cosh u sin v cos ,y = cosh u sin v sin ,z = sinh u cos v,u 0, 0 v , <

Circular plate (disPlane withcircular aperture

(9) Paraboloidal

x = 12(cosh2 + 2cos 2 cosh 2,y = 2 cosh cos sinh ,z = 2 sinh sin cosh ,, 0, <

Parabolic plate;Plane withparabolic aperture

(10) Elliptic conal

x = krsn sn ;y = (ik/k)rcn cn ,z = (1/k)rdn dn ;r 0, 2K < 2K, = K+ iu, 0

u

2K

Plane sector;Including quarterplane

= 2


36/39


Laplace & HelmholtzCoordinate system Laplace Equation Helmholtz equation

(1) Cartesian (Trivial) (Trivial)

(2) Cylinderical (Trivial) Bessel

(3) Spherical polar Associated Legender Associated Legender

(4) Parabolic cylinder (Trivial) Weber

(5) Elliptic cylinder (Trivial) Mathieu

(6) Rotation-paraboloidal Bessel Confluent hypergeometric

(7) Prolate spheroidal Associated Legender Spheroidal wave

(8) Prolate spheroidal Associated Legender Spheroidal wave

(9) Paraboloidal Mathieu Whittaker-Hill

(10) Elliptic conal Lame Spherical Bessel, Lame

(11) Ellipsoidal Lame Ellipsoidal


37/39

35

1. Associated Legender:

(5.6) (1

x2)d2y

dx2 2x

dy

dx+ n(n + 1)

m2

(1 x2)y = 02. Bessel:

(5.7) x2d2y

dx2+ x

dy

dx+ (x2 n2)y = 0

3. Spherical Bessel:

(5.8) x2d2y

dx2+ 2x

dy

dx+ (x2 n(n + 1))y = 0

4. Weber:

(5.9) d

2

ydx2 + ( 14x2)y = 05. Confluent hypergeometric:

(5.10) xd2y

dx2+ ( x) dy

dx y = 0

6. Mathieu:

(5.11)d2y

dx2+ ( 2qcos2x)y = 0

7. Spheroidal wave:

(5.12) (1 x2) d2ydx2

2x dydx

+

2(1 x2) +

2(1 x2)

y = 0

8. Lame:

(5.13)d2y

dx2+ (h n(n + 1)k2sn 2x)y = 0

9. Whittaker-Hill:

(5.14)d2y

dx2+ (a + b cos2x + c cos4x)y = 0

10. Ellipsoidal wave:(5.15)

d2y

dx2+ (a + bk2sn 2x + qk4sn 4x)y = 0

Remark 5.0.1. The spheroidal wave and the Whittaker-Hill do not belong to thehypergeometric equation regime, but to the Heun equation regime (which has fourregular singular points).


38/39


5.1 What is not known?

It is generally regarded that the Bessel functions, Weber functions, Legendre func-tions are better understood, but the remaining equations/functions are not so wellunderstood.

1. Bessel functions. OK! Still some unknowns.

2. Confluent hypergeometric equations/functions. NOT OK.

3. Spheroidal wave, Mathieu, Lame, Whittaker-Hill, Ellipsoidal wave are poorlyunderstood. Some of them are relatd to the Heun equation which has fourregular singular points. Its research has barely started despite the fact that

it has been around since 1910.4. Mathematicans are separating variables of Laplace/Helmholtz equations, but

in more complicated setting (such as in Riemannian spaces, etc)


39/39

Bibliography

[1] Milton Abramowitz and Irene A. Stegun, editors. Handbook of mathematicalfunctions, with formulas, graphs, and mathematical tables, volume 55 ofNationalBureau of Standards Applied Mathematics Series. Superintendent of Documents,

US Government Printing Office, Washington, DC, 1965. Third printing, withcorrections.

[2] G. E. Andrews, R. A. Askey, and R. Roy. Special Functions. Cambridge Uni-versity Press, Cambridge, 1999.

[3] A. Erdelyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi. Higher Transcen-dental Functions, volume 2. McGraw-Hill, New York, 1953.

[4] A. Erdelyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi. Higher Transcen-dental Functions, volume 1. McGraw-Hill, New York, 1953.

[5] T. J. Rivlin and M. W. Wilson. An optimal property of Chebyshev expansions.J. Approximation Theory, 2:312317, 1969.

[6] G. N. Watson. A Treatise on the Theory of Bessel Functions. Cambridge Uni-versity Press, Cambridge, second edition, 1944.