The Cauchy Method of Residues (Vol. 2)
The Cauchy Method of Residues (Vol. 2)
The Cauchy Method of Residues (Vol. 2)
Volume 2
Managing Editor:
M.HAZEWINKEL
Centre for Mathematics and Computer Science, Amsterdam, The Netherlands
Volume 259
by
Dragoslav S. Mitrinovic
and
Jovan D. Keckic
Dtpartment of MaJMmaticS.
Univtrsity of Btlgradt.
Btlgrade. Yugoslavia
This book is a revised and updated translation by J.D. KeCkic of the supplement
published in the second edition of D.S. Mitrinovic - 1.D. Kec"kic:
Cauchyjev raeun ostataka sa primeru:una,
Belgrade 1991. pp. 264-399. with new material added.
Table of Contents
Preface
ix
Contents of Volume 1
1.
I
1
1
Introduction
2.
Evaluation of Residues
3.
3.1.
3.2.
3.3.
3.4.
3.5.
3.6.
3.7.
4.
4.1.
4.2.
4.3.
4.4.
18
18
19
22
23
26
26
29
36
36
47
64
67
69
69
72
74
vi
Table ofCootents
6.
Applications of Calculus of Residues to Special Functions
6.1. Polygamma Functions of ArbitraIy Order
6.2. A Connection Between the Exponential and the Gamma Function
6.3. Residues of Some Functions Related to the Gamma Function
6.4. Some Integrals Involving the Gamma Function
80
80
82
84
87
93
93
95
100
101
103
108
110
7.
7.1.
7.2.
7.3.
7.4.
7.5.
7.6.
8.
(DRAGAN S. DIMITROVSKI)
110
112
114
9.
133
10.
140
147
(MIODRAG S. PETKOVlC)
11.
117
128
(GRADIMIR V. MILOvANOVlC)
11.1. Introduction
11.2. Orthogonality on the Semicircle
11.3. Existence and Representation of 1rn
11.4. Recurrence Relation
11.5. Jacobi Weight
11.6. Symmetric Weights and Gegenbauer Weights
11.7. The zeros of 1rn(z)
147
149
12.
162
150
155
156
157
159
Table of Contents
13.
vii
167
(DRAGIIA MrrROVlC)
13.1.
13.2.
13.3.
13.4.
13.S.
13.6.
13.7.
167
167
170
171
173
176
180
Name Index
187
Subject Index
190
Preface
Volume 1, i. e. the monograph The Cauchy Method of Residues - Theory and
Applications published by D. Reidel Publishing Company in 1984 is the only book that
covers all known applications of the calculus of residues. They range from the theory of
equations, theory of numbers, matrix analysis, evaluation of real definite integrals,
summation of finite and infinite series, expansions of functions into infinite series and
products, ordinary and partial differential equations, mathematical and theoretical
physics, to the calculus of finite differences and difference equations.
The appearance of Volume 1 was acknowledged by the mathematical community.
Favourable reviews and many private communications encouraged the authors to continue
their work, the result being the present book, Volume 2, a sequel to Volume 1.
We mention that Volume 1 is a revised, extended and updated translation of the book
Nau~na knjiga,
Belgrade in 1978, whereas the greater part of Volume 2 is based upon the second Serbian
edition of the mentioned book from 1991.
It seems reasonable to suppose that all those who found Volume 1 useful will also be
interested in this book. Some parts can be used for both undergraduate and graduate
courses, and like Volume 1 the book as a whole is suitable for postgraduate courses. It
ix
x
will also be of interest to researchers in complex analysis, physicists and engineers and all
those who apply complex functions.
The authors wish to express their appreciation to Kluwer Academic Publishers, well
known for their high quality productions, and particularly to the publisher Dr. D. J.
Lamer, as well as the series editor, Professor Dr. M. Hazewinkel, for their most efficient
handling of the publication of this monograph.
Belgrade,
27 January 1993
St. Sava's Day
Contents of Volume 1
1. Introduction
3. Contour Integration
24
51
108
211
253
286
308
323
343
(ISBN 90-277-1623-4)
Chapter 1
Introduction
1.1. Organization and References
1.1.1. This book is divided into chapters, each chapter into sections, and many sections
into subsections. The numeration of theorems, definitions, remarks, examples and
formulas runs continuously throughout a subsection, or a section which does not contain
subsections.
There are many cross-references in the book. If a theorem referred to belongs to the
same subsection only its number is given, while if it belongs to another, the numbers of
the chapter, section, and subsection are given. So, for example, 3.1.5 means Chapter 3,
Section 1, Subsection 5.
1.1.2. Bibliographical references are quoted after each subsection if it exists, or after each
section if it does not.
The titles of papers and books written in English, French, German, or Italian are
quoted such as they are. All the other titles are translated into English and the original
language is given in parenthesis.
There are many references to our previous book
D. S. Mitrinovic and J. D. Keckic: The Cauchy Method ofResidues. Theory
and Applications. D. Reidel Publishing Company, Dordrecht Boston - Lancaster 1984;
in further text it will shortly be referred to as Volume 1.
_ &7r(1
&7r
- -e-mR) <-,
m
Chapter 1
Ll.l. The term 2m in the first formula at the top of p. 117 should be replaced by 2n:
Hence the formula should read
e'"
21r
../2,
- -..e
J-chxd x =l+e
+<0
-eo
1r
= -trI'
cos2
1.1.3. A minor misprint is present in the formulation of Theorem 5 on p. 133: the word
its was misprinted as it. Namely, the second line of that theorem should read
was, in two places, misprinted as 2m. The correct definition of the contour C is as
n
follows
O~argz~ 21r}
IO~lzl~R,
21r
argz=-}.
1.1.5. The value of the second integral in 5.3.4.9 on p. 157 is incorrect. The correct result
is
+J"'cosax 1 d
1rcha
v.p. ocosbxl+~ x= 2chb'
1.1.6. The term e'C!l+IR+It)" in the integrand in formula (7) on p. 164 should be replaced by
e i C!l+IR+il)"/2. Also, an absolute value sign is missing. Hence, the first line offormula (7) on
p. 164 should read
Jf(z)dz
DB
= Je-.JW(R+IR+It) -1 dt.
.JW
t/(R+IR+It)"/2
1.1.7. The proof of 6.7.6 onp. 246 holds provided thatxz *ki.
1.1.8. The integrand in the formula (2) on p. 287 has a spurious factor
should be replaced by (Z2 + l); i. e. the formula (2) should read
+'"
-logF(z)=K+logz---2
y
dy.
dz
2z
0 (z2 + y2)(e2q -1)
(Z2
+ y2)2 which
Introduction
1.2.9. There is a misprint in Reference 1 on p. 349. The pages of the quoted article are:
157-160.
No
R+
R~
dgP
intr
for a closed contour r, int denotes the finite (internal) region bounded
byr
for a closed contour r, ext denotes the region bounded by
which
contains the point 00
ext
OG
v. p.
tf(z)dz
r
and
is
f(z) dz.
1.3.2. The reader is assumed to be familiar with complex analytic functions, their
derivatives and integrals, with Laurent's expansions and singularities.
Volume 1 begins with the definition of the residue. In order to make this book selfcontained we repeat some basic definitions and theorems relevant to residues.
THEOREM 1 (Cauchy)./f f is a regular function in a simply connected region G, then
tf(z)dz = 0,
r
Chapter 1
(1)
REMARK 1. In view of Theorem 1, the region G and its boundary can be replaced by any simply
connected region G I which is bounded by a closed contour
r..
Suppose that
f(z)
+'"
= LA,,(z-a)"
in the region {z 0 <Iz - al < r}. Since Laurent's series can be integrated term by tenn,
from (1) we conclude that
(2)
Resf(z)=A_1
r-"
REMARK 3. In Definition 2 the integral is taken along negatively oriented contour rsince in that
case the neighbourhood of 00, i. e. the region extr, remains left from the contour.
THEOREM 2.
We have
r='"
Resf(z) =O.
r="
REMARK 4. The above theorem does not hold for the point
fez)
=1 + 1/ z
is regular at
00,
00.
dermed by
..
'"
a (:,t 00) is a pole of order k of the function f, then
Introduction
Res/(z) =
Z=8
1
d k- 1
lim--.::l(z-at I(z).
(k-l)1 H 8 dz
f/(z)dz =2ni:t~es/(z).
1-=1
Z-G"
Chapter 2
Evaluation of Residues
2.1. Suppose that g and h are regular functions at z = a and that z =a is a zero of order m
of the function h, i.e.
(1)
L a,,(z-a)",
+~
g (z)=
,,=0
L b,,(z-a)"
+~
h(z)=(z-a)m
,,=0
(2)
Res!(z) =_1
I hIho
I
hom!
Z-Q
0
ho
hI
h2
i hm-
hm- z
... 0
0
0
hI
ao
al
Qm-I
This last formula can be written in such a way that the coefficients
from (1) do not figure in (2). Indeed, we have
(3)
Ok
g(k) (a)
=--,
k!
},(m+k) (a)
hk = - - - - - .
(m+k)!
ao,
...
(k =0, 1, ... )
and those values can be substituted in (2). This form of (2) can be found in [1], but with
an error. The table on the next page is also given in [1].
REFERENCE
(E,,=-~
z=z"
Evaluation of Residues
Type of
sin2UIaritv
Condition
fez)
Resf(z)
6=.
g(z)
h(z)
fez)
Removable
Simple pole
lim(z-a) fez)
g(z)
h(z)
g(z)
h(z)
Z-+<I
~o
=0
g(z)
~o
g(z)
h(z)
Second order
pole
fez)
z-+a
h(z)
g(z)
h(z)
g(z)
h'(a)
g(") (a)
(m+I) h(-J) (a)
Simple pole
Second order
pole
Second order
pole
g(a)
Simple pole
h(z)
g(z)
(z-a)2
Z-+<I
h"'(a)
3--.
Pole of
orderm
Pole of
orderm
Pole of
orderm
11m
Hit
h"'(a)
F(..-J)(z)
(m-I)!
F(m-I)(z)
z-+a
(m-I)!
F(z)
11m
, where
h(z)
Residue is given by
(2) and (3)
Proof Since
(2)
g(z/;.)=f(z"t;)=g(z)
Suppose that Fo is a contour around Zo such that there no other singularities in int Fo,
and let ~ be the contour obtained by rotating Fo around the origin through angle 2k1r.
n
The contour ~ encloses the pole zk. Hence
(3)
21t1r.
8
(4)
Chapter 2
=~.g(t~)~dt
2m'1
Ii
since t describes To in positive direction, as z describes lie in that direction. Applying (2)
we get
(5)
Z. Dokovi~.
l.l.lfthe function/has a pole of order k at z = a, the standard formula for the evaluation
of residue of/at z = a is, of course,
(1)
dk -
(k--l)!
Z-+lI
dZ k -
dp-I
dZp-l
where p is a positive integer> k is often more useful, particularly when / has the form
f(z)=
g(z) ,
(z-a)P
Z ~f(z)=-----
Z3
(z+ 2)
(e
z -I )
. d2 Resf(z) = -1 hm- ,
z=o
2z-+odz2
z+2
= 3, we get
Evaluation of Residues
(4)
. d(eZ-I)
Res/(z)= hm- - - - .
z--+odz
z=O
z(z+2)
1. R. C. Gupta: A New Fonnulafor Finding the Residue. Math. Educ. Sect. A 15(1981),32.
2.4. The question of practical computation of residues of some functions which have
essential singularities at z = 0 is considered in [1]. Such functions are, for instance,
Z f-+
zn exp g ( : ),
Z f-+
Z f-+
zn exp (az + :) ,
where n is a positive integer and g is an entire function. The residues are obtained by
expanding the relevant functions into series in a neighbourhood of z = 0, and they are
expressed in tenns of numerical series which can often be summed. We give a few
examples.
+>
EXAMPLE 1. If f(z)
=:~:>.z, then
.0
ResJ(z) eelz =
z=O
+OOacn + 1
2:
_n_ _
n=O (n + I)!
EXAMPLE 2. If a ~ 0, b ~ 0, then
(I)
1
Res - - = -eelz
z=O
a+bz
)n cn+'
2: ----= -b - -h e- bela;
a n=O
(n+ I)!
+00
--;
EXAMPLE 3. We have
(2)
2:
b)
+ 00 an bn
( b
ReszP-'exp ( az+-- =bp
------=- z=O
z
n=O n! (n + p)!
a
)P12 J (2V-_ab )
p
10
Chapter 2
REMARK. Using computers the desired results can be obtained in a relatively short time. For
Res
(2
Sz -
1)
= -
(3i+1)z
z=o
2 V2i I, (2 iv'!O)
---:----'-----:--
15 y'Se i i+Sy'Sei
in 5733 milliseconds.
REFERENCE
1.5. We have
Res+ (2Z2 -1) &2 (arctgz)2 = 11: (1
z=oo
-~) .
e
= -Z21
has critical singularities at z = i and z = -i. We take the branch of arctg for which
z-+'"
nnw'
2
2
3
- -W + ...
5
If Fl is the branch of F which corresponds to the chosen branch of arctg. then any
other branch ofF has the form
11
Evaluation of Residues
g (z)
Z2
Z2
so that
ZH
Z~ ( :
+ &),
Resh (z) = O.
z= co
I'
Furthermore,
Figure 2.5.
(r-
For the contour rwe choose the contour shown on Fig. 2.5. As the
radii of the circles around -; and ; (and also the distance between the
lines In and I'n~ tend to 0, we have
f-
g (z) dz -+
z2
The function
n'
12
Chapter 2
and so
Res F) (z) = Resg (z) - Res h (z)
z=oo
z=oo
z=oo
=_I_.J:g(Z)dz=-I-.2n2;(1-~)=n(1_~)
2n,j
2nl
r-
z-o
Res ez log -
z-b
z=oo
= e" - eb.
ez log ~--:.~
z-b
bas critical singularities at z = a and z = b. All the branches of F are uniform outside a
circle which encloses the points a and b. We consider the branch of the logarithm which
is 0 for z = 1. If we putz = l/w when w is in a neighbourhood of 0, then
z-o
I-ow
onw
log--=log---=
-aw (OW
1 + - +a2-w+ ... +-+
... )
z-b
(2)
I-bw
ntl
)
bnw
bw b w
+ bw ( 1 +---+-+
... +-+
....
2
n+-1
zH F.(z) + 2ik!re z ,
k eZ. Hence
all the branches of F have the same residue at z = 00, since the integral 2iknJ eZdz taken
along a circle of sufficiently large radius vanishes.
The residue of Fl at z = 00 is equal to the coefficient of w, with changed sign, in the
product ofthe series (2) and the series
I
.-
2! w2
13
Evaluation of Residues
This coefficient is
(b2!
n
an
a
-a ( 1 +-+
... +--+
... ) +b 1 +-+ ... +._-+
... )
(n + I)!
2!
(n+ I)!
=(-e"+l)+(eb-I),
and SO
ResFl (z)=et'-eb.
Z= 00
Z=ao
REMARK. The equality (1) is given in Volume 1 (2.1.2, p. 13) without a proof.
REFERENCE
1. F. Tisserand: Recueil complementaire d'exercices sur Ie calcul infinitesimal. deuxieme edition,
augementee de nouveuax exercices sur les variables imaginaires par P. Painleve. Paris 1933.
2.7. Suppose that the function / is regular inside a closed contour which encloses the
points at. ... , am and suppose that nt. ... , nm are positive integers. If the function F is
defined by
F(z)
then
(1)
2 ResF(z)=P(w),
k= I Z=Qk
=r (al ),
p(n,-l) (a 1) =f(nl-l) (a 1 ),
r (am),
P' (am) =
....
14
Chapter 2
REMARK. It was implicitly supposed that w ;t at (k = I, ... , m). If w =at for some k = I, ... , m,
then the swn on the left hand side of (I) is o.
REFERENCE
f(z+a)=
"
+ ..
k=1
k=O
2 A_kz-k + L Akzk
in the region {zl'i <lzl<r1 }. Then z=a is a pole of order n of the function / and
Res/(z) =A_I. Indeed, if we put z + a = t, this assertion becomes trivial. However, in
2=11
certain cases the expansion (1) can be useful. So, for example, for Jacobi's elliptic
functions sn, cn, dn, we have
.
1 +k2
kz
6k
sn (z + IK') = - + - - z + 0 (Z3),
cn ( Z+lK') =
2 k2 - 1
-k-;+~-Z+
0 (z,
3)
and we conclude that those function have a simple pole at z = iK', and
1
Res sn Z = -
z=iK'
Res cn Z= - -, Res dn Z = - i.
k ' z=IK'
z=IK'
REMARK 1. Jacobi's elliptic functions also have a simple pole at z = 2K + iK' and
Res
z=2K+iK'
sn z = - - ,
k
Res
z=2K+IK'
cn z = - ,
k
Res
z=2K+IK'
dn z = i.
REMARK 2. The functions sn, cn, dn and the constants k, K, K' are dermed by the so-called theta
functions St, lh., ~, 84. The constants k, K, K', are connected by the equalities
Ilk
K'
=}' V(u
I
-1) (1 _k 2 u2 )
duo
15
Evaluation of Residues
2.9. Let P",(z) =ao +a1z+ ... +a",z"', Qn(z) =bo+b1z+ ... +bnzn where a".bn
the rational function/defined by
-:f;
0 and consider
= p.. (z).
/(z)
Qn(z)
= {zllzl= R}
where R is large
IIfn>m+ 1 then
2 If n
m + 1 then
3 If n < m + 1, let T(z) be the remainder when Pm(z) is divided by Qn(z). Then
i:r.l:
.1:=1
= r<n-I)(z) .
(n -1) Ibn
Proof We have
i:r.l: =~t/(Z)dZ.
(1)
21t1 C.
.1:=1
laol+ladR+..+la",IR'"
~ -lin
bn_1 R - ... - bo ~ 0,
2 tr bn R - lin-III
as R
i>.I: =0.
.1:=1
2 Let n
fez)
m + l. Then
a", + fez) _ am
b..+1
b",+l
16
Chapter 2
b_ 1
c.
+0(1)
(R~+oo),
(n =m+l).
3 For n < m + 1, there exist polynomials Sand T such that
P,.(z) =S(z)Q,,(z) + T(z),
where T(z) =Co +CIZ+ ... +c~, C. t< 0 and s < n. Hence, by the previously established
results we have
~:rk =0
if n > s+ 1,
k=1
i:rk = T-,,-I)
(z) .
(n -1)!b"
k=1
2.10. Up to now we have been concerned with the evaluation of a residue of a given
function at a given point. We now consider, in a wayan inverse problem. Namely, instead
of determining Res/(z) for a given function/and a given point z = a, we shall now, for
r-tl
a given A
z=a isA.
Let (/J be the set of all complex analytic functions for which z =a is a regular point or
an isolated singularity, and let A be a given complex number. The general solution of the
equation in/
(1)
Res/(z)=A
z=tJ
is given by
17
Evaluation of Residues
Z-Q
Z-Q
(2)
where T
E (/J is
z=a
arbitraty.
Z=Q
=A,
which proves that (2) is a solution of (1).
Conversely, let
10 E (/J
10
is such that
Resfo(z) = A. We shall prove that a suitable choice of the arbitraty function Treduces (2)
Z=II
to fez)
=-
Z-Q
=Z-Q
+10 (z) -
Z-Q
Resfo (z)
z=a
A
+10 (z) --=/o(z).
Z-Q
and this shows that (2) contains all the solutions of (1).
REMARK. This is a special case of a more general result, proved in [1].
REFERENCE
1. J. D. K~kic: Reproductivity of Some Equations of Analysis II. Publ. Inst. Math. (Beograd)
33(47X1983),109-118.
Chapter 3
1=_I_I./'(z)dz
2m ~ /(z)
,
where c= {zllz-al= r}. Ifz = a is a zero of order n of the regular function/and if there
are no other zeros inside C, then 1= n. Similarly, if z = a is a pole of order p of the
function/which has no other poles or zeros inside C, then 1= -po
By analogy with this, B. Gavrilovic [1) defined the order of an essential singularity as
follows: If z = a is an essential singularity of the function / which has no other
singularities or zeros inside C, then the number I, defined by (1), is called the order of
that essential singularity. In connection with this definition we mention that Weierstrass
proved that the value of I is an integer, even if z = a is an essential singularity.
Gavrilovic proved the following result: If the following conditions are fulfilled: 10 / is
2/is analytic in the region int rwhere it
regular and does not vanish on the contour
can have a finite number of poles and essential singularities, then
r;
(2)
where N, P, E are the sums of orders of all zeros, poles and essential singularities off,
respectively. Formula (2) generalizes the well-known Cauchy formula for the so-called
principle of the argument.
REMARK. There exist other generalizations of Cauchy's fonnula (and also of Jensen's fonnula) of
similar kind. See, for example, Volwne I; 4.2.3, pp. 58-69.
Gavrilovic also proved in [1) the following theorem. Suppose that the function / is
regular in the region int
where r is a closed contour, and let the function g have in
int r isolated singularities (poles or essential singularities) ah ... an. Suppose that in a
r.
18
19
+'"
g(z) =G(z)+
= LA..t(z-akY;
1'=0
1'=1
A..t.,
(z-ak )
+'"
-.t/(z)g(z)dz= LLA..I,kB"".
2m r
1'=1 k=1
Formula (3) contains, as a special case, formula (2) and also the well known Cauchy
residue theorem.
REFERENCE
1. B. Gavrilovic: On Residues o/Uniform Functions (Serbian). Rad JAZU 139 (1899), 29-39.
and let
for
l~k~n
is an interpolation polynomial for the function/, and in certain cases it can be taken as an
approximation for f
If the polynomial Pn of degree n is defined by
It is often taken for granted (particularly when this formula is used in physics) that
n -+ +00 implies that Pn(x) -+ I(x) uniformly on I. However, this need not be true, not
even for functions which are analytic in a region Dee which contains the real interval I.
This fact is known as Runge's phenomenon.
Indeed, if a> 0, the function/defined by
20
Cbapter3
1
f(z)=-,
Z2+
a2
is meromorphic in the complex plane and has two simple poles at z = ;a and z = -;a.
Suppose that R is a positive number such that
for
l~k~n,
dz
l'
where x:l: ak (1 ~ k ~ n) and x:l: ia. It is clear that the absolute value of this integral is
not greater thanAIRn+2 whereA is a constant. If the functiong is defined by
g(z)=
/(z)Pn(x) ,
(z -x) p" (z)
then
Resg(z)=
Resg(z)=f(x),
Z=X
z=ak
Resg(z)= . . Pn(x)
2 /Q (Ia - x) Pn (Ia)
z=;a
(l ~k~n),
/(ak)Pn(X)
(ak-X)Pn'(ak)
Resg(z)=
z=-ia
..
Pn(x)
2 la (Ia + x) Pn ( -Ia)
Res g (z) +
z-x
k=t Z=Dk
z=;a
z=-Ia
i.e.
f(x)-P (x)+
(1)
Pn(x)
2;a(;a-x)Pn (;a)
Pn(x)
2;a(ia+X)Pn(-;a)
O.
2k+l
2m
--
and SO
(O~ k~n-l)
.
for the pomts 01' ...
an'
Then
21
Pn (x)
2ia(ia-x)Pn (ia)
2ia(ia+x)Pn(-ia)
~~----+----~~----
Pn (x)
2 iaPn (ia)
(I
I)
ia - ; + ia + x
=-
Pn (x)
x 2 + a 2 Pn (ia) .
and we shall show that for sufficiently small a the difference 1 (1) - Pn( 1) does not tend to
1 34m - 1 _ r ( 2
p (1)=----.. ' - - " ' V 2 -
2m2m
2m
)n ,
J
I
2
10gb < log- = log 2 -1. From (3), (4) and (5) follows that
(c' *- 0 constant)
+00.
22
Chapter 3
REFERENCES
Proof Let C be a circle with centre at a and radius r, such that there are no other
zeros of fin the disc {zllzl < r}, except the simple zero at a, and let
m = minlf(z)l.
zEC
Now, let K be a circle concentric with C such that If (z)1 < m for z E K, and hence for
z E int K. If u E int K is arbitrary, then If (z)1 > If (u)1 for all z E C, and so by Rouche's
theorem the functions z ~ f(z) and z ~ f(z) - f(u) have equal number of zeros in the
disc int C. This means that the equation
f(z) - f(u)
=0
has only one simple root, namely z = u in the disc int C. Therefore, by the residue
theorem
fg(z)
f'(z)
dz = 2niResg(z)
f'(z)
,
f(z)- f(u)
%=11
f(z)- f(u)
implying that
( 1)
g(u)
1 i.g(z)
=2ni~
f'(z)
dz
f(z)-f(u)
However,
g(z)
!'(z)
f(z) - f(u)
= g(z)!'(z)
f(z)
1
1- f(u)
f(z)
(2)
g(u) =IAk(f(u)Y,
.1:=0
where
= g(z)!'(z) f(f(U)k
f(z)
k=O f(z)
23
A =_I_I.g(z)!'(z) dz
k
2m ~ (f(Z),,+l
(3)
(k =0, I, 2, ... )
Formulas (2) and (3) define the expansion of g into a series with respect to the powers
of f. It is usually called Bfirmann's series.
For k > 0 the expression for the coefficientAk can be transformed as follows:
2A ( I g(Z)
I I. g'(z) d
m k = k (f(z)" c + k ~ (f(z)" z
_!I. g'(z) dz
- k ~ (f(z)" ,
and putting
!(z) = z-a
F(z)
Ak
REMARK. It is not difficult to see that BOrmann's expansion is equivalent to Langrange's expansion
(see Volume 1; 4.4.2, pp. n-79). BOrmann's result was published in 1799.
REFERENCE
1. L. Cakatov: An Introduction to the Theory ofAnalytic Functions (Bulgarian), Sofia 1972, pp.
244-245.
{zllzl~r,
that it has simple zeros 'ie/It (k = I, ... ,n) in the region D bounded by the semicircles
R>r,
and
the parts of the imaginary axis joining them, and that it has no zeros on the contour. Then
(I)
L"(I)
-\ cost =trR1- JIO~!(Rell)lcostdt
r R
/2
k-1
-_12
+~
j(~-~)logl!(;Y)!(-;Y)ldY+O(I).
2tr Y
R
r
24
Chapter 3
=~
.((~+...!....)IOgf(Z)dZ.
2m '! Z2 R2
aD
This integral breaks into four integrals (along semicircles c and C, and two parts of the
imaginary axis). The integral along the small semicircle is bounded. On the negative
imaginary axis we have Z = -iy and its contribution is
1 R( --1
1) logf(-iy)dy.
-J
2n y2 R2
r
1 1112 (-2it
1)
1 1112
-.
J
~+-2 iReltlogf(Re,t)dt=- J(cost) 10gf(Relt)dt.
2m -1112 R
R
nR
-tr/2
1 R( --1
1 ) logf(iy)dy.
-J
2n y2 R2
r
Adding up those integrals, and taking the real part of I, we obtain the right hand side of
(1).
aD
2m aD Z
f(z)
As we describe the contour, logf(z) increases by 2mn. The integrated term is therefore
purely imaginary. Since
~(rk!'~ -r~~ )
and taking real parts, formula (1) follows.
This is Carleman's theorem [1]. See also [2].
REMARK l. Ifwe replace the region D, i.e. the left half-annulus by the upper half-annulus
G=
{z I r ~ Izl ~ R,
0 ~ argz :s;
~"(_l
-~)IOgf(Z)dZ,
2m L R2 Z2
tr}
25
i_I
1 "
1 R( 2"--2
1
1 ) 10gl[(x)[(-x)ldx+O(I).
=-Jlogl[(Re")lsintdt+-J
nR 0
2n, x
R
1log [(z) dz
IoFR
I'
121<
ft
R
..
R
-Jlog[(Re u ) dt =logl[(O)1 + 1)og-11- LIOg-1
2n 0
t=1
at
i=1
bt
where ak' bk are zeros and poles, respectively, of the function[such that 0 < Ia~ <R, 0 < Ibki <R.
See Volume 1; 4.2.3, pp. 58-59.
If this disc is replaced by the left or the upper half-annulus, then starting with the appropriate
integrals (2) or (4), we obtain Carleman's fonnulas.
Instead of the disc or the half-annulus (left or upper), Mitrovic [3] took the closed half-disc
H = {zllzlS R, Irnz ~ O} where the function[is meromorphic and by considering the integral
R for z
(RIl)
lR
-J(P(x)logl[(x~
+F'(-x)logl[( -x~)dx- Irn - J e"P(Relt)log[(Re")dt
2no
2no
= Irn(tatF(at )- tPtF(bt)}
and
REFERENCES
1. T. Carleman: Ober die Approximation analytischer Funktionen durch lineare Aggregate von
vorgegebenen Potenzen. Arkiv tOr Mat. Astr. Fys. 17(1922), No.9.
2. E. C. Titchmarsh: The Theory o[Functions. London 1964.
3. D. Mitrovic: On an Equality [or Integrals (Serbian). Glasnik Mat.-Fiz. Astr. (2) 6 (1951), 193-
200.
26
Chapter 3
_1 I(t)dt
2m c t-z
(1)
and let F+ and F- be the functions represented by the integral (1) in D+ and D- ,
respectively.
The function F- + 1 is the analytic continuation of F+ in the region D- through Co'
Indeed, let C* be the contour formed from C when Co is replaced by a part of the
boundary of D (in such a way that the orientation is preserved). See Fig. 3.5. Clearly, this
part of the boundary of D is also a part of the boundary of D-. Let L - = fJD- . Then, for
any zED such that z II: Co we have
l(t) dt
2m c.t-z
(2)
=~
I(t) dt
2m c t-z
+ ~ l(t) dt.
2m r t-z
Denote the integral of the left hand side of (2) by F(z). Since C* is not in the region
D, the function F is regular in D.
If z e D+, we have F(z) = F+(z), because the first
integral on the right hand side of (2) is, by definition,
equal to F+(z), and the second integral vanishes by
Cauchy's theorem.
If z e D-, we have F(z) =F-(z) +I(z), because
the first integral in (2) is, by definition, equal to
F-(z), and the second is equal to I(z), by Cauchy's
integral formula.
The function F is therefore the
continuation of F+ onto the whole region D.
analytic
Figure 3.5.
REFERENCE
1. M. A. Evgrafov: Analytic Functions (Russian). Moscow 1968, pp. 245-246.
G=
27
e"
E(z) = - . J-dt
2m c t-z
(C= aCT).
We shall determine the analytic continuation of E onto the whole complex plane and we
shall establish an asymptotic formula for z H E(z) as z ~ 00.
Since the function E is defined by a Cauchy type integral, we obtain an analytic
continuation of E onto the whole plane by the method of 3.5. If we denote by J(z) the
value of the Cauchy integral (which defines the function E for z Ii!: G u C) the analytic
continuation of the function E in G is given by
E(z)=J(z)-e"
(z eG).
1
t-z
1
z
t2
--=----+--::--Z2
Z2(t-Z)'
we easily get
Co
c1
1
=--+-+-F(z),
z Z2 Z2
J(z)
(1)
where
Co
J 'dt,
2m c
1
= --.
e'
C1
Jte" dt,
2m c
e"
t2
2m c t-z
1
-~
F(z) = - . J-dt.
We shall first show that F(z) = 0(1) as z -+ 00. Let C 1 be the part of the contour C
which lies inside the circle It - zl < 1 and let C2 be the remaining part of C. We then have
F(z)
where
t2
e"
t2
e"
1';(z)=-. J-dt.
2m c; t-z
First of all,
In order to find an upper bound for F I , denote by w the point of C 1 which is nearest to z
and then
28
eO-
But
(1
1J
F.(z)=-. J-dt+J
2m C, t-z
C,
e;-o- dt .
t-z
l -11 ~ M
It-zl<I),
e;-ot-z
and
-1
(if
Chapter 3
2td c, t-z
z-b'
2td
where a and b are the end points of C l (Le. they are the points of intersection of C l and
the circle It - zl = 1). Since
dtl ~ 1.
I~ J-l
t-z
2m
lee-I ~ 1.
=0(1)
(z
- Co
00).
= 1. Indeed, we have
1f;
2m c
2m
1</+<0
=-. Jeo'dt =- .
11;
-1</+<0
_1<;
since the function t H eo' is periodic with period 2td. The first and the third integral on
the right hand side annul each other. In the second integral we introduce the substitution
et = W and we obtain
1
-Co = - 2td
e
e
J -dw=-Res-=-l.
Iwl=l
w=0 W
= 1,
F(z) = 0(1)
(z
I(Z)=;+O(~)
(z~oo).
00)
29
E(Z)=;+0(:2)
1 . (1)
E(z) =--e
+0 -Z2
Z
REMARK.
(z ..... oo,
z~G)
(z ..... 00,
z e G).
Co
c.
c. + ~ -12)
I (z ) =-+2"+
... +-.
zz
where
z'"
f.
c = -1- t eold t
2m c
0).
Indeed, if
J(z+ru')=J(z)
implies
~J(z+ru')=~J(z), i. e.
dz
dz
l(z+ru')=/(z).
Hence, the period 0)' of J must be of the form mO) (m e Z). However, since the derivative
of z ~ J(z + 0) - J(z) is zero, we conclude that
J(z+O)-J(z) = 0
(0
constant)
30
Chapter 3
and SO
0)
Since u = eikz, all those values of z which correspond to one value of u differ from each
other by an integral multiple of 0). If J is a uniform and a periodic function in z, then the
function J 1 is defined by
~ F(u)
When J is uniform and periodic, it is a rational function in eikz, and this will take
place if and only if all the residues of u ~ F(u)/u are zero.
We now apply this result to the case when
a+bv
u
c
u3
R(u, v)=--+-+d,
i.e.
f()
z
c
d
= a+bcosz +--+
3
sinz
sin z
(a, b, c, dER).
.
'b +-1C=
1.
0
la+,
2
31
. b - - 1.
-la+l
IC= 0
2
sm 3 z
J(Z)=J(_l _ _
2 )dz=~.
sin z
sin 3 z
sin2 z
REFERENCE
oP
oQ
oy
ox
-=-
holds if and only if ResP(x,y) does not depend on y for any pole a of the function
~=tI
xI4P(x,y).
This result was proved in [l). The following interesting application is also given there.
If the function P is defined by
+F
P ( x,Y) -_ Ax2 + 2 Bxy + Cy2 + 2 Dx + 2 Ey
--,
x2-2xy+ 1
2 (a-y)
and, on putting a =y + ~y2 -1, this residue must not involve y. Equivalently, this residue
must not depend on a when we set
Chapter 3
32
02
(3)
+1
Y=--
20
20
40 3 -40
8D+4E
4
4E
-4'
4B+2C+4F=O,
C=O,
implying
A= -B=F
C=o,
D= -E.
x2-2xy+l
Q (x, y) = 2 -2Dx
2
1'
x - xy+
and
P (x, y) dx+ Q (x, y) dy=d(Ax+D log (X2 - 2 xy + 1.
REFERENCE
3.7.3. Let a l' ... , an be relatively prime positive integers and let M> 0. Denote by N(M)
the number ofintegral positive solutions of the equation in xl' ... , Xn:
i.e.
33
Proof This problem is easily reduced to the problem of showing that the number
C(M) of integral nonnegative solutions of (I) is of order (2). Indeed. if M> 0 is large
enough and if we put
then the number of integral positive solutions of the equation (I) is equal to the number of
integral nonnegative solutions of the equation inYI' ... ,Yn :
(3)
The equation (3) is obtained from the equation (I) by means of the substitution xk =Yk+1
(k= I, ... , n).
The number C(M) of integral nonnegative solutions of (I) is equal to the coefficient of
zM in the power series expansion off, where
"
/(z)=n
(/ z
-0
k=11-zk
/< 1).
"
IT
(I
+zOk
+Z2ok + ... )=
k=l
+00
+ ..
kl=l
k n =l
2 ... 2
zO,k,+ .. +a"kn
2:n; i j
zM+1
g (z)=
f(z) = _ _ _ __
ZM+l
ZM+l
IT (l_z0k)
k=l
The function g is analytic and has a finite number of poles on the circle {z Ilzl= I}.
The poles are
34
Chapter 3
Resg (z)= 0
(forM~
z=oa
1)
(4)
:r
2: Res g (z).
k,J z=zkj
Since the positive integers aI' ... , an are, by supposition, relatively prime, the poles at
'* 1 are simple, and so
= Zkj
(5)
dn -
z=1
(6)
= ___
1 ___1_ Mn-I + k=1
n
TIOk
Ok
2(n-2)!
(n-l)!
Mn-2 +
3(
k=1
Ok)2 -
k=1
24(n-3)!
Ok 2
Mn-3 + ....
k~1
(7)
lim
C(M)=
n
M->+oo Mn-I
(n-l)!
IT Ok
k=1
... ,
REMARK 2. This text was written according to [I). The same problem was considered in [2],
Problem 27, and the complete works of Laguerre (tom I, Paris 1898, pp. 218-220) were quoted as
references. The problem is solved in [2] without the use of residues.
REFERENCES
3.7.4. If f is a function of positive real part, regular in the unit disc D, with Taylor
expansion
+'"
f(z) = 1+ Lavz v
v=J
1.
(z ED),
35
Proof. We have
I"
=~t(2-Z" -z-")z-lf(z)dz
2maD
1 2
=- {l-coskt)f(et')dt
fro
(k
1).
(k
Hence,
(k
1)
implying
2-Rea"
J{l-coskt)Ref(et' )dt ~ 0
Ib
=-
fro
i. e.
(k
1).
(k
1)
1),
Chapter 4
Jci:xr
+00
L (n) =
dx
~-I
( _ 1)n/2 n in 2
-'----'--2: (_I)k ( kn ) (n-2k)n-l,
2n(n-I)! k=O
(1)
n even
n-I
2: (- I)k (n)
(n - 2 k)n - II
k
(_I)(n-I)/2nin-12
-'----=----2n(n-I)!
I,
k=O
odd.
(-l)"/22,,-1 sinn t=
22: (-l)n(n)COS(n-2k)t+~(-1)"/2(
k=O
n-I
2 k), n odd.
k=O
=Fn(z) + F n_3(z),
n
II'
F" (z) =
2: (- l)k ( n) ei (n-2 k) z,
k=O
n' = {
2-1,
n-I
where
n even
- 2 n odd
36
n), n even
nl2
37
and where Pn-3 is a polynomial of degree n-3, such that the function
regular in the entire plane, except at z = 0, and
z-"/,,(z) is
(n-I)!
833M
REFERENCE
4.1.2. Boas and Friedman noted in [1) (see also Volume 1; 5.2.3.7, pp.120-l2l) that the
semicircular contour used for the evaluation of integrals of the form
Je I(x)dx can be
+'"
ix
Z H
Z-1I2 iz
Jsmx dx
+'"
Z-l
(1)
where k
= 112 for
= 0 for
Je I(x)dx,
+'"
ix
38
Chapter 4
where f is a rational function with the degree of the denominator greater than the degree
of the numerator.
The integral involved in the limit is not elementary and has to be estimated by some
device. The usual one is to apply the inequality sin t ~ 2thc.
The authors of [2] believe that it is simpler to break the integral in (1) into the sum
Jo
.,3
./2
The second integral causes no difficulty, and for the first, since 2cost ~ 1 for
OStSlr/3 we have
.,3
./3
R" e- Roinl dt S
o
The last integral can be evaluated explicitly and it is easily seen that it tends to 0 as
R -+ +00.
Another approach suggested in [2] is to give up the traditional use of circular-arc
contours. For the Fresnel integrals take the contour to be the square with vertices at 0, R,
R+iR, iR (initially indented at 0). The integrals along the segments [R, R + iR] and
[iR, R + iR] are
J(R +iy)"-'e'(R+IY)dy
R
and
f (x+iR)"-'e'(,,+IR)dx
o
"-I
But R2 + y2 ~ R2 and x 2 + R2 ~ R2, k -1 < 0, and the above sum does not exceed
R
= R"-I(l-e-R)
+ R"e-R -+ 0 0
as
R -+ +00.
+>
For integrals
bt
->
REFERENCES
39
4.1.3. If a, b are positive numbers and a, p real numbers such that 0 < a - p < 21r and if
z ~ /(z) is a :function of the complex variable z = re i9 which has finite limits
/(0) =lim/(z),
/(00)
.... 0
=lim/(z),
6"'1lO
within the sector from () = a to () = p, and is regular within and on the bounding lines of
the sector except for simple poles at cI' ... , cn within the sector, then
(I)
J
o
-2nif(Z-CII)/~Z)1 '
2=e.,
1'1==1
where x is real.
j /~Z) dz
P and the circles r = aro
bR"~
fJ
"'';,
IIIIP
lIIfI"
J!dz
,..P
bril' z
21lif(z-clI) /(Z)I
and r = aR,
>=1
2=',
'
Resf(z) .
.....
40
Chapter 4
.. ~
REFERENCE
{I
c+i..
-1
2ni
c-I ..
tz
e--dz=
z
for t>O
-1
2
for t=O
for t<O
c+iR
c+iR
t<0
c-iR
c-iR
Figure 4.1.4.
O~ 3tr if t> 0
2
and tr ~
2
O~ _ tr if
2
41
J
n
eCI
nl2
e-IR sin rp
I z'
R dqJ;;;;; 2 eCI
e-2 iR 'PIn
.--;--
Iz
ed
~ - - - . - 'J'l
- t(R-c)
as R -+ +00. Hence,
Je-dz= 1,
R dqJ
neW
c+iR
.
lIm
R-.+oo
-1-
lz
2n;
c-iR
1>0,
-n/2
!
i
~xp (c + ~ei 8) t R dO
Z
nl2
i; ; ; eCI
l
J
nl2
e lR cos 8
R dO
-n/2
nl2
n/2
n/2
e 2 IR 8/n dO
= ~ eCI (1 rR
e'R )
<!!... eCI -+ 0,
tR
as R -+ +00. Since the integrand is a regular function inside the contour of integration we
conclude that
Je-dz=O,
c+iR
.
lIm
-1
R-++oo 2ni
lz
c-iR
1<0.
J' -
c+iR
-1.
2nl
c-iR
1
z
1
.
.
1
R
1 -=--,
n
1
dz=-.(log(c+IR)-log(c-,R=-arctg-_2nl
as R -+ +00.
This completes the proof.
REFERENCE
42
C~4
"
... 1
t=1
(1)
then
b)
(2)
[r,
r= {zllzl=r,
RJ.
k=1
R......",
Indeed, if z E
II
F(Z)dZI = j(Ii:sin(a..
r
0 ""I
Re
k=1
0 ""I
k=1
because
43
Iff
1'=1
k=1
Since
sinl
for
I E [ 0, ;
J.
we get
/2
(2aR) dt=-, 1r (fiR
)
s 2- trJ
exp e -1,
R" o
aR"+1
1r
implying
lim JF(z)dz =o.
R-->+ ..
I.
RJ( " sina
~F(z)dz=
r F(Z)dZ+,:a x"
+ F(z)dz+
r
R
"",I
we obtain (2).
4.1.6. We have
J Vx
+00
cosx-sinxdx= -~'.
(1 +X2)
ev''2
x)( D
. cosbj:x)~
-;-dx
by
k-1
44
This function has simple poles at i and -i with residues
'I = Resf(z)
= ____
1 -
2eVi
z=;
(l + i);
= -
,e;-;;
2v 2
(1 - i).
Let r= {zllzl= R}, r= {zllzl= r}, d = {zlImz = 0, r S Re z:S R}, where R is large
enough so that the circle r encloses all the singularities off Let C = ru d u (-d) u r.
Then
f fez) dz= 2
7e
i ('I + '2).
-2.,r;(
) dx=2m(Ij+12),
o x l+x2
21
-~( hI . hI)
Jo coSX+ISIDXd
.,r;(I+~) xS +IC
.
Ii
Therefore,
+CIO
Ii
Figure 4.1.6.
eli
REMARK. If
J
I
D=
cos x
VX(I+X2)
X,
Jr
then the integral 1= f(z)dz can be expressed in tenns of a and b, where r is the contour
J !(z)dz-P !(z)dz=2nir"
ABCDA
we conclude that
J !(z)dz=-2nir
DEFA
45
JY
t
REFERENCE
4.1.7. Letfbe a rational function which may have only simple poles aI' ... , am on the
positive real axis and let zI' ... , zn be the other poles off, none of which is at the origin. If
p is a real number such that
lim zP+ 1 f(z)
z-+-o
we have:
+~
v.p.
xPf(x)dx= -
n.e-np;
2:
Res(zp f(z-ncotgnp
I akPResf(z),
z=ak
k=1
Z. then
+CO
v.p.
x Pft x ) dx= -
In
0;;;;argz<2n.
REFERENCE
=f(x);
46
chapter 4
30 fis analytic in the upper half-plane where it can have a finite number of isolated
singularities 01' .. , an;
[0, n);
+..
SO The integral
Then
(1)
J
o
f(x)log(x 2 +1)dx =
7ri(2f~esf(z)IOg(Z+i) -AJ.
"*
.1:=1
which takes the value rol2 at z = 0, is uniform in the plane cut along the imaginary axis
from -i to -ioo. Its singularities in the upper half-plane are the singularities 01' ... , an off
If CR is contour consisting of the real segment [-R, R] and the semicircle
{zllzl= R, Imz > O}, where R is large enough so that all the singularities lie inside CR ,
then
n
(2)
ff(z)log(z+i)dz = 27ri~~esf(z)log(z+i).
c.
.1:=1
.-a.
(3)
tr
-Jf(-x)log(-x+i)dx.
R
Since the functionfis even on the real axis and since for x> 0 we have
arg(-x+i)+arg(x+i) = 7r
then clearly
R
(4)
Jf(x)log(X+i)dx- Jf(-x)log(-X+i)dx
o
+""
Furthermore,
47
I<
~ Jlf(Rell)IR(logR + n)dt -+ 0,
o
asR -++00.
Hence, if R -+ +00, from (2), (3), (4) and (5) follows (1).
ExAMPLE. If fez) = 11 (Z2 + 1), all the conditions 10 - 50 are satisfied. This function has only one
singularity in the upper half-plane at z = i and
Reslog(z+i)
..,
Z2 +1
Also,
!.r
+..
A=
1
+1 dx=f,
n n) =nlog2.
-1.......;..::-~dx=ni
log(x + 1)
(-ilog2+--2
x2 +1
REMARK. Several theorems of this type were proved is Section 5.3.2. ofVolwne 1.
REFERENCE
1. D. D. Adamovic: A Collection of Solved Problems in Mathematical Analysis. Ordinary and
Partial Differential Equations and Complex Functions (Serbian), Beograd 1959.
(1)
0=0
"
Proof Since z = 0 is a simple pole of/. in a neighbourhood of that point we have the
expansion
+co
(2)
fez)
Suppose that r > 0 is sufficiently small, so that z = 0 is the only singularity off in the
region {zllzlS: r}. Integrating (2) along the arc of the circle {zllzl= r}, defined by
O~
48
C~4
(3)
8+6
8+6
+..
..+6
i=O
..
4.2.2. If a
o (cha-cosx)2
2n
eft..
a>O
(n-cotga),
a<O
sh2a
Jo (cha-cosx)2 dx=O.
26
smnx
The above equalities were proved in [1]. In the same paper the following integral
(a
* 0, b * 0,
c *0)
is also evaluated.
REFERENCE
4.2.3. Suppose that the function z = x+iy H F(z) = u(x,y) +iv(x,y) is regular in the
disc {zllzl~ I} and that F(z) is real for real z. If -1 < t < 1, then
Jo 1-2tcosf'+t
tsmf' 2 v(cos f',sinf')df'= n(F(t)-F(O.
2..
tsinf'
1=J26o 1-2tcosf'+t
2F(cOSf'+sinf')df'
after the substitution eitp =z becomes
1=-!.r. t(z2-1)F(z)
where
dz,
t(Z2 -1)F(z)
z(z - t)(1- tz)
49
implying
tsinqJ
tsinqJ
-----!.--=-2
o l-UcosqJ+t
-----!.--=-2
o l-UcosqJ+t
.
u(cos qJ,slnqJ)dqJ =
.
v(cosqJ,slnqJ)dqJ= n(F(t)-F(O.
REFERENCE
J(a,n,k) =
cos(n - k)t ft dt
_.(1-2acost+a 2 )
r..
r'
--f
( n r dz = -ifc (z-anl-azr dz ,
2 c z-a l-az
;
J(a,n,k) =
zR+ p - 1 +Zft-p-I
Zft+p-I
zR+ p - 1
and
c (z-ar(l-azr dz
J.
zR+ p - 1
dz ,
1(
c z-a )ft( l-az
are equal, which can be seen by changing the variable z = lIw in the second integral.
The function/defined by
fez)
zR+ p - 1
ft
(z-a) (l-az)
ft
has two poles: zl = a and z2 = (11 of order n, and inside the unit circle C is either zl or z2'
depending on whether lal < 1 or lal > 1.
By the residue theorem we have
so
Chapter 4
J(a,n,k) =2trResf(z)
(lal< 1),
J(a,n,k) = 2trRe~f(z)
(Ial> 1).
6=41
,=tI
(a +t),,+p-I
(
)".
t" I-a 2 -at
I-a 2
about 1= O. We get
Resf(z)=Resf(z)=
1=0
'=G
aln- kl I(n-I+ln-kl)(n-I+V)
( 1- a2
v=0
n - 1- v
a 2v
(1- a2
n=2k
(1)
n=2k+l
where
m-I
(21' -I)!!
.=1
(2 v)!!
Sm=l+2:(-l)"
F (z) = ------===
(1 +Z2)y I - z 2
51
has two critical singularities, at 1 and -l. Let Fl be the branch of F which corresponds to
the value of the square root which is positive for z = ia (a > 0), i.e. to the value of the
.Jl +a
contour C which encloses the segment [-1, 1] of the real axis. Let F= {zllzl= R}, where
R > 0 is large enough so that the contour C is contained in the region int r. By Cauchy's
residue theorem we have
1 F I (z)dz-1 F I (z)dz=2ni(ResFI (z)+ ResFI(z).
r
c
z=;
Z=-;
+ Res FI (z).
y
P'
M'
N'
Figure 4.2.5.
FI (z) dz=J,
..-.0 MN
2"
If FI(2)dzl~J
0
yrdt=2nyr
2"
c'
and so
z=oo
c'
52
Chapter 4
lim
r.....o
JFI (z)dz=lim
JFI (z)dz=O.
r.....o c.
Z=-I
z=OO
we have
ResFl (z)=
z=;
;2n
y' ,
-2 21
(_ i)2n
z=-;
-2y2 ;
i.e.
is expanded into series around w = 0, the coefficient of wln-2, with changed sign, is the
required residue.
We first determine the sign in the expression ~ w2 -1 =; ~1- w2 which corresponds
to the branch Fl' Namely, for z=-;a we took the value +~I+a2 of the function
z~~I-z2. Since
putting w = -Ilia we should get +~1 +a 2 ; we therefore take the sign +. Therefore
53
(2 n-l)!!
(2 n)!!
(1-w2)-1/2=1+-w2+ ... +
w2n + ...
n~ 3)!! _ (2 n - 5)!!
(2n-2)!!
(2n-4)!!
and SO
V2
(2n-2)!!
(2n-4)!!
implying (1).
REMARK. Fonnula (1) was recorded in Volwne 1 (5.4.3.20, p. 205) without a proof.
REFERENCE
1. F. Tisserand: Recueil complementaire d'exercices sur Ie calcul infinitesimal, deuxieme edition,
augementee de nouveuax exercices sur les variables imaginaires par P. Painleve. Paris 1933.
J
I
(1)
x2n
fX(1-x 2 )
11:
X= y3
14 .. (3n-2)
3.6 ... (3 n)
z 2n
-,;;====
fz(1-z2)
1
FI (z) dz = - 2 n i Res F) (z).
c
z=oo
Let C be the contour shown on Fig. 4.2.6, where the segments PL, MN, N'M', L'P' are
S4
chapter 4
y
M'
N'
c' M
P'
c
As b -+ 0, r -+ 0, we have
-I
.. "
Figure 4.2.6.
IJ
rn
FI (z)dz-+--;
1 fx(l-r)
N'M'
J
-1
Fl(Z)dz-+-62 6
L'P'
I
dx= --J,
JI
x 2n
I
II
Fl (Z) dz
-+ -
~6
PL
rn
xl/l (l-r)l/!
- I
Idx=-(-J)=J,
3
xII! (I - X2)1/!
Idx = - ~6 J,
where
.'
2n
3
2n
3
8=COS-+J S t n - .
I
-I
F(z)=F (-I) =-~-,
W
W2n-I'I_w2
and expanding the function w ~ (I - w2 ll3 into series, the coefficient of wln will be the
required residue (we take that branch of the square root which has the value 1 for w = 0).
We have
(1 -
and so
W 2 )-1/3 =
1 + - w2 + . . . +
3
14 (3 n-2) 2n
w +
3.6.(3n)
...
ss
14 (3 n-2)
. ,
36 (3 n)
ResFJ (z)=a
z=co
where a = 1, &-1, &-2, depending on the chosen branch. But the integral J is a real
number, and so from the equality
2 ( 1 -~) J = - 2 n i a _1._4_._..:..(3_n_-_2....::.)
36(3n)
(1)
J ~x3(1_X2}i
~-dX=
has critical singularities at z = 0, z = 1, z = -1. Outside of the contour C (see Fig. 4.2.6)
the function f is uniform and
(2)
f f(z) dz =
2 n iRes f(z).
Z=OD
(3)
ff(z)dz= J f(x)dx+C
~1I:iJ f(x)dx
eXP-7-1
-I
4Rill(exp'Ri)! If(x)ldx
+ (exp7
S6
Chapter 4
u
+(exp -4ni)(
~ni)( exp -4ni)
JIJ(x)ldx,
exp- 777-1
IfJ denotes the integral on the left hand side of (1), then from (2) and (3) we get
J-~J+J-~J= -2niResJCz),
e
z= 00
i.e.
4n
where e=exp-.
7
W 2n - 1
exp
2n
i)
-7- (1+-w2+
2
7
W 2n - 1
... +
29 (7 n-S) w"+"',
2
)
714 .. (7n)
and hence
ResJ(z)= Z=OO
29 ... (7 n -
5)
i)
n
exp (2
-.
7
2,,;
7
2,,;'
-e --7-
(Iog(l
_zn
iz
along the unit circle, indented at z = 1, applying the standard substitution z = eit and
57
2,.
(1)
~)+ ~
(log(2sin
i(t-n)fdt=O,
(0<1<2 n).
Putting n = 1,2,3 into (1) we obtain
2"
log(2sin
J
2"
~ )dt=O,
t log ( 2 sin
~) dt = 0,
2,.
t log2 ( 2 sin
~) dt =
2"
(2)
log3(2sin
~)dt=3
2"
t 2 log (2 sin
:4 ,
~)dt.
(3)
"
flog" (2 sin
~) dt
were considered by Bowman [2] and Beumer [3], but they did not apply contour
integration. It was found that the integrals of the form (3) can be expressed in terms of the
Riemann zeta-function of integral argument. Indeed, if we substitute the Fourier
expansion
log (2 sin
~) =
cos t
cos 2 t
- ...
"
2 fIog3(2sin
~)dt=
-3n ~ k- 3= -3nC(3).
o
REFERENCES
58
Chapter 4
IIIl
4.2.9. We have
(I)
..J2 (5,r
Jo t 2Jtii sin2tdt = ~
+ 2tr- 2n-Iog2 - 4 + 410g2 - 210g2 2 J.
32
6
.12
Proof Let C be the positively oriented contour consisiting of the upper semicircle
{zllzl=l, Imz~O} and the real line segment [-I, 1]={zl-I~Rez~l, Imz=O}. By
Cauchy's theorem we have
)1/2
I
f(lOgz)2 ( -z
dz=O,
c
I+z
which after the standard parametrization z = eit reduces to
(2)
eil4 J t 2 tg-ei'dt = 2
t<
t)2
I+t
I+t
J(I~
og
dt+2tri Jlogt--dt- ~ J--dt.
JI-t
JI-t
I
The integrals on the right hand side of (2) can be expressed in terms of the beta-type
integral (i.e. in terms of the gamma function).
J_t
P
-
o~
dt =.!.
Si(P-I)
(1- s)
20
and its first and second derivatives in p, evaluated at p = 0 and p = 1. Using the known
values
I ~tgfei'dt
t2
= ~ ( - 51~ -
59
Jo cos"tcosbtdt =
(I)
r(a
I)r(!..=.!.)
sin
2"+lr(I+,,;b)
tr{b-a)
/(z)
z= I,
arg(z+i)
= t,
arg(z-i)
= -t,
arg z = 0
1 x
= arg(z+i) +arg(z-i) = 0;
-i
on its upper part we have arg z = N/2, and on its lower part:
argz = -N/2. Let Cbe the contour shown on Fig. 4.2.10. where
the small circles with the centres at 0, i, -; have radii r, and the
Figure 4.2.10.
larger circle with centre at 0 has radius l. Denote by rI' r2' r3
the arcs of the small circles and by K the arc of the larger circle which comprise the
contour C. Let e(r) be the angle between OA and the imagiruuy axis.
(2)
.(
=r/(z)dz
=-i e-!(b-,,-I)
2
Jr (
I - l )" y. P-,,-I dy
~r
~b-"-I)/
+ i e2
J 1- y2
I-r(
J/(z) dz
1-,,-1 dy
)"
J/(z) dz
r,
-trl2u(r)
+ 2i
Jt!blcos"tdt + J/(z)dz.
r,
tr/2-s(r)
(3)
lim J/(z) dz =0
r-+O
(k
= I,
2, 3).
r.
(4)
-!(b-,,)/) 1 ( ) "
2
=( e~b-,,)/
-e 2
J 1-1 1-,,-1 dy-2"i tr/2Jcos" tcosbtdt.
-tr/2
60
Chapter 4
But
- - = r(a+l)r(~)
ll=.L
JI(l -y 2)" y. .b-,,-ldy=-1JI(1 -u )" utr-d u=-I B ( a+ 1,b-a)
o
20
2
2
2r (a + 1+ 2 )
td2
(S)
The right hand sides of (1) and (S) are fonnally equivalent, but (1) is more correct, since the
expression r(.at-+ 1) which figures in (S) may not be dermed, and we must assume that it denotes
the corresponding limit. For instance, if a = 0, b =4, from (1) we get the correct result: 0, while (S)
becomes
1r
4r(-I)
REMARK 2. Fonnulas (1) and (2) hold provided that b > a > -1. Ifwe let b ~ a, from (1) follows:
"'2
Jcos"tcosatdt =~1
2"+
(6)
The same result is obtained from (S) merely by putting b = a; the limit was already taken when
(1r
?+1
)'
REMARK 3. Fonnula (6) was proved by Cauchy for a > 0 using residues; see Volume 1; S.4.3.11.,
p.197.
REFERENCES
J(cosnx)(P(cosx)Y'2 dx
I<
(1)
1=
(n eN)
-I<
where
Proof If there exists a polynomial Q of degree r such that P = Q2, then (1) becomes
61
J
K
1 = (cosnx)Q(cosx)dx,
-K
~s Q(~(Z+Z-I)}"-I.
We consider the interesting case when there is no polynomial Q of degree r such that
p=Q2.
Since
= J(sinnx)(P(cosx)Y'2 dx =0,
tr
-K
we have
J
K
1 =~fz"-Ip ( !(Z+Z-I) )
Ie
2
112
dz,
p(~(Z+Z-I)
of degree 4r. Now, the main problem is to show that there exists such a positive integer r
that the algebraic equation
r:(K(Z)1I2) = 0
which corresponds to the integral
(R rationalfucntion)
62
Chapter 4
we conclude that g =2r-l. The last equality relates the genus g of the Riemann surface S
to the degree 4r of the polynomial K.
This relation can also be obtained using the Riemann-Roch theorem, which implies
that the branching index of the Riemann surface v = 2{n+g-I), where n is the number of
sheets of the Riemann surface andg its genus. In the case considered we have v = 4rl, all
branch-points being of the first order. There are 4r such points (we made the assumption
that the polynomial K has no multiple zeros): the corresponding Riemann surface has
n = 2 sheets: hence, we have 4r = 2(g+ 1), where from we again obtain g = 2r-1.
We now show that our assumption concerning the zeros of K does not restrict the
problem. Let U be any function on the Riemann surface S taking every complex value
twice. Then there exists a function u on S, satisfying the following relation:
u2 +R,(z)u+~(z) =0,
...
(z-b,,)
Assuming that al
= a2
y=-z-a,
we obtain
y2
and all the zeros are different. Therefore we may start with the function z ~ (K{Z)Jl2
generating the Riemann surface S and having all zeros different, i.e. we shall take
K{z) = (z-e,){z-e2 )
(z-e4,)
63
Ifr> 1, to the equation F;(K(z)t '2 = 0 (which is algebraic relative to the field on the
surface S) corresponds a hyperelliptic field, i.e. the corresponding integral
=~frr-'(K(Z)r2 dz
Ie
is a hyperelliptic integral.
If r
In the case when the integrals (1) are hyperelliptic, they have the following form:
The differentials of the first order have the form
L(z)dz
(K(Z)J12 '
where L is a polynomial whose degree does not exceed g - 1. Their derivatives 2g - 1
linearly independent quadratic differentials of the form
M(z)dz 2
K(z)
where the degree of the polynomial M does not exceed 2g - 2. All other quadratic
differentials (which do not have this form) have the form
N(z)dz 2
(K(Z)J12
where the degree of the polynomial N does not exceed g - 3.
It is interesting to find in which cases the residue at the pole z = 0 of the function
{o
(n ~r+l)
1= 2trB
(nSr)
where B denotes the residue at the pole z = 0 (whose order is r - n + 1) of the function
z 1--+ z"-r-I(K(z)J12.
On the other hand, if we consider the integral I
obtain
I -iH =I
=~fz-(n+r+I)(K(Z)J12 dz,
I
=I -
64
Chapter 4
I =2trB*,
Z H z-<IJ+,.+I)(K(ZJ12
at the pole z = 0 of
If we compare the last result with the result above we see that
B*=O for
n~r+l.
This last conclusion was obtained in an elementaIy way; it is more difficult to obtain the
same conclusion using Abel's differentials of the second order, which is necessary in this
case.
REMARK. This result was only briefly mentioned in Volwne I; 5.4.3.24, p. 206.
REFERENCE
(1)
a
exist. In papers [1], [2], [3] CebySev considered the problem of determining the upper and
J f(x)dx
(a~u; v~b)
REMARK
If the integral
J'
b
I(x) dx
z-x
65
f(X)
1 j f(x)dx+1 j xf(x)dx+ ... ,
j -dx=-Z
-x
Z2
i.e.
b
f(X)dx=Ao+A,+ .. . ,
z-x
Z
Z2
j
a
or
(3)
a
+ b, - a2 z+b;- -
(4)
be the m4h continuous fraction which corresponds to the right hand side of (3). Then
(5)
q~l Resgm(Z);;;;jf(x)dx;;;;
k=l z=zkhm(Z)
Resgm(z).
Assume now that zp and Zq (p < q, Zy < Zq) are two roots of the equation hm(z) = O. If
we put v = zp and then v = Zq into (5), after subtracting we obtain
66
chapter 4
Therefore, Cebyiev determined the bounds of (2) under the condition that u, v are
solutions of the equation hm(z) = O.
In paper [3] Cebyiev obtained a better result. Namely, if the expression on the right
hand side of (3) is expanded as the infinite continuous fraction
1
-------
a.z+b. - - - - a2 z +b2
(6)
L Res F(z) _~
G (z)
Res
F(z) +~ F(v)
G (z)
2 G' (v) ,
II
where the sums of the residues are taken over all the zeros of G which belong to the
interval [a, v). Hence, concludes Cebyev, if we know nothing about the functionjexcept
that it is positive on [a, b) and also the numbers (1), then LReS F(z) gives the value of
G(z)
.!.
F(v) .
2 G'(v)
REMARK 2. In [3] Cebyev examined the value of the expression F(v) and some consequences of
G'(v)
the obtained fonnula (6), but he did not use residues.
REMARK 3. In papers [2], [3] in which he ~pplied residues, Cebyev did not quote the relevant
literature (Cauchy, for instance). He only mentioned in (1] a result of Bienaime [4] which inspired
him to work on this problem.
67
REFERENCES
1. P. Cebyiev: Sur les valeurs limities des integrales. J. Math. Pures Appl. (2) 19(1874), 157160. E Collected Papers 0/ P. L. Gbylev (Russian), t. 3(1948), Moscow - Leningrad, pp. 6365.
2. P. Cebyiev: Sur la representation des valeurs limites des integrales par les residua integraux.
Acta Math. 9(1886), 35-56. == Collected Papers o/P. L. Gbylev (Russian), t. 3(1948), Moscow
- Leningrad, pp. 172-190.
3. P. Cebyiev: Sur les residus integraux qui donnent des valeurs approchees des integrales. Acta
Math. 12(1889),287-322... Collected Papers o/P. L. Gbylev (Russian), t. 3(1948), MoscowLeningrad, pp. 191-225.
4. J. Bienaime: Considerations Ii l'appui de la decouverte de Laplace sur la loi de probabilite
dans la methode des moindres carres. J. Math. Pures Appl. (2) 12(1867).
Julkevi~
and he informed us
Ostrogradski lived in Paris from 1822 till 1827 and attended Cauchy's lectures
delivered at the College de France. Cauchy talked about the results of his research even
before their publication, and naturally he lectured about integrals and residues.
Inspired by Cauchy's lectures, Ostrogradski wrote two articles:
1. Memoire sur la difficulte qui se rencontre dans Ie calcul des integrales dejinies,
lorsque la /onction a integrer est discontinue entre les Iimites de !'integration (paper is
dated: 24 July 1824);
2. Remarques sur les integrales dejinies (the paper is dated: 7 August 1824).
J/(x)dx,
b
where the function/is not continuous on [a, h). In the first article he applied a procedure
similar to Cauchy's principal value, while in the second he explicitly used residues,
although the term residue (i.e. residu) was not mentioned. It should be noted, however,
that Cauchy used word residu for the first time in his paper [1) from 1826, that is to say
after the two papers of Ostrogradski were written.
In fact, Ostrogradski in certain cases used the expression
68
(1)
Chapter 4
d"- 1
- - l i m -(z-a)II/(z),
(n-l)!
Z-+Q
dz"- 1
i.e. the residue off at the pole z = a of order n. Ostrogradski applied his results to the
integral
(2)
J'"
xsinx
cos t-cosx
d x,
~dx
J tgx
J
",/2
and
Ostrogradski submitted his papers to the Academie des Sciences de Paris. In the
opinion of Professor Julkevic, Ostrogradski was not satisfied by the obtained results and
so his manuscripts remained in the archives of the Academy. Julkevic found them there
and analysed them in [I}.
Formula (1) is really due to Cauchy, but Ostrogradski proved it in the second paper.
In the period 1824-1825 Bouniakowski also attended Cauchy's lectures and already in
1825 he successfully applied the calculus of residues in his Doctor's dissertation: Sur Ie
mouvement de rotation dans Ie milieu resistant d'un systeme de plans d'une epaisseur
constante et d'un contour determine, autour d'un axe incline, par rapport a I'horison,
suM de la determination du rayon vecteur dans Ie mouvement elliptique des planetes.
Paris 1825. A translation of this dissertation into Russian is published in Istor.-Mat.
Issled. 28(1985), 247-260.
Bouniakowski skilfully applied the calculus of residues, but he made no original
contribution. In fact, he was the first to apply residues to a specific problem of celestial
mechanics, which was considered among other questions in his dissertation. See [2} and
[3}.
REFERENCES
Issled. 16(1965),7-48.
2. V. S. Kirsanov: V. J. Buniakowski's Dissertation and A. Cauchy's Theory of Residues
(Russian). Istor.-Mat. Issled. 28(1985),261-266.
3. N. S. lennolaeva: On Doctor's Dissertation of V. J. Buniakowski (Russian). Istor.-Mat. Issled.
29(1985),241-255.
Chapter 5
Tk=exp ( -n- .
The sums
are called Gauss' sums. In Volume 1, at the end of 6.4.2 (p.232) the equality
_ 1_;n _ /S --yn
(1)
1 _;
was proved under the condition that n is an odd number. We now determine the formula
Sn' where n is any positive integer.
y
At----_"
S"=2U,,.
where
x
(n even)
(n odd)
(2)
2niz2 )
exp ( --;;-
(3)
Figure 5.1.
zHf(z)=-~-
exp (2niz) - 1
69
chapter,
70
along the contour r shown on Fig. 5.1 where AH = nl2, OA = OD = R and BC, FG are
semicircles of radius r.
Since the function (3) has poles at z = 1, 2, ... , n' and
Tk
2:n;/
Res!(z)=-.
z=k
we conclude that
l'
=Un -~T.
2 0 --~TI
2 n 2
1 2:n; iTo
2 2:n;/
1
2
BC
J!(z)
dz = -
~ Tn/ 2
FG
Hence, the sum Un is equal to the sum of the integrals along line segments HA, AB, CD,
DE,EF,GH.
Along the segment HA we have z = x + iR (x varies from nl2 to 0). The function f
becomes
2:n; i
)
exp ( -n-<x+iR)2
!(z) =---
-.
exp(2:n;i(x+iR)-I -
and so
n/2
4:n;RX)
exp ( - ----
HA
1 - exp ( - 2:n;R)
dx 'F _n__ .
4:n; R
J !(Z)dzl~-.!!-
I,DE
;
4:n;R
Along the y-axis we have z = iy, and so
CD
R-,
( 2:n;; )
exp
idy
exp ( - 2 :n; y) - 1
_--;;-y2
71
.(1
J
,
R
= -
e-2 7r y - 1
2:;
J,
R
exp ( -
I) (2
+ - - - exp - - n; y 2) d y
e2
7r
Y- 1
y2) dy.
Along the segments EF and GH we have z = n/2 + iy, and after some calculations we
find
EF
GH
J(z)dz+
Hence, as r ~ 0, R
J(z)dz=i 3 n+1
~ +00
ex p ( -
2:;
y2)d y .
we get
Sn = 2 Un = 2 (i + j3 n+ I)
+ ..
exp ( _
2:;
y2) dy,
v-;;y =x,
.f2;
or putting
+'"'
13)
v~ n
+ ..
e-ix2 dx,
and so
Je+00
2 (1 + i) =J2in
ix2
dx,
chapter s
72
implying
J
+00
e-1x2 dx =
v;n
(I - i).
Substituting the obtained value into (4) we obtain the Gauss sum
_ (1 + iln) (1 + i) _
(5)
Sn -
r-
n.
REMARK. If n is an odd number, it is easily verified that fonnulas (4) and (5) coincide.
REFERENCE
1. C. Jordan: Cours d'Analyse de I'Ecole Polytechnique, t.2; troisieme ed. Paris 1913, pp.334-339.
(z)
=L--;
+'"
,,=1
(Re z > 1)
+'"
- - = l--z+ L_"'2_lr_Z2k
~-1
2
lr=d2k)!
(lzl<21r),
then,
(2m)
(m=l, 2, ... ).
Proof. Let
qJ(z)=--,
eZ -1
f(z)
V/(z) =
z
e
2,,/.
I'
/(z) = V/(z)
z2m
(m
= 1,
2, ... ).
1
= 2m
Z-2m-IqJ(2mz)
lr=1
(2k)!
The function/has a pole of order 2m+ 1 at z =0, and using (2) we conclude that the
residue is
(3)
.-0
73
(4)
2m e2ldz2111
-0
,.=:1::11
(n
=2m n 2..
=I,
2, ... ) .
(5)
.-=1
z=x{n+~).
en we have
-(n+~J~x~(n+~).
and so,
z
1V'()I
1
11 1 11 e2m + 11 1
e2 1 -1 +-+-=2 2 2 e21d -1 +2
z=(n+fJ+iY, -(n+iJsys(n+~).
and so
Therefore,
IV'(z)1
and hence
(6)
74
ChapterS
sincem ~ 1.
From (5) and (6) follows
o=2f-I-+(-l)-2:z",,r- B:z",
...1
(2m)!'
V-
and finally
(m= I, 2, ... ).
(7)
"I(n+l,.-B"I
_
L
Bo =1,
(n = 1, 2, ... )
or
(8)
(B+lyo+l =B..1
Bo=l,
(n = 1, 2, ... ).
Bernoulli's nwnbers with odd indices, exceptB1, are zero, and from (8) we get
I I I
B. =- 30'
B2 = 42 , ...
B2 ='6'
... 1
,r
L"2=-'
n
6
....1
REMARK 2. Bernoulli's nwnbers are sometimes defined as the coefficients in the expansion
z = 1 --z+...
1
~(I)"'1
B. 211
-~.
e' -1
2"1
(2m)!
In this case we have the so-called Hennite's Connula
(2n)2M
f; k 2M = 2(2m)!B.
+..
L;-n.n)2 =-Je
+'" (
n=O
2Ir
2n
2vCOS1dt.
75
However, since
1
_ew(6+1I2)
+>
"
+>
L." I
z,,=on.
..
L." I
,,=on.
'
we conclude that
1 w(6+1I2)
Res -e
2-0 Z
n! n!
. )2
..=0 n!
'
1)
+>
--2 - - ,
mn
k=1
trz
m
trz
n
/(z) = 2'"tg-tg-cotgtrz
L (n+k)O
k=O
nO-1
2no
~
k-I
(_l)k-I~ r(a+2k-l)
(2k)!
nO+ 2k - 1
where a > 1, 0 < () < 1, and c 1, c2' ... are Bernoulli's numbers:
76
C2
I
=-,
c3 =-,
30
42
REMARK. If 1 < a < 2 and if n is a very large nwnber, Cauchy claimed that the approximation
obtained by suppressing the term which involves () is very good and easy to apply.
REFERENCE
1. A. Cauchy: Sur un emploi legitime des series divergentes. C. R. Acad. Sci. Paris 17(1843),
370-376.
_0
holds for all complex numbers z, except those for which at least one term in the infinite
series is singular.
Proof. Start from the function w H few) defined by
R'
(weC, z eC).
The value tgnw th1Z%W is bounded on any circle e" = {wllwl= n, n ENol if z satisfies the
condition zw ".# (2k+ 1)i12, k e Z when ween' Hence,
jf(w)dw < ~,
c;,
n
lim 1. f(w)dw
n....,.+GO
c;,
=O.
Therefore, by Cauchy's residue theorem it follows that the sum of residues offenclosed by
en decreases towards 0 as n -4 +00.
Res few) =
w=2111+1
2
..J~~!I)I
2z
8
th (2m+ I)m
(2m+I)3
2
'
8z2
th (2m+ I)n
2z'
77
(3)
.=_. . (2m + 1)
3(th(2m+l)nz+z2th(2m+l)1r)_~z=0.
2z
Dividing (3) by z, and letting z -+ 0, we see that the resulting formula is also true for
z = O. It is then easily reduced to (1).
The equality (1) was proved by Grosjean [1], under the condition that z iy
(y irrational real number). The above proof by residues (which is much shorter and also
holds for z = iy, Y irrational) is due to de Doe1der [2].
A special case of (1) corresponding to z = 1, is recorded in [3].
A summation similar to (1) is given in Volume 1, p. 246. Grosjean [1] rediscovered
this formula by his own methods.
REFERENCES
1. C. C. Grosjean: Proofofa Remarkable Identity. Simon Stevin 58(1984), 219-241.
2. P. 1. de Doelder: A Note on "Proof ofa Remarkable Identity". Simon Stevin 61(1987), 193-195.
3. E. R. Hausen: A Table ofSeries and Products. Englewood Cliffs 1975, p.282.
sin:rca sin:rc b
r(c+d-a-b-l)
This formula was proved by 1. Dougall. It can be obtained by integrating the function
" 1
L-sinkO>O.
k=1
Proof. Let
F(U)
1
=L" -sinkO,
k=1
l-z"
/(z)=l-z
(z
=x+iy).
are the points e-i1J, eilJ, respectively. On the circle C we have z = eit, and so
and B
78
Chapter 5
1- e
J/(z)dz =; J--I-'
elldt
l-e
l '"
fJ
AS
-fJ
fJ
.(eifJ_e-1fJ
- 21
2;
1 e2IfJ_e-2IfJ
1 ElllfJ- e-In6J
+ ... +
2
2;
n
2;
= 2; F(O).
We now integrate the function/along the circular arc re {zllz -11 =r} where 0 < r < 1,
with end points A and B. On the arc r we have z = 1 - reiu where -p 5 U S P and p is
determined so that 1 - reip is the point A. Hence
-f 1-(I-re'''
J/(z)dz=-;r J
r
re'''
r.
II
J/(z)dz + J/(z)dz =0
r
AS
implying
i.e.
2; F(O) =;
1(1-(I-re'''r)dU
-p
=;
1(
-p
-p
and therefore
2F(0)
(2)
F(O) =
However, since II-re'''1 < 1 we have Re(I- re''')'' < 1, and the integrand in (2) is
79
Lk
1
1+ -cosk6>
0,
(O<O<n)
t.1
provided that
a o ~al
~ ... ~a. ~O
and
(1 S k
s!!.).
2
REFERENCE
1. L. Vietoris: Ober das Vorzeichen gewisser trigonometrisher Summen. Sitz. Ber. Ost. Ale. Wiss.
167(1958),125-135.
Chapter 6
d )11+1
=( dx
10gI'{x).
+..
They have various applications~ for example, the sum LR(k) can be expressed in terms
k=O
lJ'(.)
r (x),
(2)
For fixed x > 0 this integral converges and defines an analytic function of p only for
Re p > O. Hence, the definition (1) is correct only for Re v<-1.
Grossman [3] obtained an extension of ]P 10gI'{x) to allp. His formula reads:
]Plogr(x)=
(3)
xP
r (p + 1)
(IOg~+y+
x
r'(P+l)
r (p + 1)
r;xl
i'X)
+1
AH ..
- xP
2n;.
A-I ..
80
F(s)C(s) _n_
r(p+l+s) sinns
ds
81
where 1 < A. < 2 and integration is along a vertical line. The function pH]' 10gI'(x)
defined by (3) is an entire function in the p-plane for each x in the plane cut along the
negative real axis (so that log 1 = 0 and xl = e,I0&r). For Re p > 0, the expressions (2) and
(3) coincide.
Using residues Grossman [3J, [4J derived from (3) two interesting expansions when
the contour is deformed in one direction or the other to lie along the real axis.
In the first case, suppose that 0 < x < 1 and deform the line contour into one which
comes from +00 on the real axis, loops around +2 in the positive sense and returns to +00.
There is a sign change because of the change in orientation in the contour integration.
The function
sl4 /(s)
=r
I'(s)(s)
n
I'(p+l+s) sinns
has simple poles at s = 2,3, ... inside the described contour, and
Res/(s)=limr I'(s)(s) (~-k)=(_I)k I1k)(k) ~ .
I-+k I'(p+l+s) SlDns
r(p+l+k)
.-k
Hence,
__I_AT r I1s)(s)
n ds=(-It,(k)
11k)
2m A-/oO I'(p + 1+ s) sin ns
I1p + 1+ k)
implying
]Plogl1x)=
x'
r(p+l)
(I
rxlog-+r+ F'(p+l) - x
I1p+l) p+1
+~(-I)k'(k)B(P+I,k)~ ).
and the summation converges for Ixl < 1 and x = 1.
For x > I, deform the contour in the other direction so that it comes from -00 on the
real axis, loops around + I in the positive sense and returns to -00. This contour encloses
simple poles at s = -2k (k = 1,2, ... ) and double poles at s = 1,0, -I, -3, -5, ... If /has
the same meaning as before, we get
TY2k-p)rt2k+I).
Res/(s)= (-ly+l.I:,
':r'
slDnp
(2n)2k+1 XU
L
r--2k
Res/(s) =
....1
x
(logx- V'(p+2);
I1p+2)
(k
=I, 2, ... );
82
Resj(s) =
1
(log2nx- y- ",(p+l);
21"{p+l)
.-0
I=I-U
(2k)
(IOgnx- ",(p+2-2k) ('(2k)
1"{p+2-2k)
(2k)
(k= 1,2, ... ).
Ixl~
REFERENCES
1 ,1+/01> F(z)
1 b+ltD F(z)
F(s)=-.
-dz+-. -dz
2m ,,_/01> s-z
2m b-ItD s-z
so
-- = JxS-Z-'dx.
1 = xs-z-'d x,
-s-z 0
+ao
z-s
J
I
,1+/00
Jxs-1dx Jr(z)x-zdz.
+ao
b+/OI>
b-/OI>
83
J1\z)x- dz = J1\z)x-zdz,
a+1ot>
)'
b+1>
b-Iot>
a-lot>
1 +'"
a+1>
1\s) =-. xl-1dx 1\z)x-Zdz,
J
2m
()
-n+l
--n
a-I>
or
1\s) =
JH(x,a)r1dx,
+'"
Figure 6.2.
X 1-4 H(x,a)
Jr(z)x-Zdz
a+1>
2m
a-lot>
=-.
.
Res1\z)x- Z = hm (z + v)1\z)x- Z
z=-v
z-+-v
Since lim (z+ v)1\z+ v)
;~-v
.
= hm
=1, we get
XV
Res1\z)x- Z =
z=-v
-v(-v+1) ... (-1)
(_1)VxV
="":""--'--v!
H(a,x)
=n-++oo
lim L
v=O
n
Res1\z)x- Z
z=- v
=n-++CO
lim L
v=0
n
(_1)VxV
vi
=e- x
e- x
Jr(z)x-Zdz
a+1>
2m
a-lot>
=-.
(a> 0),
+>
(4)
r(z) = e-'xz-1dx
(Re z>O),
84
REMARK 2. If z e C, and if the fimctionfis such that the integral
-J
f(x)r-1dx is convergent, it is
Cbaptcr6
F(z) = f(x)r-1dx
o
1 ,..,.
f(x)
=-. JF(z)x"'dz.
2m,....,.
For f(x) =e-', F(z) =r(z) the fonnulas (5) and (6) reduce to (4) and (3).
REFERENCE
1. R. Campbell: Les integrales Euleriennes et leurs applications. Paris 1966, pp. 244-248.
r(l)
(-I)"
= lim (z + n)n:z) =z-+-nz(z+l)
lim
r(z+n+l)
--....:.....:;.--=-... (z+n-I) -n(-n+l) 1
n!
Z-+-IJ
6.3.2. For a positive integer n Kurepa [1] defined the left factorial !n by
(1)
The left factorial function L is defined for complex z, such that Re z > 0, by
(2)
L(z)
= Je-'-dt
+..
(z_1
(-I
11 =L(I) =1.
N. In particular, 10 = L(O) = 0,
We now define the left factorial function L by (2) and (3) for all z
Resn:z)
z=-"
=(-1)"
n!
(n
=0,
1, 2, ... )
C. Since
8S
ResL(z) =-ResF(z)
z=-I
;=0
=-1,
;=-2
(1)2
1 1
ResL(z) = ResL(z) - Resr(z) =--+-,
0=-4
0=-3
0=-3
2! 3!
and by induction we find that
ResL(z) =
o=-n
l)k-1
I---k!
n-I (
(n
k=2
= 3,
4, ... )
The left factorial function is meromorphic with simple poles at z = -1, -3, ~, ... The
points z = and z = -2 are regular (L(O) = 0, L(-2) = 1).
REFERENCES
1. D. Kurepa: On the Left Factorial Function In. Math. Balkanica 1(1971), 147-153.
2. D. Kurepa: Left Factorial Function in Complex Domain. Math. Balkanica 3(1973),297-307.
6.3.3. The function
z H F(z)I'(p-z)
(Rep>O, x>O)
Xo
n!
respectively.
and
(-lr+ 1 I'(p+n)
n!
xp+n
= 0,
(n
=0,
1, 2, ... )
f(z)
= r(z)I'(z-l)
x'
(x> 0)
has a simple pole at z = 1, and second order poles at z = 0, -1, -2, ... We easily get
1
Resf(z)=- .
=1
x
In order to find the residue at z = 0, note that
z2f(z) =z 2r(z)r(z-1)x-;
86
Cbapter6
But
(z-Irl=-I-z- ...
l1z+I)=l1l)+F'(I)z+ ...
Res/(z) = logx-2F'(I)-1.
r-O
Applying a similar procedure, we obtain the residue at any second order pole:
(2
1)
xk
Res/(z) =
10gx--F'(k+I)-r--J:
k!(k+I)!
k!
k+1
(k
=0,
1, 2, ... ).
fez) = l1z)l1z+n)
Xz
(x> 0; n
=0,
I, 2, ... )
has simple poles at z = -k (k = 0, 1, ... n-I) and second order poles at z =-n-k
(k = 0, 1, 2, ... ). For the simple poles we easily get
As in 6.3.4. we expand (z+n+ki fez) into series, and we conclude that the residue
( 1),,+J:
( It
(I + &V/(l + k +n) )---(1 + &V/(I + k)(1 + elog x) x"+J: ,
&(k+n)!
k!
i.e. that
Res /(z)
r--,,-J:
( 1)"~J:
= (k
- +n)!k! (V/(I+k)+V/(I+k+n)-logx),
&-1
in the
87
_I_Cj F(s)r(p-s) ds
(x>O;
x'
2m c-Iot>
Rep>O)
1=-"
Res /(s)
n!
I=p+n
Let L_N be an arbitrary closed contour which encloses a part of the vertical line
Re s = c, where 0 < c < Re p, and which encloses only the "left" poles s =-n;
n = 0, I, ... , N (Fig. 6.4.1.1.). Then, by the residue theorem
2m L..N
x'
n.
,,=0
y
y
-2 -I
Rep
Figure 6.4.1.1.
Figure 6.4.1.2.
We now let N ~ +00, so that the contour L_N deforms into the contour L-aJ shown on
Fig. 6.4.1.2. which is always at a distance &> 0 from the poles, and we get
(1)
=r(p)I(-It p(p+l) ... (p+n-I) x"
..=0
= F(p)
(1+x)P
n!
88
Chapter 6
n!xn
(l+xY'
_1_Cj
O<x<1
x>1
= 0, the loop of the contour L_ should be between the points s = -I and s = 0, and
-1 < c < o. Hence, we let p -+ 0 and we delete the tenn corresponding to n = 0 in the sum (l). This
gives
REMARK. If p
I <+J/oo TT s) TT s)
+..
TTn)
+..
x"
- . ~\ ~\- dS=L(-lt-~\-1 x"=L(-lt-=-log(l+x),
2m <-I..
x'
...1
n.
...1
n
where
Ixl < 1.
I (z)
n
=_1_<+J/oO r(s)r( -n 2.
m c-/oO
s) ds
(n = 0, 1, 2, ... )
(Izl< I).
Similarly, if 0 < c < 1 and n = 0, the integral (1) can be expressed for
the sum of "right" residues at s = k:
Io(z)
1 (1)
+<0
I'(k)
=L(-I)k_,-k
= -log
k=1
1+- .
k. z
If for n = 0 we take -m-I < c < -m (m = 0, I, 2, ... ) then from the expression
-Iog(l+z) we have to take the m-th section of the Taylor expansion, i.e. we have
m
Io(z)
= -log(l+z)+ L(-I)k-I ~
k=1
89
Now let n = 0, 1, 2, ... and -m-l < c < -m (m = -1, 0, 1, ... ,n). The integral (1) can
be considered, on one hand, as the sum of the residues at the left second order poles
s = -m-l, -m-2, ... , -n and the residues at the left simple poles s = -n-l, -n-2, ... , and
on the other hand as the sum of the residues at the right second order poles
s = -m, -m+ 1, ... , and the residues at right simple poles s = 1, 2, ....
If
= F(s)F(:n-s) ,
/(s)
then
Res/(s) =(-lr l V'(I+k)- V'(I+n-k)-logz Zk
r--k
k!(n-k)!
(k
=0,
1, ... , n)
where V'(t)=F'(t)/F(t);
Res/(s)=(_l)kF(k-n)zk
r-k
k!
(n+k)! z
r-k
and hence
I,,(z) = (-1)"
(k
=1,
2, ... ),
k=_1
+(-Z),,+I
kl
(_Z)k,
L
k=O (n +k +1)!
+GO
Izl<l;
and
I (z)
"
k=O
k!(n-k)!
kl
.
+'"
+ (-I)"L
(-zr k - I ,
Izl> 1.
The above expressions for I,,(z) are analytic continuations of each other.
6.4.3. Consider the integral
K (z)
"
(n =0, 1, 2, ... )
The function/defined by
/(s)
=F(s)F(s)F(-n-s)
z'
has simple poles at s = 1, 2, ... , second order poles at s =-n-l, -n-2, ... , and third order
90
chapter 6
(k
(n+k)! Zk
=1,
2, ... )
(1
( 1)lt-k-1
ck(z) = Resf(s) = -log2 z+ ('I'(l+n-k) -2 'I'(l+k)logz+ 2 'I'(l+k)2
6=-k
. k!k!(n-k)! 2
(k
=0,
1, ... ,n).
(1)
Kn(z)
=~>k(Z)+
k=O
+00
WI
k=n+1
k=1
~)k(Z)+ LQk(Z),
3n2
2m
<-/00
(-S),....I
(n
=0,
1, 2, ... ).
The integrand
SH f()
S
( )-0
= r(s)F( -s)
-z
I
(-sr-
The case n = 1 is treated in 6.4.2. Using the same procedure we conclude that
+'"
Ln(z) = LResf(s),
k=lo=-k
Izl< 1
and
+'"
Ln(z) =- LR_esf(s),
k=O o-k
Izl> l.
91
zle
(k eN)
Res/(s)=k"
.=-Ie
z-Ie
Res/(s)=--
(k eN).
(-kr
r-Ie
The residue at s = 0 is more complicated. In order to evaluate it, we write/in the form
/(s)
and we look for the coefficient of S' in the Taylor expansion of the function
s 14 (-1)" F(I +s)F(I- s)( -zr. Multiplying the series for F(l+s), F(I-s) and (-zr',
we get
Res/(s) =
r-O
~l
log"-2V(-z)
",,0
(n-2v)!
~cv--=:-"":'-"':"
where
Izl<I
(2)
and
+'"
-Ie
Izl> 1
L,,(z) = - L-z--Res/(s),
1e=1 (-k)"
r-O
implying
(3)
L,,(z) +
(-I)"L,,(.!.)
= Res/(s)
z
...
(n = 0, 1, 2, ... )
Formula (2) defines the function Ln which is called polylogarithm of order n in the
unit disc. Formula (3) defines the analytic continuation of Ln outside the unit disc.
REMARK. Polylogarithms were briefly mentioned in Volwne 1; 8.2.5, p. 301.
=0,
1, 2, ... )
and that certain definite integrals can be expressed by means of polylogarithms, e.g.
1
Jou-log"-2
1
92
Chapter 6
+JoO.!.e-'I2(~
1
(-I)"
... klk n - I
)dl = L (-z)
(n
"=1'
=1,
2, ... ;
Rez > 0)
logl
1 ( -1) --log
1
2X+-,
,r
--=--dt
= --L2
x>1
log2 t
1 ( 1) 1
,r
dt=-~ - --log x+-Iogx,
Jt(l-x)
x
x
6x
3x
x>l.
+JoO
1
t(l-x)
+00
1
2x
3x
REFERENCE
1. 0.1. Mari~v: A Method a/Evaluating Integrals a/Special Functions (Russian). Minsk 1978.
Chapter 7
REMARK.
However, Sohocki's dissertation is -not only important from the historical point of
view, but also because it contains some serious results which were unjustly overlooked in
mathematical literature. For instance, Sohocki was the first to apply residues to
Legendre's polynomials, but the procedure was wrongly ascribed to Laurent (see, for
example, [I)).
It is very difficult to obtain Sohocki's dissertation. After many years of unsuccessful
attempts and waiting, two copies finally arrived in Belgrade from Moscow. We believe
that this is the first time that Sohocki's results are reviewed in a monograph. and so we
shall expose them with some detail.
The dissertation is divided into two parts: theory and applications. The first,
theoretical, part contains two chapters:
93
94
Cbapter7
VARIANT 2.
=A.
In mathematical literature this theorem is known as the Weierstrass theorem (see, for
example, (2) or (3) where Variant 1 is given, or (4) where Variant 2 is given), and it is
said that it was proved in 1876 by Weierstrass (5). However, this theorem can be found in
Sohocki's dissertation from 1868 and also in the book (6) by the Italian mathematician
Casorati. Therefore, Sohocki and Casorati (clearly independently from each other)
published this theorem eight years before Weierstrass, and the theorem should be called
the Casorati-Sohocki theorem. Until the publication of Markuievi~'s paper (7) the
theorem was in Russian literature also ascribed to Weierstrass (even Markuievi~ did so in
his book (4)) but after the publication (7) it is called Sohocki's theorem or the SohockiWeierstrass theorem. In Western literature the theorem still carries the name of
Weierstrass.
We quote the original formulation of Sohocki from his Master's dissertation:
/fa given function f(z) at some pOint a becomes 00 of infinite order, then the junction
fez) must take all possible values at that point.
Certain explanations are necessary. First of all, the assertion "If the function fez)
becomes 00 at a, then it takes all possible values at that point" appears to be contradictory.
However, it is only a question of inadequate terminology. Namely, when Sohocki
considers a function in a neighbourhood of an isolated singularity, he uses the following
terms. If the Laurent expansion of f around a is
f(z)=
+ ..
2:
k=-"
Ak(z-a)k
95
Sohocki says that / becomes 00 of order n at z = a, and if the Laurent expansion around a
is
+00
/(z~ = LAk(z-a)k
k=-oo
he says that it becomes 00 of infinite order at z = a. In other words, the term "function/ at
z = a becomes 00", does not mean that/(a) =00 but we have two cases:
(i) if "the :function becomes 00 of order n", it means that it has a pole of order n;
(ii) if "the function becomes
singularity.
00
take any value. This can be seen from Sohocki's proof of the theorem, and also from the
example that "the :function sin
It is interesting to note that Sohocki himself did not think that this theorem is
particularly important In the preface he emphasized some other results from his
dissertation, but this theorem is not mentioned. Nevertheless, the theorem should be
known as the Casorati-Sohocki theorem.
REFERENCES
1. Encykloplkiie der mathematishen Wissenschaften, IIA 10, Leipzig 1904-1916, pp.705-706.
2. E. T. Copson: An Introduction to the Theory of Functions of a Complex Variable. 2nd ed.
96
chapter 7
treating rRes as an operator. We quote all those theorems (but with modernized notation)
where we suppose that the functions/, g, h, ... are analytic functions with finite number of
isolated singularities. We shall give proofs of some lesser known theorems.
THEOREM
1. We have
k=1 Z=Zk
THEOREM 2.
n
Resf(z) +
k=1 Z=Zk
Resg (z).
z=1 Z=Zk
If a is a constant, then
Res af(z) = a
k=1 Z=Zk
Resf(z).
k=l Z=Zk
Resf(z, w) = Res
dWz=q
df(z w)
dw
z=q
THEOREM 4.
Z=Zo
a andw
b.
then
Res(Resf(z, w=Res(Resf(z, w).
z=a
THEOREM 5.
w=b
Ifz
w=b
z=a
Res!'(z) = O.
z=a
Resf(z)g'(z) =
z=a
THEOREM 7.
Res!'(z)g (z).
Z=Q
Ifz
f(a) = Res
z=a
f(z) .
z-a
97
Jfz
/(z)
Z=Q (z - a)n+ I
THEOREM 9.
Jfz
1 (d---(z-a)n/(z)
Res/(z)=-(n -I)!
Z=Q
n-
dzn -
Z=Q
THEOREM 10. Jfz = a is a pole or an essential singularity off. then the principal part of
Laurent's expansion around z = a is:
Res /(/) .
I=Q
z-I
Proof If z = a is a pole of order n off, then the principal part of Laurent's expansion
around that point is
i~
k=1 (z-a)k
Since
A_k= Res/(t) (t-a)k-l
t=Q
we have
A_k
--=
(z-a)k
es /(/) (/_a)k- ,
t
t=Q
(z-a)k
and so
~
L.,
A_k
--=
k=1 (z-a)k
es
~ /(/)(/_a)k- 1
L.,
t=Q k=l
(z-a)k
= es
/(/) /-a)"-(z-a)n)
t=Q
(z-a)n(/-z)
implying
~
L.,
A-k
--=
k=1 (z-a)k
R es
/(/)(t-a)n
Z=Q (z-a)n(t_z)
/(/)
+ R es-.
t=Q z-I
However, the first summand at the right hand side of the above equality is equal to 0, and
the theorem is proved.
In a similar manner we conclude that the theorem is true in the case when z = a is an
essential singularity off
98
Chapter 7
THEOREM 1l.1f z
dlogf(z)
Re
s = n;
dz
z=a
z=a
+00
2:
Ak (z - a)k
k=-co
+00
2:
Ak (g (w) - g (b)k,
k=-co
i.e.
J(g(w))g'(w)=
+co
2:
Ak(g(w)_g(f!))kg'(W),
k=-oo
or
(I)
where
g'(w)
g (w) -g (b)
+ ~ G (w),
dw
99
Resf(g(wg'(w)=A_1Res
g'(w)
w-b g(w)-.r(b)
w=b
+ Res.-!. G (w).
w=b dw
Res-G (w)=O,
w=b dw
and so
Resf(g (wg'(w)
w=b
= A_I Res
g'(w)
w=b g(w)-g(b)
g'(w)
es
dlog(g(w)-g(b
dw
w=b
m,
and hence
Resf(g(wg'(w) =mA_ 1
w=b
A_I = Resf(z),
Z=Q
and therefore
Resf(g (wg'(w)
w=b
= m Resf(z),
Z=Q
14. If
Iimzf(z)=O,
~
then
Resf(z) =0.
Z=QO
THEOREM
100
(2)
Chapter 7
J(z)=
~ Res
f(~)
f(w) +Res
.
z-w
w=o w(l-wz)
k=l W=Zk
k=l w=zk
w=o
w-z
(3)
w-z
W=Z
w-z
w=o
) .
w(l-wz)
w-z
At the end of this chapter Sohocki introduced a generalization of the residue, i.e.
residues with respect to a line. Namely, if an analytic functionJis "discontinuous along a
line I", then the residue with respect to I is defined as the integral
._1. .iJ(z)dz,
2:n,J
c
where C is a closed contour which encloses the line /, whileJhas no other singularities in
int C (except those on I). Sohocki stated that the most of the previous IS theorems hold
for this residue, but he did not develop this idea further.
101
,'was proved in Volume 1 (4.4.2, pp. 77-78). In the remark on p. 79 it was stated that
Lagrange proved formula (2) and a more general formula (3) without residues. It seems
that Sohocki gave the first proof of those formulas by means of residues.
The second Lagrange formula has an important role in the theory of symmetric
functions. Sohocki claimed that there was no proof by residues of this formula.
The result in question is as follows. Let u -:t; 0 be an arbitrary constant and let x H j{x)
be a polynomial of degree n divisible by x, but not by x2 nor by x-u. Furthermore, let the
function g be defined by
B,
x
B2
x2
... Xn
Bm
+xm
f(x)-x+u=O.
Then the sum (the symmetric function)
g (XI)
+ ... + g (xn )
+ 00 1 dk-'
L - - (f(u)g'(U))
k=ok! dU k -
in powers of u.
/(z)=%+!!.2 Q 2
q,+q2
102
Chapter 7
P 2 =%Ql+a),
P n+2 =Qn+)Pn+)+an+)Pn
(n= 1, 2, ... )
Q)=I,
Q2=ql>
Qn+2=Qn+)Qn+)+an +)Qn
(n= I, 2, ... )
then
zp
is the circle
f(Z)dW=~~ .
J z-w
q, + q2+
"
q,,+
then
1. P. L. Tchebychef: Sur les fractions continues. J. Math. Pures Appl. (2) 3(1858), 289-323.=
Collected Papers ofP. L. Gbykv (Russian), t. 2, Moscow - Leningrad 1947, pp. 103-126.
2. E. Rouche: Memoire sur Ie developpement des fonctions en series ordonnes suivant les
denominateurs des reduites d'une fraction continue. Journal de lEcole Imperiale Polytechnique,
eahier 37, 21(1858), 1-34.
3. E. Heine: Mittheiling aber Kettenbrache. J. Reine Angew. Math. 67(1867), 315-326.
103
around t
= O. Hence, if
then
+00
In + I v'1 - 2 xt + (2
P k (x)
Ik-fl-l,
k=O
and so
(1)
and this was the main equality used by Sohocki in his proofs.
1 Ifwe apply the formula (see Theorem 6 from 7.2)
dv
du
dt
1=0
dt
we get
t n+ I
(~)
dt t n
1=0
1 d
= Res -
1-0 (n
--c::-----:::
- v' 1 - 2 xl + I 2.
r:.
dt
i.e.
nRes-1-y'1-2xt+/2=Res-1_
1=0
t n+ 1
1=0
t(l-x)
t n+ 1 v'1-2xt+t2 '
104
Chapter 7
which implies
Re s1- n (I +12 - 2 xI) 1=0
12 +Ix
YI-2xI+/ 2
1,,+1
= 0,
or
(2 n - 1) x Res --::-!7='=I~~
1=0/"yl-2XI+/ 2
1
=0,
1=0/" IYI_2xI+/2
+(n-I)Res
i.e.
(2)
(n= 1, 2, ... ).
(n= 1, 2, ... )
using (2) we can evaluate Pl(x), P 2(x), ... successively. Formula (2) is usually called
Bonnet's formula.
P,,(l)=Res--= 1;
1=0/,,+11+1
Res-P" (x)
x=O
Xk+1
i.e.
= Res
1=0
(_1_ Res
--;:--;-;-171~===:)
x=ox k yl-2Ix+t 2 .
1,,+1
+1
lOS
and we get
(3)
1(2kk) Res
0k=-
2k
1=0
1.
t n - k +1 (1 +t 2)
k+2
and from (3) we see that ok = 0 if n - k is odd. This leads to the important conclusion: The
=~
2k
(2kk) Res
1=0
1 (2k)
= 2k k
n=
=~
2k
(1 +t)
n-k
( m=-
'
k=O, 1, ...
n we have m = 0 and so
(2 n -I)!!
n!
~-.
tm+1
1 (2k+l)(2k+3) .. (n+k-l)
(-I)m 2m
m!
2mk!m!
=
k+t 2m +1 (I +t 2) 2
=( _l)m(n+k-I)!!
In particular, for k
Ok.
1
Pn (x) dx=Res-
1=0 t n +1
1
dx
v'1-2xt+t 2
= Res-1- (F(t)
1=0 t n+ 1
-~ vi 1-2 tx+
t
(2),
n).
106
Chapter 7
Let
1
F(/)=-,
t
Then
(4)
and we get
-~Res Z.!(~)=~Res~ dz
PII (x)dx=Res-z-=
1=0
t,1+ I
1-0
dt til
1-0
(Z2-1)"dtdz
1
(Z2-1)"
~-=-Res
-- ,
1 dz
1
Res--=-Res
-1=0
and so
t"dt
didt
2" z=x
(Z2-
1 )" .
1 Res - P" (x) dx=2"n z=.\
z-x
However
1
d" - I
ZZ - 1 )"
Res ( ---= - -..----. (x 2 -1)",
Z=X
z-x
(n-l)!dx,,-I
and hence
implying
(S)
P (x)= __I_~_(x2_1)".
"
2" n! dx n
I - 2 xt + t 2 = T,
we have
z-x
t" dt '
107
1
TI/2 '
P (x) = Res - - n
,=0 TI/2 tn+ I '
u=du
d
J
--P (x)=Res--,
<Ix n
,=0 TJ/2 fn
d2
dx2
df
d2 U
dt 2 -
-Pn(x)=Res---,
,=0 T'/
x-f
TJ/2
--=-,
fn-
1
TJ12'
(X-t)2
T'/2'
._- ----1-3--
,=0 t n dt
dt
,=0 tn+ I
=n(n+ 1) Res
,=0 t n +1 Vl-2xt+t2
i.e.
1 -d ( 12_
dU) =n(n+I)Pn(x).
Res-
(6)
,=0 t n + 1 dt
dt
+ 3 Res
= - Res
,=0 t n +I dt
df
(X-t)2
,=0 t n - I T'/2
+ 2 Res~
,=0 tn TJ/2
t 2 -2tx
dx 2 n
+2x~P (x)-2Res
dx
If
(7)
T
,=0 tn-I T'/2
dt
dx
dx2 n
dx n
108
Chapter 7
Those were the main properties of Legendre's polynomials which Sohocki proved by
the calculus of residues. Sohocki himself emphasized in the preface that the results are
not new. Nevertheless they are important since Sohocki was the first to apply the calculus
of residues to this particular class of special functions. Later on Laurent [1] (see also [2])
also applied residues to Legendre's polynomials. and noted that in this way it is much
simpler to obtain their properties. earlier proved by Heine. without the use of residues. but
there is no mention of Sohocki. Also. Jordan in his book [3] developed the theory of
Legendre's polynomials using residues. without mentioning Sohocki.
REFERENCES
1. H. Laurent: Memoire surlesfonctions de Legendre. J. Math. Pures Appl. (3) 1(1875). 373-398.
2. H. Laurent: Traite d'Analyse. t. 5, Paris 1890, pp. 188-198.
3. C. Jordan: Cours d'Analyse de l'Ecole Polytechnique, t. 2, troisieme ed. 1913 (nouveau tirage
1959), pp. 344-349.
where qo, ql' q2' ... al' a2' ... are polynomials.
Let the functions PI' Ql' P2, Q2' ... be defined by
(n=2. 3... )
and define the sequence (Rn) by
(11= 1, 2... )
The functions R l R 2 .. were studied by Cebyiev [1]. In his Master's dissertation
Sohocki considers the problem of expanding a function/into series with respect to Ql'
Q2' and also with respect to Rl R2 ...
In solving this problem the central role is played by the following formulas. claimed to
be new
:L Res
HI
z- W
(m1=n).
where en is an arbitrary polynomial with degree lower than the degree of qn and residues
109
/(x)
+00
2: Pk (x) Qk (x),
k=l
where PrJ is a polynomial with degree lower than the degree of q", then the coefficients PrJ
are given by
12 ... "
q,,(x)-q,,(z)
x-z
Q,,(z)g(z)/(z)
/ (x) =
+00
L Pk (x) Rk (x),
k=l
REFERENCE
Chapter 8
University of Skopje
l. Suppose that:
10 The function / is analytic in the whole plane where it may have a finite number 0/
singularities~
20 On the real axis the/unction/may have only simple poles cI' ... , cn which belong
to the open interval (a, b).
v.p.
dx
. =
2
n
k=l
where a k
1
(l+x)l
b-c
=__
k
ck -a
(k
l
nl ) Resf(z)+
(cg
Uk+'
=1,
Z=Ck
2 Res ( (a-b)
]ogz
(1 +Z)l
~R
(az+b))
/ --
1 +z
... , n).
We only give an outline of the proof. By the definition of the principal value we have
b
(1)
Ck+l- r
cl-r
V'P.Jf(X)dx=!!!(J f(X)dx+kt J
a
110
f(x)dx+
J
b
f(X)dx).
111
where r > o. Each integral on the right hand side of (n is proper, and we can change the
variable as follows
at+b
l+t
(2)
X=-- --
b-x
x-a
1=--.
=b - e"
e,,-a
f( 1+1
al +b ). The point a is
"
1 +z
1. Let 0 < a < c < b < +00 and consider the principal value
b
V.P.J~dx.
x-c
a
Since
b-c
log akResf(z)=c log-,
z=c
c-a
"" Res {(a-b) logz
az+b
} =b-a-c n i,
..
(1 +Z)2 (a-c)z+b-c
C',..R
we have
J
b
v.p.
x
b-c
-dx=b-a+clog-.
x-c
c-a
Figure 8.1.
112
ChapterS
EXAMPLE 2. If 0 < 0 < c) < C2 < C3 < b < +00, we easily obtain
V.p.
;=1
b-c;
dx-'
- L.,
IT (
3
l oc._Q
g-.
~-~
k=1
k#o;
REFERENCES
1. D. S. Dimitrovski, D. D. Adamovic: Sur quelques formules du colcul des residus. Mat. Vesnik
1(16XI964),113-117.
2. D. S. Dimitrovski, M. Rajovic: Sur Ie valeur principale de I'integrale impropre. IIl-ieme note.
Ann. Fac. Sci. Univ. Skopje 2S-26(1975n6), 41-52.
zl-+f(z) logZ-Q
1. Suppose that:
1 The function /
0/singularities;
/(z)
= LA,.,(Z-XJ)II
(j
=I,
... , k)
n=-CIO
v.p.
Jf(x) dx
a
~A-2P.J=O
J=1
. 'Xk
%rder n l ,
,n", respectively,
113
where
Figure 8.2.
2:
C'-...[a.
hJ
b-z
Res J(z) log finite
distance
J
8
v.p.
x2-x-l
(x-2) (x-3)3
<Ix
z-a
~ f(z)
114
Chapter 8
exists, since
1
I(x)=--+
x-2 (X-3)3
+-(x-3)
and there are no tenns with even indices. Applying the above fonnula we conclude that the
principal value is
v.p. I -log 5 -2 log 3 + 8/45.
ExAMPLE 2. The principal value
J
8
v.p.
(x-2) (X-3)2
dx
does not exist, since there are even tenns in the Laurent expansion.
EXAMPLE 3.
J
b
v.p.
(x_c)2n
dx
exists.
REFERENCES
8.3. The Principal Value in the Case when the Limits of Integration are
Singular Points
If a and b are poles off, we define the principal value of the improper integral as
b
def
b-r
v.P.!f(x)dx=Iim! f(x)dx.
a
r-+O a+r
115
THEOREM
l. Suppose that:
1 The function f is analytic in the complex plane where it can have a finite number of
singularities;
.~. p.
Jf(x)dx
b
"
2 The principal parts of the Laurent expansions around a and b have the same
coefficients, i. e. they have the form:
"
~(z_~)V
and
In that case
.!
- v. P
/(x)dx=
Z=~
z-a)
b-z
II
B
+ 2 {'+
b-a
where 21 = n
if n
3 (b-a)3
+ ... +
}
B 1
2
(21_1)(b_a)21-1
Res /(z)
.~lZ=X.
x.-a
og -.
b-x.
J
b
"
CoROLLARY 2.
If the only singularities offon the real axis are poles at z = a and z = b of
the same order n, and if other conditions of Theorem 1 are satisfied, then the principal
value exists and
116
ChapterS
b
-v.p.!I(X)dX=
a
kanat. dal}
+ 2 { _B2 +
b-a
z -b
B4
3 (b-a)3
B' +
5 (b-a)'
Z= 00
... +
z-b
B
}
21
(2/-1)(a-b)21-1
Res/(z) = Res/(z)
z=a
z=b
and if/has no other singularities inside [a, b), then the principal value exists and
b
finite '
Res
(I (z) log
z-a)
b-z
distance
.!
2
v.P
1
dx
(x-l)(x-2)
Res/(z) = -1 "eRes/(z) = + 1.
z=l
z=2
J
2
V.p.
x-3/2
dx
(x-l)(x-2)
exists, since
Res/(z) = Res fez) = 1/2,
z=l
z=2
!(
b
v.p.
-A+ -B)
- [ f(x)-f(a)-(x-a) feb) - f(a) +c ] dx
exists.
x-a x-b
b-a
b-z
117
j f(x)
v.p.
dx
exists.
J
1
v.p.
4 X 3 -4.x2+1
[eX-(e-l)x-l+c]dx
x 2 (x-l)2
exists.
EXAMPLE 6. We have
4
1
1
v.p. {[- + - - +
x-I
1]
1
1
2 2
3
dx =-]og2+-]og---arctg-.
x-4 (1 +x2)(x-2)
2
10
17 5
5
REFERENCE
%1-'11
v. jf(x)dx=
-00
lim
%3-'31
%1+'12
(jf(x)dx+ j f(x)dx+
'iJc->O Rf<-++oo
i=l, 2, ,m; k=I,2 -R1
... , Xm
f(x)dx+ ... +
%1+'12
xm-'m,l
R2
f(x) dx+
f(x)dx)=
Xm+'m.2
if this limit, which depends on 2m + 2 variables, exists. This limit should be interpreted
as the standard "e-b" definition, i.e. the value of the above integral isA if
118
ChapterS
Instead of the above definition which involves one multiple limit, we can define
various generalized values by considering iterated limits in some prescribed order.
DEFINITION
points c l ,
...
f(x)dx.
CH ri2
(L)
(Q)
'".I,.
or
or
=Z(Ii,k)'
i, j
< n.
oJ.
v. g.
J. oJ or
d~
,,-1 ci+l-ri+hl
cl-rU
a+BI
ci+rtz
b- 82
c,,+r"2
f(x)dx
f(x) dx +
l .)).
119
The above generalized value depends on the path (L, P, Q) along which ril and ri2
tend to zero, and can therefore be multiform. In physics we come across such "integrals
along trajectories" .
If the conditions (L, P, Q) reduce to r il = r,"2 = rj l = rl2 = r and if the
iterated limit from the above definition is interpreted as the standard "e-t5" definition,
then the generalized value reduces to the principal value ofthe improper integral.
DEFINITION 2.
Definition 1 is easily extended to the case when the limits are infinite.
Suppose that the function f is continuous on (--00, +00) with the same
exceptions as in Definition 1. Ifwe take a + &1 = -Rl' b - 62 = R2, then we define the
generalized value as
DEFINITION 3.
+00
J f(x) dx
v.g.
+00
= (L, P, Q)
J f(x) dx
-00
-00
def
'ik- O
'jk ..... O
+:t J
R.--~-oo
1 ci+l-'i+l.l
R 2-->+00
f(x)dx+
f(x)dx]' .. ))).
Cn+'n2
v.p.
J~=IOg Ib-c ,.
x-c
c-a
J
b
v.p.
exists.
-dx
-=
x-c (L)
J
b
f(x)dx+
-RI
R2
Ci+'i2
ai-'ll
dx
b-c
--=Iog
- I+logK
x-c
c-a
120
ChapterS
ExAMPLE 2.
+00
v.
I-x dx=Iim
;xl
......0
R-++ 00
......0
R.....+ 00
1)
lIB
v.p.
I-x
( 210g8+--dx=lim
x2
......o
(Q)
+j"I-X
j+"I-X
--dx=
--dx
x
(Q)
;xl
-lim
......o(R(8)-Cee1/.)
1/.
I-x
--dx--IoIC
x2
+..
v.
......0
R-++..
+ ..
v.p.
J
R
log x
-x- dx = lim
R-++_
logx
-x-dx=O.
I2.>
log x l .
- - dx = - lIm (lor R_Io
x
2 ......0
R=/(s)
for instance, if
(log R)2-(Iog 8)2 - 2 C
121
+ ..
v.
dx
x (1 + x 2 )
- ..
lim {
Bj-+O
R~+ ..
J
R2
-BI
I(x) dx+
-RI
I(x) dx}
B2
log
82
VI +8"2 2
does not exist.
For 8\
J
-.
+ ..
v.p.
dx
---=logI=O
x(1 +X2)
+Joo
v.g.
dx
x (1 +X2)
(Q)
J+"
dx
x (1 + X2)
logk.
ExAMPLES. We have
J
J
2
v.p.
v.p.
(x-I/2) (x-3/2)
dx=O.
x (x-I) (x-2)
2x3+x-I
dx-nI4.
x (x-I) (X2 + 1)
n/2
v.
dtp = _I_lim
a-costp
Vl-a 2 .1-+0
az-+O
n/2
v.p.
{IOgl~I+IOgl VI+a+
vr+a-vr=al}.
Vl-a
82
dtp
1
V1+a-VI-a
---log
a-COs tp V l-a2
VI +a+ VI-a
122
ChapterS
*f
n/2
v.g.
n/2
- -dtp
-=
a-coslp (Q)
dip
a-cos Ip
VI-cz2
VI +a+ VI-a
1 The function f can have a finite number ofsingularities in the region {zlImz > O} ~
2 On the real axis the function f can have poles Xl' ... 'Xm of order nl' ... , nm
respectively~
3 max \Zj{z)1
\zI = r.
v.g. ff(X)dx=
a
* ff(x)dx
(L,P,Q)
rlk-->4J
1=1,2, ..
Rk-++ ..
m; k=1.2
where BUe is the coefficient of{z - x.)-k in the Laurent expansion off around Xv which is
convergent in the annulus
In order to prove this theorem we use the contour shown on Fig. 8.4. which has jumps
along the lines {zl Rez =xk } (k =1, ... ,m) and {zl Rez =O} in order to describe
different approaches of ril and r,7. to zero.
Starting with the Laurent expansions around the poles xv' using the Cauchy residue
theorem, and leaving out the technical details, we arrive at the formula
123
V.g.
+~
-~
__
I~
Rk-+-
'I~
Rk-+_
i-l.2, ... ,m;k-l.2
(J
1m z>O
J dx +
RI
=2xi2Res f(z)+xi
XI
J
+'12
f dx + ... +
J dx) ... )
xm+'m2
2 Resf(z)-
,=1 z=x,
.=1 z=X,
R2
X2-'21
XI-'ll
'I~
~ ~
]J"k
k- 1
'1'2
Figure 8.4.
124
Chapter 8
V.g.
.J+OO
f(x)dx=2niZResf(z)+ni
1m (zo
_00
- Iim...
m
Z
Resf(z)-
._1 z=x.
rvfc-l'O
.-1,2, ... III, k-l,2
Res/(z)
-1
) ...
)
r'l
Since the principal values are the most important, we give a few results in connection
with them. The first one is well known.
THEOREM 3. If the function f is analytic in the region {zl Imz ~ O},
Xl' ... ,
xm on the real axis, and ajinite number ofsingularities in the region {zl Imz > O},
and if
maxlzf(z>l-+ 0, as z -+ 00
then the principal value of the improper integral exists, and
+00
v.p.
-GO
THEOREM 4. The principal value of the improper integral depends only on the even
coeffiCients B v. 2k co"esponding to the terms (z - X .)-2k in the Laurent expansion around
xv
THEOREM 5. Iff has poles of arbitrary order on the real axis, the principal value of the
improper integral exists if and only if
m
PI
=1
.=
.=1
cannot exist.
THEOREM 7. If the function f has only poles of order s; 5 and
125
ExAMPLE 6. We have
+GO
v.g.
eix
:If sin a
- - dx=--a2 _x2
a
lim
ri--"'
i=l. 2. 3.4
-GO
which implies
+GO
v.p.
cosx
:If sin a
--dx=--,
a 2_x2
a
-00
+00
v.p.
sin x
--dx=O.
U 2 _X2
-00
ExAMPLE 7. We have
lIm
ri--"'
i-=l, 2. 3,4
[rl
e r4]
)og-+-)og-,
i
r2
which implies
J
+00
v.p.
cos x
:If
dx=--(e-1+sinl),
x (x-l)(l +X2)
2
-00
+00
v.p.
sin x
:If
- - - - - dx= - (e- 1 + cos 1-2).
2
x (x-I)(I +X2)
-GO
J[
+GO
v.p.
(X_I)2
-GO
(x + I)
+g
(X)] dx=O.
r,
126
Chapter 8
EXAMPLE 9. We have
*J
+00
v.g.
<Ix
(x-a)(x-b) ... (x-n)
-00
Bi~O;
,= 1,
... , n
81 )A (8;:;...
[(8
2 3)B (828;~I)N] +
where
A=
1
(a-b)(a-c) .. . (a-n)
, B=
1
(b-a)(b-c) ... (b-n)
... ,N=
v.p.
<Ix
--------0.
-00
+00
v.
1 I'1m (R 22 - R 12)
xdx= lim x21R2
=RJr++oo 2 -Rl 2R2~+OO
Rl~+OO
-00
we get
*J
+00
v.g.
(Q)
-00
xdx- lim
(R2 -R1)(R2+R1)= lim _1_(2R,
Rio R2~+ 00
RI~+ 00 2 Rl
+~)-1.
Rl
127
+00
dx
n=#:1
(x-c)" '
-00
+00
v.
dx
(X-c)"
-00
lim
81->0
82->0
_1_ [ _ _
1__
I-n (-e 1),,-1
Rhr+ OO
=
(-e 1)"-1
However, if we take
(L):
ez = e1
,,-I
y':;-I-+.....
k-.e-:1":---:"1
J+oo ~=lim
v.g.
(x-c)"
8,....0 (L)
-OD
.J+oo ~=~
-oa
(X-c)"
I-n'
ExAMPLE 12. If -00 < a < -c < c < b < +00, then the value
I I' I I I
J,.
b
dx- - l o
1 g (a+c)(b-c) +1- 1m og e2 e,
-x 2 -c2 2c
(b+c)(a-c)
2c 8i->O
e1 e4
=~ =&:l = &4) is
vp
J~=-I-IOgl (a+c)(b-c) I.
x 2 -c2
,.
2c
(b+c)(a-c)
ExAMPLE 13. Let -00 < A < a < b < C < ... < n < B < +00. The principal value
J
B
v.p.
dx
128
chapterS
where
a-
..... ,.. - - - - - - - - - - - .
(a-b)(a-c) .. (a-n)
(n-a)(n-b)(n-c) .. (n-m)
REFERENCE
1. D. S. Dimitrovski: On the Generalized and Principal Value of an Impraper Integral I
(Macedonian). Ann. Fac. Sci. Univ. Skopje 21(1971),5-12.
10 The function f is analytic in the extended plane, and it can have a finite number of
isolated singularities outside the segment [a, b] ofthe real axis;
2 0 On the real open segment (a, b) the function f can have poles c l , ... 'Cn of order
v.g.
J f(x)dx=vg(L,P,Q;Jf(x)dx
2:
C"-.[a,b]
+lim
rk ..... O
.=
'2.
n.
'"
'"
+--L. L.
{2 . b
11: ,.
. (1
2.
B-k,.
k-l
[e(k-l)in/2_1
k-I
'2v
1)] + __ _l)k-l
+Zl-k _ _ _ _ _
,k-I
,k-I
1.
-k,. + -k,.
l-k
211:i.=lk=2
where
z- a
Z=OO
,k-I
1.
+ __-----:~--k-I
'I.
129
1 b_ 1 v =Resf(z)
..-C"
b-Z)
2 B_1.v=Re f(z)log-z=c.
z-a
3 rl y ' r2 y are the radii of the semicircles around c y which tend to zero depending on
the conditions (L, P, Q);
4 the function
b-z
z 1-+ log-:--:z-a
b-z
z H f(z)log-z-a
along the contour shown on Fig. 8.5.
Figure 8.5.
130
Chapter 8
f f(x)dx= ""
b
v.g.
b-z
b-z
Resf(z)log---Resf(z)log--
L.,
z=oo
Z-Q
~[a,bl
Z-Q
'VI.....o
r V2
.=i
b_ 1 ) .
1
v.p. ff(x)dx=
a
2:
C,[a,bl
Z=oo
Z-Q
-ni
Resf(z).
,=1 z-c"
f(x)dx=
""
L.,
C,[a bl
b-z
b-z
ResfCz)log---Resf(z)log-.
Z-Q
z=oo
Z-Q
4. If all the poles on (a, b) are simple, and ifwe introduce the condition
(L)
we obtain
b
""
L.,
C,[a,bl
b-z
b-z
Z=oo
Z-Q
2:n
5. Ifall the poles on (a, b) are simple, and if fez) = 0(1 / Z2) as z ~ 00, then
A
Z2
Z3
f(z)=-+-+"
and
b-z
Resf(z) log--=O
%=00
Z-Q
Resf(z).
.=1 z=c.
131
'iy
= krZy
we obtain
b
v.g.
b z
C'\..[a,b)
L" ResJ(z).
.=1 Z=C.
V.g. j J(x)dx=
a
C'\..[a,b)
z-a
z=oo
7. If all the radii are equal, i.e. r 1 v = r2 v = r for all v = 1, ... , n, then
b
C'\..[a, b)
z-a
Z=oo
m. (_l)k-l_l
+ !~ k~2 (1 -
k) rk-1
-'lE
i ResJ(z)
.=1 Z=Cv
= 1, 2, ... , n.
v.p. jJ(x)dx=
a
C'\..[a.b)
Z=oo
z-a
-'lEi
ResJ(z).
v=1 Z=Cv
If zZJ(z) ~ M as z ~ 00 and if/is regular on [a, hI, we obtain the following equality
132
chapterS
b
f l(X) dx =
L:
C'\.[a,bJ
Res/(z) log b -z .
z-a
REMARK. The results of this section are taken from the M. A. thesis of StojQ: Mitevski, which has
not been published.
Chapter 9
University of Beograd
J~dx
x2+a2
1=
-1
Ic=
f --dz,
ez
Z2+
a2
r = {zllzl= 1,
x
Figure 9.1.
133
134
Chapter 9
1C= I + 1
l'
z2+a2
z=la
2 7 t.leia
:rr: ia ,
- . =-e
21a
where
Hence,
. a-I.
I = -:rr: eia - 1:rr:
l' = cos a + 1 - :rr: SID
a
l'
Since IE R, the expression ';cos a is its approximate value, and the correction is
contained in Re I y . Notice that the pole ia is sufficiently distant from the semicircle r.
which is the essence of this residue method.
In order to evaluate the approximate value of the integral, we use the expansion
Since
j "e-16 eel6 (I -
a2 r
i26
+00
~
"IF -
(-
k=O
+i
=and so
+00
2:
k=O
I)k a 2k
+ a4 e-i46 -
eeos 6
dO
sin (sin 0 - (2 k + I) 0) dO
(_I)ka 2k j
eeos 6
+00
+00
k=O
k=O
135
(1)
The integrals Ik can be evaluated by the Gauss - Legendre method and we get
J"
e"0s 8 sin
(sin 0 - 0) dO =
J"
e"0l 8 sin
(sin 0 - 3 0) dO =
10 =
o
II =
0.9716595188,
1.582397407.
Therefore the expansion (1) has a form suitable for the numerical evaluation of I:
1
I=g(a)= J~dx
2
(2)
-1
x 2 +a
For instance, for a= 1~, we get 1=31414.9547193153, and all the figures are
correct.
Notice that (2) is the Laurent expansion of the function a H g(a) around a = o. The
function g has a simple pole at that point.
9.2. Let us evaluate an approximate value of the integral
and
J--dx
+OO eix
1 +Xl
136
Chapter 9
1=
f --dz,
eiz
1 +z'
pole
Zo
e1z
eizo
1 + Z3
3 Zo 2
z=zo
Z--"Zo
Figure 9.2.
1= J --dx+
o
J
AB
eiz
l+z3
JO e- X
--idx=2~iBI'
l-ix 3
--dz+
If R -+ +00, then it is easily shown that the integral alongAB tends to 0, and so
+co
II
(1)
cosx
. )--dx=Re
(2nIB
I
l+x'
n)
6
2
-V3/2 cos (1
=-e
--- 3
J --e-xdx
+co
X'
l+x'
J --e
+co
x,
1 +x,
-x dx.
137
J--eIS
x,
1+ x6
dx=0.171811227 ...
wherep <0.
eipz
Figure 9.3.
(2
which is not difficult to evaluate by numerical methods, owing to the presence of the
exponential term, and the calculation becomes simpler as p increases.
For instance, for p = 10 it is enough to replace the upper bound +00 by 4 and the error
is less than 10-13 . Applying the Gauss - Legendre method we obtain 110 l'II
0.01079 184327.
9.4. Let
=J
+00
n,p
cosnx
(x-l)2+ lO-P
dx
(n, p>O).
The integrand is a highly oscillating function for large n, and its pole is near to the
integration intelVal for large p. Integrating along the contour of Fig. 9.2. we get
138
~9
Notice that the integrand on the right hand side of (1) is not oscillatory and that the pole
is at a sufficient distance. Besides, it decreases more rapidly as n increases.
Numerical examples:
110 ,3 IW -60.777375; 110,6
9.5. The integral
..
a =~ JCosnx dx
It
-.
r+pz
IW
cosnz
ZH ---
z2+l
-2635994.6.
(nEN, p>O)
IW
JI
JC+;C
-JC+;C
interval.
e1ltz
Integrating the function Z H - 2 - - 2 along
z +p
the contour shown on Fig. 9.5. as c ~ +00 we
-JC
JC
Figure 9.5.
obtain
The presence of the term e-ny ensures that the numerical calculations can be
performed with high precision. For instance, using the above formula, for p =0.01 we get
alo~90.483336, aloo~36.787794, alooo~4.539952 .10-3
REFERENCES
1. E. L. Kaplan: Numerical Integration Near a Singularity. J. Math. Phys. 31(1952),1-28.
2. D. S. Mitrinovi~ and 1. D. Keai~: The Cauchy Method ofResidues. Theory and Applications.
Dordrecht-Boston-Lancaster 1984.
3. D. D. T~i~: Introduction to Numerical Analysis (Serbian). Beograd 1987.
4. R. Piessens: Gaussian Quadrature Formulas for the Integration of Oscillating Functions.
ZAMM 50(1970),698-700.
139
Chapter 10
University ofNis
:r
:r h (z)
in the presence of a round-off error and uncertain initial data. We assume that zl' ... ,zn
are the poles of the function f, that is the zeros of the function h, which lie in the
corresponding discs ZI' ... ,Zn which are all in int r.
The inclusive calculus of residues, presented in [8], is based on the calculus of
residues of a function/taking a disc in the complex plane as its argument. Therefore, a
short review of some properties of the disc arithmetic (circular arithmetic) will be
described first, including circular complex functions.
A complex circular interval (disc) Z with centre c = mid Z E C and radius r = rad Z
(~ 0), denoted by Z = {c; r}, is the closed set
(1)
={zllz-clsr}.
140
141
The set of all discs of the form (1) will be denoted by K(C).
Addition and subtraction of discs ZI
For the above operation the inclusive isotonicity property is valid ([I; namely
IfA", Bk e K(C) andA k ~ Bk (k = 1, 2), then
AI*~ ~ BI*~
where
*e
Let A (H) be the class of analytic functions in a region H of the complex plane. The
closed set
f*(Z)={f(z)lzeZ, ZeH}
where f e A (H) need not be a disc, which makes considerable difficulties in numerical
calculations. For this reason one often operates with a circular complex function F:D .....
G (D, G ~ K(C which satisfies
(6)
F(Z) 2f*(Z)
F(z) =f(z)
(inclusion)
(complex restriction)
={mp; rp}
142
Chapter 10
1. f(z)=zn;
2. f(z)=ez;
3. f(z) = log z;
log Z
4. f (z) = arctg z;
s.
= {lOg c; -
arctg Z
(IEEZ)
log
(1 - I: I)}
= {arctg c; - ~ log (1 2
(OEEZ);
.' )(1--.I.-cl
'-)} ;
I.-cl
fez) = sin z;
sinZ = {sin c; sh r (ch2 y- sin2 X)1/2 + (ch r -1) (Ch2 y- cos2 X)1/2};
6.
f(z)=cosz;
cos Z = {cos c; sh r (Ch2 y- cos 2 X)I/2 + (ch r -I) (ch2 y - sin 2 x)1/2};
7. f(z)=p(z) =
{n
~
k-O
k=O
k.
The argument of the listed circular functions is the disc Z = {c; r}, where c =x + iy.
The property of inclusive isotonicity of a circular complex function is used in canying
out the inclusive calculus of residues. A circular function F is said to be inclusive isotone
if the implication
is valid. In particular, if Z
Z, then
f(z)=F(z)EF(Z)
is always true.
Let Z H fez) be a regular function on a closed curve r and in int r except at the
poles zl' ... ,zn of orders m l , ... ,mn respectively. Consider the class of functions f which
can be represented in the form f( z) = g (z) I h (z), where the poles of the function fare,
143
at the same time, the zeros of the function h. In that case for each pole zk (k = I, ... , n) we
have
"* o. Then
2: Res/(z),
f/(z)dz=2ni
k=l Z=Zk
(8)
The notation RESf(z) is used to emphasize that RESf(z) is an interval. On the basis of
the inclusive isotonicity property it is clear that
RES j(z)ERESj(z),
Z=Zk
ZkEZk
zEZk
(10)
k=l zEZk
On the right hand side of (10) we have a circular interval of the form {w; R} which
contains the exact value of the considered line integral. The centre w of this disc can be
taken as an approximate value of the integral, whereas the radius R determines the upper
bound of the error.
144
chapter 10
The inclusive calculus of residues will now be illustrated by a few examples. We use
the operations (2) - (5) and the expressions for elementary circular functions.
ExAMPLE 1. We shall fmd the circular interval which contains the exact value of the integral
where
by
Fa = {~lzl=4}. We take that the (exact) poles zl = td and z2 = -td of the function/defmed
/(z)
e'
= (Zl + 11)1
lie in the discs ZI = {3.l4i; 0.01} and Z2 = {-3.14i; 0.01}. respectively. Note that the exact value of
the integral II is
1
i-=;0.318309 ...
n
Since ZI c int Ij and Z2 c int Ij. we calculate the residues for both poles. Applying the
operations of circular aritlunetic and the expression for the circular complex function eZ we obtain
eZ I(ZI+ni-2)
(ZI +n i)3
{6.281593 i; O.o1p
zEZI
RESf(z) =
ZEZ2
(Z2 - n
{ - 6.281593 i;
0.01 P
J: _e_z-
j (Z2 + n2)2
r)
ZE Z 2
Ifwe take the centre 0.3l8645i of the obtained disc for an approximate value of II' then the error is
not greater than 4.17'10-3.
EXAMPLE 2. We shall determine the range of values of the integral
in the case when the variations of the poles of the integrand are not greater then 1% relative to the
exact values %1 = I and %2 =-1.
145
{I; 0.01}
(ZI + 1)3
{2; 0.Ql)2
1
RESf(z) = - - - -
ZE Z 2
(z2-1)2
{0.250025; 0.005032},
{-2; O.OW
{ - 0.250025; 0.OO25065},
so that
12 = 0 E 2 n ; (RES f (z) + RES f (z) = {O; 0.0473646}.
zEZI
ZEZ2
Hence, the range of values of the integral 12 is the disc {O; 0.0473646}.
z~f(z)=-
Z-ZI
lies in the disc Z, = {1.57; 0.001}. Applying the inclusive calculus of residues and circular
arithmetic we shall detennine the value of the interval integral
13= ff(z)dz,
={zJJzJ=2}.
r3
Since Z, c int 13, applying the operations of circular arithmetic and the expression for the
circular fimction sin Z, we obtain
RES f (z) = sin 3Z, = sin { - 4.71; 0.003} = {0.999997; 0.OOOI2},
zEZI
so that
For instance, if z,
=-nil e Z, then
,isin3z
REFERENCES
1. G. Alefeld, J. Herzberger: Introduction to Interval Computation. Academic Press, New York
1983.
2. I. Gargantini, P. Henrici: Circular Arithmetic and the Determination of Polynomial Zeros.
Nwner. Math. 18(1972),305-320.
3. P. Henrici: Uniformly Convergent Algorithms for the Simultaneous Approximation of all Zeros
ofa Polynomial. Proc. Constructive Aspects of the Fundamental Theorem of Algebra. (ed. by
B. Dejon and P. Henrici). Wiley-Interscience, London 1969,77-113.
146
chapter 10
Chapter 11
11.1. Introduction
Suppose that X is a real linear space of functions, with an inner product if, g) :,X2 -+ R
such that
(a)
(Linearity)
(b)
(a/, g) = a if, g)
(Homogeneity)
(c)
if, g) = (g,j)
(Symmetry)
(d)
(positivity)
IfX is a complex linear space of functions, then the inner product if, g) maps,X2 into
C and the requirement (c) is replaced by
(Hermitian Symmetry)
(k =0, 1, ... )
and
(Pk' p",) =0 (k *m),
is called a system of (monic) orthogonal polynomials with respect to the inner product
( , ).
..
The most common type of orthogonality is with respect to the following inner product
(f, g)=
J!(t)g(t)dJ.(t),
R
where dA.{t) is a nonnegative measure on the real line R with compact or infinite support,
147
148
Cbapter 11
I'k=
J tkcU(t),
k=O, 1, ...
exist and are finite, and Jlo > O. Then the (monic) orthogonal polynomials {Pk} satisfy the
fundamental three-term recurrence relation
(1)
k= 0, 1, 2, ... ,
(k=O, 1, ...),
ak= (tPk,Pk)
(PkoPk)
bk =
(k = 1, 2, ... ).
(Pk, Pk)
(Pk-II Pk-t)
is sometimes convenient.
Typical examples of such polynomials are the classical orthogonal polynomials of
Legendre, Cebyev, Gegenbauer, Jacobi, Laguerre and Hermite.
An other type of orthogonality is the orthogonality on the unit circle with respect to
the inner product
2n
(/, g)=
J/(ef")g(ef") dp.(O),
dl'(O)~O.
These polynomials were introduced and studied by SzegO [1]. Monic orthogonal
polynomials { <l>t} on the unit circle satisfy the recurrent relation
cPk + 1 (z) =
which is not of the form (1). For more details consult Nevai [2].
Similarly, we may consider orthogonal polynomials on a rectifiable curve or an arc
lying in the complex plane (e.g. Geronimus [3], SzegO [4]). Complex orthogonal
polynomials may also be constructed by means of double integrals. Namely, introducing
the inner product by
(J, g)=
149
JJJ(z)g(z)w(z)dxdy.
B
for a suitable positive weight function w where B is a bounded region of the complex
plane, a system of orthogonal polynomials can be generated (see Carleman [5] and
Bochner [6]).
(J, g) =
where r= {z
eel z =1
8,
or
(2)
(J, g)=
Notice that this inner product does not satisfy the conditions (c' and (d). Namely, the
second factor in (1), i.e. (2), is not conjugated, so that this product does not possess
Hermitian symmetIy; instead it has property (c).
The corresponding (monic) orthogonal polynomials exist uniquely and satisfy a threeterm recurrence relation of the form (1) from 1l.1, due to the property (zJ, g) = if, zg).
The general case of orthogonality with the complex weight function w with respect to
the inner product
(J, g)=
J!(z)g(z)w(Z)(iZ)-l dz,
i.e.
or
(3)
(J, g)=
ISO
Chapter 11
(4)
[f. g) =
JI(x)g (x)w(x)dx,
-1
which is positively definite and therefore generates a unique set of real (monic)
orthogonal polynomials {Pk}:
[Pi' p",]=O for
k*,m and
k=m
On the other hand, the inner product (3) is not Hermitian; the second factor g is not
conjugated and the integration is not with respect to the measure Iw(e j ')ld6. The
existence of corresponding orthogonal polynomials, therefore, is not guaranteed.
We call a system of complex polynomials {lrk} orthogonal on the semicircle if
[lri
[11'."
11'.. ]
=0
11'.. ]*,0
for k
*' m
for k=m
and
(k, m=O, 1, 2, ... )
,.
(5)
Re(I,I)= Re lw(el9)cl6=1=O.
o
Jg(z)w(z)dz
.
c.
=(r.J+ J
C _l
J-.
J )g(z)w(z)dz= -1+.
J g (x) w (x) dx,
c. +1
151
where rs and C.!;l are the circular parts of Cs (with radii 1 and
assume that w is such that for all g
(2)
lim
&
respectively). We
(VgEP).
g(z)w(z)dz=O
._0 C I
(3)
0=
gEP.
-I
The (monic, real) polynomials {Pk}' orthogonal with respect to the inner product (4)
of 11.2, as well as the associated polynomials of the second kind
J
1
qk(z)=
Pk (Z)-Pk(X)
z-x
-I
w(x) dx
(k=O, 1, 2, ... ),
(4)
where
(5)
Y_I
=0,
Yo
=1 for
{Pk}
and Y_I
=-1,
Yo
=0
for {qk}.
Denote by mk and Pk the moments associated with the inner products (4) and (3) of 11.2,
respectively
bo ="'0.
{z E Cllzl <1, Imz> O} and such that (2) is satisfied and the integrals in (3) exist
(possibly) as improper integrals. Assume in addition that
J
If
(7)
Re(1, 1) = Re w{e(9)dO:t O.
o
Then there exists a unique system of (monic. complex) orthogonal polynomials {Irk}
relative to the inner product (3) of 11.2. Denoting by {Pk} the (monic. rea/) orthogonal
polynomials relative to the inner product (4) ofl1.2, we have
152
(8)
Chapter 11
where
(9)
Alternatively,
(10)
(n=O, 1, 2, ... );
0_1 =1-'0'
where ak' bk are the recursion coefficients in (4) and Po = (1, 1). In particular, all On are
real (in fact, positive) if an = 0 for all n ~ O. Finally,
(n= 1,2, ... ),
Proof. Assume first that the orthogonal polynomials {Nk} exist. Putting
g(z)
=~ Nn(Z)Zk-l,
I
1 s: k < n
in (3) we find
1
(1 :;;.k<n),
nn (z) = Pn (z) -
in (3) and use the first expression for g to evaluate the first integral, and the second to
evaluate the second integral in (3). This gives
(n~
1).
Since (Nn, 1) = 0, (1, 1) = Po, and using (12) with z = 0, we get (9) for n ~ 1. Note that the
denominator in (9) (and the numerator, for that matter) does not vanish, since Re Po 0
153
by (7) and Pk(O), qk(O) cannot vanish simultaneously, {Pk} and {qk} being linearly
independent solutions of (4). For n = 0, (9) yields, by virtue of (5), 0_1 = JIo.
To show the first relation in (10), replace n by n + 1 in (9) and use (4) for z
obtain
= 0, to
,11-1 -
=iall+~
(n~ 1).
(111-1
8o
Po (-ao) + iho
mo
=Iao+-'
Po
'Po
If all an = 0, then w is symmetric and we can prove that (see Theorem 1 in 1l.6)
Po
"n by (8) and (9), it follows readily from (3) that for n ~ 2,
J
1
n,,(z)zk- 1w(z)dz=i
n,,(x)xk- 1w(x)dx=O
_I
and from (13), (9) and (8) for z = n that ("n' 1) = 0 for n ~ l. Furthermore,
(nil' nil) =
1
=[
nn(z)zn-1w(z)dz=i
_I
.. i
J [PII(x)-iO
-I
=011_1
Il _
1PII_I(X)]xn- 1 w(x)dx
-1
'
IS4
Chapter 11
proving (11).
We note from (8) that
(14)
and, similarly,
1 ~k<n,
which, applied repeatedly, gives
(IS)
Here, (11) has been used in the last step. From (14) and (IS) there follows
(16)
the inversion of (8). When k = n, the empty product in (16) is to be interpreted as 1.
EXAMPLE 1. w(z) = 1 + z.
1
(2n + 1X2n + 3)
(n~O),
*' 0,
b = n(n+1)
" (2n+1)2
so that by (10)
n-4i
00 = 3 (2 - in) ,
and by (8),
01 =
3n+8i
S (4+ in)
(n~l),
155
in
==Z2_.-2
EXAMPLE 2. w(z) =
4+in
z-
3 (4+ in)
r.
Here
P-o =
J" e
2i8
dO = 0
so that (7) is violated and thus the polynomials {Irk} do not exist, even though w(x) ~ 0 on [-I, 1]
and the polynomials {Pk} do exist. It is easily seen that q~O) = 0 when k is even, so that 0n-I is
zero for n even, and undefmed for n odd. For an explanation of this example see Theorem 1 in 11.6.
Jwk )d8;t: o.
tr
(1)
Re(l, 1) = Re
= if,
k =0, 1, 2, ... ,
(2)
7'-1
(z) = 0,
1 we get
and substituting here for zPk(z) and zPk_l (z) the expressions obtained from the basic
recurrence relation (4) of 11.3, yields
[ok + i ()k -
()k-l -
k!1;.l,
k !1;.I,
k!1;. 1.
156
Chapter 11
(z) = z- jao =
PI
we find
k~l,
THEOREM
By comparing the coefficients of zk on the left and right of (2), we obtain from (5), (6)
that
(n ~ 1).
157
w(z)=~II(z)=(I-z)a(1 +Z)II,
(1)
where fractional powers are understood in terms of their principal branches. We first
obtain the existence of the corresponding orthogonal polynomials
{t4 Jl)}.
a
THEOREM 1. We have
(2)
Jw<o.
n
dO =
31:
+ i v. p.
o
and hence Re Po
J
1
II) (ei 8)
-1
*" O.
Proof. Let CE , &> 0, be the contour formed by {JD+, with a semicircle of radius r about
the origin spared out. Then, by Cauchy's theorem
J W;X)
&,
respectively). If &~ 0 in
O=I'o-i v. p.
dx-nw(O),
-1
31:(fl.O)(z)
/I) (
"
= ( - 1),,-(0.
31:"
.
where
frIO
- Z ,
denotes the polynomial frn with all coefficients conjugated, i.e. n...(z) = frll(z).
satisfies
(1)
then
If the
158
~ll
(2)
. J
1
w(x)
-;-dx-nw(O),
O=,uo-l v.p.
-1
where the Cauchy principal value integral on the right vanishes because of symmetIy of
w. This proves (2).
Under the assumption (1) we have ak = 0 for all k ~ 0 in (4) of 11.3. In this case the
relations (6), (5) and (4) of 11.4 reduce to
k~
1.
wherem = 1, 2, ...
In particular, for the Gegenbauer weight
A> -1/2,
(3)
=C" (z)
known,
b =m =
o
b -
r (i.+l/2)
lrn
r"' r(i.+l) ,
k(k+2A-l)
k- 4(k+i.)(k+i.-l) '
k~ 1
-
as is well
159
()n=
n(n+2A-l)
1
-,
4(n+A)(n+A-l) 8n - 1
n=l, 2, ""
() _
0-
r(A+l/2)
Vn r(A+l)
i.e.
n~
1.
In particular, we have:
1 The Legendre case A. = 112:
() = (
k
+ 2)/2)
rk+l)/2)
)2
k ~ 0;
PI
ia,
J=
n
p,
0
ia z
Pn-I
ian _ ,
160
chapter 11
The elements of I n are easily computed using (6), (5), (4) of 11.4 and (10) of 11.3,
since the recursion coefficients ale> bk for the orthogonal polynomials {Pk} are known.
If the weight w is symmetric, then Pk = oLI is positive, and I n can be transformed into
a real matrix. Indeed, a similarity transformation with the diagonal matrix
transforms the complex matrix I n into the real nonsymmetric tridiagonal matrix
o
with eigenvalues 1J v = -isv . Using the EISPACK subroutine HQR (see [11 we can
evaluate all the eigenvalues 1J v (v= 1, ... ,n) of An' and then all the zeros sv= i1J v
(v= 1, ... ,n) of 7rn(z).
A theorem on the distribution of zeros of 7rn(z) in the case when the weight function w
is symmetric, i.e.
w(-z)
= w(z)
is proved in [9]. We shall quote here a particular, but important, result regarding the
distribution of zeros of the polynomials n;!(z), orthogonal with respect to the complex
Gegenbauer weight
w(z) = (l-z2
r"
2,
A. > -112.
THEOREM 1.1f A. > -112, all the zeros of ~(z) are simple and for n ~ 2 they belong to
REMARK 1. The problem of distribution of zeros for general weight functions is not solved. The
case of Jacobi weights (see (1) of 11.5) is particularly interesting. Numerical experiments indicate
that all the zeros belong to the half disc D +.
REMARK 2. For the Gegenbauer case a second order linear differential equation is found in [9]
whose particular solution is the polynomial 7r:(z). Applications of those polynomials to numerical
integration and numerical differentiation of analytic functions are given in [10).
REFERENCES
1. G. SzegO: Ober Trigonometrische und Harmonische Polynome. Math. Ann. 79(1918),323-339.
161
2. P. Nevai: Geza Freud, Orthogonal Polynomials and Christoffel Functions. A Case Study. J.
Approx. Theory 48(1986),3-167.
3. Va. L. Geronimus: Polynomials Orthogonal on a Circle and Interval. Pergamon, Oxford 1960.
4. G. SzegO: Orthogonal Polynomials. 4th ed. Amer. Math. Soc. Colloq. Publ. Vol. 23, Amer.
Math. Soc., Providence, R. 1. 1975.
5. T. Carlemann: Ober die Approximation Analytischer Funktionen durch lineare Aggregate von
vorgegebenen Potenzen. Ark. Mat. Astronom. Fys. 17( 1922123), 1-30.
6. S. Bochner: Ober orthogonale Systeme Analytischer Funktionen. Math. Z. 14(1922), 180-207.
7. W. Gautschi, G. V. Milovanovic: Polynomials Orthogonal on the Semicircle. J. Approx. Theory
46(1986),230-250.
8. W. Gautschi, G. V. Milovanovic: Polynomials Orthogonal on the Semicircle. Rend. Sem. Mat.
Univ. Politec. Torino (Special Functions: Theory and Computation), 1985, 179-185.
9. W. Gautschi, H. J. Landau, G. V. Milovanovic: Polynomials Orthogonal on the Semicircle.
Constr. Approx. 3(1987), 389-404.
10. G. V. Milovanovic: Complex Orthogonality on the Semicircle with respect to Gegenbauer
Weight: Theory and Applications. In: Topics in Mathematical Analysis (Th. M. Rassias, ed.),
World Scientific Singapore 1989, pp. 695-722.
11. B. T. Smith et al.: Matrix Eigensystem Routines - EISPACK Guide. Lect. Notes Compo Sci. Vol.
6, Springer-Verlag, New York 1974.
Chapter 12
(z) = {/+(Z),
I-(z).
zELl +.
zELl-.
Observe that K c R implies immediately the equality f+(x) =f- (x) for all x
the same time the functionf(z) has a jump discontinuity j+(x) - f-(x) on K.
R \ K. At
o(J..)
as
Izl
162
={z E C : Im(z) ~ O}
163
J
+00
_1_
(1)
2n i
f+(I)d/={/+(z),ZEA+,
I-z
0 ,zEA-.
-00
To prove (1) observe first that the involved improper integral exists due to the fact that
f+ (I)
o(.!.)
as I~ ~
III
+00.
[-1", r] and the semicircle C: IzI = r. Let z be an arbitrary point inside the contour L.
According to the Cauchy integral formula for a closed contour, we have
_1_J 1+(C) d C=_I_Jf+<C> dC +_1_ ff+(I) dl =/+(z),
C-z
2ni
C-z
2ni,
I-z
2ni
where Re(~ = I. Evidently, for z e ..1- the integral along L equals zero. Letting r ~ +00
the integral along C vanishes and the representation (1) follows at once.
In a similar manner we prove
FACT 2.1fthefunctionf-(z) is analytic in ..1- and continuous in ..1-
J
+00
_1_
2n i
(2)
(z)
1-(1) dt={
t-z
0
,zEA+,
-f-(z), zEA-.
-00
o(.!.)
It I
rp(t) =
as
+00
'P(z) =_1_.
2n,
tp(t)
t-z
dt,
zE~.
-00
Then
. lim [qi(x+ie)-qi(x-ie)]=<p(x)
.....0+
164
Chapter 12
at all points x e R.
PROPOsmON
junctions in
~+
!;;
n. u
~+,
~-.
In particular, if n. = R and if (z) or f- (z) vanishes at infinity and has a pole of order
m at some point z = a, a e ~+ or a e ~-, then by this Proposition and Liouville's theorem
we have
(z) = Pm - t (z) ,
(z_a)m
THEOREM 1.
a", k
=0C~I)
1z1-+ +00
as
values
lim f(x+i8)=f+(x),
......0+
......0+
exist at all points x e R. Then for z E K and z ;It a", k = 1, 2, ... , n the function fez) has
the representation
(3)
f(z) =_1_.
2:1f I
+00
-00
where
/+(1)-/-(1)
I-z
/I
Ok
dt+ ~ ~
Ak,p ,
k=l p=l (z-a)P
165
Proof The hypothesis on the boundaIy values of fez) implies the continuity of the
functions f+(x) and f-(x) on R Put qKx) =f+(x)- rex), x e R First we shall show
that supp( f{J) =K. To do this first consider K to be a proper subset of R, K c R The
function f (z) being analytic on the open set R \ K we have
for all x e R in R \ K. Thus f{J (x) = 0 for all such x; this proves that supP( f{J) ~ K. If
supp(f{J) c K then fez) would be analytic on a set which is larger than R \K contraIy to
the hypothesis. Consequently, supp(f{J) =K. Also, it is easy to see that in the case K = R
we have supp( f{J) = R Therefore the function
F(z)=~
+00
211: I
/+(1)-/-(1)
dt
I-z
-DO
o(.!.)
Izl
as
boundaIy values F(x) for all x e R, respectively, as &-+ 0+. Hence by Proposition 1 we
may write
(4)
Now let us put H(z) =fez) - F(z). Obviously the function H(z) is defined and analytic at
least in C \ R except at the poles off(z). Further, H(z) =
o(.!.)
Izl
According to the relation (4) we have H+(x) =H-(x) for all x e R Hence, by Proposition
2 the function H (z) is analytic everywhere in C except at the poles ak> k = 1, 2, ... , n. By
virtue of the generalized Liouville's theorem, H(z) is a rational function in C which
vanishes at infinity. Consequently the partial fraction expansion of H (z) is possible and
gives the coefficients Ak,p in (3) as can additionally be seen in [1, p. 168]. The
representation (3) is proved.
CoROLLARY
poles
ak (k = 1,
2, ...
,n)
166
Chapter 12
oC!I) as 1:1-+
+00
boundaJy values
limf(x + i e) =f+(x)
e-+O+
limf(x-i e)=](x)
e-+O+
CoROLLARY
J
+-
1
..
fez) =2 _ ,.
--
I+(t)-I-(t)
dt+
t-z
"
Res I(z)
-=------
~ z=ak
L.,
k= 1
Z-Qk
o( I!I) as 1:1-+
+00
values
lim f(x + i e) = f+(x).
e-+O+
e-+O+
f(z) =_1_.
2n,
R. For z E K we have
J
+-
--
f+(t)-I-(t)
t-z
dt.
REFERENCES
Chapter 13
JIf(t)ldt < +00 for every compact K c R (Lebesgue). For example, continuous
K
168
Chapter 13
1
tp(t)={e-l-t2.ltl<l.
Itl~1.
(T, 9'1 +9'2)=(T, 9'1)+(T, 9'2) for all 9'1> 9'2 eD,
(T, c9')=c(T. fP) for all fPeD, ceC.
A linear functional Ton D is continuous if _lim(T,
fP.. ) =0 for every sequence fPn
00
which converges to zero in D as n -+ +00.
A distribution (on R) is every continuous linear functional on D =D (R). The set of all
distributions on R is denoted by D' =D' (R). The vector structure of the set D' is defined
by the formulas
tpED,
cEC.
11-++00
In this definition the existence of the limit TeD' is postulated. However, if (Tn' fP)
169
converges for every rpeD as n~+OO then rpH(T, rp)= lim(T", rp), rpeD, is a
......+'"
member of D'.
For example, the functional 8: D ~ C defined by (8, rp) = rp(O) is the famous Dirac
distribution. The second useful distribution (connected with 8) is the principal value
distribution
rp(t)
(p(!),t rp) PT
_'" t
=
dt = ~~ J rp(t) dt
III:... t
-'"
..
rp e D.
The set of distributions that are most useful are those generated by locally integrable
functions. For every fixed locally integrable function/: R ~ C the functional T D ~ C
defined by
r:
co
(1)
(TI' fP)=
fPED,
-co
is a distribution.
A distribution T is regular if there exists a locally integrable function/such that (1)
holds. All other distributions are termed singular. For example, the distributions 8 and
170
Chapter 13
pGJ) =
R.
The following result is often used: If fJ e D and TeD' are such that
supp (7) f"'I supp (fJ) = 0, then (T, fJ) = O.
= 0,
converges to
vip)
uniformly on every
k(t; z ) = - - 2ni(t-z)
=x + iy being a parameter.
For every
Z E
171
differentiable on R and hence an element of E. If TEE' then its value at k(t, z) is well
defined and depends of z E C \ R. Thus when T operates on E it generates a function of
the complex variable z for all z E C \ R
The function
T(z)=~
2:n:1
2:n:1
e(t
t-z
is called the Cauchy integral of T or the Cauchy representation of T. It can be proved that
fez) is analytic in the domain C \ supp (T).
For example, the Cauchy integral of the Dirac 0 distribution is the function
1 ( 1) = ---.
1
~(z)=-.
A
~t' -
2:n:1
t-z
2:n:IZ
1
( t)- 2:n:;(t-z)
in
GO
J f(t)rp(t)dt,
(Tf> rp)=
rpEE,
-GO
we get at once
1 ( T , ,1)
1 JII --dt
f(t)
T(z)=-.
- =-.
A
2m
t-z
2m
II
t-z
172
Chapterl3
(ii) for eachp = 0, 1,2, ... the sequence (9'': converges uniformly to zero on every
compact of R as n
-+ +00;
(iii) for each p there exists a constant M p ' which is independent of n, such that
- oo<t< 00.
Now we say that rpn -+ rp in t)a if the difference 9'n - 9' -+ 0 in
that the space t)a is complete: the limit function 9' belongs to t)a .
t)a
as n
-+ +00.
Note
The vector space t)a' = t)a' (R) is the space of all continuous linear functionals on t)a'
where continuity has the usual meaning, i. e. that for T E t)a'
(T,rpn)-+O
if
rpn-+O in
t)a
as
n-++oo.
We have D c t)a for every a E R; and if a sequence (rpn) converges to rp in D then the
sequence also converges to rp in t) a. Accordingly every continuous linear functional on
t)a is also a distribution, that is, t)a' cD'.
For example, the distribution
p(7) is an element of
t) a'
Letf(t) be a fixed continuous function on R with f(t) = oC!I) as Itl-+ +00 and let Tf
Jf(t)q;(t)dt,
(T" q;)=
a<O.
-GO
tH----
27d(t -z)
belongs to 0 a for every a ~ -l. Therefore every distribution T E t)a', a ~ -1, acting on
this function generates a function of the complex variable z E C \ R. The function
173
1(T,,-,
1)
(-z
T(z) = - .
A
2m
analytic in the domain C \ supp (T), is also called the Cauchy integral of T
For example, the Cauchy integral of
<
1
1
t1 ) (z)= 2ni
t'
-,
2z
GO
1)
t-z
a ~ -1.
E () a',
1
= 2ni
dl
1 (/-z) =
-GO
-1
-,
2z
1m (zO,
1m (z)<O.
()a
= {z E c: Im(z) > O}
and
= Li+ u
Li-
= {z E c: Im(z) < O}
Li- so that Li
= C \ R.
(1)
e ->- 0+
D(R).
Jf(xie)f{J(x)dx
->-0+ -00
174
Chapterl3
&~
The convergence in l?a' (R) implies the convergence in D' (R), but the converse is not
true. However the following result holds.
PROPOSITION
f(z) =
0(1.)
Izl
as Izl
as&~O+.
0(1.)
Izi
as
Iz I
and f-( x - i&) ~ f- ED' (R) in D' (R) as &~ 0+ such that </+ - f-, rp) = 0 for all
rp E D(R) whose support lies in some open set fl~ R Then there exists a unique
function f(z) that is equal to J+(z) in .1+, and to f-(z) in .1-, and is analytic in .1+ u flu
.1-.
PROPOSITION
3 [4]. Let T(z) be the Cauchy integral of a distribution Tel? a'(R), a ~-l.
(2)
r+ -T- = T,
(3)
+ T-
PROPOSITION
-1
~ a < 0,
=-
~ ( T * VP';).
4 [4]. Let T(z) be the Cauchy integral of a distribution Tel? a' (R),
and T(z) =
0(1.)
Izi
as
Izl
~ +00.
175
closed set
K~R
Izl -+
Let /(Z)=oC!I) as
6--+0+
Ii!: K
+00.
Let
=
!~/(X+i&)=/+
and
1(/,+ -/,-, 1) + ~
n,
/(z)=2---=-
-t-
akA
II
k=l
~ (z:~p)P'
p=l
k
where
1
dm
dzm
z-+ak
r -/-.
Proof Put T =
First we shall prove that supp (T) ~ K. To do this, first
consider K to be a closed proper subset of R. Since the functionf(z) is analytic on the
open set R \ K we have
(/+, fP)= lim (/+ (x+is), fP)= lim (/-(x-is), fP)=(/-, fP)
.....0+
.....0+
<r -/-,
(n
<r -/-)
<r -/-).
F(z)
<r -/-)
t-z
is defined and analytic in the domain C \K. It can be shown that F(z)
=o(~)
Izl
as
Izl -+ +00. In addition by the Plemelj relations (proposition 3) F(x is) converges to a
boundary value F* in D'(R), respectively, as &-+ 0+. Now let us put H(z) =f(z) -F(z).
176
chapter 13
Obviously the function H(z) is defined and analytic at least in C \ R except at the poles of
rf
..!..) as Izl-+ +00 and
"'Uzi
K~R
&--+0+
/(z)
2nl
I-z
HO+
K andz:l: ak' k
ResJ(z)
Z
_=---'-'."k_
k=l
z-ak
K~R
..--+0+
f(z)=~
2nl
Kwe have
REMARK 1. If we replace in Theorem I the convergence in the D'(R) topology with that of
177
1 (T , ,1)
-
T(z)=-.
A
2n.
t-z
.J (T"
We recall that the function f(z) is analytic in C \ supp (T) (hence it is analytic at evety
point of the open set R \ supp (T. Also we remember that
2m r t-z
for t E R
f""\
int (/).
In order to obtain the main result (Theorem 1), we begin by the following lemma.
LEMMA
LEMMA 2.
h (z) =
_1_.
2nl
h (C) dC,
C-z
be the Riemann sum of h(z). Then for each p = 0, I, 2, ... , the sequence of p-th
A
(h~P)(z)
on every compact
178
Chapter 13
2 m. ~
(I'
)P+I
._1 ':o.-Z
II
LV
':0.
each of which is analytic in C \ rfor p = 0, 1,2, ... Now let K be an arbitraIy but fixed
compact subset of C \ r, and let d = d(r, K) be the distance from
K. Denote by I(I)
the length of r and put M
inequality
rto
=maxrlh(~I.
IhC,)(Z)I=I2!... ~
2ni k
II
for p
=0,1,2,
(j,!P)(z)
r.
r.
..
2ni
II
h (C)
(C-Z.)p+l
de = _1_.
2nl
J
r
(j,!P)(z)
converges uniformly to
j,(P)(z)
h(p) (C) de
C-zp
= 0, 1, 2,
2nl
f -de,
h (C)
C-z
1 fh(C)
h(z)=O=-.
--de,
2nl
C-z
r
zEint(l),
zEext(l).
Thus we have
PROPOSITION
J
r
tp (C) dC,
C-z
zEint(I')
179
be the Riemann sum of IP(Z). Then for each p = 0,1, 2, ... , the sequence (qf"P)(z)
converges uniformly to q/p)(z) on every compact subset ofint (l) as n -+ +00 (";0)
=11').
CoROLLARY l. Assume that the conditions of Proposition 1 hold. Put 12= R 11 int (I).
Then the sequence ( IPn (x converges to lP(x) in E(12) as n -+ +00.
Proof EvelY function IPn (x) is the restriction of the function IPn (z) on 12 and hence is
an element of E(12). In addition, evelY compact in R is a compact in C.
an arbitrary function analytic in the strip S ={z eC: -a < Im(z) <a}, 0 < a s: +00.
Further, let r be an arbitrary contour contained in S enclosing the support ofT. Then
(1)
fr
t(z)fP(z)dz=_I_.
2n,
(T,.
Proof Put 12= R 11 int (l) and let fP(x) be the restriction of lP(z) on R. Then the
value of the restriction of T at an element II' e E(12) is equal to (T, 11'). Since the integral
in (1) is a Riemann integral it can be approximated by the Riemann sums:
J=!t(Z)fP(Z)dz= lim
,........ k=l
t(zk)fP(Zk)LJzk
We now exchange the integration and the application of T. Since T is linear we can write
J= lim (Tt _1_
IP(Zk)
2n i k=l t-Zk
,........
LJZk ).
(r
( If
I m-.L,.--LJZk
t.1
n .......
T" - .
2nl
2n I
IP (Zk)
k=l t-Zk
tp(z)
d z ') =(T, -fP(l)) .
t-z
180
Chapter 13
fez) dz=_1_.
21/:1
J _1_)
(T"
t-z
S, then
dz=(T, -I).
= - tp (0).
21/:1
tp(z)
dz= -tp(O).
and
continuous in the closed half plane -:1: = {z E C : Im(z) ~ a}. If F(z) ~ 0 uniformly (with
respect to arg (z), 0 ~ arg (z) ~ 1r) as Izi ~ +00, then the function F(z) has the
representation
oo+;a
F(z)=_1_
21/: i
For z
-oo+;a
F(C)
C-z
dC,
. 1. =..1: u ..1: .
181
and let F(z) in S tend uniformly to zero as \z\--. +00. Then the function F(z) has the
representation
~+~
(1)
Fz=_l_
()
211:;
-~+~
~+~
F(C)dC __1_
C-z
211:i
-~+~
F(C)dC
C-z
~-ia
_F(x)=_l_
(2)
211:;
-~-ia
-1
~+~
F(C)dC __1_
x-C
211:i
-~+Ia
F(C)dC.
x-C
Our aim is to generalize the formula (2) in the sense of distributions in tJa'(R),
a < 0. The proof of the main result (Theorem 1) is essentially based upon
S;
(i) If
A+ia
1
((JA(Z)=-
(3)
211:i
-A+~
'1'(C)
-
C-z
d"'0,
is the Riemann sum of q,A(Z). then the sequence (q,A..,(X) converges to q,A(X) in tJa(R),
-1
S;
(ii) If
(4)
~(z)=~
211:'
.. +Ia
-oo+ia
'1'<C>dC,
C-z
is the Cauchy integral of tp( (). than the sequence ( q,A (x) converges to q,( x) in 0 a(R),
-1 S; a < 0, as A. --. +00 (over a sequence (A.n ) which approaches +00 as n --. +00).
182
Chapter 13
= 0,
ll
1, 2, ...
=Li: u Li:
((p<t,,! ( x) of functions
(5)
converges uniformly to Vf:)(x) on every compact ofR as n ~ +00. On the other hand, for
every p = 0, 1,2, ... , independent of n E N, we derive from (5) an estimate of the form
(x) I~ B". p
I ;,(p)
TA.n
1+lxl'
- oo<x<oo.
a(R) for
On the other hand, the application of these operations to the integral (3) leads to the
formula
~<t)(Z)=_I_.
2nl
(6)
A+la
-A+ia
tp<P)(C)dC
C-z
~+w
2ni
-oo+ia
C- Z
A-..oo
2n i
A+w
fPCP) (C) dC.
-A+w
C- z
183
~ (z) - ~A (z) =
-oo+ia
J'
9' (C)
C-z
dC + _1_.
-A
=_1_
211:i
J
-J
oo+ia
-A+/o
_1_.
211:1
9'(l+ia)
(I+ia)-z
211:1
9' (C)
A+ia
dC
00
dt+-1
211:i
-00
C-z
9'(t+;a)
(t+ia)-z
dt.
Now the relation 19'( f+ia) I s: Aol tl a, a < 0, implies the inequality
JI 1-I+adt+~J I
-A
00
211:
-GO
rIt
GO
tl-l+adt=Ao
:rJ
1- 1 +odt.
Since the last integral can be made arbitrarily small for sufficiently large A, independent
of zELIa' we conclude that ~,t (z) converges uniformly to ~(z) on every compact of Lla
as..t-. +00.
By considering the difference
(p = 1, 2, ... )
it is easy to see, by means of the formula (6), that rlt)(z) converges uniformly to 9'(P)(z)
on every compact of Lla as ..t -. +00. In particular, from this it follows that the sequence
(9'Y) (x) converges uniformly to
= 0,
for every p
1+JxJ
- oo<x< 00.
Consequently, ~,t(x) -. ~(x) in O_l(R) as ..t -. +00. This implies that ~A(X) -. ~(x) in
oa(R) for -1 s: a < 0. The proof is complete.
Now we can formulate the main result.
THEOREM
with T(z) =
o(~)
Izl
as
-b < Im(z) < b satisfying for every p = 0, 1, 2, ... the relation 9'(P)(z) = O(lzl~, -1 s: a
< 0, as Izl -. +00. lfO < a < b, then
184
(7)
Chapter 13
~
2n,
--/11
-_-III
2ni
_+/11
-_+/11
Proof. First of all observe that the operation (T, 'P) is well defined for a given a,
-1 S a<O. Put
-+/11
J1 =
--+/11
-+/11
J2 =
T(z)tp(z)dz,
---/11
First consider the integral J 1. It converges since If(z)'P(Z)1 =O(lzrl+a) on the straight
line (~ia,+OO+ia) as Izl
A+II/
(8)
~ +00.
-A+ia
_00
L T (Zk) tp (Zk) L1 Zk
k=l
n
(9)
l'1m
= 11m
.1...... 00 n-+oo
(T
qJ (Zk) A )
t' -1. ., --LJZk
2:n I k=l t - Zk
A+ia
J l = lim (Tt,
A-+oo
_1_.
2:n I
-.1+ia
oo+ia
qJ(Z)dZ)=(Tt,
t-z
_1_.
2:n I
-oo+ia
qJ(Z)dZ).
t-z
185
=(Tt, 2nl
~
oo-ia
-oo-ia
rp(z)
t-z
dZ).
Consequently,
J=J2 -J1 =
(Tt, 2nl
_1_.
""-ia
dz __2nl1_.
t-z
rp(z)
-",,-ia
",,+ia
-",,+ia
rp(z)
t-z
dz) .
Taking into account Lemma 1 (formula (2 we obtain (with a change of the variables),
that is
J= (T,
-,,(t,
",,-ia
oo+ia
T(z)tp(z)dz-
-oo-ia
T(z)tp(z)dz=(T, -tp(t.
-",,+ia
Name Index
The names are followed by the year of birth or the years of birth and death, if known to
the authors.
176,177,178,179,180,181,183
Cebyev, P. L. (1821-1894) 64, 66, 67, 93,
101,102,108,109,148,159
Coates, R. (1682-1716) 133
Copson, E. T. (1901-1980) 95
Curtiss, J. H. 185
Abel, N. H. (1802-1829) 64
Abramovitz, M. (1915-1958) 82
Adamovic, D. D. (1928) 35,47, 60, 75, 112
A~iev, K. (1949) 114
Ahlfors, L. V. (1907) 38
Albrecht, R. 49
Alefeld, G. (1941) 145
Aramanovich, I. G. (1918) 45,60
Bateman, H. (1882-1946) 77
Behnke, H. (1898-1979) 114
Bernoulli, 1. (1654-1705) 72,74,75
Bessel,F. W. (1784-1846) 9
Beumer, M. G. 57, 58
Bienaime, J. (179~1878) 66,67
Boas,M. 38
Boas, R. P. (1912-1992) 37,38
Bochner, S. (1899-1982) 149, 161
Bonnet, O. P. (1819-1892) 104
Bouniakowski, V. 1. (1804-1889) 67,68
Bo~,F. 57,58
Bremermann, H. 1. 185
B11rmann, H. 22, 23
Elliot, E. B. 39,40
Erdelyi, A (1908-1977) 77
Euler, L. (1707-1783) 21, 58, 83
Evgrafov, M. A (1926) 26,29
~ov,L.(188~1963) 23
Campbell, R. 84
Carleman, T. (1892-1949) 23, 24, 25, 149,
161
Casorati, F. (1835-1890) 94,95
Casteren, 1. van 58
Cauchy, A-L. (1789-1857) 3, 5, 18, 26, 27,
33, 51, 57, 58, 66, 67, 68, 75, 76, 78, 82,
100, 110, 112, 122, 133, 134, 136, 150,
157, 158, 163, 169, 171, 173, 174, 175,
Fejer, L. (1880-1959) 77
Fourier, J. (1768-1830) 57,138,185
Fresnel, A (1788-1827) 37,38
Friedmwn,E. 37,38
Frullani, G. (1795-1834) 39
Gakhov, F. D. (1~1980) 166
Gargantini, I. 141,145
Gauss, K. F. (1777-1855) 69,133, 135,136,
137
187
188
Name Index
149.161
B. (1864-1947) 18. 19
Gegenbauer. L. (1849-1903) 148. 157. 158.
159.160
Gel'fond. A. O. (1906-1968) 34
Geronimus. Ya. L. (1898) 148. 161
Gerretsen. J. (1911) 41
Gilles. J. 58
Gregory. M. B. 75
Grosjean. C. C. 77
GroSSlWlll. N. 80.81.82
Gupta. R. C. 9
Gau~.VV.(1927)
Gavrilovi~.
n'kov. L. V. 34
Jackson. D. (1888-1946) 77
Jacobi. K. G. J. (1804-1851) 14. 148. 156.
160
Jensen. J. L. VV. V. (1859-1925) 18.25
Jermolaeva. N. S. 68
Jordan. C. (1838-1922) 72. 108
Julia. G. (1893-1978) 14. 32
JuJkevi~. A. P. (1906) 67.68
Kaplan. E. L. 138
Keati~. J. D. (1945) 1.17.95.138.146
Kirsanov. V. S. 68
Kurepa. D. (1907) 56. 84, 85
101
Laguecre,E.(1834-1886) 34,148
Landau, H. J. 149. 161
Laurent, H. (1841-1908) 93. 108
Laurent, P. A. (1813-1854) 3, 4. 94. 95. 97.
98, 101, 108, 112, 114, 115, 122, 124,
135
Lebesgue, H. (1875-1941) 167
Legendre. A. M (1752-1833) 93. 94. 103.
107,108,133,135,136,137,148,159
Liouville, J. (1809-1882) 80,164,165,176
Littlewood, R. K. 139
Lunts, G. (1911) 45,60
Lutzky, M. 57
Maclaurin, C. (1698-1746) 21
Mari~v, O. I. 92
MarkuJevi~.A.1. (1908-1979) 94.95
Marsden. J. E. 6
Marsh, D. C. B. (1926) 50
Mellin. Hj. (1854-1933) 84
Metzger, J. M 75
Milovanovi~, G. V. (1948) 147. 149, 161
Mitevski, S. 132
Mitrinovi~, D. S. (1908) 1, 50.64, 138. 146
Mitrovi~, D. (1922) 25, 56, 162, 167, 185
Mittag-LeIDer, G. (1846-1927) 26
Moore, R. E. 146
Muskhelishvili, N.I. (1891-1976) 166
Nevai, P. 147, 161
Newman, D. J. (1930) 21,22
Newton. I. (1643-1727) 133
Norfolk, T. S. 10
Sansone, G. (1888-1979) 41
Schempp, VV. 22
Schoenberg, I. J. (1903Schwartz, L. (1915) 185
) 21,22
189
Name Index
Smith. B. T. 161
Sohocki, J. V. (1842-1927) 64, 67, 93, 94,
95,100,101,102,103,107,108
Sommer, F. (1912) 114
Stegun, I. A. (1919) 82
Stewart, I. 38
Stieltjes, T. J. (1856-1894) 177
Stirling, 1. (1696-1770) 21
Sveshnikov, A. (1911) 185
Szego, G. (1895-1985) 34, 148, 160, 161
Tall,D. 38
Taylor,B. (1685-1731) 88,91,103,142
Turin. P. (1910-1976) 79
Vietoris, L. (1891) 79
Vitali, G. (1875-1932) 177
Wang, E. T. H 37
Watson, G. N. (1886-1962) 14
Weber, H R. 48
Wei~,lC.(1815-1897) 18,94,95,177
Whittaker, E. T. (1873-1956) 14
Wong,R. 139
Young, R. L. 37
Zakian, V. 139
Zuser, lC. 49
Zwier,P. 75
Subject Index
division of discs 141
Additionofdiscs 141
lIII8lytic continuation 88, 90
- of Cauchy type integrals 26
- of Mittag-leIDer's function 27
- principle 164
asymptotic formula for Mittag-leIDer's
function 29
Elliptic integral 59
essential singularity, order of 14
Euler-Maclaurin's formula 21
Euler's constant 57
Euler's formula 82
evaluation of residues
- at a pole 4, 5, 8, 9
- at a regular point 4
- at a removable singularity 4
- at infinity 4
F~~-]~sinequality
76
190
191
Subject IDdex
Mittag-Leffier's function 26
- analytic continuation 27
- asymptotic formula 29
multiplication of discs 141
Newton-Coates' formulas 133
numerical evaluation of definite
integrals 133
Order of essential singularity 18
orthogonal polynomials 147
orthogonality 147
- on a curve 148
- on the semicircle 149
- on the unit circle 148
392 pp.
ISBN 90-277-1376-6
D.S. Mitrinovic and JD. Keckic: The Cauchy Method 0/ Residues. Theory and
ISBN 90-277-1623-4
R.A. Askey, T.H. Koomwinder and W. Schempp (eds.): Special Functions: Group
Theoretical Aspects and Applications. 1984,352 pp.
ISBN 90-277-1822-9
R. Belbnan and G. Adomian: Partial Differential Equations. New Methods/or their
Treatment and Solution. 1984,308 pp.
ISBN 90-277-1681-1
S. Rolewicz: Metric linear Spaces. 1985,472 pp.
ISBN 9O-2TI-1480-O
ISBN 90-277-2300-1
ISBN 90-277-1917-9
N.V. Krylov: Nonlinear Elliptic and Parabolic Equations of the Second -Order.
1987, 480.pp.
ISBN 90-277-2289-7
W.I. Fushchich and A.G. N'ddtin: Symmetries of Maxwell's Equations. 1987.228
pp.
ISBN 90-277-2320-6
P.S. Bullen, D.S. Mitrinovic and P.M. Vasic (eds.): Means and their Inequalities.
1987,480 pp.
ISBN 90-277-2629-9
V.A. Marchenko: Nonlinear Equations and Operator Algebras. 1987, 176 pp.
ISBN 9O-277-2654-X
YuL. Rodin: The RiemLlnn Boundary Problem on RiemLlnn SUrfaces. 1988,216 pp.
ISBN 90-277-2653-1
A. Cuyt (ed.): Nonlinear Numerical Methods and Rational Approximation. 1988,
480 pp.
ISBN 90-277-2669-8
D. Przeworska-Rolewicz: Algebraic Analysis. 1988,640 pp. ISBN 90-277-2443-1
V.S. Vladimirov. YU.N. Drozzinov and B.I. Zavialov: Tauberian Theorems for
Generalized Functions. 1988.312 pp.
ISBN 90-277-2383-4
G. Morosanu: Nonlinear Evolution Equations and Applications. 1988. 352 pp.
ISBN 90-277-2486-5
A.W. Leung: Systems of Nonlinear PDEs: Applications to Biology and EngineerISBN 0-7923-0138-2
ing. 1989.424 pp.
N.E. Hurt: Phase Retrieval and Zero Crossings: Mathematical Methods in Image
ISBN 0-7923-0210-9
Reconstruction. 1989.320 pp.
V.L Fabrikant: Applications of Potential Theory in Mechanics. A Selection of New
Results.1989.484pp.
ISBN 0-7923-0173-0
R. Feistel and W. Ebeling: Evolution of Complex Systems. Sel!organization.
Entropy at:ld Development. 1989.248 pp.
ISBN 90-277-2666-3
S.M. Ermakov. V.V. Nekrutkin and A.S. Sipin: Random Processes for Classical
Equations of Mathematical PhySics. 1989.304 pp.
ISBN 0-7923-0036-X
B.A. Plamenevskii: Algebras of Pseudodifferential Operators. 1989. 304 pp.
ISBN 0-7923-0231-1
N. Bakhvalov and G. Panasenko: Homogenisation: Averaging Processes in
Periodic Media. Mathematical Problems in the Mechanics of Composite MaterialS.
1989.404 pp.
ISBN 0-7923-0049-1
ISBN 0-7923-0442-X
M.I. Rabinovich and 0.1. Trubetskov: Oscillations and Waves in Linear and
ISBN 0-7923-0445-4
D.S. Mittinovic. I.E. Pecaric and A.M. Fink: Classical and New Inequalities in
Analysis. 1992. 740 pp.
ISBN 0-7923-2064-6
H.M. Hapaev: Averaging in Stability Theory. 1992.280 pp. ISBN 0-7923-1581-2
S. Gindinkin and L.R. Volevich: The Method
ISBN 0-7923-2037-9
Yu.A. Mitri>polsky. AM Samoilenko and DJ. Martinyuk: Systems of Evolution
Equations with Periodic and Quasiperiodic Coefficients. 1992. 280 pp.
ISBN 0-7923-2054-9
LT. Kigmadze and T.A. Chantmia: Asymptotic Properties of Solutions of Nonautonomous Ordinary Differential Equations. 1992.332 pp. ISBN 0-7923-2059-X