MATH 8430
Fundamental Theory of Ordinary Differential Equations
Lecture Notes
Julien Arino Department of Mathematics University of Manitoba
Fall 2006
Contents
1 General theory of ODEs 3 1.1 ODEs, IVPs, solutions ...... 3 1.1.1 Ordinary differential equation, initial value problem ...... 3 1.1.2 Solutions to an ODE ...... 4 1.1.3 Geometric interpretation ...... 8 1.2 Existence and uniqueness theorems ...... 9 1.2.1 Successive approximations ...... 9 1.2.2 Local existence and uniqueness – Proof by fixed point ...... 10 1.2.3 Local existence and uniqueness – Proof by successive approximations 13 1.2.4 Local existence (non Lipschitz case) ...... 16 1.2.5 Some examples of existence and uniqueness ...... 21 1.3 Continuation of solutions ...... 24 1.3.1 Maximal interval of existence ...... 27 1.3.2 Maximal and global solutions ...... 28 1.4 Continuous dependence on initial data, on parameters ...... 29 1.5 Generality of first order systems ...... 32 1.6 Generality of autonomous systems ...... 34 1.7 Suggested reading, Further problems ...... 34
2 Linear systems 35 2.1 Existence and uniqueness of solutions ...... 35 2.2 Linear systems ...... 36 2.2.1 The vector space of solutions ...... 37 2.2.2 Fundamental matrix solution ...... 38 2.2.3 Resolvent matrix ...... 41 2.2.4 Wronskian ...... 43 2.2.5 Autonomous linear systems ...... 43 2.3 Affine systems ...... 46 2.3.1 The space of solutions ...... 46 2.3.2 Construction of solutions ...... 46 2.3.3 Affine systems with constant coefficients ...... 47 2.4 Systems with periodic coefficients ...... 48 2.4.1 Linear systems: Floquet theory ...... 48
i Fund. Theory ODE Lecture Notes – J. Arino ii CONTENTS
2.4.2 Affine systems: the Fredholm alternative ...... 51 2.5 Further developments, bibliographical notes ...... 54 2.5.1 A variation of constants formula for a nonlinear system with a linear component ...... 54
3 Stability of linear systems 57 3.1 Stability at fixed points ...... 57 3.2 Affine systems with small coefficients ...... 57
4 Linearization 63 4.1 Some linear stability theory ...... 63 4.2 The stable manifold theorem ...... 65 4.3 The Hartman-Grobman theorem ...... 69 4.4 Example of application ...... 74 4.4.1 A chemostat model ...... 74 4.4.2 A second example ...... 75
5 Exponential dichotomy 79 5.1 Exponential dichotomy ...... 79 5.2 Existence of exponential dichotomy ...... 81 5.3 First approximate theory ...... 82 5.4 Stability of exponential dichotomy ...... 84 5.5 Generality of exponential dichotomy ...... 84
References 85
A Definitions and results 87 A.1 Vector spaces, norms ...... 87 A.1.1 Norm ...... 87 A.1.2 Matrix norms ...... 87 A.1.3 Supremum (or operator) norm ...... 87 A.2 An inequality involving norms and integrals ...... 88 A.3 Types of convergences ...... 88 A.4 Asymptotic Notations ...... 89 A.5 Types of continuities ...... 89 A.6 Lipschitz function ...... 90 A.7 Gronwall’s lemma ...... 91 A.8 Fixed point theorems ...... 93 A.9 Jordan normal form ...... 94 A.10 Matrix exponentials ...... 96 A.11 Matrix logarithms ...... 97 A.12 Spectral theorems ...... 99 B Problem sheets 101 Homework sheet 1 – 2003 ...... 103 Homework sheet 2 – 2003 ...... 113 Homework sheet 3 – 2003 ...... 117 Homework sheet 4 – 2003 ...... 125 Final examination – 2003 ...... 137 Homework sheet 1 – 2006 ...... 145
iii
Introduction
This course deals with the elementary theory of ordinary differential equations. The word elementary should not be understood as simple. The underlying assumption here is that, to understand the more advanced topics in the analysis of nonlinear systems, it is important to have a good understanding of how solutions to differential equations are constructed. If you are taking this course, you most likely know how to analyze systems of nonlinear ordinary differential equations. You know, for example, that in order for solutions to a system to exist and be unique, the system must have a C1 vector field. What you do not necessarily know is why that is. This is the object of Chapter 1, where we consider the general theory of existence and uniqueness of solutions. We also consider the continuation of solutions as well as continuous dependence on initial data and on parameters. In Chapter 2, we explore linear systems. We first consider homogeneous linear systems, then linear systems in full generality. Homogeneous linear systems are linked to the theory for nonlinear systems by means of linearization, which we study in Chapter 4, in which we show that the behavior of nonlinear systems can be approximated, in the vicinity of a hyperbolic equilibrium point, by a homogeneous linear system. As for autonomous sys- tems, nonautonomous nonlinear systems are linked to a linearized form, this time through exponential dichotomy, which is explained in Chapter 5.
1
Chapter 1
General theory of ODEs
We begin with the general theory of ordinary differential equations (ODEs). First, we define ODEs, initial value problems (IVPs) and solutions to ODEs and IVPs in Section 1.1. In Section 1.2, we discuss existence and uniqueness of solutions to IVPs.
1.1 ODEs, IVPs, solutions
1.1.1 Ordinary differential equation, initial value problem Definition 1.1.1 (ODE). An nth order ordinary differential equation (ODE) is a functional relationship taking the form d d2 dn F t, x(t), x(t), x(t),..., x(t) = 0, dt dt2 dtn that involves an independent variable t ∈ I ⊂ R, an unknown function x(t) ∈ D ⊂ Rn of the independent variable, its derivative and derivatives of order up to n. For simplicity, the time dependence of x is often omitted, and we in general write equations as F t, x, x0, x00, . . . , x(n) = 0, (1.1) where x(n) denotes the nth order derivative of x. An equation such as (1.1) is said to be in general (or implicit) form. An equation is said to be in normal (or explicit) form when it is written as x(n) = f t, x, x0, x00, . . . , x(n−1) . Note that it is not always possible to write a differential equation in normal form, as it can be impossible to solve F (t, x, . . . , x(n)) = 0 in terms of x(n). Definition 1.1.2 (First-order ODE). In the following, we consider for simplicity the more restrictive case of a first-order ordinary differential equation in normal form x0 = f(t, x). (1.2)
3 Fund. Theory ODE Lecture Notes – J. Arino 4 1. General theory of ODEs
Note that the theory developed here holds usually for nth order equations; see Section 1.5. The function f is assumed continuous and real valued on a set U ⊂ R × Rn. Definition 1.1.3 (Initial value problem). An initial value problem (IVP) for equation (1.2) is given by x0 = f(t, x) (1.3) x(t0) = x0, n where f is continuous and real valued on a set U ⊂ R × R , with (t0, x0) ∈ U.
Remark – The assumption that f be continuous can be relaxed, piecewise continuity only is needed. However, this leads in general to much more complicated problems and is beyond the scope of this course. Hence, unless otherwise stated, we assume that f is at least continuous. The function f could also be complex valued, but this too is beyond the scope of this course. ◦
Remark – An IVP for an nth order differential equation takes the form x(n) = f(t, x, x0, . . . , x(n−1)) 0 0 (n−1) (n−1) x(t0) = x0, x (t0) = x0, . . . , x (t0) = x0 , i.e., initial conditions have to be given for derivatives up to order n − 1. ◦ We have already seen that the order of an ODE is the order of the highest derivative involved in the equation. An equation is then classified as a function of its linearity. A linear equation is one in which the vector field f takes the form f(t, x) = a(t)x(t) + b(t). If b(t) = 0 for all t, the equation is linear homogeneous; otherwise it is linear nonhomogeneous. If the vector field f depends only on x, i.e., f(t, x) = f(x) for all t, then the equation is autonomous; otherwise, it is nonautonomous. Thus, a linear equation is autonomous if a(t) = a and b(t) = b for all t. Nonlinear equations are those that are not linear. They too, can be autonomous or nonautonomous. Other types of classifications exist for ODEs, which we shall not deal with here, the previous ones being the only one we will need.
1.1.2 Solutions to an ODE Definition 1.1.4 (Solution). A function φ(t) (or φ, for short) is a solution to the ODE (1.2) if it satisfies this equation, that is, if φ0(t) = f(t, φ(t)), for all t ∈ I ⊂ R, an open interval such that (t, φ(t)) ∈ U for all t ∈ I. The notations φ and x are used indifferently for the solution. However, in this chapter, to emphasize the difference between the equation and its solution, we will try as much as possible to use the notation x for the unknown and φ for the solution. Fund. Theory ODE Lecture Notes – J. Arino 1.1. ODEs, IVPs, solutions 5
Definition 1.1.5 (Integral form of the solution). The function
Z t φ(t) = x0 + f(s, φ(s))ds (1.4) t0 is called the integral form of the solution to the IVP (1.3).
Let R = R((t0, x0), a, b) be the domain defined, for a > 0 and b > 0, by
R = {(t, x): |t − t0| ≤ a, kx − x0k ≤ b} , where k k is any appropriate norm of Rn. This domain is illustrated in Figures 1.1 and 1.2; it is sometimes called a security system, i.e., the union of a security interval (for the independent variable) and a security domain (for the dependent variables) [19]. Suppose
x
x +b 0
(t ,x ) x0 0 0
x0 −b
t −a t +a t 0 t 0 0
Figure 1.1: The domain R for D ⊂ R.
y
y 0
t t 0
x 0
x
2 Figure 1.2: The domain R for D ⊂ R : “security tube”. Fund. Theory ODE Lecture Notes – J. Arino 6 1. General theory of ODEs
that f is continuous on R, and let M = maxR kf(t, x)k, which exists since f is continuous on the compact set R. In the following, existence of solutions will be obtained generally in relation to the domain R by considering a subset of the time interval |t − t0| ≤ a defined by |t − t0| ≤ α, with a if M = 0 α = b min(a, M ) if M > 0. This choice of α = min(a, b/M) is natural. We endow f with specific properties (continuity, Lipschitz, etc.) on the domain R. Thus, in order to be able to use the definition of φ(t) as 0 the solution of x = f(t, x), we must be working in R. So we require that |t − t0| ≤ a and kx−x0k ≤ b. In order to satisfy the first of these conditions, choosing α ≤ a and working on |t − t0| ≤ α implies of course that |t − t0| ≤ a. The requirement that α ≤ b/M comes from the following argument. If we assume that φ(t) is a solution of (1.3) defined on [t0, t0 + α], then we have, for t ∈ [t0, t0 + α], Z t
kφ(t) − x0k = f(s, φ(s))ds t0 Z t ≤ kf(s, φ(s))k ds t0 Z t ≤ M ds t0
= M(t − t0), where the first inequality is a consequence of the definition of the integrals by Riemann sums (Lemma A.2.1 in Appendix A.2). Similarly, we have kφ(t) − x0k ≤ −M(t − t0) for all t ∈ [t0 − α, t0]. Thus, for |t − t0| ≤ α, kφ(t) − x0k ≤ M|t − t0|. Suppose now that α ≤ b/M. It follows that kφ − x0k ≤ M|t − t0| ≤ Mb/M = b. Taking α = min(a, b/M) then ensures that both |t − t0| ≤ a and kφ − x0k ≤ b hold simultaneously. The following two theorems deal with the localization of the solutions to an IVP. They make more precise the previous discussion. Note that for the moment, the existence of a solution is only assumed. First, we establish that the security system described above performs properly, in the sense that a solution on a smaller time interval stays within the security domain.
Theorem 1.1.6. If φ(t) is a solution of the IVP (1.3) in an interval |t − t0| < α˜ ≤ α, then kφ(t) − x0k < b in |t − t0| < α˜, i.e., (t, φ(t)) ∈ R((t0, x0), α,˜ b) for |t − t0| < α˜.
Proof. Assume that φ is a solution with (t, φ(t)) 6∈ R((t0, x0), α,˜ b). Since φ is continuous, it follows that there exists 0 < β < α˜ such that kφ(t) − x0k < b for |t − t0| < β and kφ(t0 + β) − x0k = b or kφ(t0 − β) − x0k = b , (1.5) i.e., the solution escapes the security domain at t = t0 ± β. Sinceα ˜ ≤ α ≤ a, β < a. Thus
(t, φ(t)) ∈ R for |t − t0| ≤ β. Fund. Theory ODE Lecture Notes – J. Arino 1.1. ODEs, IVPs, solutions 7
0 Thus kf(t, φ(t))k ≤ M for |t − t0| ≤ β. Since φ is a solution, we have that φ (t) = f(t, φ(t)) and φ(t0) = x0. Thus
Z t φ(t) = x0 + f(s, φ(s))ds for |t − t0| ≤ β. t0 Hence Z t
kφ(t) − x0k = f(s, φ(s))ds for |t − t0| ≤ β t0
≤ M|t − t0| for |t − t0| ≤ β.
As a consequence, b kφ(t) − x k ≤ Mβ < Mα˜ ≤ Mα ≤ M = b for |t − t | ≤ β. 0 M 0
In particular, kφ(t0 ± β) − x0k < b. Hence contradiction with (1.5). The following theorem is proved using the same sort of technique as in the proof of Theorem 1.1.6. It links the variation of the solution to the nature of the vector field.
Theorem 1.1.7. If φ(t) is a solution of the IVP (1.3) in an interval |t − t0| < α˜ ≤ α, then kφ(t1) − φ(t2)k ≤ M|t1 − t2| whenever t1, t2 are in the interval |t − t0| < α˜.
Proof. Let us begin by considering t ≥ t0. On t0 ≤ t ≤ t0 +α ˜,
Z t1 Z t2 φ(t1) − φ(t2) = x0 + f(s, φ(s))ds − x0 − f(s, φ(s))ds t0 t0 Z t2 = − f(s, φ(s))ds if t2 > t1 t1 Z t2 f(s, φ(s))ds if t1 > t2 t1
Now we can see formally what is needed for a solution.
n Theorem 1.1.8. Suppose f is continuous on an open set U ⊂ R × R . Let (t0, x0) ∈ U, and φ be a function defined on an open set I of R such that t0 ∈ I. Then φ is a solution of the IVP (1.3) if, and only if,
i) ∀t ∈ I, (t, φ(t)) ∈ U.
ii) φ is continuous on I. R t iii) ∀t ∈ I, φ(t) = x0 + f(s, φ(s))ds. t0 Fund. Theory ODE Lecture Notes – J. Arino 8 1. General theory of ODEs
0 Proof. (⇒) Let us suppose that φ = f(t, φ) for all t ∈ I and that φ(t0) = x0. Then for all t ∈ I,(t, φ(t)) ∈ U (i). Also, φ is differentiable and thus continuous on I (ii). Finally,
φ0(s) = f(s, φ(s)) so integrating both sides from t0 to t,
Z t φ(t) − φ(t0) = f(s, φ(s))ds t0 and thus
Z t φ(t) = x0 + f(s, φ(s))ds t0 hence (iii). (⇐) Assume i), ii) and iii). Then φ is differentiable on I and φ0(t) = f(t, φ(t)) for all t ∈ I. R t0 From (3), φ(t0) = x0 + f(s, φ(s))ds = x0. t0 Note that Theorem 1.1.8 states that φ should be continuous, whereas the solution should of course be C1, for its derivative needs to be continuous. However, this is implied by point iii). In fact, more generally, the following result holds about the regularity of solutions.
Theorem 1.1.9 (Regularity). Let f : U → Rn, with U an open set of R × Rn. Suppose that f ∈ Ck. Then all solutions of (1.2) are of class Ck+1.
Proof. The proof is obvious, since a solution φ is such that φ0 ∈ Ck.
1.1.3 Geometric interpretation The function f is the vector field of the equation. At every point in (t, x) space, a solution φ is tangent to the value of the vector field at that point. A particular consequence of this fact is the following theorem.
Theorem 1.1.10. Let x0 = f(x) be a scalar autonomous differential equation. Then the solutions of this equation are monotone.
Proof. The direction field of an autonomous scalar differential equation consists of vectors that are parallel for all t (since f(t, x) = f(x) for all t). Suppose that a solution φ of 0 x = f(x) is non monotone. Then this means that, given an initial point (t0, x0), one the following two occurs, as illustrated in Figure 1.3.
i) f(x0) 6= 0 and there exists t1 such that φ(t1) = x0.
ii) f(x0) = 0 and there exists t1 such that φ(t1) 6= x0. Fund. Theory ODE Lecture Notes – J. Arino 1.2. Existence and uniqueness theorems 9
tt t t tt 0 2 1 0 2 1 Figure 1.3: Situations that would lead to a scalar autonomous differential equation having nonmonotone solutions.
Suppose we are in case i), and assume we are in the case f(x0) > 0. Thus, the solution curve 0 φ is increasing at (t0, x0), i.e., φ (t0) > 0. As φ is continuous, i) implies that there exists t2 ∈ (t0, t1) such that φ(t2) is a maximum, with φ increasing for t ∈ [t0, t2) and φ decreasing 0 0 for t ∈ (t2, t1]. It follows that φ (t1) < 0, which is a contradiction with φ (t0) > 0. Now assume that we are in case ii). Then there exists t2 ∈ (t0, t1) with φ(t2) = x0 but 0 such that φ (t2) < 0. This is a contradiction.
Remark – If we have uniqueness of solutions, it follows from this theorem that if φ1 and φ2 are 0 two solutions of the scalar autonomous differential equation x = f(x), then φ1(t0) < φ2(t0) implies that φ1(t) < φ2(t) for all t. ◦
Remark – Be careful: Theorem 1.1.10 is only true for scalar equations. ◦
1.2 Existence and uniqueness theorems
Several approaches can be used to show existence and/or uniqueness of solutions. In Sec- tions 1.2.2 and 1.2.3, we take a direct path: using either a fixed point method (Section 1.2.2) or an iterative approach (Section 1.2.3), we obtain existence and uniqueness of solutions under the assumption that the vector field is Lipschitz. In Section 1.2.4, the Lipschitz assumption is dropped and therefore a different approach must be used, namely that of approximate solutions, with which only existence can be established.
1.2.1 Successive approximations Picard’s successive approximation method consists in using the integral form (1.4) of the solution to the IVP (1.3) to construct a sequence of approximation of the solution, that converges to the solution. The steps followed in constructing this approximating sequence are the following. Fund. Theory ODE Lecture Notes – J. Arino 10 1. General theory of ODEs
Step 1. Start with an initial estimate of the solution, say, the constant function φ0(t) = φ0 = x0, for |t − t0| ≤ h. Evidently, this function satisfies the IVP. Step 2. Use φ0 in (1.4) to define the second element in the sequence: Z t φ1(t) = x0 + f(s, φ0(s))ds. t0
Step 3. Use φ1 in (1.4) to define the third element in the sequence: Z t φ2(t) = x0 + f(s, φ1(s))ds. t0 ... Step n. Use φn−1 in (1.4) to define the nth element in the sequence: Z t φn(t) = x0 + f(s, φn−1(s))ds. t0 At this stage, there are two major ways to tackle the problem, which use the same idea: if we can prove that the sequence {φn} converges, and that the limit happens to satisfy the differential equation, then we have the solution to the IVP (1.3). The first method (Section 1.2.2) uses a fixed point approach. The second method (Section 1.2.3) studies explicitly the limit.
1.2.2 Local existence and uniqueness – Proof by fixed point Here are two slightly different formulations of the same theorem, which establishes that if the vector field is continuous and Lipschitz, then the solutions exist and are unique. We prove the result in the second case. For the definition of a Lipschitz function, see Section A.6 in the Appendix.
Theorem 1.2.1 (Picard local existence and uniqueness). Assume f : U ⊂ R × Rn → D ⊂ Rn is continuous, and that f(t, x) satisfies a Lipschitz condition in U with respect to x. Then, given any point (t0, x0) ∈ U, there exists a unique solution of (1.3) on some interval containing t0 in its interior. Theorem 1.2.2 (Picard local existence and uniqueness). Consider the IVP (1.3), and as- sume f is (piecewise) continuous in t and satisfies the Lipschitz condition
kf(t, x1) − f(t, x2)k ≤ Lkx1 − x2k for all x1, x2 ∈ D = {x : kx − x0k ≤ b} and all t such that |t − t0| ≤ a. Then there exists b 0 < δ ≤ α = min a, M such that (1.3) has a unique solution in |t − t0| ≤ δ. To set up the proof, we proceed as follows. Define the operator F by Z t F : x 7→ x0 + f(s, x(s))ds. t0 Fund. Theory ODE Lecture Notes – J. Arino 1.2. Existence and uniqueness theorems 11
Note that the function (F φ)(t) is a continuous function of t. Then Picard’s successives 2 2 approximations take the form φ1 = F φ0, φ2 = F φ1 = F φ0, where F represents F ◦ F . Iterating, the general term is given for k = 0,... by
k φk = F φ0.
Therefore, finding the limit limk→∞ φk is equivalent to finding the function φ, solution of the fixed point problem x = F x, with x a continuously differentiable function. Thus, a solution of (1.3) is a fixed point of F , and we aim to use the contraction mapping principle to verify the existence (and uniqueness) of such a fixed point. We follow the proof of [14, p. 56-58].
Proof. We show the result on the interval t − t0 ≤ δ. The proof for the interval t0 − t ≤ δ is similar. Let X be the space of continuous functions defined on the interval [t0, t0 + δ], X = C([t0, t0 + δ]), that we endow with the sup norm, i.e., for x ∈ X,
kxkc = max kx(t)k t∈[t0,t0+δ] Recall that this norm is the norm of uniform convergence. Let then
S = {x ∈ X : kx − x0kc ≤ b}
Of course, S ⊂ X. Furthermore, S is closed, and X with the sup norm is a complete metric space. Note that we have transformed the problem into a problem involving the space of continuous functions; hence we are now in an infinite dimensional case. The proof proceeds in 3 steps. Step 1. We begin by showing that F : S → S. From (1.4),
Z t (F φ)(t) − x0 = f(s, φ(s))ds t0 Z t = f(s, φ(s)) − f(s, x0) + f(s, x0)ds t0 Therefore, by the triangle inequality, Z t kF φ − x0k ≤ kf(s, φ(s)) − f(t, x0)k + kf(t, x0)kds t0
As f is (piecewise) continuous, it is bounded on [t0, t1] and there exists M = maxt∈[t0,t1] kf(t, x0)k. Thus Z t kF φ − x0k ≤ kf(s, φ(s)) − f(t, x0)k + Mds t0 Z t ≤ Lkφ(s) − x0k + Mds, t0 Fund. Theory ODE Lecture Notes – J. Arino 12 1. General theory of ODEs
since f is Lipschitz. Since φ ∈ S for all kφ − x0k ≤ b, we have that for all φ ∈ S,
Z t kF φ − x0k ≤ Lb + Mds t0
≤ (t − t0)(Lb + M)
As t ∈ [t0, t0 + δ], (t − t0) ≤ δ, and thus
kF φ − x0kc = max kF φ − x0k ≤ (Lb + M)δ [t0,t0+δ]
Choose then δ such that δ ≤ b/(Lb + M), i.e., t sufficiently close to t0. Then we have
kF φ − x0kc ≤ b
This implies that for φ ∈ S, F φ ∈ S, i.e., F : S → S. Step 2. We now show that F is a contraction. Let φ1, φ2 ∈ S, Z t
k(F φ1)(t) − (F φ2)(t)k = f(s, φ1(s)) − f(s, φ2(s))ds t0 Z t ≤ kf(s, φ1(s)) − f(s, φ2(s))kds t0 Z t ≤ Lkφ1(s) − φ2(s)kds t0 Z t ≤ Lkφ1 − φ2kc ds t0 and thus ρ kF φ − F φ k ≤ Lδkφ − φ k ≤ ρkφ − φ k for δ ≤ 1 2 c 1 2 c 1 2 c L Thus, choosing ρ < 1 and δ ≤ ρ/L, F is a contraction. Since, by Step 1, F : S → S, the contraction mapping principle (Theorem A.11) implies that F has a unique fixed point in S, and (1.3) has a unique solution in S. Step 3. It remains to be shown that any solution in X is in fact in S (since it is on X that we want to show the result). Considering a solution starting at x0 at time t0, the solution leaves S if there exists a t > t0 such that kφ(t)−x0k = b, i.e., the solution crosses the border of D. Let τ > t0 be the first of such t’s. For all t0 ≤ t ≤ τ, Z t kφ(t) − x0k ≤ kf(s, φ(s)) − f(s, x0)k + kf(s, x0)kds t0 Z t ≤ Lkφ(s) − x0k + Mds t0 Z t ≤ Lb + Mds t0 Fund. Theory ODE Lecture Notes – J. Arino 1.2. Existence and uniqueness theorems 13
As a consequence, b = kφ(τ) − x0k ≤ (Lb + M)(τ − t0)
As τ = t0 + µ, for some µ > 0, it follows that if b µ > Lb + M then the solution φ is confined to D.
Note that the condition x1, x2 ∈ D = {x : kx − x0k ≤ b} in the statement of the theorem refers to a local Lipschitz condition. If the function f is Lipschitz, then the following theorem holds. Theorem 1.2.3 (Global existence). Suppose that f is piecewise continuous in t and is Lipschitz on U = I × D. Then (1.3) admits a unique solution on I.
1.2.3 Local existence and uniqueness – Proof by successive ap- proximations Using the method of successive approximations, we can prove the following theorem. Theorem 1.2.4. Suppose that f is continuous on a domain R of the (t, x)-plane defined, for a, b > 0, by R = {(t, x): |t − t0| ≤ a, kx − x0k ≤ b}, and that f is locally Lipschitz in x on R. Let then, as previously defined,
M = sup kf(t, x)k < ∞ (t,x)∈R and b α = min(a, ) M Then the sequence defined by
φ0 = x0, |t − t0| ≤ α Z t φi(t) = x0 + f(s, φi−1(s))ds, i ≥ 1, |t − t0| ≤ α t0 converges uniformly on the interval |t − t0| ≤ α to φ, unique solution of (1.3). Proof. We follow [20, p. 3-6]. Existence. Suppose that |t − t0| ≤ α. Then Z t
kφ1 − φ0k = f(s, φ0(s))ds t0
≤ M|t − t0| ≤ αM ≤ b Fund. Theory ODE Lecture Notes – J. Arino 14 1. General theory of ODEs
R t from the definitions of M and α, and thus kφ1 − φ0k ≤ b. So f(s, φ1(s))ds is defined for t0 |t − t0| ≤ α, and, for |t − t0| ≤ α, Z t Z t
kφ2(t) − φ0k = f(s, φ1(s))ds ≤ k kf(s, φ1(s))kds ≤ αM ≤ b. t0 t0
All subsequent terms in the sequence can be similarly defined, and, by induction, for |t−t0| ≤ α, kφk(t) − φ0k ≤ αM ≤ b, k = 1, . . . , n.
Now, for |t − t0| ≤ α, Z t Z t
kφk+1(t) − φk(t)k = x0 + f(s, φk(s))ds − x0 − f(s, φk−1(s))ds t0 t0 Z t
= f(s, φk(s)) − f(s, φk−1(s)) ds t0 Z t ≤ L kφk(s) − φk−1(s)kds, t0 where the inequality results of the fact that f is locally Lipschitz in x on R. We now prove that, for all k,
(L|t − t |)k kφ − φ k ≤ b 0 for |t − t | ≤ α (1.6) k+1 k k! 0 Indeed, (1.6) holds for k = 1, as previously established. Assume that (1.6) holds for k = n. Then Z t
kφn+2 − φn+1k = f(s, φn+1(s)) − f(s, φn(s))ds t0 Z t ≤ Lkφn+1(s) − φn(s)kds t0 Z t n (L|s − t0|) ≤ Lb ds for |t − t0| ≤ α t0 n! n+1 n+1 s=t L |t − t0| ≤ b n! n + 1 s=t0 (L|t − t |)n+1 ≤ b 0 (n + 1)! and thus (1.6) holds for k = 1,.... Thus, for N > n we have
N−1 N−1 k N−1 k X X (L|t − t0|) X (Lα) kφ (t) − φ (t)k ≤ kφ (t) − φ (t)k ≤ b ≤ b N n k+1 k k! k! k=n k=n k=n Fund. Theory ODE Lecture Notes – J. Arino 1.2. Existence and uniqueness theorems 15
The rightmost term in this expression tends to zero as n → ∞. Therefore, {φk(t)} converges uniformly to a function φ(t) on the interval |t − t0| ≤ α. As the convergence is uniform, the PN limit function is continuous. Moreover φ(t0) = x0. Indeed, φN (t) = φ0(t) + k=1(φk(t) − P∞ φk−1(t)), so φ(t) = φ0(t) + k=1(φk(t) − φk−1(t)). The fact that φ is a solution of (1.3) follows from the following result. If a sequence of functions {φk(t)} converges uniformly and that the φk(t) are continuous on the interval |t − t0| ≤ α, then Z t Z t lim φn(s)ds = lim φn(s)ds n→∞ n→∞ t0 t0 Hence,
φ(t) = lim φn(t) n→∞ Z t = x0 + lim f(s, φn−1(s))ds n→∞ t0 Z t = x0 + lim f(s, φn−1(s))ds n→∞ t0 Z t = x0 + f(s, φ(s))ds, t0 which is to say that Z t φ(t) = x0 + f(s, φ(s))ds for |t − t0| ≤ α. t0 As the integrand f(t, φ) is a continuous function, φ is differentiable (with respect to t), and φ0(t) = f(t, φ(t)), so φ is a solution to the IVP (1.3). Uniqueness. Let φ and ψ be two solutions of (1.3), i.e., for |t − t0| ≤ α, Z t φ(t) = x0 + f(s, φ(s))ds t0 Z t ψ(t) = x0 + f(s, ψ(s))ds. t0
Then, for |t − t0| ≤ α, Z t
kφ(t) − ψ(t)k = f(s, φ(s)) − f(s, ψ(s))ds t0 Z t ≤ L kφ(s) − ψ(s)kds. (1.7) t0 We now apply Gronwall’s Lemma A.7) to this inequality, using K = 0 and g(t) = kφ(t) − ψ(t)k. First, applying the lemma for t0 ≤ t ≤ t0 + α, we get 0 ≤ kφ(t) − ψ(t)k ≤ 0, that is, kφ(t) − ψ(t)k = 0, Fund. Theory ODE Lecture Notes – J. Arino 16 1. General theory of ODEs
and thus φ(t) = ψ(t) for t0 ≤ t ≤ t0 + α. Similarly, for t0 − α ≤ t ≤ t0, kφ(t) − ψ(t)k = 0. Therefore, φ(t) = ψ(t) on |t − t0| ≤ α.
0 Example – Let us consider the IVP x = −x, x(0) = x0 = c, c ∈ R. For initial solution, we choose φ0(t) = c. Then Z t φ1(t) = x0 + f(s, φ0(s))ds 0 Z t = c + −φ0(s)ds 0 Z t = c − c ds 0 = c − ct.
To find φ2, we use φ1 in (1.4).
Z t φ2(t) = x0 + f(s, φ1(s))ds 0 Z t = c − (c − cs)ds 0 t2 = c − ct + c . 2 Continuing this method, we find a general term of the form
n X (−1)iti φ (t) = c . n i! i=0
−t −t This is the power series expansion of ce , so φn → φ = ce (and the approximation is valid on R), which is the solution of the initial value problem. Note that the method of successive approximations is a very general method that can be used in a much more general context; see [8, p. 264-269].
1.2.4 Local existence (non Lipschitz case) The following theorem is often called Peano’s existence theorem. Because the vector field is not assumed to be Lipschitz, something is lost, namely, uniqueness.
Theorem 1.2.5 (Peano). Suppose that f is continuous on some region
R = {(t, x): |t − t0| ≤ a, kx − x0k ≤ b}, with a, b > 0, and let M = maxR kf(t, x)k. Then there exists a continuous function φ(t), differentiable on R, such that
i) φ(t0) = x0, Fund. Theory ODE Lecture Notes – J. Arino 1.2. Existence and uniqueness theorems 17
0 ii) φ (t) = f(t, φ) on |t − t0| ≤ α, where
a if M = 0 α = b min a, M if M > 0.
Before we can prove this result, we need a certain number of preliminary notations and results. The definition of equicontinuity and a statement of the Ascoli lemma are given in Section A.5. To construct a solution without the Lipschitz condition, we approximate the differential equation by another one that does satisfy the Lipschitz condition. The unique solution of such an approximate problem is an ε-approximate solution. It is formally defined as follows [8, p. 285].
Definition 1.2.6 (ε-approximate solution). A differentiable mapping u of an open ball J ∈ I into U is an approximate solution of x0 = f(t, x) with approximation ε (or an ε-approximate solution) if we have ku0(t) − f(t, u(t))k ≤ ε, for any t ∈ J.
Lemma 1.2.7. Suppose that f(t, x) is continuous on a region
R = {(t, x): |t − t0| ≤ a, kx − x0k ≤ b}.
Then, for every positive number ε, there exists a function Fε(t, x) such that
i) Fε is continuous for |t − t0| ≤ a and all x,
ii) Fε has continuous partial derivatives of all orders with respect to x1, . . . , xn for |t−t0| ≤ a and all x,
iii) kFε(t, x)k ≤ maxR kf(t, x)k = M for |t − t0| ≤ a and all x,
iv) kFε(t, x) − f(t, x)k ≤ ε on R. See a proof in [12, p. 10-12]; note that in this proof, the property that f defines a differential equation is not used. Hence Lemma 1.2.7 can be used in a more general context than that of differential equations. We now prove Theorem 1.2.5. Proof of Theorem 1.2.5. The proof takes four steps. 1. We construct, for every positive number ε, a function Fε(t, x) that satisfies the re- quirements given in Lemma 1.2.7. Using an existence-uniqueness result in the Lipschitz case (such as Theorem 1.2.2), we construct a function φε(t) such that
(P1) φε(t0) = x0,
0 (P2) φε(t) = Fε(t, φε(t)) on |t − t0| < α.
(P3) (t, φε(t)) ∈ R on |t − t0| ≤ α. Fund. Theory ODE Lecture Notes – J. Arino 18 1. General theory of ODEs
2. The set F = {φε : ε > 0} is bounded and equicontinuous on |t − t0| ≤ α. Indeed, property (P3) of φε implies that F is bounded on |t − t0| ≤ α and that kFε(t, φε(t))k ≤ M on |t − t0| ≤ α. Hence property (P2) of φε implies that
kφε(t1) − φε(t2)k ≤ M|t1 − t2|, if |t1 − t0| ≤ α and |t2 − t0| ≤ α (this is a consequence of Theorem 1.1.7). Therefore, for a given positive number µ, we have kφε(t1) − φε(t2)k ≤ µ whenever |t1 − t0| ≤ α, |t2 − t0| ≤ α, and |t1 − t2| ≤ µ/M. 3. Using Lemma A.5, choose a sequence {εk : k = 1, 2,...} of positive numbers such that limk→∞ εk = 0 and that the sequence {φεk : k = 1, 2,...} converges uniformly on |t−t0| ≤ α as k → ∞. Then set φ(t) = lim φε (t) on |t − t0| ≤ α. k→∞ k 4. Observe that Z t φε(t) = x0 + Fε(s, φε(s))ds t0 Z t Z t = x0 + f(s, φε(s))ds + Fε(s, φε(s)) − f(s, φε(s))ds, t0 t0 and that it follows from iv) in Lemma 1.2.7 that Z t
Fε(s, φε(s)) − f(s, φε(s))ds ≤ ε|t − t0| on |t − t0| ≤ α. t0 This is true for all ε ≥ 0, so letting ε → 0, we obtain Z t φ(t) = x0 + f(s, φ(s))ds, t0 which completes the proof. See [12, p. 13-14] for the outline of two other proofs of this result. A proof, by Hartman [11, p. 10-11], now follows.
1 Proof. Let δ > 0 and φ`(t) be a C n-dimensional vector-valued function on [t0 − δ, t0] 0 0 satisfying φ`(t0) = x0, φ`(t0) = f(t0, x0) and kφ`(t) − x0k ≤ b, kφ`(t)k ≤ M. For 0 < ε ≤ δ, define a function φε(t) on [t0 − δ, t0 + α] by putting φε(t) = φ`(t) on [t0 − δ, t0] and Z t φε(t) = x0 + f(s, φε(s − ε))ds on [t0, t0 + α]. (1.8) t0
The function φε can indeed be thus defined on [t0 − δ, t0 + α]. To see this, remark first that this formula is meaningful and defines φε(t) for t0 ≤ t ≤ t0 + α1, α1 = min(α, ε), so that 1 φε(t) is C on [t0 − δ, t0 + α1] and, on this interval,
kφε(t) − x0k ≤ b, kφε(t) − φε(s)k ≤ M|t − s|. (1.9) Fund. Theory ODE Lecture Notes – J. Arino 1.2. Existence and uniqueness theorems 19
1 It then follows that (1.8) can be used to extend φε(t) as a C function over [t0 − δ, t0 + α2], where α2 = min(α, 2ε), satisfying relation (1.9). Continuing in this fashion, (1.8) serves to 0 define φε(t) over [t0, t0 +α] so that φε(t) is a C function on [t0 −δ, t0 +α], satisfying relation (1.9). 0 Since kφε(t)k ≤ M, M can be used as a Lipschitz constant for φε, giving uniform con- tinuity of φε. It follows that the family of functions, φε(t), 0 < ε ≤ δ, is equicontinuous. Thus, using Ascoli’s Lemma (Lemma A.5), there exists a sequence ε(1) > ε(2) > . . ., such that ε(n) → 0 as n → ∞ and
φ(t) = lim φε(n)(t) exists uniformly n→∞ on [t0 − δ, t0 + α]. The continuity of f implies that f(t, φε(n)(t − ε(n)) tends uniformly to f(t, φ(t)) as n → ∞; thus term-by-term integration of (1.8) where ε = ε(n) gives Z t φ(t) = x0 + f(s, φ(s))ds t0 and thus φ(t) is a solution of (1.3). An important corollary follows. Corollary 1.2.8. Let f(t, x) be continuous on an open set E and satisfy kf(t, x)k ≤ M. Let E0 be a compact subset of E. Then there exists an α = α(E,E0,M) > 0 with the property that if (t0, x0) ∈ E0, then the IVP (1.3) has a solution and every solution exists on |t − t0| ≤ α. In fact, hypotheses can be relaxed a little. Coddington and Levinson [5] define an ε- approximate solution as Definition 1.2.9. An ε-approximate solution of the differential equation (1.2), where f is continuous, on a t interval I is a function φ ∈ C on I such that i) (t, φ(t)) ∈ U for t ∈ I; ii) φ ∈ C1 on I, except possibly for a finite set of points S on I, where φ0 may have simple discontinuities (g has finite discontinuities at c if the left and right limits of g at c exist but are not equal); iii) kφ0(t) − f(t, φ(t))k ≤ ε for t ∈ I \ S. Hence it is assumed that φ has a piecewise continuous derivative on I, which is denoted 1 by φ ∈ Cp (I). Theorem 1.2.10. Let f ∈ C on the rectangle
R = {(t, x): |t − t0| ≤ a, kx − x0k ≤ b}.
Given any ε > 0, there exists an ε-approximate solution φ of (1.3) on |t − t0| ≤ α such that φ(t0) = x0. Fund. Theory ODE Lecture Notes – J. Arino 20 1. General theory of ODEs
Proof. Let ε > 0 be given. We construct an ε-approximate solution on the interval [t0, t0 +ε]; the construction works in a similar way for [t0 − α, t0]. The ε-approximate solution that we construct is a polygonal path starting at (t0, x0). Since f ∈ C on R, it is uniformly continuous on R, and therefore for the given value of ε, there exists δε > 0 such that
kf(t, φ) − f(t,˜ φ˜)k ≤ ε (1.10) if ˜ ˜ (t, φ) ∈ R, (t,˜ φ) ∈ R and |t − t˜| ≤ δε kφ − φk ≤ δε.
Now divide the interval [t0, t0 + α] into n parts t0 < t1 < ··· < tn = t0 + α, in such a way that δ max |t − t | ≤ min δ , ε . (1.11) k k−1 ε M
From (t0, x0), construct a line segment with slope f(t0, x0) intercepting the line t = t1 at (t1, x1). From the definition of α and M, it is clear that this line segment lies inside the triangular region T bounded by the lines segments with slopes ±M from (t0, x0) to their intercept at t = t0 + α, and the line t = t0 + α. In particular, (t1, x1) ∈ T . At the point (t1, x1), construct a line segment with slope f(t1, x1) until the line t = t2, obtaining the point (t2, x2). Continuing similarly, a polygonal path φ is constructed that meets the line t = t0 + α in a finite number of steps, and lies entirely in T . The function φ, which can be expressed as
φ(t ) = x 0 0 (1.12) φ(t) = φ(tk−1) + f(tk−1, φ(tk−1))(t − tk−1), t ∈ [tk−1, tk], k = 1, . . . , n,
1 is the ε-approximate solution that we seek. Clearly, φ ∈ Cp ([t0, t0 + α]) and
kφ(t) − φ(t˜)k ≤ M|t − t˜| for t, t˜∈ [t0, t0 + α]. (1.13)
If t ∈ [tk−1, tk], then (1.13) together with (1.11) imply that kφ(t) − φ(tk−1)k ≤ δε. But from (1.12) and (1.10),
0 kφ (t) − f(t, φ(t))k = kf(tk−1, φ(tk−1)) − f(t, φ(t))k ≤ ε.
Therefore, φ is an ε-approximation.
We can now turn to their proof of Theorem 1.2.5.
Proof. Let {εn} be a monotone decreasing sequence of positive real numbers with εn → 0 as n → ∞. By Theorem 1.2.10, for each εn, there exists an εn-approximate solution φn of (1.3) on |t − t0| ≤ α such that φn(t0) = x0. Choose one such solution φn for each εn. From (1.13), it follows that kφn(t) − φn(t˜)k ≤ M|t − t˜|. (1.14) Fund. Theory ODE Lecture Notes – J. Arino 1.2. Existence and uniqueness theorems 21
Applying (1.14) to t˜ = t0, it is clear that the sequence {φn} is uniformly bounded by kx0k + b, since |t − t0| ≤ b/M. Moreover, (1.14) implies that {φn} is an equicontinuous set. By Ascoli’s lemma (Lemma A.5), there exists a subsequence {φnk }, k = 1,..., of {φn}, converging uniformly on [t0 −α, t0 +α] to a limit function φ, which must be continuous since each φn is continuous. This limit function φ is a solution to (1.3) which meets the required specifications. To see this, write Z t φn(t) = x0 + f(s, φn(s)) + ∆n(s)ds, (1.15) t0 0 0 where ∆n(t) = φ (t) − f(t, φn(t)) at those points where φn exists, and ∆n(t) = 0 otherwise. Because φn is an εn-approximate solution, k∆n(t)k ≤ εn. Since f is uniformly continuous on
R, and φnk → φ uniformly on [t0 − α, t0 + α] as k → ∞, it follows that f(t, φnk ) → f(t, φ(t)) uniformly on [t0 − α, t0 + α] as k → ∞. Replacing n by nk in (1.15) and letting k → ∞ gives Z t φ(t) = x0 + f(s, φ(s))ds. (1.16) t0 0 Clearly, φ(t0) = 0, when evaluated using (1.16), and also φ (t) = f(t, φ(t)) since f is contin- uous. Thus φ as defined by (1.16) is a solution to (1.3) on |t − t0| ≤ α.
1.2.5 Some examples of existence and uniqueness
Example – Consider the IVP 0 2 x = 3|x| 3 (1.17) x(t0) = x0 Here, Theorem 1.2.5 applies, since f(t, x) = 3x2/3 is continuous. However, Theorem 1.2.2 does not apply, since f(t, x) is not locally Lipschitz in x = 0 (or, f is not Lipschitz on any interval containing 0). This means that we have existence of solutions to this IVP, but not uniqueness of the solution. The fact that f is not Lipschitz on any interval containing 0 is established using the following argument. Suppose that f is Lipschitz on an interval I = (−ε, ε), with ε > 0. Then, there exists L > 0 such that for all x1, x2 ∈ I,
kf(t, x1) − f(t, x2)k ≤ L|x1 − x2| that is, 2 2 3 3 3 |x1| − |x2| ≤ L|x1 − x2|
Since this has to hold true for all x1, x2 ∈ I, it must hold true in particular for x2 = 0. Thus 2 3|x1| 3 ≤ L|x1| 1 1 Given an ε > 0, it is possible to find Nε > 0 such that n < ε for all n ≥ Nε. Let x1 = n . Then for n ≥ Nε, if f is Lipschitz there must hold 2 1 3 L 3 ≤ n n Fund. Theory ODE Lecture Notes – J. Arino 22 1. General theory of ODEs
So, for all n ≥ Nε, 1 L n 3 ≤ 3 1/3 This is a contradiction, since limn→∞ n = ∞, and so f is not Lipschitz on I. Let us consider the set E = {t ∈ R : x(t) = 0} The set E can be have several forms, depending on the situation.
1) E = ∅,
2) E = [a, b], (closed since x is continuous and thus reaches its bounds),
3) E = (−∞, b),
4) E = (a, +∞),
5) E = R.
Note that case 2) includes the case of a single intersection point, when a = b, giving E = {a}. Let us now consider the nature of x in these different situations. Recall that from Theorem 1.1.10, since (1.17) is defined by a scalar autonomous equation, its solutions are monotone. For simplicity, we consider here the case of monotone increasing solutions. The case of monotone decreasing solutions can be treated in a similar fashion.
1) Here, there is no intersection with the x = 0 axis. Thus if follows that > 0, if x > 0 x(t) is 0 < 0, if x0 < 0
2) In this case, < 0, if t < a x(t) is = 0, if t ∈ [a, b] > 0, if t > b
3) Here, = 0, if t < b x(t) is > 0, if t > b
4) In this case, < 0, if t < a x(t) is = 0, if t > a
5) In this last case, x(t) = 0 for all t ∈ R. Now, depending on the sign of x, we can integrate the equation. First, if x > 0, then |x| = x and so x0 = 3x2/3 1 ⇔ x−2/3x0 = 1 3 1/3 ⇔ x = t + k1 3 ⇔ x(t) = (t + k1) Fund. Theory ODE Lecture Notes – J. Arino 1.2. Existence and uniqueness theorems 23
for k1 ∈ R. Then, if x < 0, then |x| = −x, and
x0 = 3 (−x)2/3 1 ⇔ (−x)−2/3(−x0) = −1 3 1/3 ⇔ (−x) = −t + k2 3 ⇔ x(t) = −(−t + k2) for k2 ∈ R. We can now use these computations with the different cases that were discussed earlier, depending on the value of t0 and x0. We begin with the case of t0 > 0 and x0 > 0.
1) The case E = ∅ is impossible, for all initial conditions (t0, x0). Indeed, as x0 > 0, we have 3 3 x(t) = (t + k1) . Using the initial condition, we find that x(t0) = x0 = (t0 + k1) , i.e., 1/3 1/3 3 k1 = x0 − t0, and x(t) = (t + x0 − t0) .
2) If E = [a, b], then 3 −(−t + k2) if t < a x(t) = 0 if t ∈ [a, b] 3 (t + k1) if t > b 3 Since x0 > 0, we have to be in the t > b region, so t0 > b, and (t0 + k1) = x0, which implies 1/3 that k1 = x0 − t0. Thus
−(−t + k )3 if t < a 2 x(t) = 0 if t ∈ [a, b] 1/3 3 (t + x0 − t0) if t > b Since x is continuous, 1/3 3 lim (t + x0 − t0) = 0 t→b,t>b and 3 lim −(−t + k2) = 0 t→a,t 1/3 This implies that b = t0 − x0 and k2 = a. So finally, −(−t + a)3 if t < a 1 1 3 3 x(t) = 0 if t ∈ [a, t0 − x0 ](a ≤ t0 − x0 ) 1 1/3 3 3 (t + x0 − t0) if t > t0 − x0 1/3 Thus, choosing a ≤ t0 − x0 , we have solutions of the form shown in Figure 1.4. Indeed, any ai satisfying this property yields a solution. 3) The case [a, +∞) is impossible. Indeed, there does not exist a solution through (t0, x0) such that x(t) = 0 for all t ∈ [a, +∞); since we are in the case of monotone increasing functions, if x0 > 0 then x(t) ≥ x0 for all t ≥ t0. 4) E = R is also impossible, for the same reason. Fund. Theory ODE Lecture Notes – J. Arino 24 1. General theory of ODEs Figure 1.4: Case t0, x0 > 0, subcase 2, in the resolution of (1.17). 5) For the case E = (−∞, b], we have 0 if t ∈ (−∞, b] x(t) = 3 (t + k1) if t > b 1/3 1/3 Since x(t0) = x0, k1 = x0 − t0, and since x is continuous, b = −k1 = t0 − x0 . So, ( 0 if t ∈ (−∞, t − x1/3] x(t) = 0 0 1/3 3 1/3 (t + x0 − t0) if t > t0 − x0 The other cases are left as an exercise. Example – Consider the IVP x0 = 2tx2 (1.18) x(0) = 0 Here, we have existence and uniqueness of the solutions to (1.18). Indeed, f(t, x) = 2tx2 is contin- uous and locally Lipschitz on R. 1.3 Continuation of solutions The results we have seen so far deal with the local existence (and uniqueness) of solutions to an IVP, in the sense that solutions are shown to exist in a neighborhood of the initial data. The continuation of solutions consists in studying criteria which allow to define solutions on possibly larger intervals. Consider the IVP x0 = f(t, x) (1.19) x(t0) = x0, Fund. Theory ODE Lecture Notes – J. Arino 1.3. Continuation of solutions 25 with f continuous on a domain U of the (t, x) space, and the initial point (t0, x0) ∈ U. Lemma 1.3.1. Let the function f(t, x) be continuous in an open set U in (t, x)-space, and assume that a function φ(t) satisfies the condition φ0(t) = f(t, φ(t)) and (t, φ(t)) ∈ U, in an open interval I = {t1 < t < t2}. Under this assumption, if limj→∞(τj, φ(τj)) = (t1, η) ∈ U for some sequence {τj : j = 1, 2,...} of points in the interval I, then limτ→t1 (τ, φ(τ)) = (t1, η). Similarly, if limj→∞(τj, φ(τj)) = (t2, η) ∈ U for some sequence {τj : j = 1, 2,...} of points in the interval I, then limτ→t2 (τ, φ(τ)) = (t2, η). Proof. Let W be an open neighborhood of (t1, η). Then (t, φ(t)) ∈ W in an interval τ1 < t < τ(W) for some τ(W) determined by W. Indeed, assume that the closure of W, W¯ ⊂ U, and that |f(t, x)| ≤ M in W for some positive number M. For every positive integer j and every positive number ε, consider a rectangular region Rj(ε) = {(t, x): |t − tj| ≤ ε, kx − φ(tj)k ≤ Mε}