Picard–Lindelöf theorem

Differential equations
Scope
Fields Natural sciences Engineering Astronomy Physics Chemistry Biology Geology Applied mathematics Continuum mechanics Chaos theory Dynamical systems Social sciences Economics Population dynamics List of named differential equations
Classification
Types Ordinary Partial Differential-algebraic Integro-differential Fractional Linear Non-linear By variable type Dependent and independent variables Autonomous Coupled / Decoupled Exact Nonhomogeneous Features Order Operator Notation
Relation to processes Difference (discrete analogue) Stochastic Stochastic partial Delay
Solution
Existence and uniqueness Picard–Lindelöf theorem Peano existence theorem Carathéodory's existence theorem Cauchy–Kowalevski theorem
General topics Initial conditions Boundary values Dirichlet Neumann Robin Cauchy problem Wronskian Phase portrait Asymptotic / Exponential stability Rate of convergence Series / Integral solutions Numerical integration Dirac delta function
Solution methods Inspection Method of characteristics Euler Exponential response formula Finite difference (Crank–Nicolson) Finite element Infinite element Finite volume Galerkin Petrov–Galerkin Green's function Integrating factor Integral transforms Perturbation theory Runge–Kutta Separation of variables Undetermined coefficients Variation of parameters
People
List Isaac Newton Gottfried Leibniz Jacob Bernoulli Leonhard Euler Józef Maria Hoene-Wroński Joseph Fourier Augustin-Louis Cauchy George Green Carl David Tolmé Runge Martin Kutta Rudolf Lipschitz Ernst Lindelöf Émile Picard Phyllis Nicolson John Crank
v t e

In mathematics, specifically the study of differential equations, the Picard–Lindelöf theorem gives a set of conditions under which an initial value problem has a unique solution. It is also known as Picard's existence theorem, the Cauchy–Lipschitz theorem, or the existence and uniqueness theorem.

The theorem is named after

.

Theorem

Let ${\displaystyle D\subseteq \mathbb {R} \times \mathbb {R} ^{n}}$ be a closed rectangle with ${\displaystyle (t_{0},y_{0})\in \operatorname {int} D}$, the interior of ${\displaystyle D}$. Let ${\displaystyle f:D\to \mathbb {R} ^{n}}$ be a function that is continuous in ${\displaystyle t}$ and Lipschitz continuous in ${\displaystyle y}$. Then, there exists some ε > 0 such that the initial value problem

${\displaystyle y'(t)=f(t,y(t)),\qquad y(t_{0})=y_{0}.}$

has a unique solution ${\displaystyle y(t)}$ on the interval ${\displaystyle [t_{0}-\varepsilon ,t_{0}+\varepsilon ]}$.^[1][2]

Proof sketch

The proof relies on transforming the differential equation, and applying the Banach fixed-point theorem. By integrating both sides, any function satisfying the differential equation must also satisfy the integral equation

${\displaystyle y(t)-y(t_{0})=\int _{t_{0}}^{t}f(s,y(s))\,ds.}$

A simple proof of existence of the solution is obtained by successive approximations. In this context, the method is known as

.

Set

${\displaystyle \varphi _{0}(t)=y_{0}}$

and

${\displaystyle \varphi _{k+1}(t)=y_{0}+\int _{t_{0}}^{t}f(s,\varphi _{k}(s))\,ds.}$

It can then be shown, by using the

to |φ(t) − ψ(t)|, where φ and ψ are two solutions, shows that φ(t) = ψ(t), thus proving the global uniqueness (the local uniqueness is a consequence of the uniqueness of the Banach fixed point).

See Newton's method of successive approximation for instruction.

Example of Picard iteration

Let ${\displaystyle y(t)=\tan(t),}$ the solution to the equation ${\displaystyle y'(t)=1+y(t)^{2}}$ with initial condition ${\displaystyle y(t_{0})=y_{0}=0,t_{0}=0.}$ Starting with ${\displaystyle \varphi _{0}(t)=0,}$ we iterate

${\displaystyle \varphi _{k+1}(t)=\int _{0}^{t}(1+(\varphi _{k}(s))^{2})\,ds}$

so that ${\displaystyle \varphi _{n}(t)\to y(t)}$:

${\displaystyle \varphi _{1}(t)=\int _{0}^{t}(1+0^{2})\,ds=t}$

${\displaystyle \varphi _{2}(t)=\int _{0}^{t}(1+s^{2})\,ds=t+{\frac {t^{3}}{3}}}$

${\displaystyle \varphi _{3}(t)=\int _{0}^{t}\left(1+\left(s+{\frac {s^{3}}{3}}\right)^{2}\right)\,ds=t+{\frac {t^{3}}{3}}+{\frac {2t^{5}}{15}}+{\frac {t^{7}}{63}}}$

and so on. Evidently, the functions are computing the Taylor series expansion of our known solution ${\displaystyle y=\tan(t).}$ Since ${\displaystyle \tan }$ has poles at ${\displaystyle \pm {\tfrac {\pi }{2}},}$ this converges toward a local solution only for ${\displaystyle |t|<{\tfrac {\pi }{2}},}$ not on all of ${\displaystyle \mathbb {R} }$.

Example of non-uniqueness

To understand uniqueness of solutions, consider the following examples.^[3] A differential equation can possess a stationary point. For example, for the equation dy/dt = ay (${\displaystyle a<0}$), the stationary solution is y(t) = 0, which is obtained for the initial condition y(0) = 0. Beginning with another initial condition y(0) = y₀ ≠ 0, the solution y(t) tends toward the stationary point, but reaches it only at the limit of infinite time, so the uniqueness of solutions (over all finite times) is guaranteed.

However, for an equation in which the stationary solution is reached after a finite time, the uniqueness fails. This happens for example for the equation dy/dt = ay^2/3, which has at least two solutions corresponding to the initial condition y(0) = 0 such as: y(t) = 0 or

${\displaystyle y(t)={\begin{cases}\left({\tfrac {at}{3}}\right)^{3}&t<0\\\ \ \ \ 0&t\geq 0,\end{cases}}}$

so the previous state of the system is not uniquely determined by its state after t = 0. The uniqueness theorem does not apply because the function  f (y) = y^2/3 has an infinite slope at y = 0 and therefore is not Lipschitz continuous, violating the hypothesis of the theorem.

Detailed proof

Let

${\displaystyle C_{a,b}={\overline {I_{a}(t_{0})}}\times {\overline {B_{b}(y_{0})}}}$

where:

${\displaystyle {\begin{aligned}{\overline {I_{a}(t_{0})}}&=[t_{0}-a,t_{0}+a]\\{\overline {B_{b}(y_{0})}}&=[y_{0}-b,y_{0}+b].\end{aligned}}}$

This is the compact cylinder where  f  is defined. Let

${\displaystyle M=\sup _{C_{a,b}}\|f\|,}$

this is, the

of) the slopes of the function. Finally, let L be the Lipschitz constant of  f  with respect to the second variable.

We will proceed to apply the

metric

on ${\displaystyle {\mathcal {C}}(I_{a}(t_{0}),B_{b}(y_{0}))}$ induced by the uniform norm

${\displaystyle \|\varphi \|_{\infty }=\sup _{t\in I_{a}}|\varphi (t)|.}$

We define an operator between two function spaces of continuous functions, Picard's operator, as follows:

${\displaystyle \Gamma :{\mathcal {C}}(I_{a}(t_{0}),B_{b}(y_{0}))\longrightarrow {\mathcal {C}}(I_{a}(t_{0}),B_{b}(y_{0}))}$

defined by:

${\displaystyle \Gamma \varphi (t)=y_{0}+\int _{t_{0}}^{t}f(s,\varphi (s))\,ds.}$

We must show that this operator maps a complete non-empty metric space X into itself and also is a contraction mapping.

We first show that, given certain restrictions on ${\displaystyle a}$, ${\displaystyle \Gamma }$ takes ${\displaystyle {\overline {B_{b}(y_{0})}}}$ into itself in the space of continuous functions with the uniform norm. Here, ${\displaystyle {\overline {B_{b}(y_{0})}}}$ is a closed ball in the space of continuous (and bounded) functions "centered" at the constant function ${\displaystyle y_{0}}$. Hence we need to show that

${\displaystyle \|\varphi -y_{0}\|_{\infty }\leq b}$

implies

${\displaystyle \left\|\Gamma \varphi (t)-y_{0}\right\|=\left\|\int _{t_{0}}^{t}f(s,\varphi (s))\,ds\right\|\leq \int _{t_{0}}^{t'}\left\|f(s,\varphi (s))\right\|ds\leq \int _{t_{0}}^{t'}M\,ds=M\left|t'-t_{0}\right|\leq Ma\leq b}$

where ${\displaystyle t'}$ is some number in ${\displaystyle [t_{0}-a,t_{0}+a]}$ where the maximum is achieved. The last inequality in the chain is true if we impose the requirement ${\displaystyle a<{\frac {b}{M}}}$.

Now let's prove that this operator is a contraction mapping.

Given two functions ${\displaystyle \varphi _{1},\varphi _{2}\in {\mathcal {C}}(I_{a}(t_{0}),B_{b}(y_{0}))}$, in order to apply the Banach fixed-point theorem we require

${\displaystyle \left\|\Gamma \varphi _{1}-\Gamma \varphi _{2}\right\|_{\infty }\leq q\left\|\varphi _{1}-\varphi _{2}\right\|_{\infty },}$

for some ${\displaystyle 0\leq q<1}$. So let ${\displaystyle t}$ be such that

${\displaystyle \|\Gamma \varphi _{1}-\Gamma \varphi _{2}\|_{\infty }=\left\|\left(\Gamma \varphi _{1}-\Gamma \varphi _{2}\right)(t)\right\|.}$

Then using the definition of ${\displaystyle \Gamma }$,

${\displaystyle {\begin{aligned}\left\|\left(\Gamma \varphi _{1}-\Gamma \varphi _{2}\right)(t)\right\|&=\left\|\int _{t_{0}}^{t}\left(f(s,\varphi _{1}(s))-f(s,\varphi _{2}(s))\right)ds\right\|\\&\leq \int _{t_{0}}^{t}\left\|f\left(s,\varphi _{1}(s)\right)-f\left(s,\varphi _{2}(s)\right)\right\|ds\\&\leq L\int _{t_{0}}^{t}\left\|\varphi _{1}(s)-\varphi _{2}(s)\right\|ds&&{\text{since }}f{\text{ is Lipschitz-continuous}}\\&\leq L\int _{t_{0}}^{t}\left\|\varphi _{1}-\varphi _{2}\right\|_{\infty }\,ds\\&\leq La\left\|\varphi _{1}-\varphi _{2}\right\|_{\infty }\end{aligned}}}$

This is a contraction if ${\displaystyle a<{\tfrac {1}{L}}.}$

We have established that the Picard's operator is a contraction on the Banach spaces with the metric induced by the uniform norm. This allows us to apply the Banach fixed-point theorem to conclude that the operator has a unique fixed point. In particular, there is a unique function

${\displaystyle \varphi \in {\mathcal {C}}(I_{a}(t_{0}),B_{b}(y_{0}))}$

such that Γφ = φ. This function is the unique solution of the initial value problem, valid on the interval I_a where a satisfies the condition

${\displaystyle a<\min \left\{{\tfrac {b}{M}},{\tfrac {1}{L}}\right\}.}$

Optimization of the solution's interval

We wish to remove the dependence of the interval I_a on L. To this end, there is a corollary of the Banach fixed-point theorem: if an operator Tⁿ is a contraction for some n in N, then T has a unique fixed point. Before applying this theorem to the Picard operator, recall the following:

Proof. Induction on m. For the base of the induction (m = 1) we have already seen this, so suppose the inequality holds for m − 1, then we have:

By taking a supremum over ${\displaystyle t\in [t_{0}-\alpha ,t_{0}+\alpha ]}$ we see that ${\displaystyle \left\|\Gamma ^{m}\varphi _{1}-\Gamma ^{m}\varphi _{2}\right\|\leq {\frac {L^{m}\alpha ^{m}}{m!}}\left\|\varphi _{1}-\varphi _{2}\right\|}$.

This inequality assures that for some large m,

and hence Γ^m will be a contraction. So by the previous corollary Γ will have a unique fixed point. Finally, we have been able to optimize the interval of the solution by taking α = min{a, b/M}.

In the end, this result shows the interval of definition of the solution does not depend on the Lipschitz constant of the field, but only on the interval of definition of the field and its maximum absolute value.

Other existence theorems

The Picard–Lindelöf theorem shows that the solution exists and that it is unique. The

${\displaystyle y(t)=0,\qquad y(t)=\pm \left({\tfrac {2}{3}}t\right)^{\frac {3}{2}}}$.

Even more general is Carathéodory's existence theorem, which proves existence (in a more general sense) under weaker conditions on  f . Although these conditions are only sufficient, there also exist necessary and sufficient conditions for the solution of an initial value problem to be unique, such as Okamura's theorem.^[5]

Global existence of solution

The Picard–Lindelöf theorem ensures that solutions to initial value problems exist uniquely within a local interval ${\displaystyle [t_{0}-\varepsilon ,t_{0}+\varepsilon ]}$, possibly dependent on each solution. The behavior of solutions beyond this local interval can vary depending on the properties of  f  and the domain over which  f  is defined. When  f  is globally Lipschitz, solutions can be extended and are defined over the entire real line. However, if  f  is only smooth, solutions may not be defined for certain values of t. For instance, the equation dy/dt = y² with initial condition y(0) = 1 has the solution y(t) = 1/(1-t), which is not defined at t = 1. Nevertheless, if  f  is a smooth

See also

• Complete vector fields
• Frobenius theorem (differential topology)
• Integrability conditions for differential systems
• Newton's method
• Euler method
• Trapezoidal rule

Notes

References

• ISBN 9780070992566
  .
• Lindelöf, E. (1894). "Sur l'application de la méthode des approximations successives aux équations différentielles ordinaires du premier ordre". Comptes rendus hebdomadaires des séances de l'Académie des sciences. 118: 454–7. (In that article Lindelöf discusses a generalization of an earlier approach by Picard.)
• Zbl 1263.34002
  .

External links

• "Cauchy-Lipschitz theorem". Encyclopedia of Mathematics.
• Fixed Points and the Picard Algorithm, recovered from http://www.krellinst.org/UCES/archive/classes/CNA/dir2.6/uces2.6.html.
• Grant, Christopher (1999). "Lecture 4: Picard-Lindelöf Theorem" (PDF). Math 634: Theory of Ordinary Differential Equations. Department of Mathematics, Brigham Young University.