Processing math: 100%
Skip to main content

Section 1.6 Existence and Uniqueness of Solutions

If xβ€²=f(t,x) and x(t0)=x0 is a linear differential equation, we have already shown that a solution exists and is unique. We will now take up the question of existence and uniqueness of solutions for all first-order differential equations. The existence and uniqueness of solutions will prove to be very importantβ€”even when we consider applications of differential equations.

Subsection 1.6.1 The Existence and Uniqueness Theorem

The following theorem tells us that solutions to first-order differential equations exist and are unique under certain reasonable conditions.

Let us examine some consequences of the existence and uniqueness of solutions.

Example 1.6.2.

Consider the initial value problem yβ€²=y1/3 with y(0)=0 and tβ‰₯0. Separating the variables,

yβˆ’1/3yβ€²=dt.

Thus,

32y2/3=t+C

or

y=(23(t+C))3/2.

If C=0, the initial condition is satisfied and

y=(23t)3/2

is a solution for tβ‰₯0. However, we can find two additional solutions for tβ‰₯0:

y=βˆ’(23t)3/2,y≑0.

This is especially troubling if we are looking for equilibrium solutions. Although yβ€²=y1/3 is an autonomous differential equation, there is no equilibrium solution at y=0. The problem is that

βˆ‚βˆ‚yy1/3=13yβˆ’2/3

is not defined at y=0.

Example 1.6.3.

Suppose that yβ€²=y2 with y(0)=1. Since f(t,y)=y2 and βˆ‚f/βˆ‚y=2y are continuous everywhere, a unique solution exists near t=0. Separating the variables,

1y2dy=dt,

we see that

y=βˆ’1t+C

or

y=11βˆ’t.

Therefore, a solution also exists on (βˆ’βˆž,1) if y(0)=βˆ’1. In the case that y(0)=βˆ’1, the solution is

y=βˆ’1t+1,

and a solution exists on (βˆ’1,∞). Solutions are only guaranteed to exist on an open interval containing the initial value and are very dependent on the initial condition.

Remark 1.6.4. Solutions Curves Cannot Cross.

The Existence and Uniqueness Theorem tells us that the integral curves of any differential equation satisfying the appropriate hypothesis, cannot cross. If the curves did cross, we could take the point of intersection as the initial value for the differential equation. In this case, we would no longer guaranteed unique solutions to a differential equation.

Subsection 1.6.2 Picard Iteration

It was Emile Picard (1856–1941) who developed the method of successive approximations to show the existence of solutions of ordinary differential equations. He proved that it is possible to construct a sequence of functions that converges to a solution of the differential equation. One of the first steps towards understanding Picard iteration is to realize that an initial value problem can be recast in terms of an integral equation.

Suppose that \(u = u(t)\) is a solution to

\begin{align*} x' & = f(t, x)\\ x(t_0) & = x_0, \end{align*}

on some interval \(I\) containing \(t_0\text{.}\) Since \(u\) is continuous on \(I\) and \(f\) is continuous on \(R\text{,}\) the function \(F(t) = f(t, u(t))\) is also continuous on \(I\text{.}\) Integrating both sides of \(u'(t) = f(t, u(t))\) and applying the Fundamental Theorem of Calculus, we obtain

\begin{equation*} u(t) - u(t_0) = \int_{t_0}^t u'(s) \, ds = \int_{t_0}^t f(s, u(s)) \, ds \end{equation*}

Since \(u(t_0) = x_0\text{,}\) the function \(u\) is a solution of the integral equation.

Conversely, assume that

\begin{equation*} u(t) = x_0 + \int_{t_0}^t f(s, u(s)) \, ds. \end{equation*}

If we differentiate both sides of this equation, we obtain \(u'(t) = f(t, u(t))\text{.}\) Since

\begin{equation*} u(t_0) = x_0 + \int_{t_0}^{t_0} f(s, u(s)) \, ds = x_0, \end{equation*}

the initial condition is fulfilled.

To show the existence of a solution to the initial value problem

xβ€²=f(t,x)x(t0)=x0,

we will construct a sequence of functions, {un(t)}, that will converge to a function u(t) that is a solution to the integral equation

x(t)=x0+∫tt0f(s,x(s))ds.

We define the first function of the sequence using the initial condition,

u0(t)=x0.

We derive the next function in our sequence using the right-hand side of the integral equation,

u1(t)=x0+∫tt0f(s,u0(s))ds.

Subsequent terms in the sequence can be defined recursively,

un+1=x0+∫tt0f(s,un(s))ds.

Our goal is to show that un(t)β†’u(t) as nβ†’βˆž. Furthermore, we need to show that u is the continuous, unique solution to our initial value problem. We will leave the proof of Picard's Theorem to a series of exercises, but let us see how this works by developing an example.

Example 1.6.6.

Consider the exponential growth equation,

dxdt=kxx(0)=1.

We already know that the solution is x(t)=ekt. We define the first few terms of our sequence {un(t)} as follows:

u0(t)=1,u1(t)=1+∫t0ku0(s)ds=1+∫t0kds=1+kt,u2(t)=1+∫t0ku1(s)ds=1+∫t0k(1+ks)ds=1+kt+(kt)22.

The next term in the sequence is

u3(t)=1+kt+(kt)22+(kt)32β‹…3,

and the nth term is

un(t)=1+1+∫t0kunβˆ’1(s)ds=1+∫t0k(1+ks(kt)22!+(kt)33!+β‹―+(kt)nβˆ’1(nβˆ’1)!)ds=1+kt+(kt)22!+(kt)33!+β‹―+(kt)nn!.

However, this is just the nth partial sum for the power series for u(t)=ekt, which is what we expected.

Subsection 1.6.3 Important Lessons

  • Existence and uniqueness of solutions of differential equations has important implications. Let xβ€²=f(t,x) have the initial condition x(t0)=x0. If f and βˆ‚f/βˆ‚x are continuous functions on the rectangle

    R={(t,x):0≀|tβˆ’t0|<a,0≀|xβˆ’x0|<b},

    there exists a unique solution u=u(t) for xβ€²=f(t,x) and x(t0)=x0 on some interval |tβˆ’t0|<h contained in the interval |tβˆ’t0|<a. In particular,

    • Solutions are only guaranteed to exist locally.
    • Uniqueness is especially important when it comes to finding equilibrium solutions.
    • Uniqueness of solutions tells us that the integral curves for a differential equation cannot cross.
  • The function u=u(t) is a solution to the initial value problem
    xβ€²=f(t,x)x(t0)=x0,
    if and only if u is a solution to the integral equation
    x(t)=x0+∫tt0f(s,x(s))ds.
  • Existence and uniqueness of solutions is proved by Picard iteration. This is of particular interest since the proof actually tells us how to construct a sequence of functions that converge to our solution.

Reading Questions 1.6.4 Reading Questions

2.

What important rule from differential calculus do we use when solving a first-order differential equation?

Exercises 1.6.5 Exercises

1.

Which of the following initial value problems are guaranteed to have a unique solution by the Existence and Uniqueness Theorem (Theorem 1.6.1)? In each case, justify your conclusion.

  1. yβ€²=4+y3, y(0)=1
  2. yβ€²=√y, y(1)=0
  3. yβ€²=√y, y(1)=1
  4. xβ€²=txβˆ’2, x(0)=2
  5. xβ€²=txβˆ’2, x(2)=0
  6. yβ€²=xtany, y(0)=0
  7. yβ€²=1ty+2t, y(0)=1
Hint
  1. There exists a unique solution to \(y' = 4 + y^3\text{,}\) \(y(0) = 1\text{,}\) since \(f(t, y) = 4 + y^3\) and \(\partial f(t, y)/\partial y = 3y^2\) are continuous at the point \((0, 1)\text{.}\)
  2. The Existence and Uniqueness Theorem does not apply to \(y' = \sqrt{y}\text{,}\) \(y(1) = 0\text{,}\) since \(f(t, y) = \sqrt{y}\) is not continuous at \((1, 0)\text{.}\)
  3. There exists a unique solution to \(y' = \sqrt{y}\text{,}\) \(y(1) = 1\text{,}\) since \(f(t, y) = \sqrt{y}\) and \(\partial f(t, y)/\partial y = 1/(2 \sqrt{y}\, )\) are both continuous at the point \((1, 1)\text{.}\)
  4. The Existence and Uniqueness Theorem does not apply to \(x' = t/(x - 2)\text{,}\) \(x(0) = 2\text{,}\) since \(f(t, x) = t/(x - 2)\) is not continuous at \((0, 2)\text{.}\)
  5. There exists a unique solution to \(x' = t/(x - 2)\text{,}\) \(x(2) = 0\text{,}\) since \(f(t, x) = t/(x - 2)\) and \(\partial f(t, x)/\partial x = - t/(x-2)^2\) are both continuous at the point \((2, 0)\text{.}\)
  6. There exists a unique solution to \(y' = x \tan y\text{,}\) \(y(0) = 0\text{,}\) since \(f(x, y) = x \tan y\) and \(\partial f(x, y)/\partial y = x \sec^2 y\) are both continuous at the point \((0, 0)\text{.}\)
  7. The Existence and Uniqueness Theorem does not apply to \(y' = y/t + 2t\text{,}\) \(y(0) = 1\text{,}\) since \(f(t, y) = y/t + 2t\) is not continuous at \((0, 1)\text{.}\)
2.

Find an explicit solution to the initial value problem

yβ€²=1(tβˆ’1)(y+1)y(0)=1.

Use your solution to determine the interval of existence.

3.

Consider the initial value problem

yβ€²=3y2/3y(0)=0.
  1. Show that the constant function, y(t)≑0, is a solution to the initial value problem.
  2. Show that
    y(t)={0,t≀t0(tβˆ’t0)3,t>t0
    is a solution for the initial value problem, where t0 is any real number. Hence, there exists an infinite number of solutions to the initial value problem. [Hint: Make sure that the derivative of y(t) exists at t=t0. ]
  3. Explain why this example does not contradict the Existence and Uniqueness Theorem.
4.

Let Ο•n(x)=xn for 0≀x≀1 and show that

limnβ†’βˆžΟ•n(x)={0,0≀x<11,x=1.

This is an example of a sequence of continuous functions that does not converge to a continuous function, which helps explain the need for uniform continuity in the proof of the Existence and Uniqueness Theorem.

5.

Consider the initial value problem

yβ€²=2ty+ty(0)=1.
  1. Use the fact that yβ€²=2ty+t is a first-order linear differential equation to find a solution to the initial value problem.
  2. Let Ο•0(t)=1 and use Picard iteration to find Ο•n(t).
  3. Show that the sequence {Ο•n(t)} converges to the exact solution that you found in part (a) as nβ†’βˆž.
Proof of the Existence and Uniqueness Theorem

In Exercise Group 1.6.5.6–14, prove the Existence and Uniqueness Theorem for first-order differential equations.

6.

Use the Fundamental Theorem of Calculus to show that the function u=u(t) is a solution to the initial value problem

xβ€²=f(t,x)x(t0)=x0,

if and only if u is a solution to the integral equation

x(t)=x0+∫tt0f(s,x(s))ds.
7.

If βˆ‚f/βˆ‚x is continuous on the rectangle

R={(t,x):0≀|tβˆ’t0|<a,0≀|xβˆ’x0|<b},

prove that there exists a K>0 such that

|f(t,x1)βˆ’f(t,x2)|≀K|x1βˆ’x2|

for all (t,x1) and (t,x2) in R.

8.

Define the sequence {un} by

u0(t)=x0,un+1=x0+∫tt0f(s,un(s))ds,n=1,2,….

Use the result of the previous exercise to show that

|f(t,un(t))βˆ’f(t,unβˆ’1(t))|≀K|un(t)βˆ’unβˆ’1(t)|.
9.

Show that there exists an M>0 such that

|u1(t)βˆ’x0|≀M|tβˆ’t0|.
10.

Show that

|u2(t)βˆ’u1(t)|≀KM|tβˆ’t0|22.
11.

Use mathematical induction to show that

|un(t)βˆ’unβˆ’1(t)|≀Knβˆ’1M|tβˆ’t0|nn!.
12.

Since

un(t)=u1(t)+[u2(t)βˆ’u1(t)]+β‹―+[un(t)βˆ’unβˆ’1(t)],

we can view un(t) as a partial sum for the series

u0(t)+βˆžβˆ‘n=1[un(t)βˆ’unβˆ’1(t)].

If we can show that this series converges absolutely, then our sequence will converge to a function u(t). Show that

βˆžβˆ‘n=1|un(t)βˆ’unβˆ’1(t)|≀MKβˆžβˆ‘n=1(K|tβˆ’t0|)nn!≀MK(eK|h|βˆ’1),

where h is the maximum distance between (t0,x0) and the boundary of the rectangle R. Since |un(t)βˆ’unβˆ’1(t)|β†’0, we know that un(t) converges to a continuous function u(t) that solves our equation. 9 We must a theorem from advanced calculus here to ensure uniform continuity. Any sequence of functions that converges uniformly, must converge to a continuous function.

13.

To show uniqueness, assume that u(t) and v(t) are both solutions to

x(t)=x0+∫tt0f(s,x(s))ds.

Show that

|u(t)βˆ’v(t)|≀K∫tt0|u(s)βˆ’v(s)|ds.
14.
  1. Define 10 A similar argument will work for \(t \leq t_0\text{.}\)
    Ο•(t)=∫tt0|u(s)βˆ’v(s)|ds,
    then Ο•(t0)=0 and Ο•(t)β‰₯0 for tβ‰₯t0. Show that
    Ο•β€²(t)=|u(t)βˆ’v(t)|.
  2. Since
    |u(t)βˆ’v(t)|βˆ’K∫tt0|u(s)βˆ’v(s)|ds≀0,
    we know that
    Ο•β€²(t)βˆ’KΟ•(t)≀0.
    Use this fact to show that
    ddt[eβˆ’KtΟ•(t)]≀0.
    Conclude that
    Ο•(t)=∫tt0|u(s)βˆ’v(s)|ds=0
    for tβ‰₯t0 or for all tβ‰₯t0 and u(t)=v(t).

Subsection 1.6.6 Projectβ€”Picard Iteration