ConceptComplete

The Picard-Lindelof Framework

The existence and uniqueness theory for ordinary differential equations establishes when initial value problems have solutions, how many solutions they have, and how solutions depend on initial data.


Initial Value Problems

Definition8.1Initial value problem (IVP)

An initial value problem for a first-order ODE system consists of a differential equation and an initial condition:

x=F(t,x),x(t0)=x0,\mathbf{x}' = \mathbf{F}(t, \mathbf{x}), \quad \mathbf{x}(t_0) = \mathbf{x}_0,

where F:ΩRn\mathbf{F}: \Omega \to \mathbb{R}^n is defined on an open set ΩR×Rn\Omega \subset \mathbb{R} \times \mathbb{R}^n and (t0,x0)Ω(t_0, \mathbf{x}_0) \in \Omega.

A solution is a differentiable function x:IRn\mathbf{x}: I \to \mathbb{R}^n defined on an interval II containing t0t_0 such that x(t)=F(t,x(t))\mathbf{x}'(t) = \mathbf{F}(t, \mathbf{x}(t)) for all tIt \in I and x(t0)=x0\mathbf{x}(t_0) = \mathbf{x}_0.

Definition8.2Lipschitz condition

A function F(t,x)\mathbf{F}(t, \mathbf{x}) satisfies a Lipschitz condition in x\mathbf{x} on a set SR×RnS \subset \mathbb{R} \times \mathbb{R}^n if there exists a constant L>0L > 0 (the Lipschitz constant) such that:

F(t,x)F(t,y)Lxy\|\mathbf{F}(t, \mathbf{x}) - \mathbf{F}(t, \mathbf{y})\| \leq L\|\mathbf{x} - \mathbf{y}\|

for all (t,x),(t,y)S(t, \mathbf{x}), (t, \mathbf{y}) \in S. If F\mathbf{F} has continuous partial derivatives Fi/xj\partial F_i / \partial x_j, then F\mathbf{F} is locally Lipschitz with L=supF/xL = \sup \|\partial \mathbf{F}/\partial \mathbf{x}\| on compact subsets.


The Picard-Lindelof Theorem

Theorem8.1Picard-Lindelof existence and uniqueness

Let F:ΩRn\mathbf{F}: \Omega \to \mathbb{R}^n be continuous on an open set ΩR×Rn\Omega \subset \mathbb{R} \times \mathbb{R}^n and Lipschitz in x\mathbf{x} on Ω\Omega. Then for every (t0,x0)Ω(t_0, \mathbf{x}_0) \in \Omega, the IVP x=F(t,x)\mathbf{x}' = \mathbf{F}(t, \mathbf{x}), x(t0)=x0\mathbf{x}(t_0) = \mathbf{x}_0 has a unique solution on some interval [t0δ,t0+δ][t_0 - \delta, t_0 + \delta].

More precisely, if R=[t0a,t0+a]×Bb(x0)ΩR = [t_0 - a, t_0 + a] \times \overline{B_b(\mathbf{x}_0)} \subset \Omega, M=supRFM = \sup_R \|\mathbf{F}\|, and δ=min(a,b/M)\delta = \min(a, b/M), then there exists a unique solution on [t0δ,t0+δ][t_0 - \delta, t_0 + \delta].

ExampleVerifying Lipschitz conditions

Consider x=x2x' = x^2, x(0)=1x(0) = 1. Here F(t,x)=x2F(t,x) = x^2, F/x=2x\partial F/\partial x = 2x. On any compact set xK|x| \leq K, the Lipschitz constant is L=2KL = 2K. So Picard-Lindelof guarantees local existence and uniqueness. The solution x(t)=1/(1t)x(t) = 1/(1-t) exists on (,1)(-\infty, 1) and blows up at t=1t = 1.

Compare with x=xx' = \sqrt{|x|}, x(0)=0x(0) = 0. Here F(t,x)=xF(t,x) = \sqrt{|x|} is not Lipschitz at x=0x = 0 (F=1/(2x)|F'| = 1/(2\sqrt{|x|}) \to \infty). Indeed, both x(t)=0x(t) = 0 and x(t)=t2/4x(t) = t^2/4 (for t0t \geq 0) are solutions, so uniqueness fails.


Integral Equation Formulation

RemarkEquivalence with integral equations

The IVP x=F(t,x)\mathbf{x}' = \mathbf{F}(t, \mathbf{x}), x(t0)=x0\mathbf{x}(t_0) = \mathbf{x}_0 is equivalent to the Volterra integral equation:

x(t)=x0+t0tF(s,x(s))ds.\mathbf{x}(t) = \mathbf{x}_0 + \int_{t_0}^{t} \mathbf{F}(s, \mathbf{x}(s))\,ds.

This reformulation is essential because:

  1. It incorporates the initial condition directly.
  2. The integral operator has better regularity properties than differentiation.
  3. It provides the foundation for the Picard iteration scheme.
ExampleIntegral equation for a linear ODE

For x=2txx' = -2tx, x(0)=1x(0) = 1, the integral equation is:

x(t)=1+0t(2sx(s))ds.x(t) = 1 + \int_0^t (-2s \cdot x(s))\,ds.

Starting with x0(t)=1x_0(t) = 1: x1(t)=1t2x_1(t) = 1 - t^2, x2(t)=1t2+t4/2x_2(t) = 1 - t^2 + t^4/2. The pattern gives x(t)=et2x(t) = e^{-t^2}, confirming the exact solution.