Initial Value Problem

Work in progress. This note is still being written and incomplete.

Also known as IVP. A first-order ODE together with the solution value at the left endpoint.

y′(t)=f(t,y),a≤t≤b,y(a)=αy'(t) = f(t, y), \quad a \leq t \leq b, \quad y(a) = \alpha

Numerical methods approximate yy at discrete points of [a,b][a, b].

Lipschitz Condition

f(t,y)f(t, y) satisfies a Lipschitz condition in yy on D⊆R2D \subseteq \mathbb{R}^2 iff a constant L>0L > 0 exists with

∣f(t,y1)−f(t,y2)∣≤L∣y1−y2∣|f(t, y_1) - f(t, y_2)| \leq L|y_1 - y_2|

for all (t,y1),(t,y2)∈D(t, y_1), (t, y_2) \in D. LL is a Lipschitz constant for ff.

Partial Derivative Test

If DD is convex and a constant L>0L > 0 exists with ∣∂f∂y(t,y)∣≤L\left| \dfrac{\partial f}{\partial y}(t, y) \right| \leq L for all (t,y)∈D(t, y) \in D, then ff satisfies a Lipschitz condition in yy on DD with Lipschitz constant LL.

Existence and Uniqueness

Let D=[a,b]×RD = [a, b] \times \mathbb{R}. If ff is continuous on DD and satisfies a Lipschitz condition in yy on DD, then the IVP has a unique solution y(t)y(t) for a≤t≤ba \leq t \leq b.

Stable Solution

Let δ(t)\delta(t) be a continuous perturbation of the equation and δ0\delta_0 a perturbation of the initial value. The perturbed problem is

z′(t)=f(t,z)+δ(t),a≤t≤b,z(a)=α+δ0z'(t) = f(t, z) + \delta(t), \quad a \leq t \leq b, \quad z(a) = \alpha + \delta_0

The solution y(t)y(t) is stable iff constants ε0>0\varepsilon_0 > 0 and k>0k > 0 exist such that for every ε∈(0,ε0)\varepsilon \in (0, \varepsilon_0):

  • whenever ∣δ(t)∣<ε|\delta(t)| < \varepsilon for all t∈[a,b]t \in [a, b] and ∣δ0∣<ε|\delta_0| < \varepsilon
  • the perturbed problem has a unique solution z(t)z(t)
  • ∣z(t)−y(t)∣<kε|z(t) - y(t)| < k\varepsilon for all t∈[a,b]t \in [a, b]

Here:

t∈[a,b]t \in [a, b]
Here:
kk
stability constant, independent of ε\varepsilon

Well-Posed Problem

The IVP is well-posed iff:

  • a unique solution y(t)y(t) exists
  • the solution is stable

Well-Posedness from Lipschitz

Let D=[a,b]×RD = [a, b] \times \mathbb{R}. If ff is continuous and satisfies a Lipschitz condition in yy on DD, then the IVP is well-posed.

If the Lipschitz condition fails, check the definition directly.

Worked Examples

Example 1

Find a Lipschitz constant in yy for f(t,y)=t2−3yf(t,y) = t^2 - 3y on D=[0,2]×RD = [0,2] \times \mathbb{R}.

Solution Reveal Hide
∂f∂y(t,y)=−3\frac{\partial f}{\partial y}(t,y) = -3∣∂f∂y(t,y)∣=3\left| \frac{\partial f}{\partial y}(t,y) \right| = 3

The partial derivative test gives L=3L = 3.

Example 2

Show that the IVP y′=1+tsin⁡(ty)y' = 1 + t\sin(ty), 0≤t≤20 \leq t \leq 2, y(0)=0y(0) = 0 has a unique solution.

Solution Reveal Hide
∣∂f∂y(t,y)∣=∣t2cos⁡(ty)∣≤4\left| \dfrac{\partial f}{\partial y}(t, y) \right| = |t^2 \cos(ty)| \leq 4

on [0,2]×R[0, 2] \times \mathbb{R}. Since f(t,y)=1+tsin⁡(ty)f(t, y) = 1 + t\sin(ty) is continuous, ff is Lipschitz in yy with L=4L = 4. The IVP has a unique solution on [0,2][0, 2].

Example 3

Use the Lipschitz condition to decide whether f(t,y)=sin⁡(t)+4yf(t,y) = \sin(t) + 4y has a unique solution for every initial value on [0,1][0,1].

Solution Reveal Hide

ff is continuous on [0,1]×R[0,1] \times \mathbb{R}.

∂f∂y(t,y)=4\frac{\partial f}{\partial y}(t,y) = 4∣∂f∂y(t,y)∣=4\left| \frac{\partial f}{\partial y}(t,y) \right| = 4

ff is Lipschitz in yy with L=4L = 4. Every IVP with this right-hand side has a unique solution on [0,1][0,1].

Example 4

Show that f(t,y)=∣y∣f(t,y) = \sqrt{|y|} does not satisfy a Lipschitz condition in yy on D=[0,1]×[0,1]D = [0,1] \times [0,1].

Solution Reveal Hide

Take y1=h2y_1 = h^2 and y2=0y_2 = 0, where 0<h≤10 < h \leq 1. Then

∣f(t,y1)−f(t,y2)∣=∣h−0∣=h|f(t,y_1) - f(t,y_2)| = |h - 0| = h∣y1−y2∣=h2|y_1-y_2| = h^2

If a Lipschitz constant LL existed, then h≤Lh2h \leq Lh^2. Hence 1≤Lh1 \leq Lh.

This fails when h<1Lh < \dfrac{1}{L}. No Lipschitz constant exists on DD.

Example 5

Show the IVP y′(t)=y−t2+1y'(t) = y - t^2 + 1, 0≤t≤20 \leq t \leq 2, y(0)=0.5y(0) = 0.5 is well-posed.

Solution Reveal Hide

f(t,y)=y−t2+1f(t, y) = y - t^2 + 1 is continuous on D=[0,2]×RD = [0, 2] \times \mathbb{R}, and ∣∂f∂y∣=1\left| \dfrac{\partial f}{\partial y} \right| = 1. So ff is Lipschitz in yy with L=1L = 1, and the IVP is well-posed.

To see the stability bound directly, take a constant equation perturbation δ\delta and an initial perturbation δ0\delta_0

z′(t)=z−t2+1+δ,z(0)=0.5+δ0z'(t) = z - t^2 + 1 + \delta, \quad z(0) = 0.5 + \delta_0

Both problems have closed-form solutions

y(t)=(t+1)2−0.5et,z(t)=(t+1)2−δ+(δ0+δ−0.5) ety(t) = (t+1)^2 - 0.5 e^t, \qquad z(t) = (t+1)^2 - \delta + (\delta_0 + \delta - 0.5)\, e^t∣z(t)−y(t)∣=∣(δ0+δ) et−δ∣≤(2e2+1) ε|z(t) - y(t)| = \left| (\delta_0 + \delta)\, e^t - \delta \right| \leq (2 e^2 + 1)\, \varepsilon

for ∣δ∣<ε|\delta| < \varepsilon and ∣δ0∣<ε|\delta_0| < \varepsilon, using et≤e2e^t \leq e^2 on [0,2][0, 2]. The bound holds with k=2e2+1k = 2 e^2 + 1 for every ε>0\varepsilon > 0.

Example 6

Show the IVP y′(t)=te3t−2yy'(t) = t e^{3t} - 2y, 0≤t≤10 \leq t \leq 1, y(0)=0y(0) = 0 is well-posed.

Solution Reveal Hide

f(t,y)=te3t−2yf(t, y) = t e^{3t} - 2y is continuous on D=[0,1]×RD = [0, 1] \times \mathbb{R}.

∣∂f∂y∣=∣−2∣=2\left| \frac{\partial f}{\partial y} \right| = |-2| = 2

ff is Lipschitz in yy with L=2L = 2, so the IVP is well-posed.

Wolfram Mathematica notebooks

A companion notebook for this note. Open it in Mathematica and evaluate the cells to reproduce the worked example.

The summary notebook defines every function used in this module, with documentation and no examples, as a starting point for your own problems.

Written by November 3, 2026 3 min read
Was this helpful?