Nonlinear Functional Analysis: Gerald Teschl
Nonlinear Functional Analysis: Gerald Teschl
Nonlinear Functional Analysis: Gerald Teschl
Gerald Teschl
Gerald Teschl
Fakultät für Mathematik
Nordbergstraße 15
Universität Wien
1090 Wien, Austria
E-mail address: Gerald.Teschl@univie.ac.at
URL: http://www.mat.univie.ac.at/~gerald/
The present manuscript was written for my course Nonlinear Functional Analysis
held at the University of Vienna in Summer 1998 and 2001. It is supposed to give
a brief introduction to the field of Nonlinear Functional Analysis with emphasis
on applications and examples. The material covered is highly selective and many
important and interesting topics are not covered.
It is available from
http://www.mat.univie.ac.at/~gerald/ftp/book-nlfa/
Acknowledgments
I’d like to thank Volker Enß for making his lecture notes available to me and
Matthias Hammerl for pointing out errors in previous versions.
Gerald Teschl
Vienna, Austria
February 2001
iii
iv Preface
Contents
Preface iii
v
vi Contents
5 Monotone operators 53
5.1 Monotone operators . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
5.2 The nonlinear Lax–Milgram theorem . . . . . . . . . . . . . . . . . 55
5.3 The main theorem of monotone operators . . . . . . . . . . . . . . 57
Bibliography 61
Glossary of notations 63
Index 65
Chapter 1
where o, O are the Landau symbols. The linear map dF (x) is called derivative of
F at x. If F is differentiable for all x ∈ U we call F differentiable. In this case we
get a map
dF : U → L(X, Y )
. (1.2)
x 7→ dF (x)
1
Qm we call F continuously differentiable and write F ∈ C (U, Y ).
If dF is continuous,
Let Y = j=1 Yj and let F : X → Y be given by F = (F1 , . . . , Fm ) with
Fj : X → Yi . Then F ∈ C 1 (X, Y ) if and only if Fj ∈QC 1 (X, Yj ), 1 ≤ j ≤ m, and
in this case dF = (dF1 , . . . , dFm ). Similarly, if X = m
i=1 Xi , then one can define
the partial derivative ∂i F ∈ L(Xi , Y ), which is the derivative of F considered as
a function
P of the i-th variable alone (the other variables being fixed). We have
dF v = ni=1 ∂i F vi , v = (v1 , . . . , vn ) ∈ X, and F ∈ C 1 (X, Y ) if and only if all
partial derivatives exist and are continuous.
1
2 Chapter 1. Analysis in Banach spaces
The set of all r times continuously differentiable functions for which this norm is
finite forms a Banach space which is denoted by Cbr (U, Y ).
If F is bijective and F , F −1 are both of class C r , r ≥ 1, then F is called a
diffeomorphism of class C r .
Note that if F ∈ L(X, Y ), then dF (x) = F (independent of x) and dr F (x) = 0,
r > 1.
For the composition of mappings we note the following result (which is easy to
prove).
Lemma 1.1 (Chain rule) Let F ∈ C r (X, Y ) and G ∈ C r (Y, Z), r ≥ 1. Then
G ◦ F ∈ C r (X, Z) and
then
sup |dF (x)| ≤ M. (1.7)
x∈U
1.1. Differentiation and integration in Banach spaces 3
since we can assume |o(ε)| < εδ for ε > 0 small enough, a contradiction. 2
As an immediate consequence we obtain
|F | = sup
Qm
|F (x1 , . . . , xm )| < ∞. (1.11)
x: i=1 |xi |=1
Rt (1.18)
This even shows that F (t) = F (a) + a
(Ḟ (s))ds for any F ∈ C 1 (I, X).
shows that x is a fixed point and the estimate (1.20) follows after taking the limit
n → ∞ in (1.22). 2
Next, we want to investigate how fixed points of contractions vary with respect
to a parameter. Let U ⊆ X, V ⊆ Y be open and consider F : U × V → U . The
mapping F is called a uniform contraction if there is a θ ∈ [0, 1) such that
we infer
1
|x(y + v) − x(y)| ≤ |F (x(y), y + v) − F (x(y), y)| (1.26)
1−θ
and hence x(y) ∈ C(V, U ). Now let r = 1 and let us formally differentiate x(y) =
F (x(y), y) with respect to y,
|∂x F (x(y), y)| ≤ θ by Theorem 1.2. Hence we get a unique continuous solution
x0 (y). It remains to show
x(y + v) − x(y) − x0 (y)v = o(v). (1.28)
Let us abbreviate u = x(y + v) − x(y), then using (1.27) and the fixed point
property of x(y) we see
(1 − ∂x F (x(y), y))(u − x0 (y)v) =
= F (x(y) + u, y + v) − F (x(y), y) − ∂x F (x(y), y)u − ∂y F (x(y), y)v
= o(u) + o(v) (1.29)
since F ∈ C 1 (U ×V, U ) by assumption. Moreover, |(1−∂x F (x(y), y))−1 | ≤ (1−θ)−1
and u = O(v) (by (1.26)) implying u − x0 (y)v = o(v) as desired.
Finally, suppose that the result holds for some r − 1 ≥ 1. Thus, if F is
C r , then x(y) is at least C r−1 and the fact that d x(y) satisfies (1.27) implies
x(y) ∈ C r (V, U ). 2
As an important consequence we obtain the implicit function theorem.
Theorem 1.6 (Implicit function) Let X, Y , and Z be Banach spaces and let
U , V be open subsets of X, Y , respectively. Let F ∈ C r (U × V, Z), r ≥ 1, and fix
(x0 , y0 ) ∈ U × V . Suppose ∂x F (x0 , y0 ) ∈ L(X, Z) is an isomorphism. Then there
exists an open neighborhood U1 × V1 ⊆ U × V of (x0 , y0 ) such that for each y ∈ V1
there exists a unique point (ξ(y), y) ∈ U1 × V1 satisfying F (ξ(y), y) = F (x0 , y0 ).
Moreover, the map ξ is in C r (V1 , Z) and fulfills
dξ(y) = −(∂x F (ξ(y), y))−1 ◦ ∂y F (ξ(y), y). (1.30)
Proof. Using the shift F → F − F (x0 , y0 ) we can assume F (x0 , y0 ) = 0.
Next, the fixed points of G(x, y) = x − (∂x F (x0 , y0 ))−1 F (x, y) are the solutions
of F (x, y) = 0. The function G has the same smoothness properties as F and
since |∂x G(x0 , y0 )| = 0, we can find balls U1 and V1 around x0 and y0 such that
|∂x G(x, y)| ≤ θ < 1. Thus G(., y) is a uniform contraction and in particular,
G(U1 , y) ⊂ U1 , that is, G : U1 × V1 → U1 . The rest follows from the uniform
contraction principle. Formula (1.30) follows from differentiating F (ξ(y), y) = 0
using the chain rule. 2
Note that our proof is constructive, since it shows that the solution ξ(y) can
be obtained by iterating x − (∂x F (x0 , y0 ))−1 F (x, y).
Moreover, as a corollary of the implicit function theorem we also obtain the
inverse function theorem.
8 Chapter 1. Analysis in Banach spaces
Proof. Fix x0 ∈ Cb (I, U ) and ε > 0. For each t ∈ I we have a δ(t) > 0 such
that |f (x) − f (x0 (t))| ≤ ε/2 for all x ∈ U with |x − x0 (t)| ≤ 2δ(t). The balls
Bδ(t) (x0 (t)), t ∈ I, cover the set {x0 (t)}t∈I and since I is compact, there is a finite
subcover Bδ(tj ) (x0 (tj )), 1 ≤ j ≤ n. Let |x − x0 | ≤ δ = min1≤j≤n δ(tj ). Then
for each t ∈ I there is ti such that |x0 (t) − x0 (tj )| ≤ δ(tj ) and hence |f (x(t)) −
f (x0 (t))| ≤ |f (x(t)) − f (x0 (tj ))| + |f (x0 (tj )) − f (x0 (t))| ≤ ε since |x(t) − x0 (tj )| ≤
|x(t) − x0 (t)| + |x0 (t) − x0 (tj )| ≤ 2δ(tj ). This settles the case r = 0.
Next let us turn to r = 1. We claim that df∗ is given by (df∗ (x0 )x)(t) =
df (x0 (t))x(t). Hence we need to show that for each ε > 0 we can find a δ > 0 such
that
sup |f∗ (x0 (t) + x(t)) − f∗ (x0 (t)) − df (x0 (t))x(t)| ≤ εδ (1.32)
t∈I
|f∗ (x0 (t) + x(t)) − f∗ (x0 (t)) − df (x0 (t))x(t)| ≤ εδ(t) (1.33)
whenever |x(t) − x0 (t)| ≤ δ(t). Now argue as before. It remains to show that df∗
is continuous. To see this we use the linear map
such that we know the solution for ε = 0. The implicit function theorem will show
that solutions still exist as long as ε remains small. At first sight this doesn’t seem
to be good enough for us since our original problem corresponds to ε = 1. But
since ε corresponds to a scaling t → εt, the solution for one ε > 0 suffices. Now
let us turn to the details.
Our problem (1.37) is equivalent to looking for zeros of the function
Lemma 1.8 ensures that this function is C r . Now fix λ0 , thenR G(0, λ0 , 0) = 0
t
and ∂x G(0, λ0 , 0) = T , where T x = ẋ. Since (T −1 x)(t) = 0 x(s)ds we can
apply the implicit function theorem to conclude that there is a unique solution
x(λ, ε) ∈ C r (Λ1 × (−ε0 , ε0 ), Dr+1 ). In particular, the map (λ, t) 7→ x(λ, ε)(t/ε) is
in C r (Λ1 , C r+1 ((−ε, ε), X)) ,→ C r (Λ × (−ε, ε), X). Hence it is the desired solution
of our original problem. 2
Chapter 2
2.1 Introduction
Many applications lead to the problem of finding all zeros of a mapping f : U ⊆
X → X, where X is some (real) Banach space. That is, we are interested in the
solutions of
f (x) = 0, x ∈ U. (2.1)
In most cases it turns out that this is too much to ask for, since determining the
zeros analytically is in general impossible.
Hence one has to ask some weaker questions and hope to find answers for them.
One such question would be ”Are there any solutions, respectively, how many are
there?”. Luckily, this questions allows some progress.
To see how, lets consider the case f ∈ H(C), where H(C) denotes the set of
holomorphic functions on a domain U ⊂ C. Recall the concept of the winding
number from complex analysis. The winding number of a path γ : [0, 1] → C
around a point z0 ∈ C is defined by
Z
1 dz
n(γ, z0 ) = ∈ Z. (2.2)
2πi γ z − z0
It gives the number of times γ encircles z0 taking orientation into account. That
is, encirclings in opposite directions are counted with opposite signs.
In particular, if we pick f ∈ H(C) one computes (assuming 0 6∈ f (γ))
Z 0
1 f (z) X
n(f (γ), 0) = dz = n(γ, zk )αk , (2.3)
2πi γ f (z) k
11
12 Chapter 2. The Brouwer mapping degree
2. The number of zeros can change if two zeros with opposite sign change (i.e.,
opposite signs of f 0 ) run into each other.
3. The number of zeros can change if a zero drops over the boundary.
Hence we see that we cannot expect too much from our degree. In addition, since
2.2. Definition of the mapping degree and the determinant formula 13
it is unclear how it should be defined, we will first require some basic properties a
degree should have and then we will look for functions satisfying these properties.
Its complement CV(f ) = Rn \RV(f ) is called the set of critical values. We set
C r (U , Rn ) = {f ∈ C r (U, Rn )|dj f ∈ C(U , Rn ), 0 ≤ j ≤ r} and
(D3). If U1,2 are open, disjoint subsets of U such that y 6∈ f (U \(U1 ∪ U2 )), then
deg(f, U, y) = deg(f, U1 , y) + deg(f, U2 , y) (additivity).
Before we draw some first conclusions form this definition, let us discuss the
properties (D1)–(D4) first. (D1) is natural since deg(f, U, y) should have something
to do with the solutions of f (x) = y, x ∈ U , which is the same as the solutions
14 Chapter 2. The Brouwer mapping degree
Theorem 2.1 Suppose deg satisfies (D1)–(D4) and let f, g ∈ Dy (U , Rn ), then the
following statements hold.
(ii). If y 6∈ f (U ), then deg(f, U, y) = 0 (but not the other way round). Equiva-
lently, if deg(f, U, y) 6= 0, then y ∈ f (U ).
(iii). If |f (x) − g(x)| < dist(y, f (∂U )), x ∈ ∂U , then deg(f, U, y) = deg(g, U, y).
In particular, this is true if f (x) = g(x) for x ∈ ∂U .
Proof. For the first part of (i) use (D3) with U1 = U and U2 = ∅. For the
second part use U2 = ∅ in (D3) if i = 1 and the rest follows from induction. For
(ii) use i = 1 and U1 = ∅ in (ii). For (iii) note that H(t, x) = (1 − t)f (x) + t g(x)
satisfies |H(t, x) − y| ≥ dist(y, f (∂U )) − |f (x) − g(x)| for x on the boundary. 2
Next we show that (D.4) implies several at first sight much stronger looking
facts.
Theorem 2.2 We have that deg(., U, y) and deg(f, U, .) are both continuous. In
fact, we even have
(i). deg(., U, y) is constant on each component of Dy (U , Rn ).
(ii). deg(f, U, .) is constant on each component of Rn \f (∂U ).
Moreover, if H : [0, 1] × U → Rn and y : [0, 1] → Rn are both continuous such
that H(t) ∈ Dy(t) (U, Rn ), t ∈ [0, 1], then deg(H(0), U, y(0)) = deg(H(1), U, y(1)).
N
X
deg(f, U, 0) = deg(f, U (xi ), 0). (2.7)
i=1
It suffices to consider one of the zeros, say x1 . Moreover, we can even assume
x1 = 0 and U (x1 ) = Bδ (0). Next we replace f by its linear approximation around
0. By the definition of the derivative we have
Lemma 2.3 Two nonsingular matrices M1,2 ∈ GL(n) are homotopic in GL(n) if
and only if sgn det M1 = sgn det M2 .
16 Chapter 2. The Brouwer mapping degree
shows that diag(±1, 1) and diag(∓1, −1) are homotopic. Now we apply this result
to all two by two subblocks as follows. For each i starting from n and going down
to 2 transform the subblock diag(mi−1 , mi ) into diag(1, 1) respectively diag(−1, 1).
The result is the desired form for M .
To conclude the proof note that a continuous deformation within GL(n) cannot
change the sign of the determinant since otherwise the determinant would have to
vanish somewhere in between (i.e., we would leave GL(n)). 2
Using this lemma we can now show the main result of this section.
Lemma 2.5 (Sard) Suppose f ∈ C 1 (U, Rn ), then the Lebesgue measure of CV(f )
is zero.
Proof. Since the claim is easy for linear mappings our strategy is as follows.
We divide U into sufficiently small subsets. Then we replace f by its linear ap-
proximation in each subset and estimate the error.
Let CP(f ) = {x ∈ U |Jf (x) = 0} be the set of critical points of f . We first
pass to cubes which are easier to divide. Let {Qi }i∈N be a countable cover for
U consisting of open cubes such that Qi ⊂ U . Then it suffices S to prove that
f (CP(f ) ∩ Qi ) has zero measure since CV(f ) = f (CP(f )) = i f (CP(f ) ∩ Qi )
(the Qi ’s are a cover).
Let Q be any of these cubes and denote by ρ the length of its edges. Fix ε > 0
and divide Q into N n cubes Qi of length ρ/N . Since f 0 (x) is uniformly continuous
on Q we can find an N (independent of i) such that
Z 1
0 ερ
|f (x) − f (x̃) − f (x̃)(x − x̃)| ≤ |f 0 (x̃ + t(x − x̃)) − f 0 (x̃)||x̃ − x|dt ≤ (2.12)
0 N
18 Chapter 2. The Brouwer mapping degree
for x̃, x ∈ Qi . Now pick a Qi which contains a critical point x̃i ∈ CP(f ). Without
restriction we assume x̃i = 0, f (x̃i ) = 0 and set M = f 0 (x̃i ). By det M = 0 there
is an orthonormal basis {bi }1≤i≤n of Rn such that bn is orthogonal to the image of
M . In addition, there is a constant C1 such that Qi ⊆ { n−1 i ρ
P
i=1 i | |λi | ≤ C1 N }
λ b
(e.g., C1 = n2(n/2) ) and hence there is a second constant (again independent of i)
such that
n−1
X ρ
M Qi ⊆ { λi bi | |λi | ≤ C2 } (2.13)
i=1
N
(e.g., C2 = nC1 maxx∈Q |f 0 (x)|). Next, by our estimate (2.12) we even have
n
X ρ ρ
f (Qi ) ⊆ { λi bi | |λi | ≤ (C2 + ε) , |λn | ≤ ε } (2.14)
i=1
N N
and hence the measure of f (Qi ) is smaller than CN3nε . Since there are at most N n
such Qi ’s, we see that the measure of f (Q) is smaller than C3 ε. 2
Having this result out of the way we can come to step one and two from above.
Step 1: Admitting critical values
Z N Z
X
φε (f (x) − y)Jf (x)dx = φε (f (x) − y)Jf (x)dx
U i=1 U (xi )
N
X Z
= sgn(Jf (x)) φε (x̃)dx̃ = deg(f, U, y), (2.16)
i=1 Bε0 (0)
where we have used the change of variables x̃ = f (x) in the second step. 2
Our new integral representation makes sense even for critical values. But since
ε depends on y, continuity with respect to y is not clear. This will be shown next
at the expense of requiring f ∈ C 2 rather than f ∈ C 1 .
The key idea is to rewrite deg(f, U, y 2 ) − deg(f, U, y 1 ) as an integral over a
divergence (here we will need f ∈ C 2 ) supported in U and then apply Stokes
theorem. For this purpose the following result will be used.
where Df (u)j is the determinant of the matrix obtained from f 0 by replacing the
j-th column by u(f ).
Proof. We compute
n
X n
X
divDf (u) = ∂xj Df (u)j = Df (u)j,k , (2.18)
j=1 j,k=1
where Df (u)j,k is the determinant of the matrix obtained from the matrix associ-
ated with Df (u)j by applying ∂xj to the k-th column. Since ∂xj ∂xk f = ∂xk ∂xj f we
infer Df (u)j,k = −Df (u)k,j , j 6= k, by exchanging the k-th and the j-th column.
Hence n
X
divDf (u) = Df (u)i,i . (2.19)
i=1
20 Chapter 2. The Brouwer mapping degree
(i,j) (i,j)
Now let Jf (x) denote the (i, j) minor of f 0 (x) and recall ni=1 Jf ∂xi fk = δj,k Jf .
P
Using this to expand the determinant Df (u)i,i along the i-th column shows
n
X n
X n
X
(i,j) (i,j)
divDf (u) = Jf ∂xi uj (f ) = Jf (∂xk uj )(f )∂xi fk
i,j=1 i,j=1 k=1
X n n
X Xn
(i,j)
= (∂xk uj )(f ) Jf ∂xj fk = (∂xj uj )(f )Jf (2.20)
j,k=1 i=1 j=1
as required. 2
Now we can prove
Proof. Fix ỹ ∈ Rn \f (∂U ) and consider the largest ball Bρ (ỹ), ρ = dist(ỹ, f (∂U ))
around ỹ contained in Rn \f (∂U ). Pick y i ∈ Bρ (ỹ) ∩ RV(f ) and consider
Z
deg(f, U, y ) − deg(f, U, y ) = (φε (f (x) − y 2 ) − φε (f (x) − y 1 ))Jf (x)dx (2.21)
2 1
U
where
Z 1
u(y) = z φ(y + t z)dt, φ(y) = φε (y − y 1 ), z = y 2 − y 1 , (2.23)
0
R
and apply the previous lemma to rewrite the integral as U divDf (u)dx. Since the
integrand vanishes in a neighborhood of ∂U it is no restriction to assume
R that ∂U is
smooth
R such that we can apply Stokes theorem. Hence we have U divDf (u)dx =
∂U
D f (u)dF = 0 since u is supported inside Bρ (ỹ) provided ε is small enough
(e.g., ε < ρ − max{|y i − ỹ|}i=1,2 ). 2
2.3. Extension of the determinant formula 21
Remark 2.9 Let me remark a different approach due to Kronecker. For U with
sufficiently smooth boundary we have
Z Z
1 1 1 ˜ = f , (2.25)
deg(f, U, 0) = n−1 Df˜(x)dF = n Df (x)dF, f
|S | ∂U |S | ∂U |f |n |f |
Since f˜ : ∂U → S n−1 the integrand can also be written as the pull back f˜∗ dS of
the canonical surface element dS on S n−1 .
This coincides with the boundary value approach for complex functions (note
that holomorphic functions are orientation preserving).
Our final step is to remove the condition f ∈ C 2 . As before we want the degree
to be constant in each ball contained in Dy (U , Rn ). For example, fix f ∈ Dy (U , Rn )
and set ρ = dist(y, f (∂U )) > 0. Choose f i ∈ C 2 (U , Rn ) such that |f i − f | < ρ,
implying f i ∈ Dy (U , Rn ). Then H(t, x) = (1 − t)f 1 (x) + tf 2 (x) ∈ Dy (U , Rn ) ∩
C 2 (U, Rn ), t ∈ [0, 1], and |H(t) − f | < ρ. If we can show that deg(H(t), U, y) is
locally constant with respect to t, then it is continuous with respect to t and hence
constant (since [0, 1] is connected). Consequently we can define
Proof. If f −1 (y) = ∅ the same is true for f + t g if |t| < dist(y, f (U ))/|g|.
Hence we can exclude this case. For the remaining case we use our usual strategy
of considering y ∈ RV(f ) first and then approximating general y by regular ones.
Suppose y ∈ RV(f ) and let f −1 (y) = {xi }N j=1 . By the implicit function theorem
i
we can find disjoint neighborhoods U (x ) such that there exists a unique solution
xi (t) ∈ U (xi ) of (f + t g)(x) = y for |t| < ε1 . By reducing U (xi ) if necessary, we
can even assume that the sign of Jf +t g is constant on U (xi ). Finally, let ε2 =
dist(y, f (U \ N i
S
i=1 U (x )))/|g|. Then |f + t g| > 0 for |t| < ε2 and ε = min(ε1 , ε2 ) is
the quantity we are looking for.
It remains to consider the case y ∈ CV(f ). pick a regular value ỹ ∈ Bρ/3 (y),
where ρ = dist(y, f (∂U )), implying deg(f, U, y) = deg(f, U, ỹ). Then we can
find an ε̃ > 0 such that deg(f, U, ỹ) = deg(f + t g, U, ỹ) for |t| < ε̃. Setting
ε = min(ε̃, ρ/(3|g|)) we infer ỹ − (f + t g)(x) ≥ ρ/3 for x ∈ ∂U , that is |ỹ − y| <
dist(ỹ, (f + t g)(∂U )), and thus deg(f + t g, U, ỹ) = deg(f + t g, U, y). Putting it
all together implies deg(f, U, y) = deg(f + t g, U, y) for |t| < ε as required. 2
Now we can finally prove our main theorem.
Proof. Our previous considerations show that deg is well-defined and locally
constant with respect to the first argument by construction. Hence deg(., U, y) :
Dy (U , Rn ) → Z is continuous and thus necessarily constant on components since
Z is discrete. (D2) is clear and (D1) is satisfied since it holds for f˜ by construction.
Similarly, taking U1,2 as in (D3) we can require |f − f˜| < dist(y, f (U \(U1 ∪ U2 )).
Then (D3) is satisfied since it also holds for f˜ by construction. Finally, (D4) is a
consequence of continuity. 2
To conclude this section, let us give a few simple examples illustrating the use
of the Brouwer degree.
First, let’s investigate the zeros of
Theorem 2.15 Let X and Y be Banach spaces and let K be a closed subset of X.
Then F ∈ C(K, Y ) has a continuous extension F ∈ C(X, Y ) such that F (X) ⊆
conv(F (K)).
Proof. Consider the open cover {Bρ(x) (x)}x∈X\K for X\K, where ρ(x) =
dist(x, X\K)/2. Choose a (locally finite) partition of unity {φλ }λ∈Λ subordinate
to this cover and set
X
F (x) = φλ (x)F (xλ ) for x ∈ X\K, (2.34)
λ∈Λ
2.5. Kakutani’s fixed-point theorem and applications to game theory 25
fixed-point theorem to set valued functions. This generalization will be more suit-
able for our purpose.
Denote by CS(K) the set of all nonempty convex subsets of K.
K = hv1 , . . . , vm i, m ≤ n, (2.38)
unless x1 is one of the vertices, this doesn’t help us too much. So lets choose a
better function as follows. Consider the k-th barycentric subdivision and for each
vertex vi in this subdivision pick an element yi ∈ f (vi ). Now define f k (vi ) = yi
and extend f k to the interior of each subsimplex as before. Hence f k ∈ C(K, K)
and there is a fixed point
m
X m
X
k
x = λki vik = λki yik , yik = f k (vik ), (2.40)
i=1 i=1
of the i-th player. We will consider the case where the game is repeated a large
number of times and where in each step the players choose their action according to
a fixed strategy. Here a strategy si for the i-th player is a probability
Pmi k distribution
1 mi k
on Φi , that is, si = (si , . . . , si ) such that si ≥ 0 and k=1 si = 1. The set
of all possible strategies for the i-th player is denoted by Si . The number ski is
the probability forQthe k-th pure strategy to be chosen. Consequently, if s =
(s1 , . . . , sn ) ∈ S = ni=1 Si is a collection of strategies, then the probability that a
given collection of pure strategies gets chosen is
n
Y
s(ϕ) = si (ϕ), si (ϕ) = ski i , ϕ = (k1 , . . . , kn ) ∈ Φ (2.43)
i=1
(assuming all players make their choice independently) and the expected payoff
for player i is
X
Ri (s) = s(ϕ)Ri (ϕ). (2.44)
ϕ∈Φ
Theorem 2.17 (Nash) Every n-person game has at least one Nash equilibrium.
Theorem 2.19 (Product formula) Let U ⊆ Rn be a bounded and open set and
denote by Gj the connected components of Rn \f (∂U ). If g ◦ f ∈ Dy (U, Rn ), then
X
deg(g ◦ f, U, y) = deg(f, U, Gj ) deg(g, Gj , y), (2.52)
j
m
X X
deg(g ◦ f, U, y) = sgn(Jg (z)) deg(f, U, z)
j=1 z∈g −1 (y)∩Gj
m
X X
= deg(f, U, Gj ) sgn(Jg (z)) (2.53)
j=1 z∈g −1 (y)∩Gj
m
X
= deg(f, U, Gj ) deg(g, Gj , y). (2.54)
j=1
Moreover, this formula still holds for y ∈ CV(g ◦ f ) and for g ∈ C by construction
of the Brouwer degree. However, the case f ∈ C will need a closer investigation
since the sets Gj depend on f . To overcome this problem we will introduce the
sets
Ll = {z ∈ Rn \f (∂U )| deg(f, U, z) = l}. (2.55)
Observe that Ll , l > 0, must be a union of some sets of {Gj }m
j=1 .
Now choose f˜ ∈ C such that |f (x) − f˜(x)| < 2 dist(g (y), f (∂U )) for x ∈
1 −1 −1
U and define K̃j , L̃l accordingly. Then we have Ul ∩ g −1 (y) = Ũl ∩ g −1 (y) by
2.7. The Jordan curve theorem 31
and hence fore every i we have Ki ⊆ Gl for some l since components are maximal
connected sets. Let Nl = {i|Ki ⊆ Gl } and observe that we have deg(k, Gl , y) =
32 Chapter 2. The Brouwer mapping degree
P
i∈Nl deg(k, Ki , y) and deg(h, Hj , Gl ) = deg(h, Hj , Ki ) for every i ∈ Nl . There-
fore,
XX X
1= deg(h, Hj , Ki ) deg(k, Ki , y) = deg(h, Hj , Ki ) deg(k, Ki , Hj )
l i∈Nl i
(2.59)
By reversing the role of C1 and C2 , the same formula holds with Hj and Ki
interchanged.
Hence X XX X
1= deg(h, Hj , Ki ) deg(k, Ki , Hj ) = 1 (2.60)
i i j j
n n
shows that if the number of components of R \C1 or R \C2 is finite, then so is
the other and both are equal. Otherwise there is nothing to prove. 2
Chapter 3
33
34 Chapter 3. The Leray–Schauder mapping degree
Our next aim is to tackle the infinite dimensional case. The general idea is to
approximate F by finite dimensional operators (in the same spirit as we approx-
imated continuous f by smooth functions). To do this we need to know which
operators can be approximated by finite dimensional operators. Hence we have to
recall some basic facts first.
then
n
X n
X
|Pε (x) − x| = | φi (x)x − φi (x)xi | (3.4)
i=1 i=1
n
X
≤ φi (x)|x − xi | ≤ ε.
i=1
2
This lemma enables us to prove the following important result.
3.3. The Leray–Schauder mapping degree 35
Theorem 3.2 Let U be bounded, then the closure of F(U, Y ) in C(U, Y ) is C(U, Y ).
|FN (xn ) − FN (xm )| ≥ |F (xn ) − F (xm )| − |FN (xn ) − F (xn )| − |FN (xm ) − F (xm )|
ρ ρ
≥ ρ−2 = (3.5)
4 2
Theorem 3.4 Let U be a bounded open subset of a (real) Banach space X and let
F ∈ Dy (U , X), y ∈ X. Then the following hold true.
(ii). deg(1l, U, y) = 1 if y ∈ U .
(iii). If U1,2 are open, disjoint subsets of U such that y 6∈ f (U \(U1 ∪ U2 )), then
deg(1l + F, U, y) = deg(1l + F, U1 , y) + deg(1l + F, U2 , y).
(iv). If H : [0, 1]×U → X and y : [0, 1] → X are both continuous such that H(t) ∈
Dy(t) (U, Rn ), t ∈ [0, 1], then deg(1l + H(0), U, y(0)) = deg(1l + H(1), U, y(1)).
Proof. Except for (iv) all statements follow easily from the definition of the
degree and the corresponding property for the degree in finite dimensional spaces.
Considering H(t, x) − y(t), we can assume y(t) = 0 by (i). Since H([0, 1], ∂U )
is compact, we have ρ = dist(y, H([0, 1], ∂U ) > 0. By Theorem 3.2 we can pick
H1 ∈ F([0, 1] × U, X) such that |H(t) − H1 (t)| < ρ, t ∈ [0, 1]. this implies
deg(1l+H(t), U, 0) = deg(1l+H1 (t), U, 0) and the rest follows from Theorem 2.2. 2
In addition, Theorem 2.1 and Theorem 2.2 hold for the new situation as well
(no changes are needed in the proofs).
3.4. The Leray–Schauder principle and the Schauder fixed-point theorem 37
Theorem 3.5 Let F, G ∈ Dy (U, X), then the following statements hold.
(ii). If y 6∈ (1l + F )(U ), then deg(1l + F, U, y) = 0 (but not the other way round).
Equivalently, if deg(1l + F, U, y) 6= 0, then y ∈ (1l + F )(U ).
(iii). If |f (x) − g(x)| < dist(y, f (∂U )), x ∈ ∂U , then deg(f, U, y) = deg(g, U, y).
In particular, this is true if f (x) = g(x) for x ∈ ∂U .
Theorem 3.7 (Schauder fixed point) Let K be a closed, convex, and bounded
subset of a Banach space X. If F ∈ C(K, K), then F has at least one fixed
point. The result remains valid if K is only homeomorphic to a closed, convex,
and bounded subset.
38 Chapter 3. The Leray–Schauder mapping degree
Theorem 3.8 Let U ⊂ X be open and bounded and let F ∈ C(U , X). Suppose
there is an x0 ∈ U such that
Corollary 3.9 Let U ⊂ X be open and bounded and let F ∈ C(U , X). Then F
has a fixed point if one of the following conditions holds.
Proof. (1). F (∂U ) ⊆ U and F (x) = αx for |x| = ρ implies |α|ρ ≤ ρ and hence
(3.9) holds. (2). F (x) = αx for |x| = ρ implies (α − 1)2 ρ2 ≥ (α2 − 1)ρ2 and hence
α ≤ 0. (3). Special case of (2) since |F (x) − x|2 = |F (x)|2 − 2hF (x), xi + |x|2 . 2
3.5. Applications to integral and differential equations 39
for |x−x0 | < δ. In other words, F is continuous. Next we note that if U ⊂ C(I, Rn )
is bounded, say |U | < ρ, then
Z
τ (t)
|F (U )| ≤ sup f (t, s, x(s))ds ≤ (b − a)M, (3.12)
x∈U a
f (t0 , s, x)| ≤ ε1 and |τ (t)−τ (t0 )| ≤ ε2 for |t−t0 | < δ. Hence we infer for |t−t0 | < δ
Z
τ (t) Z τ (t0 )
|F (x)(t) − F (x)(t0 )| = f (t, s, x(s))ds − f (t0 , s, x(s))ds
a a
Z τ (t0 )
Z τ (t)
≤ |f (t, s, x(s)) − f (t0 , s, x(s))|ds + |f (t, s, x(s))|ds
a τ (t0 )
≤ (b − a)ε1 + ε2 M = ε. (3.13)
Theorem 3.11 Let F be as in the previous lemma. Then the integral equation
has at least one solution x ∈ C(I, Rn ) if |λ| ≤ ρ/M (ρ), where M (ρ) = (b −
a) max(s,t,x)∈I×I×Bρ (0) |f (s, t, x − y(s))| and ρ > 0 is arbitrary.
Proof. Note that, by our assumption on λ, λF maps Bρ (y) into itself. Now
apply the Schauder fixed-point theorem. 2
This result immediately gives the Peano theorem for ordinary differential equa-
tions.
and the first part follows from our previous theorem. To show the second, fix ε > 0
and assume M (ε, ρ) ≤ M̃ (ε)(1 + ρ). Then
Z t Z t
|x(t)| ≤ |f (s, x(s))|ds ≤ M̃ (ε) (1 + |x(s)|)ds (3.17)
0 0
43
44 Chapter 4. The stationary Navier–Stokes equation
η∆V ∂i ∂i vj . (4.2)
d
ρ∆V vj (t, x(t)) = η∆V ∂i ∂i vj (t, x(t)) − (∂j p(t, x(t)) + ∆V Kj (t, x(t)), (4.3)
dt
where ρ > 0 is the density of the fluid. Dividing by ∆V and using the chain rule
yields the Navier–Stokes equation
Taking the completion with respect to the associated norm we obtain the Sobolev
space H 1 (U, R). Similarly, taking the completion of C01 (U, R) with respect to the
same norm, we obtain the Sobolev space H01 (U, R). Here C0r (U, Y ) denotes the
set of functions in C r (U, Y ) with compact support. This construction of H 1 (U, R)
implies that a sequence uk in C 1 (U, R) converges to u ∈ H 1 (U, R) if and only if uk
and all its first order derivatives ∂j uk converge in L2 (U, R). Hence we can assign
each u ∈ H 1 (U, R) its first order derivatives ∂j u by taking the limits from above. In
order to show that this is a useful generalization of the ordinary derivative, we need
to show that the derivative depends only on the limiting function u ∈ L2 (U, R).
To see this we need the following lemma.
Lemma 4.1 (Integration by parts) Suppose u ∈ H01 (U, R) and v ∈ H 1 (U, R),
then Z Z
u(∂j v)dx = − (∂j u)v dx. (4.9)
U U
And since C0∞ (U, R) is dense in L2 (U, R), the derivatives are uniquely determined
by u ∈ L2 (U, R) alone. Moreover, if u ∈ C 1 (U, R), then the derivative in the
46 Chapter 4. The stationary Navier–Stokes equation
Proof. Again we can assume u ∈ C01 (U, R) and we assume j = 1 for notational
convenience. Replace U by a set K = [a, b] × K̃ containing U and extend u to K
by setting it 0 outside U . Then we have
Z x1 2
2
u(x1 , x2 , . . . , xn ) = 1 · (∂1 u)(ξ, x2 , . . . , xn )dξ
a
Z b
≤ (b − a) (∂1 u)2 (ξ, x2 , . . . , xn )dξ, (4.12)
a
where we have used the Cauchy-Schwarz inequality. Integrating this result over
[a, b] gives
Z b Z b
2 2
u (ξ, x2 , . . . , xn )dξ ≤ (b − a) (∂1 u)2 (ξ, x2 , . . . , xn )dξ (4.13)
a a
This scalar product will be more convenient for our purpose and hence we will use
it from now on. (However, all results stated will hold in either case.) The norm
corresponding to this scalar product will be denoted by |.|.
4.2. An insert on Sobolev spaces 47
Next, we want to consider the embedding H01 (U, R) ,→ L2 (U, R) a little closer.
This embedding is clearly continuous since by the Poincaré-Friedrichs inequality
we have
d(U )
|u|2 ≤ √ |u|, d(U ) = sup{|x − y| |x, y ∈ U }. (4.15)
n
Moreover, by a famous result of Rellich, it is even compact. To see this we first
prove the following inequality.
where xi = (x̃1 , . . . , x̃i , xi+1 , . . . , xn ). Squaring this equation and using Cauchy–
Schwarz on the right hand side we obtain
n Z
!2 n Z 2
X 1 X 1
2 2
u(x) − 2u(x)u(x̃) + u(x̃) ≤ |∂i u|dxi ≤n |∂i u|dxi
i=1 0 i=1 0
Xn Z 1
≤ n (∂i u)2 dxi . (4.18)
i=1 0
is compact.
Proof. Pick a cube Q (with edge length ρ) containing U and a bounded sequence
u ∈ H01 (U, R). Since bounded sets are weakly compact, it is no restriction to
k
The last term can be made arbitrarily small by picking N large. The first term
converges to 0 since uk converges weakly and each summand contains the L2 scalar
product of uk − u` and χQi (the characteristic function of Qi ). 2
In addition to this result we will also need the following interpolation inequality.
Proof. We first prove the case where u ∈ C01 (U, R). The key idea is to start with
U ⊂ R1 and then work ones way up to U ⊂ R2 and U ⊂ R3 .
If U ⊂ R1 we have
Z x Z
2 2
u(x) = ∂1 u (x1 )dx1 ≤ 2 |u∂1 u|dx1 (4.24)
and hence Z
2
max u(x) ≤ 2 |u∂1 u|dx1 . (4.25)
x∈U
4.2. An insert on Sobolev spaces 49
Here, if an integration limit is missing, it means that the integral is taken over the
whole support of the function.
If U ⊂ R2 we have
ZZ Z Z
u dx1 dx2 ≤ max u(x, x2 ) dx2 max u(x1 , y)2 dx1
4 2
x y
ZZ ZZ
≤4 |u∂1 u|dx1 dx2 |u∂2 u|dx1 dx2
ZZ 2/2 ZZ 1/2 ZZ 1/2
2 2 2
≤4 u dx1 dx2 (∂1 u) dx1 dx2 (∂2 u) dx1 dx2
ZZ ZZ
≤4 u2 dx1 dx2 ((∂1 u)2 + (∂2 u)2 )dx1 dx2 (4.26)
and applying Cauchy–Schwarz finishes the proof for u ∈ C01 (U, R).
If u ∈ H01 (U, R) pick a sequence uk in C01 (U, R) which converges to u in H01 (U, R)
and hence in L2 (U, R). By our inequality, this sequencepis CauchyR in L4 (U, R)
and
qR converges to a limit v ∈ L4 (U, R). Since |u|2 ≤ 4 |U ||u|4 ( 1 · u2 dx ≤
R
1 dx u4 dx), uk converges to v in L2 (U, R) as well and hence u = v. Now take
the limit in the inequality for uk . 2
As a consequence we obtain
1/4
8d(U )
|u|4 ≤ √ |u|, U ⊂ R3 , (4.28)
3
and
is compact.
50 Chapter 4. The stationary Navier–Stokes equation
Proof. Let uk be a bounded sequence in H01 (U, R). By Rellich’s theorem there
is a subsequence converging in L2 (U, R). By the Ladyzhenskaya inequality this
subsequence converges in L4 (U, R). 2
Our analysis clearly extends to functions with values in Rn since H01 (U, Rn ) =
⊕nj=1 H01 (U, R).
where Z
a(u, v, w) = uk vj (∂k wj ) dx. (4.35)
U
In other words, (4.34) represents a necessary solubility condition for the Navier-
Stokes equations. A solution of (4.34) will also be called a weak solution of the
Navier-Stokes equations. If we can show that a weak solution is in C 2 , then we can
read our argument backwards and it will be also a classical solution. However, in
4.3. Existence and uniqueness of solutions 51
general this might not be true and it will only solve the Navier-Stokes equations in
the sense of distributions. But let us try to show existence of solutions for (4.34)
first.
For later use we note
Z Z
1
a(v, v, v) = vk vj (∂k vj ) dx = vk ∂k (vj vj ) dx
U 2 U
Z
1
= − (vj vj )∂k vk dx = 0, v ∈ X. (4.36)
2 U
R
We proceed by studying (4.34). Let K ∈ L2 (U, R3 ), then U Kw dx is a linear
functional on X and hence there is a K̃ ∈ X such that
Z
Kw dx = hK̃, wi, w ∈ X. (4.37)
U
Moreover, the same is true for the map a(u, v, .), u, v ∈ X , and hence there is an
element B(u, v) ∈ X such that
and hence
ηv − B(v, v) = K̃. (4.40)
So in order to apply the theory from our previous chapter, we need a Banach space
Y such that X ,→ Y is compact.
Let us pick Y = L4 (U, R3 ). Then, applying the Cauchy-Schwarz inequality
twice to each summand in a(u, v, w) we see
XZ 1/2 Z 1/2
2 2
|a(u, v, w)| ≤ (uk vj ) dx (∂k wj ) dx
j,k U U
XZ 1/4 Z 1/4
4
≤ |w| (uk ) dx (vj )4 dx = |u|4 |v|4 |w|. (4.41)
j,k U U
Theorem 4.7 Let X be a Hilbert space, Y a Banach space, and suppose there is
a compact embedding X ,→ Y . In particular, |u|Y ≤ β|u|. Let a : X 3 → R be a
multilinear form such that
|a(u, v, w)| ≤ α|u|Y |v|Y |w| (4.42)
and a(v, v, v) = 0. Then for any K̃ ∈ X , η > 0 we have a solution v ∈ X to the
problem
ηhv, wi − a(v, v, w) = hK̃, wi, w ∈ X. (4.43)
Moreover, if 2αβ|K̃| < η 2 this solution is unique.
Proof. It is no loss to set η = 1. Arguing as before we see that our equation is
equivalent to
v − B(v, v) + K̃ = 0, (4.44)
where our assumption (4.42) implies
|B(u, v)| ≤ α|u|Y |v|Y ≤ αβ 2 |u||v| (4.45)
Here the second equality follows since the embedding X ,→ Y is continuous.
Abbreviate F (v) = B(v, v). Observe that F is locally Lipschitz continuous
since if |u|, |v| ≤ ρ we have
|F (u) − F (v)| = |B(u − v, u) − B(v, u − v)| ≤ 2α ρ |u − v|Y ≤ 2αβ 2 ρ|u − v|. (4.46)
Moreover, let vn be a bounded sequence in X . After passing to a subsequence we
can assume that vn is Cauchy in Y and hence F (vn ) is Cauchy in X by |F (u) −
F (v)| ≤ 2α ρ|u − v|Y . Thus F : X → X is compact.
Hence all we need to apply the Leray-Schauder principle is an a priori estimate.
Suppose v solves v = tF (v) + tK̃, t ∈ [0, 1], then
hv, vi = t a(v, v, v) + thK̃, vi = thK̃, vi. (4.47)
Hence |v| ≤ |K̃| is the desired estimate and the Leray-Schauder principle yields
existence of a solution.
Now suppose there are two solutions vi , i = 1, 2. By our estimate they satisfy
|vi | ≤ |K̃| and hence |v1 − v2 | = |F (v1 ) − F (v2 )| ≤ 2αβ 2 |K̃||v1 − v2 | which is a
contradiction if 2αβ 2 |K̃| < 1. 2
Hence we have found a solution v to the generalized problem (4.34). This
solution is unique if 2( 2d(U )
√ )3/2 |K|2 < η 2 . Under suitable additional conditions on
3
the outer forces and the domain, it can be shown that weak solutions are C 2 and
thus also classical solutions. However, this is beyond the scope of this introductory
text.
Chapter 5
Monotone operators
53
54 Chapter 5. Monotone operators
strictly monotone if
Note that the same definitions can be made if X is a Banach space and F :
X → X ∗.
Observe that if F is strongly monotone, then it automatically satisfies
hF (x), xi
lim = ∞. (5.7)
|x|→∞ |x|
(Just take y = 0 in the definition of strong monotonicity.) Hence the following
result is not surprising.
F (x) = y (5.8)
Proof. Set
G(x) = x − t(F (x) − y), t > 0, (5.9)
then F (x) = y is equivalent to the fixed point equation
G(x) = x. (5.10)
where L is a Lipschitz constant for F (i.e., |F (x) − F (x̃)| ≤ L|x − x̃|). Thus, if t ∈
(0, 2C
L
), G is a contraction and the rest follows from the contraction principle. 2
Again observe that our proof is constructive. In fact, the best choice for t is
clearly t = CL such that the contraction constant θ = 1 − ( CL )2 is minimal. Then
the sequence
C
xn+1 = xn − (1 − ( )2 )(F (xn ) − y), x0 = x, (5.12)
L
converges to the solution.
and
|a(x, z) − a(y, z)| ≤ L|z||x − y|. (5.15)
Then there is a unique x ∈ X such that (5.13) holds.
Proof. By the Riez theorem there are elements F (x) ∈ X and z ∈ X such that
a(x, y) = b(y) is equivalent to hF (x) − z, yi = 0, y ∈ X, and hence to
F (x) = z. (5.16)
as underlying Hilbert space. Next we multiply (5.19) by v ∈ H01 and integrate over
U
Z Z
∂i Aij (x)∂j u(x) + bj (x)∂j u(x) + c(x)u(x) v(x) dx = f (x)v(x) dx. (5.22)
U U
where
Z
a(v, u) = ∂i v(x)Aij (x)∂j u(x) + bj (x)v(x)∂j u(x) + c(x)v(x)u(x) dx
ZU
f (v) = f (x)v(x) dx, (5.24)
U
where −b0 = inf b(x), c0 = inf c(x) and we need to control the middle term. If
b0 ≤ 0 there is nothing to do and it suffices to require c0 ≥ 0.
If b0 > 0 we distribute the middle term by means of the elementary inequality
ε 1
|u||∂j u| ≤ |u|2 + |∂j u|2 (5.26)
2 2ε
which gives Z
b0 2 εb0 2
a(u, u) ≥ (a0 − )|∂j u| + (c0 − )|u| . (5.27)
U 2ε 2
b0
Since we need a0 − 2ε > 0 and c0 − εb20 ≥ 0, or equivalently 2cb00 ≥ ε > 2a
b0
0
, we see
2
that we can apply the Lax–Milgram theorem if 4a0 c0 > b0 . In summary, we have
proven
Theorem 5.4 The elliptic Dirichlet problem (5.19) has a unique weak solution
u ∈ H01 (U, R) if a0 > 0, b0 ≤ 0, c0 ≥ 0 or 4a0 c0 > b20 .
hF (x), xi
lim = ∞. (5.28)
|x|→∞ |x|
whenever xn → x.
58 Chapter 5. Monotone operators
Fn (xn ) = yn , xn , yn ∈ Xn . (5.30)
More precisely, let Pn be the (linear) projection onto Xn and set Fn (xn ) = Pn F (xn ),
yn = Pn y (verify that Fn is continuous and monotone!).
Now Lemma 5.1 ensures that there exists a solution S un . Now chose the sub-
spaces Xn such that Xn → X (i.e., Xn ⊂ Xn+1 and ∞ n=1 Xn is dense). Then our
hope is that un converges to a solution u.
This approach is quite common when solving equations in infinite dimensional
spaces and is known as Galerkin approximation. It can often be used for
numerical computations and the right choice of the spaces Xn will have a significant
impact on the quality of the approximation.
So how should we show that xn converges? First of all observe that our con-
struction of xn shows that xn lies in some ball with radius Rn , which is chosen
such that
hFn (x), xi > |yn ||x|, |x| ≥ Rn , x ∈ Xn . (5.31)
Since hFn (x), xi = hPn F (x), xi = hF (x), Pn xi = hF (x), xi for x ∈ Xn we can
drop all n’s to obtain a constant R which works for all n. So the sequence xn is
uniformly bounded
|xn | ≤ R. (5.32)
Now by a well-known result there exists a weakly convergent subsequence. That is,
after dropping some terms, we can assume that there is some x such that xn * x,
that is,
hxn , zi → hx, zi, for every z ∈ X. (5.33)
And it remains to show that x is indeed a solution. This follows from
implies F (x) = y.
t 7→ hF (x + ty), zi (5.37)
is continuous for t ∈ [0, 1] and all x, y, z ∈ X, since this condition together with
monotonicity can be shown to imply weak continuity.
60 Chapter 5. Monotone operators
Bibliography
61
62 Bibliography
63
64 Glossary of notations
Glossary of notations
66
Index 67
stationary, 44
n-person game, 27
Payoff, 27
Peano theorem, 40
Poincaré inequality, 47
Poincaré-Friedrichs inequality, 46
Prisoners dilemma, 28
Proper, 35
Reduction property, 29
Regular values, 13
Regulated function, 4
Rellich’s compactness theorem, 48
Rouchés theorem, 12
Sard’s theorem, 17
Simple function, 4
Stokes theorem, 19
Strategy, 27
Symmetric multilinear function, 3