Arsenin - Basic Equations and Special Functions of Mathematical Physics
Arsenin - Basic Equations and Special Functions of Mathematical Physics
Arsenin
Basic Equations
and
Special Functions
of
Mathematical Physics
Jl
Basic Equations and Special
Functions of Mathematical Physics
T
BASIC EQUATIONS AND SPECIAL
FUNCTIONS OF MATHEMATICAL PHYSICS
V. Ya. Arsenin
Translated by
S. Chomet, King’s College, London
Printed in Poland
Contents
Preface 7
PART 1
PART 2
Chapter 10 Gamma Function 215
Chapter 11 Cylindrical Functions 223
Chapter 12 Spherical Harmonics 268
Chapter 13 Chebyshev-Hermite and Chebyshev-Laguerre
Polynomials 291
Appendix Definition of Generalised Functions. The S
Function 303
Answers and Solutions 321
Bibliography 355
Index 357
5 A
Preface
V. Ya. Arsenin
7
PART ONE
1
1
4
Linear Equations
Substituting these derivatives into (la) and collecting up terms
involving the same derivatives, we obtain the transformed equation
aiiM« + 2 a 12 i/^ + a22M„+-F1(M^, un, u, £, rj) = 0 (4)
where
an = an <p2x-\-2al2<px<py+ci22<p2y
£(<?; vOl 2
*12 alia22 (^12 all 022) ( 6)
un, u, f, V) = 0 (7)
This is the canonical form.
We shall now summarise the procedure for finding y{x, ^) and
yt(x, y) without discussing the conditions for their existence:
a. If an — a22 — 0 in D, then an ^ 0. Dividing both sides of
(1) by 2ax2 we obtain the canonical form (7).
b. Suppose that *11+ 022 ^ 0 in D. We shall confine our attention
to the case when at least one of the coefficients au , a22 is not iden
tically zero in a part Dx of D. Suppose an has this property.
We shall take <p(x, y) and ip(x, y) in (3) to be functions which
ensure that the coefficients a n and a22 in the transformed equation (4)
5
Mathematical Physics
will vanish, i.e. we shall suppose that they are the solutions of the
following equations:
a\\(P2xJr'^a\2(Px(PyJr a22(pl = 0
auW l+^uVxW y+^wl = 0
These yield
<Px _ _ — a i2 ± V V 2 Wx _ — fli2=fcr'V2
(fy Ify Q.\ J
Consequently, each of the equations in (8) can be split into the
following:
<Px+*i(x, y)<py= 0, Wx+h(x, y)v>y = 0 (9)
where
a \2 - l /] V 2 an~\~X V2
h ( x , y) - ?i(x, y) = ( 10)
mi mi
The two equations in (9) are equivalent to
dy dy
S i = i , ( x ’ y)• no
P= y), <* = ^ 2 J~ = v ( x , t )
so that (13) and, consequently, (la) will assume the required canon
ical form (in the new variables)
upp+iiaa+ F 2(iip, ua, u, p, a) = 0
3. If (la) is parabolic in D, there exist in D functions cp(x, y)
and y)(x, y) such that the transformation given by (3) reduces
Equation (la) to the canonical form
u ^ + F^Ui, un, u, f , >]) = 0 (14)
The procedure for finding the functions 9?(.v, _y) and ip(x, y)
can be summarised as follows. To begin with, we must find the
function <p(x, y) which will ensure that the coefficient a n in the
transformed equation will vanish, i.e. it will be a solution of the
equation
an cpl+2a{2cpxcpy+ a 22(p2y = 0 (15)
As in the case of hyperbolic equations we are assuming that an
is not identically zero in any part of D. We then solve (17) for
fxhfy, but in contrast to the hyperbolic case (see (9)), we obtain
only one equation, namely,
<Px+Kx,y)rPy = 0 (16)
where
Mathematical Physics
Any solution of (16) which is not identically equal to a constant
can be taken for <y(x, y). The coefficient a 12 in the transformed
equation will then also vanish, which follows from the condition
that (la) is a parabolic equation, and from the identity given by
(6). For y>(x, y) we can take any twice continuously differentiable
function which does not lead to a22 being zero. Dividing the equation
transformed in this way by a22, we obtain the required canonical
form. An equation of the parabolic type has only the single family
of characteristic curves
8
Linear Equations
Similarly, an elliptic equation will reduce to the form
Substituting v = —~ /J 2, ju = —y ) we obtain
( 21 )
, 1 1 / ,
“«i + T TS L _•(-(“< - u-i) = 0
9
Mathematical Physics
Fig. 1.1
Mpp^rHao 0
a
dx
i/Z
|/ x
(or —
\
L= —
\ /y
7=^)
| x /
The general integral of this equation is
/ * + Vy = c
We shall therefore substitute
£ = 1Sx + \/y
and >i can be set equal to any function <p(x, y) which does not lead
to «22 being zero in the transformed equation. Let
'/ = I -v
10
Linear Equations
The canonical form of the equation is then
h ,,— - = 0
7]
PROBLEMS
1. Reduce the following equations to canonical form:
(a) = 0
x z u x x — y 2Uyy
string (Fig. 2.1). We shall also assume that the vibrations are small,
i.e. the square of the displacement can be neglected, and that they
take place in the (x, u) plane.
The assumption that the vibrations are small ensures that the
tension T in the string is independent of time t. To show this,
consider a segment (xl3x2) °f the undisturbed string. Its initial
length is x 2—x 1, whereas the length at time t is
Xl *2
^ ]/1 T ul dx « jj 1 •dx = x2—Xi
X\ x t
It follows that, to within terms of the second order in ux, the length
13
Mathematical Physics
of a fixed segment of the string is independent of time, i.e. the
segment is not extended. Hence, in view of Hooke’s law, it follows
that the tension T is also independent of time (to within terms
of the second order in ux). Consequently, T can be a function
of only x :
T = T(x)
Since we are considering the transverse vibrations of the string,
we shall only be interested in the component of the tension along
the it axis, which will be denoted by Tu. It is clear that
Xi
12 >Z*2
= J [7’(x2)wv(x2, t ) — 7’(x1)w,(xi , t)] d r + jj $ / ( I , t) d f dr (1)
h <i xi
This is the integral form of the equation for small transverse vibra
tions of a segment of a string.
If w(x, t ) has continuous second-order derivatives and T(x)
has a continuous first-order derivative, then using the mean value
theorem for the integrals in (1) we obtain
14
Partial Differential Equations in Physics
where £l512>£3 e [*u -v2], Ti, rz> ri e [*u C]- Dividing both sides of
(2) by V2r V2.x and proceeding to the limit as V2r -> 0 and V2x' -> 0,
we obtain the following differential equation for the small transverse
vibrations of a string
If we now apply Newton’s second law to the segment (a,, x2) of the
rod in a time interval V1t = (2- ( i ,w e obtain
*2
*1
n 12x2
= $ {^(.v2)A:(.X2) Ux(x2, t ) —S (* i)£ (* i) ux(x ,, r)} d r + J $/(£, r ) d |d r
n n*i
where /(x , /) is the density of the external force acting on a cross-
section at a distance x from the origin. This is the integral form of the
equation for the small longitudinal vibrations of the rod. Assuming
that the function u(x, t) has continuous second-order derivatives and
that k(x) and iS(x) have continuous first-order derivatives, we can
readily show that the differential equation for the small longitudinal
vibrations of a rod is
When >S(x), k(x) and p(x) are constants. Equation (5) assumes the
form
a2i'xx+F {x, t) = u„
where
«2 = f f o = - /( .v . o
Equations (3) and (5) are identical in form and differ only in the
notation (Sk instead of T, and pS instead of p). They are hyperbolic
at all points since, by definition, T(x), S(x) and k(x) are positive.
cos a .^ cosy = 1
- 1
V1 +
Consequently, cos a « 1 and hence Tc « T.
The second consequence of the fact that the vibrations are small
is that the tension T is independent of the time t. To show this,
consider a part S of the undisturbed membrane. This area is equal
to dx d^. At time t this area is equal to
s
B, B2
yt I
I
Fig. 2.2 1
x, X
y2
= $ [T{X2 , v)- T{x!, .)■)] d r = 0 (6)
y\
Xi
Since the intervals (x,, x 2) and (j^, y 2) are arbitrary, it follows from
(6) and (7) that 7’(x1, j ’) = T(x2, y) and 7’(.v,>’1) = T(x, y2), and
so on.
Let S be a segment of the membrane at time t, bounded by the
contour C (Fig. 2.3). Let S x and C1 be the projections of S and C
on to the (x, >0 plane.18
18
Partial Differential Equations in Physics
will cut the (x, y) plane along the normal n to C1 at the point M x
(Fig. 2.3). Fig. 2.4 shows the curve L along which this plane cuts
the surface S. It is clear that
tana
Tu — T sin a T
|/ l + t a n 2a
Consequently,
r^dU ,, \ x du d /l
Pu $r„d/ = T-^— dl
on £ dn cos/?
where /? is the angle between the elements dl and d/t .
Since /? < y (as in the condition a < y discussed above)
cos /? ^ cos y = 1/ j/ l +m2+iPy w 1. Therefore,
P “ = T lC1 T n i l '
^ jj pn ds
s’
where n is a unit vector in the direction of the inward normal to S.
Gauss’s theorem yields
^ pn ds = —jj jj Vp dr
s ‘ b'
where Vp is the pressure gradient. In the absence of external forces,
the equation of motion for the volume D can be written in the form
20
Partial Differential Equations in Physics
p~+Vp = 0 (9)
sss — +div(pz;) dr = 0
-^y-+div(pa) = 0 CIO)
o= p = p 0(l + a) (12)
Po
and confine our attention to small vibrations for which we can
neglect second order terms involving the ratio a, the velocity, and
the gradients of velocity and pressure, then equations (9) and (11)
21
Mathematical Physics
can be considerably simplified (this is referred to as linearisation).
In fact, subject to these assumptions, we have
1
~~ j 7 ~~ ~~~ 0 O~\~0Z ...) "(1 o')
P Po 1-1-0' Po Po
Po(l+o)y ss Po(l+yo) (13)
v (x , y, z, t) = v(x, y, z, 0) —a2 Vo dr
6
u = . f ( x , y , z)H n2 \<7dr
o
u, = a2a
a» = a2at
1 dH
curl E =
c dt
div E = 0
(18)
div H = 0
1 dE
curl H =
c dt
23
Mathematical Physics
dn
where k is the thermal conductivity. The thermal conductivity may
be a function of temperature, position and time
k = k(u, M, t)
Consider a part D of the body bounded by the surface S. Let f { M , t)
be the density of heat sources and let us evaluate the heat balance
for D in a small time interval V2U
2 , = 5 J J / ( M , r)d rV 2/
D
s
where Qy represents the heat gain from the sources and Qz the heat
lost from D. The derivative dujdn is evaluated in the direction of the
outward normal to the surface S. Finally,
Qz = {jjjjjcpu, dxV2t
D
24
Partial Differential Equations in Physics
Application of Green’s theorem to the first integral yields
For the two- and three-dimensional cases the above three types of
boundary condition are of the form
u 's = ft (M , t) (type I) (24)
-k^~ =
dn s
we have the boundary condition of type III.
Other types of boundary condition will be discussed later. The
boundary conditions considered above are linear because the
unknown function or its derivatives enter into them in linear forms.
They are called homogeneous if the right-hand sides {pi, v, fi) are
identically zero, and inhomogeneous in all other cases. It is evident
that similar boundary conditions are encountered in problems
leading to elliptic equations. Physical interpretation of each presents
no special difficulties.
We shall now formulate the corresponding three types of boun
dary condition for equations of the form
div(fcVw) —qu+ f{M , 0 = puu (27)
and
div(fcVu) —qu-rf(M, t) = pu, (28)
where k, q and p are functions of the coordinates of M. All the
equations considered above belong to this type with q = 0.
| HM)
PROBLEMS
*0
Substituting these functions into (6), we obtain d’Alembert’s formula
^ <r(x-aiyv<p(x-r at) , 1 f , ^
u(A, t ) = ------------ 2 J ? '0 ) dz (8)
x — at
3* 35
Mathematical Physics
lead to
|Ui(-Y, t ) - U z(x, t) | < £, -CO < .V < CO, t < /,
+ - y [(pi(x+at)-<p2(x+at)]
x+ a t
+ 2~ W i ( z ) - y 2(z)] dz
x — at
and, consequently,
+ -y l <Pi{x+at)-cp2{x+at)\
x+ at
+ \
x — at
Idz
x+ at
1 1 I f
< T 4 + T 4 + -£r 1
x — at
idz
= <5+<5t < 5 (1 + 0
If we let 5 = e/(\ + b), the inequality
| itX(x, t)—u2(x, t) | < e
will be satisfied for all oo < .v < oo, t ^ /,, which was to be proved.
B 6
I <Pn(x)—<pn+k(x) | < - and | y>„+k(x)—ipn(x) \ <
for all —00 < jc < 00. In view of the theorem proved above, we
have the following inequalities for all t < tx and —00 < x <00
These hold for any n > N and any positive integral k. However,
this means that the sequence of solutions {u„(x, 0} converges
uniformly in the above range of xc, t to some function u(x, t). This
function is called the generalised solution of Cauchy’s problem
(2)-(3). It is given by
x + al
x — at
37
Mathematical Physics
or
This function and its derivative ut(x, t) will assume the prescribed
values cp(x) and y(.v). It follows that in the above case, the d’Alem
bert formula will also yield the (generalised) solution of Cauchy’s
problem. The problem can also be solved in another way by using
the generalised functions and their convolution (see Appendix,
Section A.l).
Definition The fundamental solution G(x, t) of the wave equation
a2uxx = u,t is defined as the solution of the Cauchy problem
a2Gxx = G,t , G(x, 0) = 0, Gt(x, 0) = <5(x)
It is readily verified by direct substitution that
G(x, t) = — [y]{x+at)-r]{x—.at)\
function. In point of
fact (see Appendix, Section A.l),
Gx(x, t) = — [b{x+at)—b{x~at)\
Gxx(x, t) = — [b'(x+at)—b'{x—at))
Consequently,
a2Gxx(x, t) = G tt(x, t)
G( x , 0 ) = 2a- h ( x ) - V(x)] = 0
38
The Method o f Characteristics
The solution of Cauchy’s problem
a2vxx = vtt, v(x, 0) = 0, v,(x, 0) = y>(x)
will be written in the form of the convolution
v(x, t) = G (x, t ) * y(x) (9)
In fact, by evaluating the derivatives of the convolution (see
Appendix, Section A.l), we obtain
vxx = Gxx * V, v,t = G„ * rp
and, consequently,
a2vxx—v„ = (a2Gxx- G tt) * tp = Q*ip = 0
v (x, 0) — G(x, 0) * ip (.x) = 0 * ip (.x) = 0
v t(x, 0) = Gt(x, 0) *y>(x) = d(x) * y ( x ) = y(x)
The convolution G * ip can also be written in the form
oo Qt
° 2R XX = R ,t
Next,
w (.v, 0) = R,(x, 0) = cp(x)
wt(x, 0) = Rtt(x, 0) = G„(x, 0) * <p(x) = 0 * y(x) = 0
If the function cp(z) is continuous, then w(x, t) can be written
in the form
x + at
v ,{ x,t ) = R,(X, t ) = + l ^ \
I X -Q t ]
<p(x—at)+<f (x + a t)
= 2
i.e.
, . cp(x— a t)+ < p { x + a t)
w ( x , t) = --------------2 -------------
or
x+at
/ .x tyix—at)+cp(x+at) 1 f / , ,
« < * , / ) = -------------2------------ + 2a ) w(z)dZ
x — at
zv(M,t) = \ w ( M , t , T ) d r (17)
o
In fact,
Consequently,
pwt t = f { M , t ) + \ p W n dr (19)
o
Let us evaluate the quantity div (kVw)—qw. The operation
div (A:V) can, of course, be carried out under the integral sign.
We have
3.3.2 Let us now use the above method to solve Equation (1).
It is required to solve the Cauchy problem
a ^ x x + f i x , t) = utt
u(x, 0) = <p(x), ut(x, 0) = y>(x)
a2vxx =
fl(-v, o) = <7(a), V,(x, 0) = y>(x)
42
The Method o f Characteristics
zv(x, t) = ^ W(x, t, t) dr
o
where W(x, t, r) is a solution of Cauchy’s problem
a2Wxx = Wtt,
w \ t=x= 0 , w t\ =z= f ( x , T)
and, consequently, can be written in accordance with the d’Alembert
formula
x X a { t — x)
W(x, f, t) = - ^ 5 / 0 »T) dz
x — a t t — x)
Therefore,
t x - \ - a ( t — t)
0 = 2^ 5 S f{z,x)dzdx
0 x -a (t-x )
— at
The function
x+at
44
The Method o f Characteristics
will then be a solution of the boundary-value problem. In fact, it
satisfies the homogeneous wave equation because it is a super
position of forward and reverse waves. It satisfies the boundary
condition in view of Lemma 1. Let us verify the initial conditions:
„(Xj0) = ^ w + ^ j Vi(z)d2
= viW = v W . -v > o
It follows that the initial conditions are also satisfied.
The boundary-value problem
Q Ux x ■ Ut t ,
Fig. 3.3
problems defined by (21) and (22) are prescribed only on the semi
infinite straight line for which x > 0. However, from the boundary
condition for the problem defined by (21) it follows that
0 ( _ z) = - F(Z)
Consequently, 0 ( —xy) = —i 7(x1). It follows that instead of the
forward wave arriving from the point —Xj we can consider the
reverse wave which has left the symmetric point x: at time t = 0.
This reverse wave will reach the point x = 0 in a time t1. From
t = t x onwards, it must be replaced by the forward wave which
has left the point x = 0 at time t = t x and has propagated with
displacement —0 ( —xi). It follows that when the boundary con
dition »(0, /) = 0 is satisfied at x = 0 we have the phenomenon
of reflection in which the magnitude of the displacement is conserved
but its sign changes.
It can be shown in a similar way that when the boundary con
dition ux(Q, t) — 0 is satisfied at x = 0, we have the phenomenon
of reflection with the conservation of the magnitude and of the
sign of the displacement.
The method of characteristics can also be used to construct
solutions of homogeneous boundary-value problems on a finite
segment subject to boundary conditions of types I and II. To be
specific, consider the first boundary value problem
a2uxx = utt
u(x, 0) = <p(x), wf(x, 0) = y>(x), 0 < x < / (23)
u(0, t) = 0> r/(/, t) = 0, t> 0
46
The Method o f Characteristics
To obtain the solution let us continue the initial functions <p(x)
and y(x) along the entire line and assume that they are odd with
respect to x = 0 and x = /. Let 9?2(x) and V>z(.x) represent the
functions continued in this way. The function
x+ at
47
Mathematical Physics
We already know how to solve the problem for v(x, t). The
propagation of the boundary value will now be considered. Since
the only reason for the appearance of a perturbation is the situation
at the boundaries, we shall seek the solution in the form of the
forward wave
zv(x, t) = 0>{x—at)
From the initial conditions we find that
w(x, 0) - 0(x) = 0, x > 0
It is evident that wt(x, 0) = —a0'(x) = 0 will also be satisfied
for x > 0. From the boundary condition we find
0 ( —at) = I> 0
and therefore
0 z > 0
z < 0
«(G 0 = 4 ^ 2 - (26)
SM
where P denotes the variable integration point.
Let da) be an element of solid angle subtended by the area
do- at M so that dcr = r2 dco. It follows that we can also write
Lemma
V2m =V2(m)
where on the left-hand side the Laplacian V2« is evaluated with
respect to the coordinates of M, and the right-hand side is evaluated
with respect to the coordinates of point r. Henceforth, we shall
omit the subscript r.
4 49
Mathematical Physics
d 'm
(29)
ut(r, 0) = y(r) = ^ ^ 2 J J V(p ) da
Let
a-v„ = v„
v(r, 0) = <px(r), vt(r, 0) = y,(r)
t;(0, t) = 0
^ 0 - » ( f + a ) + * (f- a 0 + T J V2(z)dz
r— at
Consequently,
r+ ai
_ <Pi{r-\-at)+<p2{r—at) , 1
u(r, 0 = ~ ^ "v + J V2(z)l
r
r — at
+ {aiy(at)+aty)(—at)}
Since the functions y{z) and Tp{z) are even, whereas y'(z) is odd,
it follows that
Using (28), (29) and (30), we obtain Poisson’s formula for the
required solution of Cauchy’s problem (24)-(25):
1 d f f y(P)d<r i f f v C f ) da
u(M, t) = (31)
An a dt J J at Ana J J at
cat
«i/
Thus, if we assume the existence of the solution of Cauchy’s
problem for the three-dimensional space, we find that it must be
given by (31). Consequently, this is the only solution.
3.7.2 We can now solve Cauchy’s problem in three-dimensional
space for the inhomogeneous wave equation
a2V2u + f { M , t) = utt
u(M, 0) = <p{M), ut(M, 0) - rp(M)
We can split it into two subsidiary problems:
1. a2V2v — vtt
v (M, 0) = <p(M) , v,(M, 0) — ip (A/)
The solution of this problem is given by Poisson’s formula.
2. a2V2w-\-f(M, t) = wtt
w ( M , 0) = 0, wt( M , 0) = 0
It is evident that u = v+ w .
Problem 2 can be solved by the method described in Section 3.3.
In particular, we must first solve the auxiliary problem
a2V2W = Wtt
w \ , ^ = 0, n a = t= / ( M , t)
52
The Method o f Characteristics
and then from Poisson’s formula we have
= J W (M ,t, r)d r
o
or
at
r
do- dr
0 cr
where r is the distance between the integration point P and the point
M, i.e. r = rMP. This integral can, clearly, be written as an integral
over the region D% bounded by the sphere S$:
(32)
/ l/--)(5 (P ,A /o )
u - dv
r
53
Mathematical Physics
Using the fundamental property of the 8 function, and remembering
that for any continuous function <p(M)
0 if M0 $ D
S $ 5 ? ( p ) W M 0)dV= |
D (p(Mo) if Mo g D
we obtain
0, where at < rMMo
u(M, l) 1_______ 1 /./ r MM (33)
- , where at > rMMo
.4 Jta2 rMMoJ \ a
<
p(P\) d a
uppers^ a at cosy
54
The Method of Characteristics
where y is the angle between the normals to the (x, y ) plane and
to the sphere 5 $ at the point P. It is evident that
\P P u IPtMl2
cosy
\MP\ at
\ / ( a t f —( x - £ ) 2- ( y - r ] ) 2
at
where £, rj are the coordinates of the point P r which is the projection
of P or, to the (x, y) plane, and (x, y) are the coordinates of the
point of observation M. Therefore,
and, similarly,
+ f f y(£. V ) & d r )
(35)
Ina j J |/ a2t2— {x — £)2—(y — rj)2
SM
We are now in a position to solve Cauchy’s problem in two-
dimensional space and for the inhomogeneous equation. It can be
reduced to the problem just considered and to Cauchy’s problem
for the inhomogeneous equation with zero initial values. The latter
problem can be solved by the method described in Section 3.3.
We shall not repeat all the steps leading to this solution and shall
merely quote the final result
55
Mathematical Physics
by the curve S (Fig. 3.6). For t < t, the circle Y m does not contain
points belonging to D, and therefore the integrals in (35) are zero so
that u(x, y, t) = 0. For any t > tx the circle Y m contains the domain
D or part of it and, therefore, u(x, y, t) ^ 0 for these values of t.
PROBLEM S
( 2)
and
u{M, 0) = u, (M, 0) = <p,(A0
[respectively u(M , 0) = (3)
59
Mathematical Physics
60
Separation o f Variables
We shall not reproduce the proof of this theorem here (the one
dimensional case will be discussed in Chapter 9).61
61
Mathematical Physics
The Fourier series for the function f ( M ) in terms of the functions
00
{0n(M)} is defined as the series ^ Cn0 n(M), in which the coefficients
n= 1
C„ are given by
C„ = jj/(M )p(M ) 0„(M) dr
3 P ^ndr J
D
62
Separation o f Variables
Cn ~ \ p ( P ) < r ( P ) i ’,(P) i z r
Dn
II «!
M 2 = \ p ( P ) 0 2n(P) drP
D
( d0
In view of Green’s theorem, the last integral is equal to j /cw — dr,
and therefore 5
R[w, 0] = ^ k(7w- V 0 ) 6 t + jj qw0> d r — ^ k z v ^ —da ( 8)
d d s
If we are dealing with the first or second boundary-value problem,
f d0 ld0
we have \ k w ——d r = 0, so that w = 0 I —— = 01. In all
i dn s \ d% 63
63
Mathematical Physics
W
(M,O= rt=
X^"(')0n
-1 (M) (9)
where
y n(t) = j f ^ \ p ( P ) u ( p , t ) # , ( p ) A T F (io)
Using Equation (5) for 0„(M) we can transform the last formula
to the form
i?[w,0„]
^ (0 = U® nf
P U,t d r
64
Separation o f Variables
and hence, if we compare this result with ( 10) we have
11 n11 D
a, = = 1 [p<,P)r i (P)<fn(.P) dr
l /, I 4 ll«'„lr ^
which was to be established.
By assuming the existence of the solution of (l)-(3), we have
come to the conclusion that it can be represented by the series (7)
and, therefore, it is a unique solution.
4.1.3 Let us now consider the following properties of eigenfunc
tions and eigenvalues.
Property 1 If 0 is an eigenfunction corresponding to an eigenvalue
A, then C 0 (where C is a constant) is also an eigenfunction corre
sponding to the same eigenvalue.
Property 2 If and <f>2 are eigenfunctions corresponding to the
eigenvalue A, any linear combination such as C1$ 1+ C2$ 2 is also
an eigenfunction corresponding to the same eigenvalue A. The
validity of these statements is obvious.
Property 3 Eigenfunctions (P, and <f>2 corresponding to different
eigenvalues A, and A2 (Ax ^ A2) are orthogonal in the domain D
with a weight p(M), i.e.
] p ( P ) 0 , ( P ) 0 2(P) dr = 0
D 65
5 65
Mathematical Physics
- R [ ® 2, , 0 2] = (X2- X y ) \ p 0 i0 2dx
D
jj p 0 \ dX-\-By ^p0y02 dr = 0
D D
\ p 0 y 0 2 d x f B 32 p 0 l d T + B 33 ^ p 0 20 3dx = 0
D D D
and so on.
66
Separation o f Variables
Continuing this orthogonalisation process, we can consider r
eigenfunctions @l} 0 2, ■ d>r corresponding to the given eigenvalue
A and these will now be orthogonal in pairs. Assuming further that
this orthogonalisation process has been performed, we conclude
that any two linearly independent eigenfunctions of the bound
ary value problem (5)-(6) are orthogonal in the domain D with
a weight p.
L [ 0 1—i 0 2]+(ot—iP)p(0l—i 0 2) = 0
Therefore, A = a —ifi and 0 = 0 x—i 0 2 are the eigenvalue and
eigenfunction of the problem, respectively. Using Property 3 we
have
which is impossible.
\ 0 nL[*n] d r + A ^ p ^ d r ^ O
D D
and hence
ll^ ll2
Since R[0„, 0 „] ^ 0, we have A„ ^ 0.
5* 67
Mathematical Physics
Remark For the first and third boundary-value problems all the
eigenvalues are positive. For the second eigenvalue problem with
q(M) = 0, we have A = 0 as the eigenvalue and 0 = 1 as the corre
sponding eigenfunction.
J. = — ( » = 1 ,2 ,3 ,...)
„ n2 4.x2 n2n269
69
Mathematical Physics
are the eigenfunctions. For each we can find the following eigen
functions
Vn{t)= C„e-“2^ {n = 0, 1, 2, ...)
The required solution of the problem will, in view of Theorem 2
(Section 4.1.2) be the function
CO
u(x, t) = Cn e - ^ ^ c o s — x
n—0
where
/ t
C0 = I ^ , C„ = - J ^ ) c o S y f d f ( n = 1 ,2 ,3 ...)
6 o
i
ll^oli2 = /, I M 2 = jjcos2— f df = y (« = 1, 2, ...)
0
71
Mathematical Physics
The factor B {\///,) will be included in yIJ(t). Substituting (22) into
(20), we obtain the equation for the eigenvalues:
cos /i
1 w 2\
7(/?i+/i2) t* I
where ju = ]/Xl. ptx, jx2, ... are the positive roots of this
equation. The eigenvalues will then be the numbers
Xn Pi
0„(x) = y c o s y ^ + ^ s i n - J-x
W{X, 0 = 2 C nZ -^'^nix)
n= 1
The coefficients C„ can be found from the initial condition using the
orthogonality of the eigenfunctions 0 „(x):
/
u{\, 0
)=cp(x), ut(x, 0
) - <Pl (x)
w(0,r) = 0, Tux(l,t) -- m0utt(l, t)
Among the class of functions 0(x)W(t) we shall seek solutions which
satisfy only the boundary conditions. Separating the variables we
have
W "+ a22V = 0 (23)
&"+?.& = 0 (24)
0(0) = 0 (25)
The boundary condition on the right can be written in the form
T 0 ,(t)'P(t)-mo0(I)'P"(i) = 0
Replacing W"(t) in (23) by 1F(t) and dividing both sides of the
equation by ^ ( 0 , we obtain
0'(t)+hX0(I) = 0 (26)
where h = o2m0IT.
The solutions of (24) which satisfy (25) are of the form 0(x)
= sin \ flx . From (26) we find the equation for the eigenvalues
K > 0:
-1
tan fi = (A — } / 2 I
IjT ’
The corresponding eigenfunctions are
= siny *
i
Dn = ---- \L~ii2 \\P0<P1 (0 0 n(f) d^+ Wo-PlCO (O)
aPnWn\\ L> >
l |0 J 2 = P o $ ^ (£ ) d l+ m 0^ 2(/)
o
Let us now return to the properties of the eigenvalues and eigen
functions. We note, to begin with, that since for the functions 0
belonging to class A we have R[0,<P] ^ 0,
R[0 , 0]
inf = /i ^ o
<PeA im i2
= ~ \ { f L [ 0 ] + 0 L [ f ] } d r - 2 p \ P 0 f dr
D D
= - 2 \ f { L [ 0 ] + p p 0 ) dr
D
\ f { L [ 0 ] + p p 0 } dr = 0
D
For the case qdM) > qi(M) the proof is practically the same as
before.
The function 0 " which is equal to 0 ' in D' and zero in D " —D'
(shaded region) will clearly belong to class A " . If we perform this
operation for each of the functions in A', we obtain the new class
77
Mathematical Physics
and
^ p 0 2dr = ^ p 0 2 dr
D" D'
— [k0']~ q 0 + X p 0 = 0
o
u(x, 0, 0 = 3 (33)
>1
II
= <p(x, y). ut( x , y , 0) = tpi{x,y) (34)
Among functions of the form <P(x, y ) IF(r) we shall seek solu
tions of (31) which satisfy only the boundary conditions (32), (33).
Substituting this function into (31), (32) and (33) and separating
the variables, we obtain the following Sturm-Liouville problem
V2<£+A0 = 0 (35)
0 (O ,y)= 0 (l,y) = O (36)
0 (x,O ) = 0 ( x , l ) = O (37)
This problem can also be solved by the method of separation of
variables. We shall seek solutions of the form &(x, y ) = A(x) B(y).
Substituting this function into (35) and separating the variables
we obtain
A= 0
A " + p A = 0 (38)
B"+ (X—p)B = 0 where B"+xB = 0 (39)
From (36) and (37) we find that
A (0) = A (I) — 0 (40)
B ( 0) = B (I) = 0 (41)
79
Mathematical Physics
A„(x ) = sin -y -x
and
n2k 2
ak p ( * = 1 ,2 ,...)
_ . . . nk
Bk(y) = sin — y
The eigenvalues K m and K,n will clearly be equal and the cor
responding eigenfunctions
. .n n nk , , . nk . nn
<P„fk = sin — x sin—y y and n ~ sin —-j—x sin —-—
80
Separation o f Variables
oo
i i
(42)
I
where Ck are the Fourier coefficients of the expansion of f{M )
in terms of the functions {0fc}.
A sufficient test for the completeness of the system {#„} If a function
F(M) is continuous in D, and for any e > 0 there exists a linear
combination Sn = ...+a„tf>n for which $ p(F—S„)2 dr < e,
then the set {&„} is complete. D
Consider a fixed e > 0. For any function f ( M ) which is square
integrable in D, we can find a function cp(M) which is continuous
in D and is such that \p(J—?>)2dr < e/4. This requires proof but
D
we shall not give it here. For the function <p(M) and for the
chosen e, it is possible by hypothesis to find a linear combination
S n = ot.l0 1Jr ...Jr oi.n<hn, for which
$ P (<P-Sn)2 d t <
D ^
\p (f- ^ c k®k)2 dr
D
6 81
Mathematical Physics
in which Ck0 k is the partial sum of the Fourier series for /(A/).
k~l
It is clear that
n 2 n n
$ „ (/-£ a * , d t = \ p f 2i r - 2 Y J Ck J p / ^ d r H - ^ Q I I ^ I i 2
D X k=\ ' D k=\ D k —l
n 2 n
\p lf~ Y j dT = $ p / 2d r - ^ a | i 0 t ||2
D ' k^l ' i) k =\
n n
2 dr = ^ a | | 0 t ||2
D k=I
\f{ M ) dr = ] T c n 5 $„(M) dr
D' n=\ D'
n
Proof Consider the difference d„ = jj / d r — ^ Ck jj 0 k dr
D' k^l D’
n
^ = !!!/-E c^ )dT
D' k= 1
n n
< b'Sf - E
k^i
dr < [ \ f - Y JCk0 k I dr
D k =\
dr
\ f ~ ^ C k0 k dr Ck&k
> 1 VP
D k= 1
<
$p/2d r - ^ Cf\d>k\\2
D k= I
< 2 \ p { f - g f d r + 2 \ p [ g - ^ C fc0fc)2dr
D D k= 1
^ p ( / —^ C k^ k) dr < dr < £
D k= 1 D
= (48)
n= l
where
Using Equation (5) for p@n under the integral sign in (49), we
obtain
*P > „,«] l f
L[ m] dr
1
^ (0 = ^ sin 10:n ( t - 0 ) f n(0) d0
o
r ' ‘In + H = °
We shall seek the solution in the form of the sum
u(M, t) = v (M , t ) Jr w { M , t)
86
Separation o f Variables
(^ H r °
where
f i ( M , t) = f ( M , 0+£[*>i]-/Wim
<P(M) = <p(M)—v1(M, 0), , 0)
We have already considered this problem in Section 4.3.1. The
form of the function v x(M, t) is either guessed or it is found by
Duhamel’s method (see Section 4.3.3).
87
Mathematical Physics
Example 6 It is required to solve the problem
a2uxx = utt
u(x, 0) = rp^x) , ut(x, 0) = cp2(x)
w(0, 0 = /«i(0, i<(/, 0 = M O
We shall assume that the functions p.x(t) and ju2(t) are twice dif
ferentiable and will take
V {( X , t) = — y —- j« i ( 0 + ' y / u 2 ( 0
zv(0, t ) = w(l, t) = 0
The function w(x, t) will also be sought in the form of the sum
w = R(x, t)+Q(x, t), where R(x, t) is the solution of the homo
geneous boundary-value problem
a2Rxx = Ra
R ( x y 0) = cpi (.v), = q>2{x), R ( 0, 0 = K(0 0 = 0
and is given by (see Example 1)
co
<?Qxx+f(X, t) = Qu
Q(x, 0) = <2((*, 0) = 0
Q(P,t) = Q ( l , t ) =0
w h e re/O , 0 = - y ^ z 'O ) -
According to Section 4.3.1,
00
2 0 . 0 = ^ ¥„(/) sin— x
n=i
The functions are given by
r
*„( 0 sin]/A„(t-d)f,(6)d8
where
/,(<>) = t o d f = — K - o v r o o - ^ 'w ]
o
2 ^ 1 + P,l2 ’ ~ 2a2/+/32/2
II
A
90
Separation o f Variables
L [ u ] = p ( M ) \ U,t (59)
\u,
(a4 nJ r ^ U) s = ^ (6l)
This involves the following steps.
1. First we solve the problem defined by (59)-(61) with a sta
tionary inhomogeneity on the boundary condition, i.e. with a
boundary condition of the form
| a— , t)
w!t=t = ut\t=z= o
- t —x, x)i](t—x)\ dr
ot
From this expression and from the formulae for the derivative
t
= ^— r)v2(t—r)} dr
0
t
- ^ — {p(M, r)r](t—T)} dr
0
t
= ^ (M , t) — 7](t-r)dr
0
t
= ^ p (M , r)d(t—r) dr = p ( M , f)
o
Let us substitute (62a) into (59) and use Equations (63) and (64).
For internal points in D we have, in view of (62),
and using the continuity of wrr in the region (M e Z), t > 0), we
find that (for t —r -> 0)
L[w(M , 0, t)] = p(M)wtt{M , 0 ,t) = 0
since w(M, 0, t) = 0 for internal points of D and, consequently,
L[zo(M, 0, 0] = 0.93
93
Mathematical Physics
Therefore,
un(M, t) = \ t t , t) dr
6
so that
t
L[u) -putt = jj d -{L[w\~pwt,} dr = 0
o
since L[w\ = pwtt. This means that (62a) is the solution of (59)-(61).
Consider the special case /u(M, t) = Q(t). Proceeding by analogy
with the above discussion, the solution can be obtained as follows.
1. We first solve Equation (59) subject to the boundary condition
(“In + ^ U) s =
where t is a fixed number.
2. The function Q(x) R(M, t —r) >?(/ —t) will be the solution
of Equation (59) subject to the boundary and initial conditions of
the form
(« +P“j s = 2 ( t )
H't = x = Ui\, = x = 0
We note that in view of the initial conditions, for all internal points
of the domain D we have
R(M, 0) ee Rt(M, 0) = 0
3. The function
Q(t)[R(M, t - T ) r / ( t - z ) - R ( M , t - T - d : ) j j ( l - r - d r ) ]
94
Separation o f Variables
»r =T= Wf!t =t = o
4. The solution of the original boundary-value problem will
be the function
U(M, 0 = ^ Q ( r ) ~ R ( M , t - t) dr
o
This can be verified by direct substitution as in the previous case.
It is thus sufficient, in this case, to find the solution R(M, t) for
the problem with the very simple (stationary) inhomogeneity on
the boundary condition Q(t) = 1.
95
Mathematical Physics
Therefore, CO
or
§ Q (r)-^ R (x , t - r ) d r + Q ( t ) R ( x , 0 )
o
The function R(x, 0) is zero at all internal points within the range
[0, /] and for x = 0, while for x = / we have R(x, 0) = 1.
If it is required to solve the problem
a2uxx = ut , u(x, 0) = 0, u(0, /) = £?i0), u(l,t) = Qz{t)
the solution can be sought in the form of the sum u = v-\-w, where
for?; and w we have the following two problems:
v: a2vxx = vt, v(x, 0) = 0, v(0, t) = v(l, 0 = 0
zu: a2wxx = w,, w (x, 0) = 0, w(0, t) = 0, zo(l, t) — Q2(0
Each of these problems can be solved by Duhamel’s method,
as shown in the above example. This method can also be used
for the solution of boundary-value problems on a semi-infinite
straight line.
t)=pu„
j x (kvx) = p v tt (65)
P: k { l ) v \ l , t ) + JPi
± . k ( 0)^2(0, 0
2a 2aj
98
Separation o f Variables
Integrating this relationship between 0 and t, we find that
o)]
— [kux] = p u t (70)
where a > 0 and a < (MD—Mr )/2T. The function ^(.v, t) is con
tinuous in D and, consequently, it reaches its maximum value in
D at some point (x1; tj). It is evident that v(x 1, t 1) ' ^ M D since
zj(a, /) > u(x, t) in D.
The point (xq, /,) cannot lie on any of the three boundaries
/ = 0, x = 0, x = /. In fact,
j ^ ( k u x) + f ( x , t) =pu,
j ^ ( k u x) + f ( x , t) = pu,
PROBLEMS
a2vxx = vt (3)
v(x, 0) = 0 (4)
and is continuous in Dlt where \v(x, t ) \ ^ 2 M throughout D{.
107
Mathematical Physics
AM l b 2 \ AM
\ v ( ± b , t ) \ < 2M < ® ( ± O /) -
108
Green’s Functions for Parabolic Equations
Let us derive this function. Consider, to begin with, the following
special Cauchy problem
a2uxx = u, (7)
(0, .v < 0
u (x, 0) = tp (.y) = v (x) = x>0 (8)
This function satisfies the first condition in (10). From the second
condition we obtain a relationship for C:
GO
|/4aM
u(x, t ) = f z~yl dy
(t t H ; S
109
Mathematical Physics
0 j /4 a lt
w (x,/) = ~ { [ dy
\/n J \n J
or
u(x, 0 = y 1+ 0
\ l/4a2t
Q
u (.v, 0) = iv ( x - x l)~ r]( x - x 2)}
CP ( x 2 -Yj)
110
Green’s Functions for Parabolic Equations
—
cp
G { x S u 0+—
cp
G (x -h , 0
^ Uxx = Mt
u(x, 0) = d(x—x 0) —6 (x + x 0)
or, correspondingly,
G**(x, x0, t) = G (x —x0, 0 + G ( a '+ a 0, t)
( * - * o ) 2~
1 4a2/
e 4a2‘ + e
\/ 4na2t
d^ G ( x - £ , t) = <p{£)G{x-£, 1) d£
cp
This also applies to any other element of length d | along the line
t = 0. In view of the remark in Section 5.2, it is natural to suppose
that the temperature due to the action of all such elements, i.e. the
temperature due to the initial temperature distribution u(x, 0) = y(x)
is given by
oo
U ( x , t ) = $ <r(£)G(x-S,t)d£ (15)
— CO
If this is correct, the function given by (15) will, in fact, be the solu
tion of Cauchy’s problem (13)~(14). To verify this statement, it is
sufficient to show that this function satisfies Equation (13) for all a
in the range — oo < x < oo and t > 0 and the initial condition
given by (14).12
112
Green’s Functions for Parabolic Equations
8 113
Mathematical Physics
Remark It follows from (16) that if the initial values <fi(x) and
<p2(x) differ by less than e, i.e. \<P\(x)—f 2(x)\ < E f°r ah then the
corresponding solutions of Cauchy’s problem, ux(x, l ) and u2(x, t),
will also differ from each other by less than e, i.e.
\ 9>(£)G,d£ = - J — G M ,/) d |
9 ( * —£)2
— 00
\ — — 9?(.v+2act]/1 ) e a*da
^ 2 I/ n t
00 *
Remark 1 It follows from (15) that the heat propagates along the
rod instantaneously. In fact, suppose that the initial temperature
<p{x) is positive over a finite segment (xl5 x 2) of an infinite rod and
is zero outside the segment. The temperature at an arbitrary point x
is then given by
•*2
u ( x , t ) = 5 <p(£)G(x—£ , t ) d £
It is evident that this function is positive for any x and any small
t > 0. We have reached this conclusion because the physical assump
tions which we have used in the formulation of Cauchy’s problem
were not strictly accurate [for example in writing Equation (13)].
1 -J iL
G ( x , 0 = - / = , - e 4a2'
X Ana11
with the initial function y(x), i.e.
u(x, t ) = G(x, t) * <p(x)
If the initial function rp(x) in this formula is taken to be an arbitrary
finite generalised function, then u(x, t) will also be a solution of
Cauchy’s problem.
8* 115
Mathematical Physics
y 4a2t
u{xt t ) = - —
\n
[ e“®2d a-----i-
J 1/tt
f
t
e“®2 da
y 4a2t
0 ]/ 4a2, 0 —oc
= 5 ) - )/?t
* U +$
p4a2r
M|+ M2 , »2- i l i f / -V
| 74a2r
m2mx* = , !/(*, 0) = A e- *2
u ( x , t ) — A jj e~i2C(x—£, t) d£
116
Green’s Functions for Parabolic Equations
u ( x , t ) = - ^ jj Q - i x + l a a ^ t )2e - “2 d a
\ ' 71
f J00
— — 00
00
Ae~*2 — ^ g - t a a ^ l - ( 4 a 2f + l ) a 2 ^
]//7l
4 x 2a 7t oo I W T /-—r-r-
, _,2 ------- r- * - I — ' + 1/ 1 + 4 a 2< g
/I e l + 4 a 2< f \ j / i + 4 0 2r
Y 71 \ da
^ G ( x - £ , t - t) = / ( f , i ) G ( x - f , / - r ) d |d r
cp
In view of the remark in Section 5.2, it is natural to suppose that
the temperature distribution due to all such sources distributed
along the line in the time interval between 0 and t will be given by
t OO
w{x, t) = 5 5 /(£> r)G(x—£, t —t) d£ dr (22)
0 —oo
This formula is, in fact, the solution of the problem (19), (21). The
condition given by (21) will, of course, be satisfied. To verify that
the function given by (22) does, in fact, satisfy Equation (19) can
be confirmed by analogy with the analysis given in connection with
Equation (15).
If the heat source acts only at the point £o, but is a function of
time, the function f(x, t) in (19) is of the form
f ( x , t) = f ( t ) d ( x —£0)
The solution of the problem (19), (21) is then of the form
t 00
w(x, t) = $ d(£—£o)f(j)G(x—£, t —T) dc dr
0 —00
/
= S./It)<j(.y —fo, t ~ r) dT
0
where we have used the properties of the (5-function.
118
Green's Functions for Parabolic Equations
Consider now the Cauchy problem for the equation of heat conduc
tion in two- and three-dimensional space. To begin with, consider
the homogeneous equation
a2 V2u = u, (23)
— 00 — 00 — 00
0 - 00—00—00
(26)
120
Green’s Functions for Parabolic Equations
— . 2 . n n . Tin . Tiant
G(x, x0, t) = -j- y sin ~j~x sin— x0sin ——
n= 1
The situation is quite similar for the second and third boundary-
value problems. However, this method is not ideal for bounded
regions and is not usually employed.
PROBLEMS
All the results which are derived in this chapter can be deduced
from a small number of formulae and relationships. These will
first be derived.
6.1.1 Suppose that the functions u(M) and v{M) have the following
properties:
1. These functions and their first-order derivatives are continuous
throughout a closed domain D, bounded by the surface S, except,
possibly, for a finite number of points.
2. The functions and their first-order derivatives are integrable
in D.
3. The functions have second-order partial derivatives which are
integrable in D.
Subject to these conditions, we have
R [ u , v ] = — J vL[u] dr
D
D D s
124
Green's Functions for Elliptical Equations
(3)
5 ^ d,i=W )dr
For a doubly-connected region D', bounded by two concentric
spheres SR and S Ri(Rx < R), centred on the point M0, Equation (2)
can be written in the form
\^{vL[u\ — uL[v]} dr -- ^ d—
-U—v 1\ dcr
,
dr
d' sR \
du ,
k \ v ----- u ——1dcr (4)
dr dr
125
Mathematical Physics
6.1.2 Mean value theorem If the function u(M) is harmonic in
a spherical domain DR and is continuous, together with its first-
order partial derivatives, in DR = DR-\-SR, then its value at the
centre M 0 of DR is equal to the arithmetic mean of its values on
the sphere S R, i.e.
u (M °) = $ u^ d<7
Sr
^2 \ u (M ) da = 4nR]u(M% M* e SRl
Sr 1
If we now go to the limit i!, -+ 0 we obtain Equation (6).
We have considered harmonic functions in three-dimensional
space. In two-dimensional space (plane) the mean value theorem
is given by
v (M ) < H S+ H^ Hs
Proof Suppose that the two functions w, and u2 are the solutions
of this problem. Their difference u = ul—u2 is a function which
is harmonic in D, continuous in D and zero on S. In view of the
above theorem, the maximum and minimum values of u are zero
and, consequently, u = ux—u2 = 0 throughout D.
Theorem If a sequence of functions u1, u2, ..., u„, ..., which are
continuous in a closed and bounded domain D and are harmonic
in D, converges uniformly on the boundary of the domain, it will
128
Green’s Functions for Elliptical Equations
also converge uniformly in D. The reader is recommended to
prove this theorem using the Cauchy criterion for the convergence
of the sequence.
6.2.2 Theorem All the solutions of the second internal boun
dary-value problem
L[w\ = 0 ,~ | = 0
on s
Using (1) for u = w and v — w we obtain
dr = 0
9
129
Mathematical Physics
Fig. 6.1
(ID
131
Mathematical Physics
where
oq = ot^M) , <x2 = ol2(M) , a\ , a2 > 0 and ocj-f-a2 # 0
One of these methods is the method of separation of variables.
We shall illustrate it by the following two examples.
r i - (?&’)-).& = 0 (15)
dr
W"+XW=0 (16)
Using (14) we find that
V(<P+7ji ) = 'I'{cp) (17)
For X < 0, Equation (16) does not have solutions satisfying (17).
Consequently, X ^ 0. For X > 0 we find that
Wipp) — ylsin \/Xcp-\-Bcos ]/Xtp
From (17) we find that X 2n = Inn and hence Xn = n2, where n is
an arbitrary non-negative integer. Therefore, the eigenvalues of
the problem (16)-(17) are
X„ = n2 (n = 0, 1 ,2 , ...)
132
Green’s Functions for Elliptical Equations
and the corresponding eigenfunctions are
1, sin<p, cosq>, ..., sinncp, cosn<p, ...
For X = 0 the general solution of (16) is
W{(p) = A ay + B 0
It is only for A 0 = 0 that this will satisfy the condition given by
(17). Therefore, X = 0 corresponds to the eigenfunction W0(<p) = 1.
Consider now Equation (15). For X = n2 we have
r2& " + r & ' - n 2® = 0
The general solution of this equation is of the form
<Z>o(0 = C o + D o ln - ( « = 0) (19)
The coefficients A n and Bn can be found from (13) using the ortho
gonality of the eigenfunctions within the range [0, 2n] with the
weight p — 1:
CO
Ao= — J m m , = m cos « f d^
o 6
2n
whose coefficients (A0, B0, C„A„, C„B„, D„An and D„Bn) can be
determined from the boundary conditions
p (n 1 ,2 ,...)
. nn
/ , m = 2 j C„ sin——x
n= 1
and
eo
. nn
fi(x) = y , ( c ncosh ~ b f - D n sinh — 6 sin - f x
n= 1
Flence,
Mathematical Physics
and
i
„ , jt/2 , „ . , n n , 2 f „ .^ t
C„ cosh — 6 + Z)nsinh — b = - j ] /2(f) sin di
0
In Sections 11.1 and 12.4 we shall quote other examples of the
application of the method of separation of variables for equations
of the elliptic type, which will require the use of special functions.
Another method of solving the boundary-value problems for
elliptic equations is the method of Green’s functions. It consists
of the following. To begin with, one finds the solution of the prob
lem (10)-(11) for special values of the functions f ( M ) and y{M).
In particular, one solves the problem
L(G) = — d ( M , P) (29)
( - c + “ £ L - 0 (3o)
This solution is called the Green’s function for the problem (10)-(11).
We shall require that the required function G(M, P) be con
tinuous (together with its first-order partial derivatives if a2 ^ 0)
everywhere in a closed domain D except, perhaps, for the point
P at which G may have a singularity.
If Green’s function has been found, it can be used to find the
solution of the original problem (10)-(11). This can be done with
the aid of Green’s formula as follows
136
Green’s Functions for Elliptical Equations
u ( P ) = - \ k < p { M ) ~ d o M- ^ G { M , P ) f { M ) d z M (33)
S D
\ {GjL[G2]—G2L[G,]} d z M = \ k |g , - C 2S daM
D S ' n 1137
137
Mathematical Physics
The left-hand side is equal to
do-A/
6' '
r dG2 dGx
- ^ - G xG2- { - ^ G 2Gx = 0
dn a2 a2
G(M, P) = P)
\ V2WdrM = - 1
and, therefore,
= -1
ip(R) =
4nR
It follows that Green’s function G(M, P) is given by
* G(M , P) - ~ ~ + v ( M , P) (36)
mp
G(M, P) = — In + v ( M , P) (360
In ’Mr
The function v(M, P) is defined as the solution of the problem
i dv\ «i 1
V2w = 0, «2
_ r MP dn S 471
t ( P ) = ^ k ( M ) G ( M , P)tp(M) d<rM
s
140
Green’s Functions for Elliptical Equations
or
where
C = const da M
G{M, P) = — ^— + v
47trMp
G(M, P) = --------------------
4 7irMP 47irMPl
Fig. 6.2
G(M ’ p) = Tr~ln
1.71 T~~
t Mp _ i1.71 l n 7‘A/^P! “ i l n 7'MP2
^ + ^2 l l n rr_\tP3
"
In fact, in this case, the functions is equal to
1 iIn ----------
— -— 1 1— In
i -------
1 b -=—
,1 m
i ------
1
2^ i'mp1 2.71 rMp2 27i rM
It is harmonic within the right-angle D (as a function of the point
M) and is equal to —(l/27r)[In(l/rWP)] on its sides. This follows from
the fact that if M e lx then rMP = rMPl, rMPi = rMPi, if M e l2,
then rMp — rMPl, rMPl = >'mp3-
Example 3 Determine Green’s function for the first (internal)
boundary-value problem for the circle
G (A /,P ) = ^ - l n — + ®
2ai rAfp
142
Green’s Functions for Elliptical Equations
The problem reduces to the determination of the function v, which
is harmonic within the circle and is equal to —(\/27i)[\n(\jrMP)] on
its boundary.
Let P be the singular point of Green’s function, and let P x be
the image of P with respect to the boundary of the region (the circle
C). Pl is the image of P with respect to the circle C if both these
points lie on the same ray drawn through the centre of the circle
since the triangles OMP\ and OMP are similar. Therefore, the
function
2n PfMPy
is the required function. Consequently, Green’s function for
the first (internal) boundary-value problem for a circle is of the
form
^T>2TI '..2MPi~~P\
dG 1 1 [ R 2+rl1P- p - „2
dn !IC 2n 2 R r MP
2 2Rr2MPl
Substituting for rMPi from (37) and for p { from the formula
p, = R2/p, we obtain
dG 1 R 2- p 2
dn 2siR rlMP
From the triangle OPM we find that rjfP = R2jr p 2—2Rpcos(0—yi)
and therefore
dG 1 i?2—p2
dn 2tiR R2+ p 2- 2 R p c o s ( 0 - y )
Substituting this value into (39) we obtain the Poisson integral
2rt
(R2- p 2M 6 ) d6
(40)
- s- S R2jrp 2—2Rpcos(Q—ip)
where (p, yf) are the polar coordinates of P, and (R , 0) are the polar
coordinates of the point M on C.
PROBLEMS
144
Green’s Functions for Elliptical Equations
4. Derive Green’s function for the first internal boundary-value
problem for a plane layer 0 < z < /* (L[//] = V2//).
5. Using the principle of the maximum and minimum of harmonic
functions, show that Green’s function for the first boundary-value
problem for the domain D is positive in D {L[ii\ = V2w).
6. Using the Poisson integral prove the following theorems.
(a) Any harmonic function which is positive on a plane is equal
to a constant.
(b) Harnak’s theorem. Let {iti(x, >>)}, i = 1,2, ..., be harmonic
functions in a finite domain D bounded by the contour r , which
00
Potentials
fMP
where rMP is the distance between M and P.
If there are charges elf e2, ..., e„ at the points P lt P2, Pn,
the electrostatic potential due to these charges is given by
e\ <?2
u(M) = + + 0)
f‘M P l r MP2 r M Pn
146
Potentials
contain the amount of charge p(P) drP. The potential due to this
charge is approximately equal to
p{P) dr,
I'm p
u ( M) = [ ^ l drp (2)
g r MP
u (M )= $p(P)ln(— W (3)
u ( M ) = \ f ( M , P ) drP (4)
D
where
Let <5, = min {5, <$,}. For all points M, such that M M 0 < d2, we
then have inequalities (5) and (6) and, consequently, the inequality
u ( M ) —u(M 0) 1 < e
This proves the theorem.
It is evident that
! ij p(p) d rP ^ e \p \ drP < A \ dtp < A drP
fMP ^ r MP f MP
r2S ^mp
to
JMo WM0
°Mo "Mo 'M
O
where T$ is a spherical region centred on M and of radius 2d
(.D 5Mo c= T lf) . Transforming to polar coordinates, we have
2d 7t 271
f _^Tp_ _ ^ C C f /-sin 6 dr dddtp = 8And1 (r = rMp)
J
t 2<5
rMp J J
0 0 0
J
1M
du
dz \ (L r MP
l z2 P (P) drP (7)
sign and the formulae given by (7) are valid for du/dx, du/dy and
dujdz. To be specific, consider the integral
(i-x)p(P)
d r.
r MP
It is evident that
(t-x)p(P)
r3 drp < e
r MP
'M0
it is sufficient to take 6 < e/S tiA .
V2it = - 4 n p ( M ) (8)
We shall suppose that p(P) has bounded and integrable first-
order partial derivatives dp/dH, dp/dr], dp/dt,.
150
Potentials
Proof The second-order derivatives d2u/dx2, d2u/dy2, d2u/dz2 can
not be obtained by differentiating the right-hand side of (7) under
the integral sign since this yields divergent integrals, for example:
3 ( Z - x ) 2p(P)
—j — di> + d rP
I' m p r MP
1l x7
\
J
TM 0
p(P)f t ± - W
d£ \ r\fp I
=
1 Mo
1 - * dT,-
rMP
ji(
1 Mo
P -W
rMP I
dux \ P (P)
« d r , - cosa dfTr (9)
dx rMP i r MP
'M0 ^ Af0
dp_
( £ — x 0)
d2ul d£ POP) ( I —x0) cos a da>
dri
dx2 M=Mq rpi 5 rM0P - ! rM0P
1 Mo 3M0
However, (£—x0)/rM P = cos a, and therefore
dp_
d2Ui (£ —*o)
dr p P{P) cos2ada> (10)
'dx2 - s \ r2
M = Mo -,5 I' m q P PMqP
1M q *M
d2Ui
dy2 m =m 0
S
[A/n
(.v-yo) O
tP-
Yj
M qP
d tp -
A
\
3Mo
f M0 P
cos2/?da> (11)
dp
d2u i ( C - z 0)
dC
d T p - ( 4 — cosVddp (12)
dz2 M = M q *5 'M qP A rMoP
‘ M q 3M q
where /? and y are the angles between the normal to S\/o and the y
and z axes, respectively.
Combining (10), (11) and (12), we have
dp_ dp_
d£ dij
V2w|m =m0 — V2iti M = Mq J v2 co sad rP + \ —— cos/? drP
T„
<5 n t 0p
- -3 '" o p
J M 0 1 M q
(* ^P_
, \ d£ f P do-p
+ —2 cosy dtp — \ —2 (13)
j i Pm q P •> r M 0p
' M q 'M q
15l
Potentials
V2w| m = m 0 — 4jcp(M o)
follows that A = 0. From the condition that the potential and its
first-order derivatives must be continuous, we find that
2 , C, 4 C,
— ~ n R 2p + B = -jj- and — —n R p =
R2
and hence Cx — (4/3) ti R 3p , B = 2 n R 2p. Therefore,
P (B )
I' m p
j j
it = V21— *p = V21— * p = —47tS (M) * p(M) — —4np(M)
v( M ) = [ ^ - d o P (17)
s rMP
This integral is called the potential due to a simple layer. We shall
suppose henceforth that the function p(P) is bounded, \p(P)\ ^ H,
and the surface S' is a Lyapunov surface. The surface S is defined
as a Lyapunov surface if it has the following properties.154
154
Potentials
[ d(TP
' r MP
da p
< H
1/ ( x —£)2- r ( y - rj)2
'Mo
ds
co sy ]/(x —l)2+ (y —v)2
-M0
From the third property of Lyapunov surfaces, <5 can be taken to
be so small that for points P e S 5Mo we have cos y > 1/2. We there
fore have
ds
V ( x —£)2- r ( y —v)2
JM0
< 2H
s ds
}r( . x - t y + < y - v ?
In order that the integral |u1(A/)| be less than the given number e,
it is sufficient to take 5 < 1/8 tiH.
w(M) — e
rMP2 rMPi
or
zv(M) = <?/?— (—^
on \ rMP p=p*
157
Mathematical Physics
(19)
Fig. 7.2
158
Potentials
along the positive direction of the normal at each point. The geo
metric locus of the end-points of these lines is a surface S l which
is separated from S by the distance h. Suppose that on 5 there are.
negative charges distributed with density v(P)/h and on S x there are
positive charges with the same density (Fig. 7.2).
We thus have a double layer of charges of opposite signs, which
can be regarded as a set of dipoles distributed over S and S { with
the density v(P)/h. The potential due to the dipole on an element
da of S and S x is v(P)(d/dn)(\/rMP)da. The potential due to all the
dipoles is therefore given by
F ig . 7 .3
Mathematical Physics
The double-layer potential can, therefore, also be written in the
form
w(M) = ^ i'(-P) d(oMP (21)
s
7.3.3 It follows from (21) that the potential due to a double layer
is defined also at the points M on the surface S. We therefore have
the following two properties.
Property 1 The potential due to a double layer is defined every
where.
w( M) = v ( P ) ^ - \ \ da,,
rMP \P=P* g
St
( dcojvfp+ ^ dcoMP — 0
Si s2
The algebraic sum of all solid angles dcoMP will therefore also be
equal to 4n.
It follows that, in this case again, w(M) = 4nv0. If the point M
lies outside S, the solid angle dcoMP corresponding to elements d(TP
on .S', (Fig. 7.5) will be negative, whereas the solid angles da)MP
corresponding to elements do> on S2 will be positive. Therefore,
162
Potentials
The function w2(M) is continuous at M0. Therefore, for arbitrary
e > 0 the quantity \w2{ M )~ w 2(M0)\ will be less than e/3 if M M 0
is sufficiently small. Next,
and, consequently,
and
|«7,(Af0)l < y
Therefore,
\ w ( M ) - w { M 0)\ < \ w l{M)\ + \wl{MQ)\Ar \w2{M)—w2{MQ)\ < e
where w(M0) and w(M0) are the values of w{M) and w(M) at M0
on S. From these equations we have
163
Mathematical Physics
V2»2 = 0 (31)
u2\s = <p(M)-u1(M)\s = F(M) (32)
We shall seek the solution of this problem in the form of a double-
layer potential
w{M )= U p f - ^ d o p
i rMp
with suitably chosen v(P). For any v(P) this potential is harmonic
in D (in view of Property 2, Section 7.3). To satisfy the boundary
condition (32), it is necessary that the relationship weKt(M) = F(M) 164
164
Potentials
u d£
c
Using (34), the boundary condition leads to an integral equation
for v(s)\
\ ^ v ( S ) 6 ^ n v { s ) = F(s) (35)
c
We shall seek the solution of this equation in the form
v(s) = — F ( s) ± A
71
A d£+F(s)+JiA = F(s)
165
Mathematical Physics
and hence
JL $ f« )d f+ 2 ^ = 0 and
c c
It follows that
u(M) =
c c c
F($) df
71
or
1 f 2rRcos(p—r:
F(i) di
“ = L71 .)\2 Rr2
166
Potentials
PROBLEMS
Integral Equations
<p(M)—X J K ( M yP M P ) d r P = f ( M )
D (M )
Let
/ '( * ) = <P(x) (3)
We then have
X
^ W = J’o + S / © df (5)
o
Substituting (4) into (5) we obtain
X |
y ( X) = + J ( X - S ) e p ( s ) dS (6)
6
We have thus expressed the function y(x) and its derivatives
y'(x) and y"(X) in terms of the function tp(x) defined by (3), (4) and
(6). Substituting these values into (1) and taking the functions a(x)
and b(x) under the integral sign, we obtain the following Volterra
integral equation of the second kind:
£M = = -/(.v )
0 0
( 8)
Using the properties of the <5 function, the integral on the left-
hand side of Equation (8) is equal to G(s1} s2)—G(s2, ^i), whereas
the right-hand side is zero. For the first and second boundary-value
problems this follows directly from the fact that the functions
Gi and Gz, or the derivatives dGJdx and dGJdx, vanish at the
ends of the segment (x = 0 and x = /). For the third boundary-
value problem we can express the derivatives dGi/dx and dG2/dx
at the end of the segment in terms of Gl and G2:
dGi F1n(r. „ dG2
= h— C (0 , j ,). — —■(7(0, s2)
dx x= 0 x=0 0=1
dGi dGi 2
= -~Pl
^ G(l, s^ , G(l, s2)
dx x=I &-2 dx x= l 0=2
or
j-f e
k(x) Gx(x, s) ^ q(x) G(x, s) dx = —1
<7*0+0, s)—Gx(s—0, s)
k(s)
w M -P o m = 0 ( }
These relationships can be regarded as a system of equations
for a, and Since at least one of the numbers and /?, is not
zero, the determinant of the system (10) is zero:
Vi'(0) y i (0)
= 0
H(0) J i (0)
This determinant is the value of the Wronskian at x = 0 for the
solutions j,(x ) and >q(x). It is known that the Wronskian made
up of the solutions of a given linear homogeneous equation either
is identically zero, or is everywhere non-zero. Since in our case
w(0) = 0, the Wronskian for _y,(x) and y,(x) is identically zero.
Hence, it follows that jj(x) = Cy^x).
j C iyi(x), x s
G(x, s) on
{C2y 2(x), x > s
173
Mathematical Physics
C 1
k{s) k{s)
and hence C = 1. Green’s function is therefore of the form
* < s
(13)
{X’ S ) ~ \ y i(s)y2(x), x> s
From (9) and (13) it follows directly that
or
( <
— J G(x, s)f(x) d x + ^ ( x ) d ( x —s) dx = &(x)[(7(x, £)./(■*)
o o
—y(x)Gx(x, s)]£
From the boundary conditions (16) for y (x) and G(x, s) it follows
that the left-hand side of this equation is zero. Consequently,
i
J G(x, s) f ( x ) dx = y(s)
o
Changing the symbol used for the integration variable and using
the symmetry of Green’s function, we obtain the formula stated
in the theorem.
It follows that _y(x) satisfies (15) and this proves the theorem.
= 0, L2/ ( l ) + p 2y(l) = 0
0 (22)
L[y] = —f i x )
L[y] = —?>{x—s)+<P0(x)<P0(s)
« ./(0 )-/? ,j(0 ) = 0, «2/ ( / ) ] - p 2y ( l ) = 0
176
Integral Equations
which is continuous in [0, /] and is orthogonal to the eigenfunction
<P0(x) corresponding to the eigenvalue I = 0, i.e. such that
i
^ <P0(x)G(x, s) dx = 0
o
The corresponding properties of Green’s function, in this case,
can be established in a similar way.
we obtain
12 177
Mathematical Physics
where
b b
ctij = J ai(s)bj(s) dj, pi = dj
a a
b
<Pn(x) = SK(X, j)9Jn_i(j) ds (30)
a
178
Integral Equations
If we change the order of integration in the last integral, we obtain
b b
<Pi{x) = / ( * ) + '- $ Ki(x, s ) f ( s ) dS + / } 5 K2(x , t)<p0(t) d t
where
b b
S s ) f ( s ) d s + A n J K n(x, 0 ^ , ( 0 d t
\ k „(x , s ) f ( s ) d^
It is evident that
b
[A"3(jc, j)| < jj s )K2(s , f)| d j < A 3(b—a)2
a
b
\Kn(x, j)| < S \Kx{x, s)K„_y(s, 0! ds < A \ b - a y - 179
12* 179
Mathematical Physics
and, therefore,
b b
I $ Kn(x, s)f(s) ds | < An(b—d)n~1$ 1/0)1 ds < AnB { b - a ) n
a a
^ BAn\A\n(b—a)n (32)
n= 1
is the majorant series for (31). If |A| < 1/A(b—a), then the series
(32) will converge. Consequently, for such A the series (31), and
together with it the sequence of functions y„(x), will converge
uniformly to the function <jp(x). This function is a solution of (25).
In the limit, as n -> oo in (30), we obtain
b
q>(x) = X^Kix, s)(p(s) ds —/ (x)
and, therefore,
V(x) = V(x)
8.4.2 Since the series given by (31) converges for |A| < \/A(b—a),
00
it follows that for these values of A, the series ^ An\X\n~x{b — a)n~x
n= 1
will also converge. However, this series is the majorant series for180
180
Integral Equations
00 00
2 A"-1ATn(x, j) so that ]>] Xn~lK„{x, s) converges uniformly. The
n = -1 n= l
series given by (31) can therefore be written in the form
b oo
v(x) = / 0 ) + ;- ‘y)}/(‘y) dy
a n= 1
00
where the function R(x, s, X) = /,-n_1 Kn(x, s) is called the resol-
71 = 1
vent of (25).
It follows that if we know the resolvent of (25), we can use (33)
to obtain its solution. We have defined the resolvent only for small
values of \X\. However, it can also be defined for any finite domain
in the complex plane of the variable X by analytical continuation
(except, perhaps, for a finite number of singular points in this
domain). If this has been done, then (33) yields the solution of
(25) for any X except for these singular points. We shall not pause
to consider this in any detail.
<Pn(x) = /(* ) + * SK { X , S ) tp n^ i s ) dj
181
Mathematical Physics
The sequence converges uniformly in the range [a, b] for any
values of In point of fact, it is clear that
In general,
unique solution.
K ( X , S) = ^ tf ;( .v ) 6 ;( /)
f= l
182
I n te g r a l E q u a tio n s
The solution of the equation with the kernel K{x, s) will be the
approximate solution of the original equation.
A third method consists of the following. The ranges [a, b\
of the variables x and s are divided into a grid of n identical parts
b
<p(x) ^
K (x >sj)Vi ^ S j + f i x )
j'=i
Assuming in this expression that x is equal to x {(i = 1, 2, ..., n),
let us consider the set of equations
n
Solving this system for gj,-, we obtain the values of the approxi
mate solution at the grid points. We shall not pause to consider
these approximate methods in any detail and refer the reader to
specialist literature.
8 .6 F R E D H O L M ’S THEOREM
In this section and in the next chapter we shall confine our attention
to Fredholm’s integral equations of the second kind
b
has the trivial solution g?(x) = 0 for any value of the parameter
However, for some values of 2, it may also have non-trivial solutions.
and, hence,
b
(<Pi—<P2) —?-<
\ K(x,s)((pl—<p2)ds = 0
a
( m = 1, 2, ..., n)
= Po(*)+ 2 C‘I(M X)
<7=1
where <p0(x) is a solution of (36), <pq(x) are eigenfunctions of the
kernel K(x, j) corresponding to the eigenvalue 7. and Cq are arbitrary
constants.
Remark In Section 8.4 it was shown that the inhomogeneous
Volterra equation has a unique solution for any values of the
parameter /. Consequently, in accordance with Fredholm’s theorems,
the Volterra equation has no eigenvalues.
9
S<pi(x)vz(x) dx = 0
a
Let us multiply the first of these by <pi{x) and the second by <pi00
and subtract one from the other. If we now integrate the resulting
identity with respect to x within the range [a, b\, we obtain
1 j \b b b
K ------y H Pi 0 0 ^ 2 0 0 dx = jj jj K(x >s)<pi(s)<p2(x) d j d.v
' 1 Z' a a a
bb
' ~ \ \ K (x »J) P2(*) P i00 d s d.v
aa
189
M a th e m a tic a l P h y s ic s
b b
= ^ K(x, s)<p\(s)cp1{x) ds dx
a a
Consequently,
b
| ^ dx = 0
b b
190
Integral Equations with Symmetric Kernels
Let us multiply the second of these identities by i and subtract
the results from the first identity. We have
b
V>i(.y)-/> 2(.y) = (a— j)[vi(j)—/>2(j)]dJ
a
Therefore
/. = a —//? and <p(x) = (x)—iy2(x)
are, respectively, the eigenvalue and the corresponding eigenfunction.
Since / / /. (/? # 0), it follows from Property 3 that the functions
<p(x) and <p(.v) are orthogonal, i.e.
b b
5 < P ( x ) l p ( x ) dx = jj {y>?(x) Vwl{x)} = 0
a a
Integrating this inequality with respect to x within the range [a, b\,
we obtain
1 " n a
191
Mathematical Physics
Since all the Xn lie in the finite range [A0, B0], it follows that all the
numbers X\ are less than B2, X2 < B2, where B2 = max {A%, BQ.
p j
If we replace all the I 2 in the sum ^^'yY by the larger number B2,
n=l *n
we obtain for any integer p
p b b
^ ^ K 2(x, j) d j dx
n=1 a a
Integrating this inequality with respect to x within the range [a, b],
and bearing in mind that the eigenfunctions are normalised, we
have for any integer p > 0
b b b l>
Jf jj jj ^ 2(x, s) d.v dx for p < X2 jj ^ K 2(x, s) ds dx
n = 1 a a a a
192
Integral Equations with Symmetric Kernels
which is impossible. It follows that to each eigenvalue X there
corresponds only a finite number of eigenfunctions.
It follows from Section 8.3 that a degenerate symmetric kernel
can only have a finite spectrum. In fact, in order that a homogeneous
system of linear equations ((27) for /?, = (); Chapter 8) for the
coefficients Cf should have a non-zero solution, it is necessary
that the determinant of this system D(X) be zero. The eigenvalues
can be found from this equation. It is evident that it has only
a finite number of roots. The converse is also true; if a symmetric
kernel K(x, .?) has a finite spectrum, it is degenerate. In fact, suppose
that /j, X2, ..., is the spectrum of the kernel and <Mx), <p2(x),
... ,<pn(x) is the set of all the corresponding eigenfunctions of the
kernel (complete system). Consider the symmetric continuous
function
n
bb
= ,u $ $ K (n)(x, s)<pq(x)y(s)ds dx
a a
U -- — v.
/I
13
193
Mathematical Physics
Since the functions <pp(x) are orthonormal, the last integral
is equal to
when n and cp{x) are the eigenvalue and eigenfunction of the kernel
K(x, s), since
b b n
b
= /i J K M (x, = ip(x)
a
i.e. the kernel K(x, s) is degenerate. This means that the following
theorem must be valid.
Let us substitute
b
Acp = \ K (x, s)q>(s)ds
Integral Equations with Symmetric Kernels
It follows from the definition of iterated kernels that
b
A (A<p) = A 2<p = ^ K2(x, s)cp (^)dj
a
and, in general,
b
A n<p = A(An~l<p) = \k„ (x, s)<p(s)ds
a
£sn- l
rjV = o
T =T
since £sn = 1.
13* 195
Mathematical Physics
v ( x ) = ^ ( p p( x )
P=1
or
« .* * ) = o)
7=t p
in which the series converges absolutely and uniformly within the
range a < x, s < b.
We shall show, to begin with, that the series on the right of (3)
converges absolutely and uniformly for a < x, s < b. For this let
us estimate the finite series
m+q m+q
and the fact that \ap\ tends monotonically to zero as p -> oo. From
Bessel’s inequality
oo b
z
m+q
Pj>(*) Pj.fr). D
(5)
p=m
K~2
Since |/ m| oo as m -> oo, it follows from (5) in accordance with
Cauchy’s criterion that the series (3) converge absolutely and uni
formly.
Let
OO
= J] i PpfrOPpfr)
P- 1 'v
197
Mathematical Physics
This function is continuous within the range a < s, s < b. We must
show that Kn(x, s) = <P(x, s). Let us suppose that this is not so.
The symmetric function
Q(x, s) = Kn(x, s)~&(x, s)
will by Theorem 1 have an eigenvalue fi and an eigenfunction
ip(x), i.e.
b
y>(x) = n \ Q { x , s)y{s)ds
a
b
= [i
a 1 p= 1 P )
= 11\ k „(x,s)y>(s)ds
198
Integral Equations with Symmetric Kernels
We have used the orthogonality of the function ip{x) to all the eigen
functions (pp{x). Consequently, as was shown in Section 9.2, ip(x)
should be a linear combination of the functions cpp{x). However,
this is impossible since y>(x) is orthogonal to all the functions <pp{x).
Therefore, it cannot be assumed that Q(x, 5) 0 0.
/(* ) = X fpVpW
P= 1
(7)
where
b
fp =
a
This series converges absolutely and uniformly in the range [a, b\.
To prove this theorem the following lemma will be convenient.
199
Mathematical Physics
it is necessary and sufficient that it should be orthogonal to each
eigenfunction of the kernel, i.e.
b
S0(*)Pp(.^d* = O ( p = 1, 2, . . . ) (9)
a
Proof Sufficiency:
b bb
SQ(x)<pp(x) dx = /.p 5 5 K(x, s)<rP(s)Q(x) d.<r d.v
a aa
b b
= Ip 5 <pp(s) { 5 K(x, s)Q(x) dA-} ds = 0
a a
a = J $ X! ^ d*
a a p =1 p
oo b b
= J 4 J <Pp(x)Q(x) dA J <pp(s)Q(s ) ds = 0
p= 1 ^ a a
bb b
0 = Ji = 55 IS K 2( x , t)K2(t, s) d/} Q(x)Q(s) ds
aa a
b b b
= 5{[5 02(-v)d.v] [ j Kz(l, s)Q(s) di]} i t
a a a
b b
= 5{S*2(.Y,oe(*)dA )2dt
a a02
200
Integral Equations with Symmetric Kernels
Consequently,
b
J K2(x , t)Q(x) dx = 0 (10)
a
o ')
p
Proof of the theorem We shall first show that (11) converges abso
lutely and uniformly. From the Cauchy inequality we have
n +q
< ( 12)
201
Mathematical Physics
20) = y > ,p P ) - /P )
= j ; - f' = 0
Consequently, in accordance with the above lemma, it is orthogonal
to the kernel K(x, s), i .e.
b
$tf(x, j)g (x )d x = 0 (13)
a
Next, since the functions Q(x) and <pp(x) are orthogonal,
b b 00
b
= - Se w /w <1*
202
Integral Equations with Symmetric Kernels
Q (*) = X X Vpix)—f { x ) = 0
7=i Ap
and this proves the theorem.
£(*) = pX!
—i
o6)
Substituting (15) into (14), we obtain
oo b oo
or
oo b go b
2 c P<pP(x) = s K(x’ s ) f ^ ds+ ?- 2 c p \ K(x, s)<pp(s) ds
p —1 a p=1 a
X! < > p M s £ c„ ^ l
P = 1 P= 1
and hence
r
W = —
3 +I — C
1 '-'P for
1U1 C
'“'P = ;—77—
;
'■p Ap '-p '•
?>(*)=/(-'■)+; V - r f r ^ * ) (l7>
r^
or
b
$ f(* )< P r +t(x) dx = 0 (t = 0, 1, 2, 9)
9 .6 E X P A N S IO N THEOREM
i [ / J = - ^ [ * / '] - ? / = ~F(x)
is positive (negative).
It is readily shown that the necessary and sufficient condition
for the kernel K(x, 5) to be positive (negative) definite is that all its
eigenvalues Xp be positive (negative).
206
Integral Equations with Symmetric Kernels
In fact, the function
b
f ( x ) = ^ K(x, s)h(s) d?
a
f ( x ) = jj ( x , s) h (5 ) ds = ^ — -<PP(x) (24)
a p= 1 P
K(x,s) = Y f M l M
U ^
where epp and are the eigenfunctions and eigenvalues of the kernel.
We shall not prove this theorem here.
It is important to note that all the theorems and facts referring
to the Fredholm equations and described in this chapter and in
207
Mathematical Physics
Section 8.5 are valid also for the multidimensional case including
kernels of the form
H{P, Q) d
K(P, Q)
\PQ\a ’
where H{P, Q) is a continuous kernel, \PQ\ is the distance between
the points P and Q and d is the dimensionality of the space.
Example 1
.ysin ( a w ) Cl a.x
d.9 e
a2j\-s2 2
208
Integral Equations with Symmetric Kernels
where x > 0 and a > 0. Adding and subtracting these equations,
we obtain
CO
71
$ sin (vs) P i / ,_p—as ds = 1 / —
A L' 2 a2+ s 2 l / f [>
14 209
Mathematical Physics
2(a) defined by (28), the values of 7.(0) = 1/(2 jj H(z)dz) will increase
6
monotonically and continuously in the range 0 < a < A up to the
value /.(A) = 0. Therefore, to each value 7. e [0, 7.(0)] there corre
sponds a definite value a (a > 0), which can be determined from
(28) and, consequently, the solutions of (27) are
Q eu + C 2e-“
210
Integral Equations with Symmetric Kernels
where Cl and C2 are arbitrary constants. The solution corresponding
to 2(0) is
eoor_e-a*
<p(x) = lim = X
oc-»0 2a
which is readily verified by direct substitution.
00
When 2 > 11(2 $ H(z) dz), the solution must be sought in the
6
form e±'Px. Substitution of these functions into (27) then shows
that the values 2 = 2(//?) = 2(—//?), where
correspond to the real solutions cos fix and sin fix of Equation (27).
Therefore, the spectrum of Equation (27) is continuous; all the
non-negative X are its eigenvalues. Thus, if we take the equation
00
Gamma Function
The gamma function (or the Euler integral of the second kind) is
defined by
00
r ( z ) = $ e -'f'-M f (1)
6
This function has the following properties.
Ut z-1dz = 0
c
= z $ / z- le - 'd t = : f ( z )
o
Successive application of (2) yields
_______ r { z + n + \ ) _______
(3 )
z (z + 1)... (z + n —1) (z+«)
(-1 )"
Re f ( —«)
n\
216
Gamma Function
Property 5 /'(/z + l) = n\
Direct evaluation yields r ( 1) = 1. Using (3) we have
r { n + l ) = n\ (2.)
Therefore, I \ z ) can be regarded as a generalisation of the factorial
function to arbitrary complex numbers.
Proof
r ( x ) r ( l — x) = 4 ^ u2x~Ie~lj2du jj v~2x+le~v2dv
ou
, 2x-l
u
e iu2+v2) I — | du dv
- S0i0
Substituting u = r cos 93 and v = r sin<p, we obtain
03 j i /2
r ( x ) r ( l - x ) = 4 5 \ e~r2r(cot <p)2x~xdr dy
0 0
00 n /2 n /2
111
M a th e m a tic a l P h ysics
p ex/?d/3 ji
r ( x)r ( m = ^ T tV = i n ^
— 00
A zo) = -^(rt+l) = n\ ^ 0
Hence, the supposed existence of the zero z0 leads to a contradiction.
Property 8
5 e~'tz“]d/ = 0
c
Moreover,
0 = J e“ 'tz_I dt
c
F ig . 1 0 .3
where y\, y'z are shown in Fig. 10.3, and t k = ak-\-ifik (k — 1, 2, 3, 4).
Let us proceed to the limit in this equation as a -> oo, keeping j3k
constant. The integrals jj and $ will tend to $ and —jj , respectively.
K y[ r2 J'i
If we can show that the integrals jj and jj > evaluated over the
> l l2 tjt*
ranges t { t 2 and t3t4, will then tend to zero, this will prove the lemma.
Consider jj , bearing in mind that z = x-{-iy\
tit 2
= e-“ J \t\x~le~y!lTg'dt
fl<2
Fi(z) = § e {t z ~l dt
r
(upper edge)
1
The integral jj e~‘tz~xdt is not an improper integral but is a continuous
r
00
function of z. The integral jj e~'tz~xdt converges uniformly in any
1
band —N < Re z < IV, since for all t > 1 we have |e~Uz-1| < e~T‘v_1,
00
and the integral $ c~‘tN~xdt converges. Consequently, the integral
i
00
jj t~‘tz~xdt is also a continuous function of z in any strip for which
220
Gamma Function
— Re z < N, and hence in the entire plane of the variable z.
Hence, it follows that the function F1(z)is continuous in the entire
plane of the variable z.
Moreover, using Morera’s theorem, we can establish the analy-
city of the function Ffz). The integral \ Fj(z)dz, evaluated over an
c
arbitrary piecewise-smooth closed contour C is zero since tz~ l is an
analytic function of z and
Sf.(‘)d- = CS{i+S!d
C r 1
z
1 oo
= S (S z~ttz~l dz} d f+ § (J e- f / z-1dz) dt — 0
r c 1C
The change in the order of integration is justified because the integral
CO
Proof We have
221
Mathematical Physics
'_e‘ r
t~z ldt =
27ii (e - 2jiiz- f ) S e 2ni Y 1
F ig . 1 0 .4
11
Cylindrical Functions
za\c/0a2—a0v2\-i-z17+1[a! (a + I)2—atv2\-i-za+2[a2(a+2)2
—a2v2-\-a0]+ ...J\-za+n[an(o+n)2—anv2-j-an_2]-h ... = 0
If (2) is to be a solution of (1), we must have
a0(a2—v2) = 0
tfitOH-l)2—v2] = 0
az [(cr'b2)2— = 0
fl„[(u+n2) - ^ ] - f a „ . 2 = 0,...
From the first equation we find that
a=
since a0 # 0. Let us take a = v. From the second equation, we then
have ai = 0. Moreover,
a n~Z
an = (io-\-ri)2—v2 ’ n = 2, 3,
02k —
2lk+vr ( k + v + \ ) r ( k + \ )
We therefore obtain the formal solution of Equation (1) in the form
of the generalised power series
, M_ V (~it w k*’ m
^ L ixk+ v+ i)t\k+ i)
This series is, in fact, the solution of (1) in the region in which it
converges. It is readily verified using, for example, the d’Alembert
test, that the series converges everywhere except possibly for z — 0.
224
Cylindrical Functions
(-l)V „ (z )
r(^ + i)r(^ + » + i)
j= 0
11. 1.2 The behaviour of Bessel’s functions J,(z) and J_v{z) for
non-integral v in the neighbourhood of z = 0 is a manifestation of
a property exhibited by linearly independent solutions of the fol
lowing class of equations (with singular points)
and hence
d \fi ___ C __
dx y ! k(x)y](x)
Integrating this identity within the range [x, x j, a < x < .vi, we
obtain
*1
yz(xi) _ r Cdt
yi(x) = y f x )
> ip i) J k(t)y\(t)
X
or
*i
.y2p'i) 1; / x_i_ Cthjx) ( x - a )m (/-rz )-2
>’20')
y f x i ) y i ^' ‘ 2/w?i(f1)i/i(fi)
x < < X], ^ = £ ip)
and, therefore,
j ’2(x) = ^ i O ) > 'i O ) # i P ) 0 —a)
where
>2Oi) |___________£ ______ _ _
A fx)
j ’i (X]) 2m(x\ —a)1"'v(i\)u] 0 )
-c » iO )
B fx)
2ni(p{£i)u\{$\)
226
Cylindrical Functions
Since ^i(x) and B^x) have no singularities at x — a, this proves
the theorem for the first case. In the second case, similar calculations
lead to the following result
v2(.v) = A 2(.v) y \ (.y)+ B2(x ) \n{x—a)
where
C ln (.Y i— a ) = C
A 2{x )
yi(xi) <p(t2)u2i(h) ’ 2' <p(t2)u\{£2)
x < f2 < * i, h = h(x)
Since the functions A2(x) and Bz(x) have no singularities at x = a
this proves the theorem for the second case.
The above theorem is of considerable importance in formulating
boundary-value problems for Equation (4) within the range [a, b]
where either one or both end-points of this range are singularities
of the equation. If it is required to find the solution ^i(x) which
is bounded in [a, b], then by writing the general solution in the
form
y = Qj'iCjO+C^Ok)
we can find one of the arbitrary constants from the condition
that the solution must be bounded (C2 = 0). Therefore, the con
dition that the solution must be bounded plays the role of a boun
dary condition, so that it must be formulated as one of the boundary
conditions. .
11. 1.3 The Bessel functions Jv(Xx) are orthogonal with the
weight p(x) — x. More precisely, for any v > — 1
^ x Jv x j Jv xj dx = 0 when o i^ j3 (5)
0
where both a and /? are the roots of one of the three equations
j v(y) = 0, J'v(y) = 0, yJ'v(y)-\-hJv(y) = 0
+ jj Jv{Xx)—Jp{Xx) - ^ y , ( JMx)|dx
0
= (jf2—)}) xJp(?.x)JP(px) dx
or
We shall show that for v > —1 and ). = a/l, p = fi/l, the left-hand
side of Equation (7) will vanish. We note that, using Equation (3),
we can write
where Ri(x) and R2(x) are power series. When v > —1 the last
expression will vanish when * = 0. Substituting X — ajl, / 1 = /?//
in Equation (7), we obtain
-■*"(«) = - - / ; ( « ) + (
229
Mathematical Physics
Substituting this expression for J"(v) into (10), we obtain
Jv ( 11)
We note that the square of the norm in the range [a, 6] of any
function T„(2.\;) satisfying Equation (6) can be calculated from
the formula
to.
z d Y v(z)
xY„(?.x)dx Yl(z)
~2?3 dz
- Y , ( x x ) — Y,(j,x) 1*
L
and then pass to the limit as /u -> In evaluating the limit of the
right-hand side, we must use l’Hopital’s rule in Equation (6) just
as in the calculation of
Theorem 2 All the zeros of the Bessel functions J,,(z) with real
v > —1 are real.
230
Cylindrical Functions
231
Mathematical Physics
We then have
00
CO CO
or
Jv(rei9) = r veiwp [Afr, <p)-iD,{r, <P)] ]
( 12)
Jv(re~i9) = r ve~lvp[Av(r, <p) + iDv(r, <p)] I
where
00 00
or
(13)
Remark It follows from (3) that the roots of the equation Jv{z) = 0
are arranged symmetrically with respect to z = 0 on the z plane.
Using the identities
d 7_r(z) — iQ) \z*Jv(z)\ = zvJv_ f z ) (14)
d: zv
which can be verified directly, the reader can readily prove the
following theorem.
233
Mathematical Physics
The prime will represent the derivative with respect to 2. Next,
from the first formula in (14) we have
z /'(z )-l'/,(z ) = —zJv+i(z)
and hence, using (16), we obtain
J'v(z)\z =y — Jv+\(y) (18)
Returning now to (17), we have
dy 1
- r - M z)lz-y (19)
dv Jv+i(y) ^
Differentiating the identity
i ; ( z ) |i ,( z ) - j ,( 2 ) } / ; ( 2 ) dz dz
ov or
or
dym = 2v ?
— Jl{z) dz > 0
dl' ymJv+l(.Yni) 3 z
This proves the theorem.
= .&(R) = o, [<z>(0)j/<oo
r
R(t) = C t ~ 'alt for / > 0
The general solution of the equation for <Z>(r) 7can be written
in the form
<P{v)= AJ,(]r?.r)-rBN»(\r).r)
where N0 \ X r is a solution of the equation for 0 (r) which is linearly
independent of J0(]/2r).
By the theorem in Section 11.1.2 r) is unbounded in the
neighbourhood of r = 0. Since the required solution must be
bounded, we have B = 0.
Consequently, &{r) = A J qQ/X r). It is clear that we can set
A = 1. From the boundary condition for r = R, we find the equa
tion for the eigenvalues
Join) = o, n = \'X R
235
Mathematical Physics
By Theorems 1-3 (Section 11.1.4) this equation has an infinite
number of simple real roots
Mi < M t< ••• < Mn < •••
u(r, t) = ^ C„ e R2 J0 rj
2 Mk
Ck dr
R
236
Cylindrical Functions
7 -i- i 2 \
u(r, t) ^ Ci e + C2e ^ •/<>(— ->*
Fi and fi2 can be found from tabulated values of J0(x). They are
//, = 2.4048, fxz = 5.5201
c
This formula is obtained by double integration by parts of the second
term; A and B are the end-points of the contours of integration.
Let us now take the function K{z, I) to be (l/jr)e“ ,ZBlnf and let us
takea(£) as a solution of the equation
v " Jrv2v = 0
for example e”L The contour C is chosen so that all the above
operations can be performed and so that the expression Kv'—Kfi
237
Mathematical Physics
Ci
(26)
H[2)(z) = — [ e- ‘” ln{+w df
c2
which are called Hankel functions.S
o
the original integral over Cx will converge uniformly for Re z > (5 > 0.
The uniform convergence of the integrals
$ //<X)(z)dz = J J e - i2sinf+,'^dzd£ = 0
L Ci L
J7T ?! J
0 0
Consequently, the Hankel function H^l){z) can be written in the
form
CO oo
or
00 00
= - ^ - e ~ ivn[ d/5 _ J _ C Q—zsinhfl+i'fl d ^j
• JZ J 7Z J
0 0
71
+ — ( e - izslnf+i,,|d£ (262)
n J
o
As z tends to infinity, over the sequence z„ = x0+/y„, where x 0 is
fixed, x0 > 0, y n > 0 and y n -> oo, the first two integrals in (26j)
and (262) remain bounded since
00 00
(j e-znsinh^=FV/)d^
0 0
240
Cylindrical Functions
The integral
71
^ g —^nSinfg/Arosinf-iVI d £
0 0
will also be bounded as y„ -*■ oo. Consequently, remains
bounded for z„ -+ go. The integral
= j e>,nsin1e ~ sin1+‘‘,f d £
o o
in (262) tends to infinity as y„ -> + oo, since sin £ > 0 within the
range of integration. Consequently, i/<2)(z„) -»■ oo as z„ -> co. This
proves the linear independence of Hankel functions.
In fact,
71 J
Ck
16 241
Mathematical Physics
since the substitution of the limits into the integrand yields zero.
The procedure for //*2)(z) is the same. Next,
]_
H % (z)-H f\(z) ^ e-i*slnf4W2/ Sin fd f
71
Ck
On the other hand,
Since the functions Jv(z) and H (k\ z ) (k — 1,2) are analytic for
Re z > 0, it is sufficient to show that (30) is valid for z = x > 0.
Proof Let j v(z) represent the right-hand side of (30). We then have
Fig. 11.2
Let us substitute
«-£(f)'i( k= 0 y
-
i2 k +v
(-i)i-
—- = J v{x)
k=0
(&+V-M) r ( k + \ )
where we have used Equation (7) of Chapter 10.
Therefore,
^ W = y [ ^ V ) + t f ?>(*)]
or
1 e - i z slna + i r a d a
Jv(z) = (e'™—e“iwt) ^ e -2B,nh/'-
2ti
o
16* 243
Mathematical Physics
or
00 71
; „ ( z ) = — ( e- i2Si"“+inIda (34)
2 71 J
— 71
It is readily shown that the functions Jr(z) and N v{z) are linearly
independent. If this were not so, there would exist a constant C
such that for Re z > 0 we would have Nv(z) = CJv{z). However,
from (35j) we would then have
(1). From (35J it follows that the recurrence relations (28) and (29)
of this chapter are valid also for J„(z) and iV„(z):
l c \ , c2 , , cn
co H - - - - - - - - 1— 2~ + • • •
Z Zz zf
Ck
lim zn = 0 where n = 0, 1, 2, ...
z —>CO
ze£ k=D
The following notation is used:
Cn
/(* )' z
+
^
— +
ZH ^
or
Cl
/(z ) = c0+
245
M a th e m a tic a l P h y s ic s
and hence
ci = limz {/(z)—c0}
Z —►00
n—1
\ ^ Ck
c„ = lim z"
Z-*00 ^ > - k2= 0, 7
/ ( * ) ~ c 0+ - ^ + . . • + — +
then also
If
V (z ) c,
/(* ) = 7n + \
? ( z) Z
then
C|
y>(z) - <p(z) C0+ - i- - + —
z z
246
Cylindrical Functions
0 (A) =1 - ~ \ e (2dt
- J
y it J
V 71
We have
1 . 1- 3
e*2/(-v) = ^ - - - T O +
2a 22a 3 ^ 23a 5 -
we obtain
1•3 • 5 . . . ( 2« - l )
!*»(*)! < 2" +!x 2n+l ° { ^ Y)
Consequently,
( - i r ' ( 2« - 3) ! !
+ n v2n—1
2" A
and
3
23a5
' ( - l ) ’- '( 2n - 3) ! ! ^ + ( ) |_ _ l_ jj
n . J in — 1
247
Mathematical Physics
11.3.3 One of the most widely used methods for obtaining asymp
totic representations is the method of steepest descents. This is
based on the fact that for large values of the variable x the magni
tude of the integral
/(x ) = $y(£)ex,,(f)d£
c
is determined mainly by the part C„ of the contour of integration
C, on which |e;"p(f)| = ex Rcfp^l is large in comparison with the
values of this modulus on the remainder of C. The integral over Cu
is estimated more readily the shorter this region and the steeper
the fall of the quantity x Re 9?(£). In the method of steepest des
cents, one tries to deform the contour of integration C into the
most convenient form C in the above sense. By Cauchy’s theorem
this deformation does not affect the value of the integral provided
we remain in the region in which <p(£) and y>(£) are analytic. Since
the function <p(£) = u(a, /?)+iv(a, /?), £ — a+ij8 is analytic, the
direction of the most rapid variation of u(a, /?) is the same as the
direction of the line v{<x, /S) = const. The deformed contour Cn
must contain the point £0 = a0+i'/?0, at which u(ct, /S) reaches its
maximum value (among the values of this function on C).
It is readily shown that ^'(£0) = 0. In fact, the derivative of
u(a, ft) along the line C, taken at the point £0, is equal to zero since,
at this point, the function u(a, /?) reaches its maximum value
(along C).
dv
Next, = 0, since in the neighbourhood of the point
ds f=fo
£ = £0, we have v = const along C. Therefore,
da dv
P '( f o ) +i = 0
ds !f=fo ds Is=£0
The point £0 for the surface u = i/(a, /S) is, of course, a saddle point.
Therefore, in applying the method of steepest descents to the
asymptotic form of the integral \y»(£)ew(f)d£, the path of integration
c
C must be deformed into C which passes through £0 where
9> ' ( £ o ) = 0, and in the neighbourhood of this point coincides
with the line v(a, j3) = const = v ( a 0, fi0).
11.3.4 We shall use the method to deduce the asymptotic behav
iour of the Hankel functions 7/{,J)(x) and H (v2)(x) for large values
of x (x > 0).
Let us establish the following theorem.
248
Cylindrical Functions
■ i ^ + i ^ + [l- 7 ) y = 0 (lmv= 0)
has an asymptotic representation of the form
, ’(A) = ^ s in (x + '5o)[ 1 + 0 ( ^ ) ]
where A0 and S0 are constants.
y i(x)
y(x) =
}/ x
The differential equation for ^i(x) is
/ -----:
"ill (36)
y i + \ i■ hi = o
249
Mathematical Physics
where
and hence
sin2[f+«5(f)]d! (42)
For a fixed a and b -> oo, the right-hand side of (42) tends to a limit.
Consequently, the left-hand side will also have a limit:
lim 8{b) = (50
b —►oo
We therefore have
00
s (a) = <30+ ^ ~ 2-s in 2[£+<K£)]d£
However,
CO 00
y 2-sin2(£+<5)d£
x x
X X
and therefore
S (a) — <50+ O
x
250
Cylindrical Functions
From (38) and (41) we have
and consequently
Repeating the argument given for <5(.v) and d(b), we arrive at the
conclusion that the limit
lim In A ( b ) = In A 0
b —* oo
exists and
ln,4(x) = \ n A Q+ 0 ^ —~ j
Consequently,
^ ) - ^ [ l + o (l)]
Therefore,
= ^o^in(x+<50) | l + o |- ^
and
y(x) = —^sin(.x+<50) 1+ 0
Vx d )
where <p{£) = — z sin £. The saddle points S0 can be found from the
equation
251
Mathematical Physics
where k is an integer. Since the contour Cx lies in the strip
—7t < Re £ < 0, we shall be interested only in £„ = —ti/2:
<p(|) = —i sin | = —i sin(a + i/5) = —i sinacosh /? + cos a sinh/?
Therefore, v(a,/?) = —sin a cosh/?. In the neighbourhood of the
point £0 the deformed contour C, has the equation
—sin a cosh/? = —sin a0 cosh/?0 = 1
or
sin a cosh /? = —1
The direction of this curve at the point £0 = —n/2 is determined
by the angular coefficient
_d£_| = _ J _ = _!
da WI = ~ Y" sin a a„_^3L
~ 4
Fig. 11.3
_J_ 1-23
since
a2 —
— 7^1 1 T^j o2
@
2 X ° \/ 2 x a\/2
and
-1 v —
lx-° .x~ia
cos a2 sinh /?2 = = -+ ...
x a\/2 2 1/2 -3! 1/2 ‘ 3! 2 (/ 2
~ s \ . f\ a I ***
2x?° 6! x6<T
Similarly,
f g — ix sin I + ir | ^ gX cos 3 3 sinh 03
1J I
cm
x (cos 3 sinh 0—cos 33sinh 03) VP ]d^T'i
x 1 1 e
ii3. I.)
253
Mathematical Physics
11.3.7 Let us estimate the integral over Cu . We have £ = —rt/2
-\-rei3n/4 in the range (£2j—7t/2), 0 < r < x~a, and £ = —7t/2 + re~i:z/4
in the range (—jt/2, ?3), 0 < r < x- ?. These two formulae can be
combined into one:
t 7l
./ n n \ x a —i — —i—
iIac—v ——— J r* f.x[cos(re 4 ) —l]+ /vre 4
= e \ e
where
vr r6 r l0 p 14
mv(x, r) = —r=- +.v
1/2 6 7 “ 10[+ 144
.2+4 n
,»+i +
(2+ 4 n)\
,,4n
vr r4 r8
ntz(x, r) = —x 2------ — + ... + ( _ l ) n+i--------
|/2 4! 8! ^ ; (4n) !
or
254
Cylindrical Functions
As .V -* oo, the functions mk(x, OjX a) ( k , j = 1 ,2 ) tend to zero
as x~°. Consequently,
X~a _ jcr_2
R= + S e 2 dr
-x~a
where p^.v) and p2(.v) are bounded as x -*■ oo. Substituting r\Zx/2
= y, we obtain
R — { 1 ( p i + 2‘p2)x a
— {\+Pi{x)x~oO(c-0-5xX 2<!)}
where
P(x) = pi(x)-rip2(x),
Pi{x) = p ( i ) + / 0 ( e J '5sl 2°)
and the modulus of p3(x) is bounded for large x > 0. Therefore,
for large x > O w e have
255
Mathematical Physics
Using the theorem on the asymptotic representation of real
cylindrical functions (see Section 11.3.4) and the uniqueness of the
asymptotic representation, we conclude that and a2(x)x~a
are small quantities of the order of 1jx. Consequently, p(x)x~a and
p(x)x~a are small quantities of the same order.
Therefore, we obtain the following improved versions of
(43M46):
+ Z)2sin - y 1+0
F i g . 1 1 .5
17
Mathematical Physics
two neighbouring zeros of the Bessel function Jv(z) (and of the
Neumann function) tends to n as the absolute magnitudes of the
zeros increase without bound. Plots of Bessel and Neumann func
tions are shown in Figs. 11.4 and 11.5. that
— zjw = 0 (53)
K (54)
k=a r ( k + v + \ ) r ( k + \ )
258
Fig. 11.7
17*
Mathematical Physics
Using Equations (47) and (49), it is readily shown that for large
\z\ and | arg z | < n j l —d,
(55)
(56)
(59)
(60)
=o
260
Cylindrical Functions
1
beri(z) (61)
I 2 Z nl(n-\-1)!
1
beijCzJ (62)
w n=Q
where E\y\ is the integral part of the number y.
It is clear that
2 2
net 7in
W„(z) sin ^ z, /„
<Pn(r) = C J o ^ r j + D ^ o i ^ r j
nn \
set D„ = 0. Therefore,
(—j^-r I becomes infinite at r = 0, we musi
nn
0 n(r) = C J 0
261
Mathematical Physics
. sin
sin— z
h
The general solution of (64) (for A = A„) can be written in the form
262
Cylindrical Functions
R < r < oo and the function I0 \^j~ r j *s not bounded in this range of
/•[since for large values of r it behaves as | / /;/2^r2« r e(Jin/,,)r(l + 0(l//-)],
we must set Dn = 0. Hence
<£„(/•) = C„K0
? ')
The solution of the original problem will therefore be sought in the
form
for x > 0
y& =
for x < 0
If we set
C\ — —C2 = D x = D z =
and
Ci — C2 — D\ — —Do — —
1'
'-1 /3 -A * 4-V» for x ^ 0
3
Ai(x')
l7W
J-m yl.vlw +A» f |.v H for x < 0
for a' ^ 0
for .v < 0
264
Cylindrical Functions
PROBLEMS
Spherical Harmonics
'//(1, 0 = 1 ^ 7 = i + f + •
Consequently, P„(l) = 1.
Substituting .v = —1 into (1), we obtain
n -i,n =
268
Spherical Harmonics
(4)
d T [( l“ *V l + ; ° ' = 0 W
The second solution of (5), which is linearly independent of
P„(x)> has a logarithmic singularity at the point x = ± 1 , in accord
ance with the theorem in Section 11.1.2.
X ck*k- For / = «(«-!-1) this series cuts off at the term of power n,
k=0
i.e. for / = n(n-\-1) the solution will be a polynomial of order n,
i.e. P„(x). It differs from the Legendre polynomial of order n by
a constant factor. This factor can be chosen so that P„( 1) = 1.
- [ ( l - . v 2) P i ] + I ( H l ) P t (.v) = 0
Let us multiply the first of these by Pk(x) and the second by Pn(x).
Let us then subtract one from the other and integrate the result
with respect to x within the range [—1, 1]. The result is
-i v
1
= [k(k+ !)—«(« + !)] S p n(x)Pk(x) dx
270
Spherical Harmonics
or
i l
\ - - { ( 1 -X^iP'nPk-PnPL)} d-Y= (k ~ n ) (A'+ /J+ 1) [ Pn(x) P(x)dx
J d.v ,
-i ~L
Consequently,
j, = (A;-,)(l+B+l) W-. =0
when n ^ k.
IIP-1 Sp - p ’ 6 x =
-i
\ p^ r L
-l
x K - ^ L pA ix
2n-l
jj xP„Pn^ d x
n
where we have used the orthogonality of the polynomials P„ and
Pn_2. In the last integral the product xP n can be expressed in terms
of P„+1 and P„-1 in accordance with (6). This yields
2 2/z—l 1* jrc+ 1 n
l,p”L U " - 'i 2 ^ + T n + , + 2„ + i
In- 1 f
P L i d.v
2/2+1 J
or
2 / 2-1
LPJI2 = n — 1 li (9)
2 / 2+1
Here again, we have used the orthogonality of the polynomials
P--i and Pn+t.
272
Spherical Harmonics
From (9) we have
3-iiPii;: 2
\\P k\\2 ( 10)
2k+\ 2k+l
and hence ||Pil|2 = 2/3.
Fig. 12.1
w0 \ °) = ^T l cnrnPn(cos 0) (13)
nZ0
The coefficients c„ can be determined from the second boundary
condition using the orthogonality of Legendre polynomials:
18* 275
Mathematical Physics
The last two integrals can be evaluated using Equations (7) and (6)
of this chapter. The final result is
2/7+1
cII Pn+m
2 n
In this problem we have used the following theorem on the
expansion of the function 97(f) into a Fourier series in terms of
the Legendre polynomials:
If a function 97(f) is piecewise-continuous together with its
first-order derivative then at each point at which 97(f) is
continuous, its Fourier series in terms of the Legendre polynomials
converges to this function.
We shall not prove this theorem here.
1
^ di
1
(y/ sinO)+/<y/ = 0
sinO dd
Bounded solutions of the equation for ip are obtained when
H = //(/7+1) in the form of the Legendre polynomials
y){0) = P„(cos 0)
Substituting
9 + -)
/(r) =
\ r
276
Spherical Harmonics
we obtain
>2 n
t/ \ t ■ <p = 0
<J> ^-------- 9" -t-
gip cos 0
-% Jn +l/2(p)Pn (COS6) (14)
Z
»= 0 VP
l» ± y M = .+ ± 1 $ ew/».{{)df
IP 2
Integrating n times by parts on the right, we obtain
2 c » ^n + l/2 (P ) _ 1 cIpjp ^
[e ^ (W )]L ,
2/1+1 | p iP oipf
This result is valid for any p. For large p we can replace the function
yn+1/2(p) by its asymptotic representation. The result is
2c„
]/ 2 j 1 \ 71
cos p - /! +
2« + l PI 71 I
2 |2 c „ ^
j/ir(2 /i+ l) P
eip*Pn(!d
IP
277
Mathematical Physics
This shows that
71
I *
— 2/"sinlp—n —
a2Au — utt
it follows that v(M, t ) will also be a solution of this equation.
The functions a, u0 and v will be interpreted as velocity potentials.
On the surface of the sphere we must then have = 0.
Therefore, the problem fori/ can be formulated as follows:
dv ditp
V < CO
dn r —R dn r= R
278
S p h e r ic a l H a r m o n ic s
1)1
Bm Bm = 0
r2 \
If we now substitute Bm = Dm/[/r, we obtain
D ••
mT
l k2 (15a)
The coefficient ocm can be found from the boundary condition (152)
which yields
BUR) = - a 'U R )
Hence, we find that
j'mikR) h j k r )
&(M) = ^ ®UM)
w=0 m —0 ' '
h'UkR)
d_ dfc+1
( l - ^ ) * +1 Pn(x)
dx dxft+I
Proof Substituting
<s = S
A k„ , s = J ( l - ^ ) k^ P n( x ) - ^ P s ( x ) d x
dfc ifc-i
= C1- ^ ) dxh
r dk~l , d
d.x
(1~ X ) d ^ P"(x)
f d*~! d
Ai s = - ) -j z t t W n(-V)]<
d.x d.x
-1
The second factor in the integrand can be transformed with the
aid of (19) (replacing k by k — 1). The final result is
l
< • = ( n - k + l)(n+fr) jj
2i '
or
A k„,s = (/i - k + \) { n + k ) A kn:? (22)
Using this last expression for A kf , Ak~ f and so on, we obtain
4k _ (»-r/c)! a0
Ak
CM n s ----- -
— (23)
(n-k)\ n's
Since
282
Spherical Harmonics
and
it follows that
1
A k„,s = J P*(x)Pks(x)dx = 0 for n^s
-i
and
(«+£)! 2
= ILP„fcli2 = ( n - k ) \ 2n+l
Remark It follows from (21) and (22) that the /c-th order deriva
tives of the Legendre functions are orthogonal in the range
[—1, 1] with the weight p(x) = (1—x 2)k.
(r2F')—XF = 0 (25)
dr
1 d I . . d Y \. 1 d2Y
sin0 + -XY = 0 (26)
sin0 dO d6 sin20 dtp2
Definition The solutions of (26) which are bounded for 0 < 6 < zi,
0 < <p ^ 271 and are such that Y(0, tp-\-2n) = 7(0, <p) are called
surface harmonics (or generalised spherical harmonics).
If the bounded solutions of (26) are sought in the form 0(0) 0 (99),
<P{(p-\-2n) = 0(tp), the differential equations for 0 (0) and 0 (99) are
1 d f*
(0 ' sin0) + X 0= 0 (27)
sin0 d0 sin20
0"-Yn& = 0 (28)
Since the function 0 (9?) is periodic, we find that ju — k z, where k is
an integer and, therefore,
0 (9?) = A cos^9?+ S s\nk<p
If we substitute cos 6 = £ into (27), we obtain
d 2W AW k2
0= 0 (29)
r+ 1- f 2
which is the same as (16). When X — n(n-\-1), this equation has
the solution 0 = P£(|), which is bounded in [—1, 1]. Consequently,
surface harmonics of the form 0 ( 0) 0 (99), where 0 (99+ 2:1) = 0(<p)
are given by
Y k„(0, 99) = P kn(cos0) sink <p
and
f„"l ( 0 , 9?) = P£(cos0) coskcp (30)
Y°„(0, <p) = P°(cos0) = P„(cos0)
These functions are called the fundamental surface harmonics
of order n. It is clear that the functions
yn(0,9>)= 2
k=>—n
c *y "(0>9O (31)
284
Spherical Harmonics
are also surface harmonics. They are called surface harmonics of
order n. When 2 = n(« + 1), Equation (25) has solutions of the
form
1
F\(r) — /•" and F2(r)
and, consequently,
M r , 0,<p) = rnYn(0, q>)
285
Mathematical Physics
The square of the norm is given by
271 71 2.1 1
Proof Since
n
V«(0,<p)= X C “ Y kn ( 8,p )
k = —n
dfc \ ^
= (1—£2)fc/2^ r / | aq£"~2q cos ktp
R—0
0 < r < R, 0 < 8 < 7i, 0 < (p < 2n, t > 0, which satisfies the
boundary conditions
u(R, 8, <p, t ) = 0, \u(0, 8, (p, t)| < oo
u(r, 8, cp-'r 2n, t) = u(r, 8, <p, t)
and the initial conditions u(r, <P, 0) = fir, 8, cp).
z " + - z '+ = 0
r
The general solution of this can be written in the form
z(r) = M /„+i/2(]/a r)-\-N7_„_1/2(j/a r)
and, consequently,
Jn+1/2O) = 0, n = /a R
287
Mathematical Physics
Let
MI,n >P'2,ii! •••■>Ms, n5 •••
be the positive roots of this equation. We then have aSi„ = /j 2,sJ R 2
Solutions of (37) are of the form
J
d n , s (r ) —
I r
Returning now to Equation (35) we find that
a^l,n
in '
Bs. n — C s, n 6
> ^ J „ + i ^ \ — r\[Cs.n.kYi(d,<p)
n=A _ iI /c=0
0 s— I_A *
+ Ds.n,k Y ; k(0,<p))e R2
The coefficients Cs.„.k and D3 n k can be determined from the
initial conditions using the orthogonality of the functions Yk(0, <p)
and the Bessel functions. The final result is
( « — /c)!(2/? + l )
ne(n-\-k) ! [J!hU2(/<s, n)]2R2
R 2n n
X J jj jj/‘3/2/(/-, O,9p)y„+i/2/ — !-/-Jy5(0>93)sin0drdOd95
66o ' >
n = ___ (/i—A') !(2/»~r 1)
Us" " k neR\nYk)\[J^iz(M s,n)}2
R 2.-t n
X jj J jj r3/2/■(/-, 0, y)Jn+ii?Y,s~Li \ Y„ k(0, (p) s\nO dr do dtp
00 0
288
Spherical Harmonics
PROBLEMS
1. Calculate P„(0).
2. Determine whether the k-th order derivative of the Legendre
polynomials P2n( x \ where k is fixed, are orthogonal in the range
[0, 1]. If they are, determine the weight.
3. Solve Problem 6 of Chapter 2 for arbitrary initial data,
placing the origin of the coordinates at the fixed end of the string.
4. Determine the electrostatic field inside and outside a hollow
sphere whose upper half is charged to the potential Vx, and lower
half to the potential Vz.
5. Find the expansion in terms of surface harmonics for the
spherical charges induced in a perfectly conducting earthed sphere
by point charge e placed (a) inside the sphere and (b) outside
the sphere.
6. Consider the polarisation of a dielectric sphere of radius R
placed in the field of a point charge if the permittivity is e = ex
for r < R and e = e2 for r > R.
7. Find the potential of a simple layer in the form of a
circular disc.
8. Calculate the potential at all points of a conducting sphere
of conductivity a when a current / enters the sphere at one of its
poles (0 = 0) and leaves at the pole Q = n.
9. A sphere loses heat through its surface to the surrounding
medium which is maintained at zero temperature. A point heat
source of strength Q is placed at an internal point of the sphere.
Find the steady-state distribution of temperature inside the sphere.
10. Determine the potential due to a point charge placed
between two conducting earthed concentric spheres r = Ri and
r = R2, also determine the density of surface charges.
11. Find the steady-state distribution of temperature in a sphere
of radius R, a part of whose surface ^ (0 < a) is at a temperature
u0 = const, while the remainder of the surface S z is maintained
at zero temperature.
12. A sphere of radius R is heated by a plane-parallel beam of
heat of density q incident on its surface. The sphere loses heat
to the surrounding medium maintained at zero temperature in
accordance with Newton’s law of cooling. Find the steady-state
distribution of temperature.
13. Consider the oscillations of a gas in a spherical container
due to small oscillations of the wall, assuming that the oscillations
19
289
Mathematical Physics
Chebyshev-Hermite and
Chebyshev- Laguerre
Polynomials
Ot (=0
On the other hand, the n-th order derivative dnHjdtn for t = 0
is given by
dnH 2x 1 - 0
n\ f e2x
-d t
~d?~ t=o 2ni tn+ i
19* 291
Mathematical Physics
Cl
~ ( e - * y ) + A e - * 2j ’ = 0 (5)
dx
$ H„(x)Hp(x)e-**dx = 0 if n * p ( 6)
— OD
292
Chebyshev Polynomials
Proof Consider the two identities
- f e~**H„(x) = 0
- f [ e - * ,H;(x)]+2pe-*‘Hr(x) = 0
Multiply the first by Hp(x) and the second by H n(x), subtract one
from the other and integrate the difference with respect to x between
—oo and + o o . The result is
00
— 00
00
= 2( p - n ) 5 Hn(x)Hp( x ) e - x2dx
— 00
jj - ^ { ( H pH'n- H nH'p)z-*2}dx
— CO
and, consequently,
CO
00
2(p—n) jj HnHpt 2dx = (HpH'n- H nH'p) z - x2
—00
(9)
n=l n=0
CO 00
( 10)
n= 0 n=0
Equating the coefficients of equal powers of t in (9) and (10), we
obtain the recurrence relations (7) and (8), respectively.
The identity given by (8) can be used to evaluate the integral
2 (« + l)
The identity (7) can be used to evaluate the square of the norm
OO
ll^Ji2 = s e-xlH 2n dx:
—oo
OO oo
II/-/JI2 = \ e - ^ /7 1_1(.v){T/,1+1(.v)+2;7//n_I(.v)}d.v
—00
OO
= 2n J Q~x2H„_i(x)dx
or
|7/„||2 = 2n\\H„_i ||2 (11)
294
Chebyshev Polynomials
In this we have used the orthogonality of the polynomials H„.^
and H n+1. From (11) it follows that
00
Therefore,
\\Hn\\2 = 2"n\\/n (12)
13.1.4 The Chebyshev-Hermite polynomials form a normal
system. Consequently the theorem of Section 12.6 applies to the
polynomials. All the zeros of the polynomials Hn(x) are simple
and real.
13.1.5 The Chebyshev-Hermite functions
Hn(x) - * 2'2
rp„(x) (13)
lltf.ll C
are frequently used. They vanish at infinity. These functions form
an orthogonal system with the weight p(x) = 1:
00
5 Wn(x)xpp(x)dx = 0 if n+ p (14)
— CO
IIVn II = 1 (15)
From Equation (5) for the polynomials H n(x) one can readily
obtain the differential equation for the functions y„(x)\
y>"+(?.~x2)y> = 0 (A = 2 /i+ l) (16)
Example Determine the values of E for which the Schrodinger
equation for the linear harmonic oscillator
j 2mE m2ojl \
V + | - ^ 2 ---------------------- = 0 (1?)
has a solution which is bounded in the range —oo < x < oo.
In this equation in, oj0 , E are the mass, the natural frequency and
the total energy of the oscillator, respectively, and h is Planck’s
constant divided by 2 n . _____
If we substitute z = \/oj0m/hx into (17) this equation reduces
to the form given by (16), in which I = 2E/o)0h:
d2y> 2E
z2\xp — 0 (18)
dz2 o>0fi
295
Mathematical Physics
1 f d 2u d2u d2u \
V2w -c2(a2+ /?2) 0 (20)
c2(a2+ p 2) + ~df? d?\
We shall seek the solution of this equation in the form
it = A(ot) B(fi) D(z). Separating the variables we obtain the fol
lowing differential equations:
A"-sr(]x~?.2c2ix2)A = 0 (21)
B " - { jjl+X2c2P2)B = 0 (22)
D " + / 2D = 0
where A2 and /a, are unknown parameters. Substituting £ = j/A ca
and 7] = / ] /Ac/9 into (21) and (22), we obtain
d2^
i z \A = 0 (23)
d£2 Ac
and
Let us multiply the first of these by La(x) and the second by La(x),
subtract one from the other and then integrate with respect to x
within the range (0, oo). The result is
A Ta
x“+1e -x n
dx
= ( p - n ) $ La(x)La(x)xa e -xdx
o
The left-hand side of this equation can clearly be written in the
form
\ - ^ { x a+ , e - x[La(Lay - ( L ayLa]} dx
o
298
Chebyshev Polynomials
Consequently,
4 Lan(x) = - L * t \ ( x ) (29)
dx
To do this, we must find the connection between the derivative of
the function and its partial derivatives, dLa/dx and dLa/dt. Direct
calculation shows that
and
dx
Substituting (24) into these identities, we obtain
00 __
( 1 - 2 t + t z) Y , " L n ( x ) t n- 1 = [ a + l - x - ( a + l ) f ] ^ £ ; ( x ) / B (30)
n= l n= 0
and
CO
- / tnL*+t(x) (31)
«=0 n= 0
299
Mathematical Physics
—(n- 1 + a) L“_2(x)}— dx
oo
— [ Xa c~xL*_x(x) [—x£“(x)]dx
n J
We have used the orthogonality of the polynomials and £„_2
and the orthogonality of L an and L“_,. If we use (28) again, we obtain
00
|]L“ ||2 = — \ .va e - JCL “_ 1( x ) { ( n + l ) L n
% I(x )
n oJ
— (2 /7 + 1 + o t)Z ,5 S (A ')+ (n + a )L n -i(:v:)} d x
00
= 5 -v * e -'[i;-iW ]2d.v = IHS-.lr
or
l i r a ||2 _ n + a
l+nll — -------- l+n-1 II
lira |,2 (33)
r
7 > + « + l) , = r ( n + a + 1)
n\r(a+2) 1 ( + } n\ "
300
Chebyshev Polynomials
Therefore,
Ta |,2 T(>Z+ a -H )
(34)
LA ~ n\
# K * )= m « * v w (35)
L-'n
(37)
where
e - * = 2 c nL°(x)
/i = 0
301
Mathematical Physics
c“= i k f ) = izjjiF$
Integrating n times by parts, we obtain
n\ \ _ d”-1
c’ = ? > + « + i ) | e" d y ^ (v e~'>
d n-2 'co
+ e - T:^ r y(x',+ae -JC)| + . . . +
dxr‘ )|
|0
JA-n+ae -2xdA-j
Substitution of the limits yields zero.
Since
00
A' /?
< | S|/n (0 -/» .(0 1*1 < S\ M t ) - M 0 \ dr
A-o a
F„(x)—Fm(x)\ <
3 0 - a ) *o — 6M
+ j « * -)< ■
Consequently, the sequence {Fn(x)} converges uniformly'in any
segment [a,/?] B). Property 2 is therefore satisfied and this
proves the theorem.
gn(x) =
1+e"
and the interval (—oo, co) (Fig. A.l).
and the interval (—co, co) (Fig. A.2). In any segment [a, ft] c=
(— oo,0) or [a,/3] c= (0, oo), the sequence converges uniformly to
zero. However, it is not uniformly bounded. The sequence of func-
tions {F„(x)}, where F„(x) = \ /„(/)d/ will also converge uniformly
— QO
in any segment [a, /?] within the intervals (—co, 0) or (0, oo) and
is uniformly bounded (by the number 1) in the interval (—oo, oo).
Consequently, by Theorem 3 it is a fundamental sequence. If this
is so, then by Theorem 2 the sequence {/„(/)} is also fundamental.
Example 3 Consider the sequence of functions {^n(x)}, where
<p„(x) is a piecewise-linear continuous function equal to zero outside
the interval (—1/«, 1/«) (Fig. A.3). In any segment [a, /?] within
the intervals (—oo, 0) and (0, oo) it converges uniformly to zero.
However, it is not uniformly bounded. The sequence of functions
*
{$„(x)}> 0 n(x) — jj <pn(t)dt, will also converge uniformly in any
— CO
segment [a, /3] in the intervals (—oo, 0) and (0, oo) and is uniformly
bounded (by the number 1) in the interval (—oo, oo). Consequently,
by Theorem 3 it is a fundamental sequence. If this is so, it follows
from Theorem 2 that the sequence {<?„(*)} is also a fundamental
sequence.
Example 4 Consider the sequence of functions {<f „ ( a')} where
71
306
Appendix
and the interval (—00, 00) (Fig. A.4). Since the sequence of functions
{ ^ (x )} , where !/y„(x) = \ y>„(t)dt, satisfies the conditions of Theo-
—00
rem 3, this sequence is also fundamental. Consequently, by Theorem
2 the sequence (yn(x)} is a fundamental sequence.
9'(b), x > b
V(x) = (p{x), a ^ x ^ b
tpifl), x ^ a
It is clear that
x > x0
i ] ( x - X 0)(f.(x)
“ lo, X < x0
|'/ ? J = S Pn(t)dt
— 00
We then have
00
—00
yn(-v) = J d„(t) dt
—co
is uniformly bounded, since 0 ^ y„( x ) ^ 1, and in any segment
[a, /3] which does not contain x — 0 converges uniformly to zero
if [a, p\ c ( —oo,0), and to unity if [a, {$] <=(0, co). In point of
fact, in the former case {[oc, /5] cz(—oo, 0)} we can find a number
no such that for all n > n0 we have en < \(3\. Consequently, the
segment [a, fi] will lie to the left of all intervals ( —e„, en). Therefore,
all functions d„(x) with n > n0 will be zero in [a, /?]. Hence, it follows
that y„(x) — 0 in [a, /?] for n > m . In the second case {[a, /J] <=(0, oo)}
we can find a number n0 such that for n > n0 we have e„ < a. Con
sequently, the segment [a, /3] will lie to the right of all intervals
(—e„, £»), n > n0. Therefore, for x e [a, /S] and n > n0,
x e„
309
Mathematical Physics
The reader is recommended to prove this. Consequently, any 8
sequence {<5„(x)} defines the <5 function <5(x). The <5 sequence
{(5n(x—x0)} defines the <5 function <5(x—x0).
Multi-dimensional 8 sequences are defined in a similar way.
For example, in three-dimensional space the d sequence is defined
as being of the form
{<5,,(x)<5n(y)(5n(z)}
i.e. as consisting of the products Sn(x) <5„(y) <5n(z), whereas the 8
function in three-dimensional space is defined by 6 sequences of
this form. Therefore, by definition, 8(M) = 6(x) 8>(y)d(z) and
8 ( M , M 0) = {<5„(x—x0)8u(y—y0)8u(z—z0)}
= <5(x—x0) <5(y—y 0) 8(z—z„)
where x, y, z are the coordinates of the point M and xo, Jo, zo are
the coordinates of the point M 0.
We have shown above that the sequences {A 8„(x—x0)} are
fundamental sequences. It is possible to prove a more general
statement, namely that for all functions y(x) which are continuous
in the neighbourhood O(x0) of the point x = x 0, the sequence
{ f ( x ) ^ ( .r - x 0)} is fundamental in <9(x0) and is equivalent to the
sequence {9?(x0)^„(x—x0)}.
Proof Consider the functions
X X
X
^ $ 1^ (0 —y(^o)l^»(^—x0) dt
— OC
^ 5 \<p(t)—9K-vo)l 8n(t—x 0) dt
*0+en
< £ jj (5„(/—x0) d t = e
Xo—*n
310
A p p e n d ix
i.e.
|Fn( x ) - 0 B(x)| < £ (A.l)
Since the sequence {$n(x)} converges uniformly [(0„(x)
= <p(xo)y„(x—x0); see proof of Theorem 4], it follows that the
sequence {.Fn(x)} converges uniformly in any segment [a, /?] c= 0 (xo).
From this and from the inequality given by (A.l) it follows that
the sequence {^(x) <5„(x—xo)} is fundamental in 0(xo) and is equiva
lent to the sequence {<p(xo) Sn(x2—Xo)}.
The generalised function q>(x) S(x—x0) specified by the funda
mental sequence {^(x) (5„(x—xo)} is defined as the product of the
function <p{x) which is continuous in the neighbourhood of the
point x = xo and the d function <5(x—xo).
The above statement means that
(f (x) h(x x0) = cp(xQ) <5(x —x 0)
fix)= j m dt
-- oc
A.7 The integral of the product of fi(.v—.v()) and of an arbitrary
continuous function y{.v), i.e.
b
\ ? (a-) (5(.v-.v„) d.v
312
A p p e n d ix
lim $r(-'')<5nCv—x0)dx
5 (p(x)d(x—x0) dx = 0
Let x0 e (a, b). We can then find a number n0 such that for n > n0
the intervals (x0—e„, x0+e„) will lie entirely in the segment [a, b\.
Consequently, for such n,
b x0-]-En
5<p(x) S„(x —x0) dx = 5 (fix) 6„(x—x0) dx
a xo £n
Let us apply the mean value theorem to the last integral. We obtain
Wo-Ff/i -TO+f/i
\ (f (x)bn(x — x 0) d x — </■(£„) jj A„(x—x 0) d x = (p(£„)
Xq £n Xo
where i n e [x0—e„,x0+e„]. A sn oo, we have £„ -+ ,v0. Therefore,
since the function<p(x) continues in [a, b\, weobtain
b
lim J <p(x) d„(x—x 0) dx = lim $>(£„) = y(x0)
313
Mathematical Physics
I °« x =£ 0
«K-v) = \ co, x = 0
314
A p p e n d ix
within the interval outside which <p(x) = 0. The following are the
simplest properties of the convolution.
-^ r [/(■*) * ? « ] = U * v i r) = f ( x ) * ^ ' \ x )
( p = 1 ,2 , ..., A)
This follows directly from (A.4) and from the fact that the function
<p(x) is localised. If the function f(x) has continuous derivatives
up to order A and <p(x) is integrable and finite, then
( P = 1 ,2 , ..., A) (A.7)
Proof Let [a, /?] cz (A' B'). The segment [a—b, /?—a] then belongs
to (A, B). Since {/„(x)} is a fundamental sequence, there exists
an integer A > 0 and another sequence {/^(x)} such that
Fn(k) = f„(x) and {Fn(x)} converges uniformly in [a—A, /S—a].
We have
/„(x) * (p{x) = F {nk\ x ) * <p(x)
and, using Equation (A.7), we obtain
Proof Let [a, /3] c: {A, B). For any e > 0 one can then find
n0(s) such that for n > n0 and all x in [a, fi] and all t in (—e„, e„),
= | J <p(x~t)d„(t)&t— $ <p(x)dn(t)dt
— 00 — CO
CO
^ jj |9?(.Y—r)—9?(.Y)|<5n(r)dt
— CO
= J i9?(*—o —
—£n
For n > nQ(e) and any x e \<x, /S|, the last integral is less than
en
e J d„(t) dt = e
— en
2. 8(otx) = —<3(a-)
a
5. J q ( x ) 8 ( x - x ()) dx = <p(A'o)
— CO
OC
6. jj d(x—t)d(s—t)dt = 8(x—s)
— oc
318
Appendix
7. 6 \x) = — d(x)
00
9. jj d'(x—t)d(s—t)dt = d'(x—s)
dn n\
10. ^ (5(.v) = ( - l ) n— 6(x)
d.v'1 .v
1L $ <p(x) - ^ d(x ~ s) dx =
—CO
if 9?(n)(x) is continuous at x = s.
CO
12. \ $ $ <p(M)d(M, M 0) dvM = <p(M0)
—CO
13. The Fourier transform of the 6 function 6(,y—.v0) is given by
00
« ({ )= — t <5(.v-.v„)e«'d* = —
j/2.
V 27i J }/ In
and, consequently,
-i$(x-xo)
di
r — 00 -00
= &(x—x 0)&(y—yo)&(z—z0).
Therefore
d(M, M 0) = d( x —-X0) (5(j^ y0) d(z Zq)
319
Mathematical Physics
in Cartesian coordinates.
0 — oo 0
C H A PT E R 1
C H A PT E R 2
0(0, t) = 0(7,/) = 0.
8. a2uxx+(H/cp)I(t) = u„, u(x, 0) = ut(x, 0) = 0,
w(0, t) = w(/,/) = 0, c is the velocity of light.
/), x<0
„,(* ,/), *>0 ( '= 1 . 2 ) .
w(x, 0) = <P(x), u,(x, 0) = p,(x), a2i = Et/Pi (i = 1,2),
Mi(0, t) = m2(0, 0 , 0 = E2u1x(0, /).
10. Tuxx = [p+m 0<5(x—x0)] m„, u (x, 0) = <p(x), wf(x, 0) = 9>,(x),
u(0, t) = u(l, t) = 0.
Or:
( W ^X, /), X < X0
u (x,t) = \ a 2 ( U i ) x x = (M,)rr O' = 1,2),
I U2(x, t), X > x 0
m ( , 0) = 9)( x ) ,
x wf(x, 0) = ( P i ( x ) , i/,(0, 0 = 0, u2(l, t) = 0,
0 = w2(^o, t), T [w2x(x0, t) — u l x ( x 0 , /)] = mQutt(xQ, t).
11- Tuxx+ F ( t ) d ( x —v0t) = pu„, u(x, 0) — t/f(x, 0) = 0.
12. vx-\-Li,-\-Ri = 0, ix-\-Cvt+Gv = 0, v (x, 0) = cp (x),
/ (x, 0) = cp!(x).
13. vx+ Li, = 0, /X+ C a, = 0, w(x, 0) = cp (x), / (x, 0) = <p, (x),
—v(0, t) = R0i (0, t), c0v t(l, t) = i (I, t).
14. vx-\-Lit = 0, ix+Cv, = 0, v(x, 0) = cp (x), i(x, 0) = <p,(x),
- v ( 0 , t) = Li»i,(0, t), v(l, / ) - £ ( / ) = L ^ i t(l, t).
. . \ v i (x,t), x < 0 |/i(x , 0 . *<0
15
l®2(*> i), x >0 *(*» 0 — | ^ ^ x > o
(vk)x+ Lk(ik),+Rkik = 0, (ik)x+ Ck(vx), = 0, v ( x , 0 ) = <p(x),
i (x, 0) = 9>!(x), i,(0, t) = i2(0, t),
®2«(0, 0 - ® i « ( 0 , 0 = (l/c o )/,(0 , 0 (k = 1, 2 ) .
322
Answers and Solutions
For the current i (x, t) :
(ik)xx — CkL k(ik)u + Ck Rkik, (0, t) = z2(0,t).
A=/ k=2
2c 0 c 2c 3c 0 c 2c 3c
Fig. B .l
4=/ k=2
/ N / \ . /X /X / \ / \ / \ / \
~2c ~c 0 c 2c ~3c ~2c -c 0 c 2c 3c
k--3
. y x .
-4c ~3c -2c - c 0 c 2c 3c 4c
Fig. B.2
A =4
0 c 2c 3c 4c 0 '^ g,'2 c 3c 4c Jc 3c
Fig. B.3
5. n ( x , t ) = \/2a{F(x-\-at) —F ( x —at)},
where F(z) = p/p{r]{z-xQ) - p { z - \ - x 0)}.
Hint: solve the problem: a2uxx = utt, m(0, t) = 0,
u(x, 0) = 0, tt,(x, 0) = (p/p)d(x—x 0) , 0 < a- < go .
i/ i (.y, 0 = / ( r - .v /« ,) + j J l |£ |~ f P i £j - / ( ( + .v/a|).
V P\ E\ + | IhFi
324
Answers and Solutions
^ 2
Refracted wave: u2{x,t) = — ()l 1----- f ( t —xla2).
V p i E i +X/'Pi E,
8. u (x, tx) = 0, u(x, t2) = — u(x, 0), u (x, t4) - u(x, 0).
x+af
a + v o
12 .
CHAPTER 4
where A = y ^ F 0x0(/—x0) .
CO cos-^j-x s i n a t
ip v*
4. “ ( + ' ) = — 2_,
n= 1 1+
Pn pl
where ftn are the positive roots of the equation ft tan ft hi.
00
7. (a) i?cr = 7i a I \/ 15 ;
(b) i?cr = 0;
(C) i?cr is equal to the smallest positive root of the equation
( 1 - h R ) t a n ^ - R = $ - R.
a a
8. For boundary conditions of type I:
nl n2 p2\ , . 7in .np
= + 0 n.P ( x , y ) = sin— x sin — >>
(n,p = 1, 2, ...)•
In the case of { l — k), ?.np = ji2/l2{n2-\-p2),
K i = ;-2,i = 5 n2/l2.
, , , - . 2 71 . 71 .7 1 .271
but &it2(x, y) = sin— xsm-^-.y ^ ®2 ,i = sin— x sin-y-y.
Therefore, to a single eigenvalue X = 5ti2//2 there correspond
two linearly independent eigenfunctions.
For boundary conditions of type II:
i
K ,p = n
i l n2 , P2 \ r
0 n,p(x,y)
\= nn nP
COS-j-XCOS ~ y y
(n,p = 0, 1, 2, ...);
for boundary conditions of type III:
^ 'n , p P ‘n '~ \~ tX p >
327
Mathematical Physics
tan I t a n ]/ a £ = ~
/Zj /?2— /z3/z4—a
tan ]//? nz =
hsh6- p
2 lg 2/z+J_
n . » (* ,o = ^ { 7rfl(2/7-(-1)3 C° S 21 ° (2/z+l)2
n=0 V
328
Answers and Solutions
2«+l } . 2/Z+l
rrx— - / -v .
~2I a 2
8F0/ V i : 1)" • 2 « + i
12.
(2/z+ l)2 Sln 2 / ^
/7—0
2 /7 + 1
X cos - 2 r . ^ t.
n2n2
= c «= y \ sin+ r - d-
12
6
( //l 3 V>>)
w h ere z;(x) = / l i e a + /l2 e a • T h e co efficien ts A\ and
A2 are d e ter m in ed fro m th e fo llo w in g e q u a tio n s:
/7,(«3—Mi)
A\ A2
\/ hi —a h l
}/lh a = / / 2W2,
h2-V e Jr A 21h2
a
329
Mathematical Physics
X„ = — ( 2 « + l) 2 or Q (0 = - i j j «,(0, r )d r.
/
16. 2 ( 0 = w(x, O dx ((Puxx—fiu = ut),
0
oo
where u(x, t) = v(x)-\- C„e~a2A,,,sin 71XC -P*
21
/7= 0
aun
v(x) = — cosh -V?
— (/—x ) ,
/- i/fl a
|//9 cosh----- /
2«+1
sin- -71 £ d|.
= - A f) TT
n= 0
where Xn = T ( 2 n + l f / 4 l 2, a2 = l/i?C .
sinj/An (l—x)
X —=
A [/?(/?0+ / ? / ) + ^ o ^ c o s A / ’
330
Answers and Solutions
where a2 = 1/RC, X„ are the positive roots of the equation
R tan],/X„ I = —Rq\/ X .
w
4H, - a 2Xn> • 2 /1 + 1
19. H(x, l) = H0- , .. | i ^ sin six,
7i yZ_J
z—i —2/z4
2/z+ l 21
n=0
where a2 = c2/4 tio/li, Xn — Jiz(2 n+l)z/4l2.
00
Xe —
°2V" —1;—t^tCOS———
2/z+l n x .
(2 /i+ 1)2 21
c"= wVS®(f)0”(f)df’
+/zsinj/2„ x,
2h ] / r
positive roots of the equation tan yX I = ^ 2- .
2
Qr
24. &u = 0, ur(R,(p) Q
uv (r, 0) = (r, 7i) = -- -Q- r
nkR ’ 2k R 2k R'
0(0, 0 = 0, '1,(1, l j = - k a2 = p G j p .
(h—Io/4) sin— x
Ofln
»(.Y,0 = 4 a ] T
Iokn(2pn sin2ia„) I
n=i
where I is the polar moment of inertia of the transverse cross
section of the rod, G is the shear modulus, 7) is the moment
of inertia of the rod, p is the linear density of the rod and it,,
are the positive roots of the equation tan u = I J I 0p-
, 2n-\-1
4a& l2\ ^ 0)1 ---- 1— a^it — (2/z+ \)na sin cot
26. (a) U(x,t ) = —^ r ^ (2« +1 )2[to212—(2 « +1 )27i2a2\
n=1
. 2n~\~ 1
X sin----— x;
(b) Replace sin cot by cos cot in the last formula (co ^ ktxjl,
k = 1 ,2 ,3 ,...).
, . Tin
00
col sin — a t —n a n sin cot
2F0a l S p . nn . T i n
27. u( x,t) sin — t sin — .V
71T Z -J (uo2l2—n2a2n2)n
7T2a2
28 ■ i/(x, t) = —^ s i n — e ~—1,r2 (" '- "+ 0 (t) dr
X---<* —?T2/?2tf2
+ ^C „e /2 'sin — .v,
n=l
2 (* 7Zft
where A' is the thermal conductivity, C„ = ~ \/( + ) s in — £df
, I 0 I
2/7+ 1 . 2/7+ 1
4~ 'y i Cn sin —+ Ttis/sin——— n x .
n= 0 ~2.1 21
n^no
2/7+1
where C„
na(2n-\-1)
J[w^i+)+£+<+] sin-
TT“
77^ ds,
Mathematical Physics
. . 2aA , . o)
V*{X) = 1 E S (~ l) SmV X
and similarly for F = A cos cot.
Hint: the solution should be sought for (a) in the form u = v(x)
sin ait+zo(x, /); in case (b) in the form u = ^ (x ) sin co t
+ v 2(x)t cos cot+w{x, t), where a(x) sin cot (correspondingly,
zq(x) sin cot+v2(x)t cos cot) satisfies the equation and the boun
dary conditions of the problem.
00
2 qR2 1 &
31. u(r, t) = Fi(r)+tF2(r) e r.
k xr x:
n= 1
pi cos p n
334
Answers and Solutions
• Pn
sin— x lt 0 ^ x < x0,
pn
•
sin— x0
a,
where # n(x) =
1 Pn
s in a- 2- V—*o)» x0 < x < /,
sin-^—( / - x 0)
a2
c n = p ^ j j / ( * ) p ( * ) 0 n(X)d*,
P __
*1 cot— k2
un are the positive roots of the equation — x0 =
Q-i O-i Q2
cot — (x0—/)• The eigenfunctions &n(x) are orthogonal with
a2
the weight p(x) in the range [0, /].
33. u(x, t) is the solution of the problem
a2uxx = [1 + (C0/C)(5(x—x0)]wt,
u(0, t) = uil, t) = 0, u(x, 0) = /(x ),
( u f x , t), 0 < x < x0,
or u(x, t) = \ . . ,
1 u2(x, t), x0 < x < I,
= (Ui)„ (i = 1, 2). u f 0, 0 = 0 = u2(l, 0,
Ui (x0, t) = u2(x0, t ) , u(x, 0) = f { x ) ,
ku 2xixo, t) k u lx(x0, t) Cqut(x0, t),
00
sin p„x V
00
34. u(x, l) = E0Jr 2 l E uC 'y e 12RC
A A
35. (a) u(r, <p) = ~R rcos<p = —a;
(b) u = A + - ^ y ;
(c) u = Axy;
AEB B—A 2 2
(d) it = — + - 2 r 2 -~- ( a 2 - v 2) .
(c) „ = ;
(<0 « = + 4 B) y ~ -12^ [ 3 (- v 2+ r ) i - 4 .- )] + c ,
where D is an arbitrary constant.
38. C =
1
u0 __ L
i i - i L l + i
La c e2 \ b c,
Hint: solve the problem A«j = 0 for a < r < c, Ah2 = 0
for c < r < b, u,{d) = u0, u2(b) = 0, Wl(c) = u2(c).
du
- £ , ------ , «0 is the potential difference between the
E> r=c dn r —c
plates of the capacitor.
1 I £i - £ 2 1
39. u - u0 for R < r < c.
l , £l - £2 1
/? 1" £2 C
£l 1
e2 '
It = Uq for r > c.
1 , «i — 1
R 1 £■> c
Ui—u2+0.25A(Rl—Rf) R^ .
40. (a) u = lh+ ^ . ( r 1- R ^ + ~‘- "V "" ln ^ ;
\nR7—Ini?, r
1
41- u = I ^ - R ^ - R . R ^ R . + r J — - -
.2/1+1 2n+l
CO
cosh —--— 7ix cos ——— 7iy
2b 2b
42. =
in ,2/j+l
n= 0 (2 n + 1 ) cosh - —— 7ia
2b
Hint : see Example 2 of Section 6.3.
22
337
Mathematical Physics
CHAPTER 5
RC G1
1. u(x, 0 = ^ — 0| x- | /
41
— « ( ] / x 2+ t 2+ ( z + zo)2 >0»
where u{r, /) is the solution of the last problem.
2 j\ \ } /4 a 2(t—T) J \ ] / 4 a 2(t—r)
10
7 f e ) ) l e(r)dr
P-D2 ('•+D2i
6. w(r, /) = ~ 4a2, _ e 4a2/ ^
1s ^ 5 fW)
|/ 4 a27it
1 »» t2 (r-t)1 C+f)2
g 4a2(/ —i ) __g 4a2(/ —t) d | dr.
1/471 a2 ,) ,1 ] /1—t
00 }
338
Answers and Solutions
00
8. u(r,
0 ' '
/ 00
r2+i2 rt
4a2{t-z)Jo d£ dr.
t —r 2 a2(t—t)
00
e -/ir (x-xp)2
9. G(x, x 0; t) = — -e 4a2' .
| / 47ia2r
CHAPTER 6
22 *
339
Mathematical Physics
1
3. (a) G (M, P) =
Z_J Arc \ fMP„ ntp' I ’
n= 0
where P„ are the points with coordinates (p„, 0o,<po), P'n are
the points with coordinates (<pn, 0o, tpQ) and (p'0, 0o,<po) are
the coordinates of P,
A ) for n = 2k,
*2,
e„ =
for n = 2k-\-1,
i(f)
A for n = 2k,
(£ ) Po
A A for n = 2&+1,
Ri Po
i 2k
A Po for n = 2k,
Rz
Pn [2k+2
A p0 for n = 2 k-\-1,
A
i 2k
A A for n = 2k,
Ri Po
Pn :
R 2 \ 2k R 2
for ii = 2 k + 1.
Ri Po
(b )C (« ,f) = - V I. l’"------- In
where Pn are points with coordinates ( p n,( p 0), P'n are points
with coordinates (p'n ,<Po), ( P o , T o ) are the coordinates of P;
P n > P n , e n and <?' are determined by the same formulae as in
case (a).
340
A n s w e r s a n d S o lu tio n s
1
4. G(M, P) -
r MP'n
CHAPTER 7
B, rC ,
r
Ri
B = C = D/R2+ 4 n \ ( 1 / |- 1 / R 2) | 2p(f)d£.
_ , ^ I 47zRPo> r < R’
2. wtr) = t ,
\A n R 2p0jr, r > R.
| 27ip0(R2- r 2/ 3 ) - M f r 1, r ^ R,
“( r ) ~ i M O / r - l /r O , r
where A/ = 4/3 7iR*p0, r — ]/x 24-y2-\-(z—h)2,
fi = l /^ + Z + ^ + Z z ) 2,
p0 is the charge density, (0, 0, h) are the coordinates of the
centre of the sphere of radius R.
Hint: calculate the effect of a perfectly conducting plane
z = 0. Reflect the original sphere in the plane z = 0. The
solution in this case then takes the form
M
r < R,
u(r) =
ir - 4rx / r > R,
341
M ath em atical Physics
m (4 - —In/?—
\ 2 2 R2
4. u(r) =
M In— ,
r
where M = 7iR2p0, p0 is the charge density, R is the radius
of the circle.
5. The potential of the simple layer 0 x < 1 of density p0 is
/
1_____
v(x, y) = p0 jj In d f.
o |/ ( f —x)2jr y 2
«. w( M) = v0(*] —
cos <p
djP=v0>f ,/(f d£ •
2ti
- , * , X J C (p2- ^ 2) /( 6 )d 0
7. (a) ro(p,<p) 2 jiJ / ? 2 _1_ p 2 _ 2 ^ p c o s ( o _ (p) ’
CO
1 f y /(g ) d£
(b)
* —
3OO ( i s - x f + y 2 •
CHAPTER 11
a2<In / \
1. «(/', o = 8/4 \ —3~ - — -e
Z_j a ^ iC ^ n )
/2= 1
where a„ are the positive roots of the equation J0(a) = 0.
2. « (M ) = ^ C „ e r),
342
Answers and Solutions
1
C„ = /H iJ i-tU .
/ ’( * ) + ¥ ( * ) = 0, f ( r ) = % d f - r > ) + e .* 2 j „ j q
2h tii —
k
00
a-a21
n
3. //(r, ?) - H0- 2 H 0 V _ i _ R2
1 )
' / _ J
c- = w - o f 5 d r- m = ^ lni - -
/„ are the positive roots of the equation $ '(i? 2) — 0.
5. a2&.u+Q/p = utt)u(R,t) = 0, u(r, 0) = u,(r, 0) = 0, |w |< o o ,
“ d
\ a I -1
/ to \
70l — R)
a /
T I-=rr
“n \ cos——
anfl t.
70
R R
(b) u(r, t ) = f ( r ) cos to —
r (to2R 2— a2a l ) /] (a„) ’
n^= ,1 an
00
2A
, for = R2,
Bn aanJ 2(an)
B.
2A ^ol—' > )
Raa„ Jf(an)
CO
. a„r . a a„
7. (a) w(r, 0 = /( r ) sin a>f— ^ Cn
/J=l
2,4 a>
where C„ = ~ H dr,
a)
aanRJ0\ — RU^a,,) 0
m = a
'•Hr*
a„ arc the positive roots of the equation J0(a) = 0;
00
aotn
(b) u(r, l) = f ( r ) cos to t— ^ DnJ0 -r cos
aRa„
where Dn = C„
to
344
Answers and Solutions
, 2P V ^ ^ aal„,
ru ^ L Z -. 3 .4 " > K « > ) r 7”( J? T m /! '■
[1, n 7^ 0,
where e„ = t a£n) are the positive roots
^ 2, n 0,
of the equations J„(a) = 0.
“ s in h 4 z /J -r'
R
9. u(r, z) = j/0- 2 ( F 0- K 1) y ------- --
7 . _• i &n
n=i sinh-—L/jan,/1(an)
where a„ are the positive roots of the equation / 0(a) = 0-
10. u(r, z) = —~ — ,
v } 7iR2a tiR
- /
2 7xo Z_i s
n= 1
sinh^ zy°(x*‘Wir'') + const,
cosh-—hot2J I M
R h1cosh— z
R
11. u(r, z) = --y ^ ^
27ikZ—i <x2nJf(an)
n=1 a .s in h - ^ h + J tfa cosh /i
X70( » ,
where a„ are the roots of the equation J0(oc) = 0.
Hint: the solution u(r, z) should be sought in the form tt = v(r)
-I~w{r, z).
The formulation of the problems for v(r) and w(r, z) is:
, , , , V 2 qR J° [ r r) - - - »
n= 1
qR
where F0(r) = u0■
4k
,(n ) ' 2
n,m = IIL__
ft ,f ’
are th e p o sitiv e r o o ts o f th e e q u a tio n s
yJ'v'n(y)+RhJVn(y) = 0,
v„ are th e p o sitiv e r o o ts o f th e e q u a tio n
(Jh Jrh2)
t a n v a = —5— r - - v .
v —ni h2
00 00
;(«)
Ak —
— (y P
1R
CO oc 00
X c o s «<p dr d z dtp,
R h 2n ,.
4 (*f* f* / yln) \ . Tim
z>- ‘= j j j r f( r - 9’ , V - M V
X sin n<p dr d z d<p,
y [n) an d 7nm>t are d eterm in ed as in P ro b lem 13.
71
K ± m {x) = y JV,/2 ( x ) = - / . 1;2W - |l/ / 71X
- COS X ,
-V -
347
Mathematical Physics
tf iW M = y 71X
e ,JC, H V A x) = - • ] / | e -
/ 2 p -\X
7lX
w0
21. ?;(r) = K0(fir).
K a(PR)
u2~ u l . 71n
22. u (r, z) = sin—^ z ,
h
n = l '
-a2 / 2 n2k2\
25. u(r, z, f) = v(r, z)+ ^ ^ C n>*e l a*+- * r y
n= 1 fc=l
/ a n \ . 7ik
y J A ~ R rr a i r z ’
348
Answers and Solutions
Rh
u0 2Mo V 1 (—1)" r l nn \ o- nn
v(r, z ) = X ^ + — 2_, - 7 — -T M ~ r T nT 2-
A
Hint', the solutions should be sought in the form
u(r, z, t) = v(r, z)+ w (r, z, /)>
where t>(r, z) is the solution of the problem (see Problem 22)-
V zv = 0, v (R ,z ) = 0, v(r, 0) = 0, v(r,h) = u0, \v\ < oo.
CHAPTER 12
* m v
1. ^ ( 0 ) = (- 1 )' 22k(k I)2 ’ P2f t + l ( 0) — 0.
Cn = ( 4 n - l ) ^ ( | ) P 2n - i( y ) df,
o
J3> = ____ y )d { .
a]/2n(2n—\) j]
4. E — -V h , where u(r, 0) is the potential of the field
00
' Z-j
> P2n+1(C0sfl),
2
n=0
u(r, 0) = Oj
V 1- V 2
/ d \ 2,1+2
V, 2 ^ C nU P2n+1(cosO), r p P ’
4 « + 3 (-1)"(2«)!
Cn =
2n+ 2 22n( n \f '
00
\ ' R2n+
(b) u(r, 0) = —e / I yn+l *•«+* COS0).
* * t c\ r
n=0
///«/: the resultant potential V(r, 0), due to the point charge
of the induced charges, should be sought in the form of the sum
V(r, 0) = e j r ^ u i r , 0),
00
r < R,
XjCnnr)
n__n
=0' ’ p^cos^’
where u(r, 0) =
R
v Dn i P„(cosO), r > R,
n= 0
where rx is the distance between the point (r, 0) and the point
(r0, 0) at which the charge is located. Use the expansion
350
A nsw ers an d Solutions
CO
, .. e e —e
u2(r, 8) = -------\-e —
£->ri
■y __ «
Z—i n e1-|- («
(n + l)e2 /•on+1'-^ r P '‘(COSO)'
n= 0
Hint', the expansion coefficients are found from the conditions
diii ! du?
u{{R, 8) = u2(R, 8), e2‘
1 dr [r=R z dr j,=* ’
where r, is the distance from the point (r, 8) to the point (r0, 0),
where the charge is located;
(b) If the charge lies outside the sphere (r0 < R), then
e
u‘(r’ e) = ^£i r i - +n = 0 + > + i ) eT ^ p "(cos8)'
CJJ
2/1 + 1
u2(r,, 8) = e ---- ^ P n ic O s O ) .
Z—: «e, + (« + l)£ 2 f
n=0
7.
00
^
| ^ Z_i
2 ( 7) + [£ .(O )+ -P .-2 (O )]P .(co s0 ), r > R.
( n=0
„ 1 _ .
l/ f e 0) = —
351
Mathematical Physics
/ y 4 //-f3 / r \ 2"+1
8. v(r, 0) = P 2n+1 (cos 0).
2noR 2—J 2 n + l U ;
/?=0
In view of the symmetry of the problem v(r, ti/2) — 0. There-
00
, Q Q n + l - R h rn
0rn n , ..
9. u(r, ° ) - ^ + 4- 7 2 j ' n + R h / p +f p »(cos °) >
where is the distance between the point (r, 0) and the source
(ro,0 ).
Hint: the solutions should be sought in the form of the sum
Q -v(r, 0),
u(r, 0) =
4jrA:r]
where a(r, 0) is the solution of the problem
1 Jl \
A , _ 0. * . , ( * 0)+l,v(R, 0) =
Use the expansion of l//'i into a series in terms of the
Legendre polynomials.
p V i 1 i-2,,+1_ /?2”+1
10. u(r,0) = - ----A ^
2 - i [ R 2n+l_ R f71+ I n+1
where r ( is the distance between the point (r, 0) and the charge
placed at (ro, 0).
352
Answers and Solutions
11 . «0, °) = - ~ | l - c o s a - ^ [P„+1(cosa)
n= 1
ln
—Pn-!(cos a)] I—I Pn(cos 0) |
■)
„ / rn_ ^ L 1 i r cos0
12‘ ^ 2k \2 h R "h P l + P/z
uu
(4 n + l)P 2n(O)r2nP2n(cos0)
- s (2n+/zP) (2n - 1) (2n+2) i?2n
/J=1
i/w f: the boundary condition for the problem is
f qjk cos 0, 0 ^ 0 < n/2
ur(R, 6)+hu(R, 6) = | 0; sr/2 < 0 < re .
+ Pn(C0S 0) ^ % ( 0 - |^ ^ n +l/2 ,
^•Ai+1/200 2'•^n+l/2(a) = ^ ’
t
n ( 0 = — i/"(T)sin-fl“«-((-T )d i,
0Ct. Cl t) iv
* 0
_2 j r M* V . +10( ^ » ) d i-
Ah —
nR n+1 n ( n + 1)
Jn+ 1/2 (a Jt) 1
a*
where w(r, 0, t) is the velocity potential a2Vzu = u„,
u(r, 0,0) = ut(r, 0,0) = 0, ur(R, 0, t) = Pn(cos 0)f(t), |w |<co.
14. ).n m k = l/R 2(a<
i"))2+ k 2 are the eigenvalues, and
23 353
Mathematical Physics
= ^ y"+‘'J( ^ r)F;(cos9){sintr
are the eigenfunctions, where o4n) are the positive roots of
the following equations:
(a) ^ ,+i/2(a) = 0,
(b) Jn+1/2(k) 2aL^n+ll2(a) =
(c) 2ctJ'a+1p(a)—( l —2Rh)JB+lf2(ci) = 0,
h is the constant in the condition vr{R, 6, q})-{-hv(R, 6, <p)= 0.
co co n
15. w(r, <?, <?) = X X X {Cn,m,kc° s k cPJrD n rn'ks\nk(p}
m—1 n = 0 /c = 0
a2[>}]2
1 / a(n) \ R*
X —^ / n+1/2| — rjP i(co s 0)e
This book list is included for the benefit of Western readers who
do not have access to the Russian references quoted in the original
text - Editor.
General
morse, P. m ., and feshbach , h ., Methods of theoretical physics
(2 volumes), McGraw-Hill, New York and London (1953).
courant, r ., Methods o f mathematical physics, Volume 2, Inter-
science, New York and London (1962).
Chapters 1-7
epstein , B., Partial differential equations, McGraw-Hill, New York
and London (1962).
Sommerfeld, A., Partial differential equations in physics, Academic
Press, New York and London (1949).
Chapters 8 and 9
mikhlin , s. G., Integral equations, Pergamon Press, Oxford (1957).
Chapters 10-13
Kantorovich , l . v ., and krylov , v. i., Approximate methods of
higher analysis, P. Noordhoff, Groningen (1958).
erdelyi, A., Higher transcendental functions. Volume 3, McGraw-
Hill, New York and London (1955).
Whittaker , E. T., and w atson , G. n ., Modern analysis, Cambridge
University Press, London (Fourth edition, 1927).
Appendix
lighthill , M. J., Introduction to Fourier analysis and generalised
functions, Cambridge University Press, London (1958).
355
Index
357
Index
Cylindrical functions 223-67 Elliptic equations 5, 6, 7, 9, 25
asymptotic representations 245, 29, 131
256 Green’s functions for 124-45
differential equation 279 Energy integral 98
examples 261-63 Equivalent sequence 303,307,311,
312
Error function 110
d’Alembert formula 34, 35, 36, Error integral, asymptotic representa
38, 40, 43, 44, 51 tion of 246, 255
d’Alembert test 224 Euler integral 215
Degenerate kernels 177, 178, 182, Expansion theorem 61, 97, 199-
194 203, 205-6
Degenerate symmetric kernel 193
(5-function 113, 118, 119, 307,
310-14, 317-19 Force density 14
definition 307 Fourier coefficients 201
<5-sequence 308-10, 313, 316, 318 Fourier integral 209
multi-dimensional 310 Fourier law 24
Differential equation, for cylindrical Fourier method 59-106
functions 279 Fourier series 62, 64
Legendre 270 Fourier transform 319
linear 3 Fredholm integral equations 170,
Differentiation of sequences, term-by- 176, 177, 178, 182-83, 185, 1S8,
term 318 207
Diffusion coefficient 25 definition 168, 169
Diffusion equation 25 Fredholm theorems 183-87, 203
Dipoles, distribution of 158 Fundamental sequence 303-12,
potential due to 158-59 315-18
Dirac (5-function. See (5-function Fundamental solution 108, 111,
Divergence operator 22 119-21
Duhamel’s method 87, 91-96 convolution of 115
Fundamental surface harmonics
284
Eigenfunctions 59-81,83-84,133, orthogonality 285
135, 177, 184, 186, 187, 195-97,
206, 236, 262
definition 61 Gamma functions 215-222
orthonormal 190 Gauss’s theorem 19, 20, 21, 49
properties 65, 74 Generalised Chebyshev-Laguerre
properties of set of 81-84 polynomials 297
simplest properties of 188-94 Generalised functions 38, 303-20
Eigenvalues 59-81, 84, 132, 135, Generalised spherical harmonics 284
173, 184, 186, 187, 195-97, 206, Green’s formula 124-26, 151, 175
235, 236 Green’s function 108-12, 174-77
definition 61 for boundary-value problem 170,
properties 65, 74 176, 205
r-fold degenerate 78 for elliptical equations 124-45
simple (or non-degenerate) 78 for parabolic equations 107-23
simplest properties of 188-94 method of 131-34
Electric field in vacuum 23 partial derivative 172
Electromagnetic waves, propagation Green’s theorem 25, 63
of 301
Electron motion in Coulomb field
301 Hankel functions 237-45, 251,
Electrostatic potential 146 280
358
Index
Hankel functions-cwzr/wwecf Kernel s , iterated-continued
asymptotic behaviour 248 spectrum of 194-97
definition 258 negative definite 206
linear independence of 241 positive definite 206
Harmonic functions 124—26, 128, symmetric 188-211, 206
130-32, 150, 160, 285 degenerate 193
Harmonics 62 specified in finite domain 208
Heat balance 24 specified in infinite domain
Heat propagation, along infinite 208-11
straight line 107, 112 Volterra 169, 188
in two- and three-dimensional space
119-22
Heat transfer 24 Lame coefficients 320
on semi-infinite straight line 111 Laplace equation 154, 167, 283
Heat transfer equation 25 homogeneous 3
Hilbert theorems 174, 175, 176, inhomogeneous 3
206 with constant coefficients 4
Hilbert-Schmidt theorem 199- with two independent variables
204, 206, 207 3-12
Homogeneous equations 3, 41, Laurent expansion 243
42, 119, 120, 168, 169, 183-86 Legendre differential equation 270
Homogeneous polynomial of order Legendre functions, associated
n 286 280-83
Homogeneous wave equation 32, graphs of 283
49, 54 orthogonality of 281
Hooke’s law 14, 15, 17 Legendre polynomials, examples
Huygens’ principle 57 275-80
Hydrodynamics 20 generating function 268
Hyperbolic equations 5, 7, 8, 10, integral of 272
20, 41, 59 orthogonality 270, 272, 273,
275, 277
properties of 268-80
Improper integral 147-51, 160 zeros of 273
Independent variables, transforma l’Hopital’s rule 52, 229, 230
tion of 5, 7 , Linear differential equations 3
Infinite medium, vibrations of 48 Linear harmonic oscillator 295
Inhomogeneous equations 3,117, Linear integral equations 168-69
121, 168, 169 Linearisation 22
expansion of solution 203-5 Localised function 314-16
Inhomogeneous wave equation 40, Lyapunov surface 154, 156, 162
48, 52, 55
Integral equations 146, 168-87
Fredholm’s. See Fredholm Maxwell’s equations 23
linear, classification of 168-69 Mean value theorem 126, 313
problems leading to 169-77 Membrane, transverse vibrations of
Volterra 169, 170, 181, 187
with degenerate kernels 177-78 16-20
with symmetric kernels 188-211 Method of characteristics 32-58,
Iterated kernels. See Kernels 68
Method of separation of variables
59-106, 122, 132
Kernels 182, 183, 184 principle of 59-81
classification 206-8 Method of source functions 121
degenerate 177, 178, 182, 194 Method of steepest descents 248
iterated, expansion of 197-99 Method of travelling waves 35
359
Index
Morera’s theorem 216, 220, 221 Reflection of waves at fixed and free
ends 45
Resolvent 181
Natural frequencies of vibration 62 Rod, longitudinal vibrations 15-
Negative definite kernels 206 16
Nernst’s law 25 vibrations of 26
Neumann function 244, 258
definition 258
Newton’s law 27
Newton’s second law 14, 16, 19 Schrodinger equation 295
Non-homogeneous equations 185 Semi-infinite straight line 46, 51
Norm 63 boundary-value problems 43, 47
Separation of variables, method of.
See Method of separation of
One-dimensional wave equation 15 variables
Orthogonality 190 Set of functions, definition 81
associated Legendre functions 281 test for completeness of 81
Bessel functions 227-30, 232, Smooth functions 213, 314, 316,
236 318
eigenfunctions 190 Spherical functions 286
fundamental surface harmonics of order n 285
285 Spherical harmonics 268-90
Legendre polynomials 270,272, generalised 284
273, 275, 277 Steady-state conditions, determina
surface harmonics 285 tion of 25, 28
Steady-state solution 71
Steady-state temperature distribution
Parabolic equations 5, 8, 9, 99 130
boundary-value problems 59 String, vibration of 13-15, 26, 32,
Green’s functions for 107-23 40
Partial differential equations, in Sturm-Liouville problem 61, 79,
physics 13-31 170, 176
second-order 3 Surface harmonics 283-88
Planck’s constant 295 fundamental 284
Poisson formula 52-55 orthogonality 285
physical interpretation of 56 of order n 285
Poisson integral 144, 167 orthogonality 285
Positive definite kernels 206 Symmetric kernels 188-211, 206
Potentials 146-67 degenerate 193
application to boundary-value specified in finitedomain 208
problem 164-67 specified in infinite domain 208-
due to dipoles 158-59 11
due to double layer 157-64
due to simple layer 154-57, 167
volume 146-54
Propagation of electromagnetic waves Temperature distribution 24, 112
301 along homogeneous rod 70
Propagation of heat, along infinite steady-state 130
straight line 107, 112 Term-by-term differentiation of se
in two- and three-dimensional space quences 318
119-22 Thermal conductivity 24, 130
Three-dimensional space 49, 52,
310
Rational numbers 303 Two-dimensional space 54, 55
Real numbers 303 Two-dimensional wave equation 20
360
Index
Transformation of independent Volterra integral equations 169,
variables 5, 7 170, 181, 187
Volterra kernels 169, 188
Uniqueness of solutions, for propa Volume potential 146-54
gation of heat along infinite properties of 148-54
straight line 107
of boundary-value problems 96-
102, 126-31 Wave equations 23, 26, 32, 38
Uniqueness problem 96 homogeneous 32, 49, 54
Uniqueness theorem 107, 128 inhomogeneous 40, 48, 52, 55
Unit step function 307, 312 one-dimensional 50
two- or three-dimensional 48
Waves, reflection at fixed and free
Vacuum, electric field in 23 ends 45
Vibration, natural frequencies of 62 Weierstrass theorem 311
of elastic rod 15-16 Wronskian 78, 173, 174, 226
of infinite medium 48
of membrane 16-20
of rod 26
of string 13-15, 26, 32, 40, 72 Young’s modulus 15
361