ODE Lecture 14
ODE Lecture 14
ODE Lecture 14
March-June 2023
Lecture-14
(Sections 7.1-7.3 of Differential Equations by S. L. Ross, 3rd Edition)
There are many physical problems in which more than one element or
object interact with each other. To mathematically model such
physical problems, we need a system of two or more differential
equations. For example,
• Biology: Predator-Prey Model
• Mechanics: Motion of certain Spring-Mass System (Two masses
moving on frictionless surface)
• Electrical Circuit: Parallel LRC Circuit
• Mixture Problems: Connected Mixing Tanks
Consider the population problem of two different species, say, Foxes and Rabbits
interact living in the same ecosystem. A simple model supposes that rabbits eat
only grass and foxes eat only rabbits. In other words, fox is a predator and rabbit
is a prey.
When rabbits are present, the number of interactions/ encounters between these
two species per unit time is proportional to their populations F and R , that is,
proportional to the product F R .
Thus when rabbits are present, there is a supply of food to foxes, so foxes are
added to the system at the rate of bF R, b > 0 . Therefore
dF
= −aF + bF R, a > 0, b > 0.
dt
But when foxes are present, the population of the rabbits is decreased by the rate
at which the rabbits are eaten during their encounters with the foxes, that is,
−dF R, d > 0. Therefore
dR
= cR − dF R
dt
where c > 0 and d > 0 .
Thus, this model leads to the following nonlinear system of ODEs.
Lotka-Volterra Predator-Prey Model: F -Foxes (Predator) and R -Rabbits (Prey)
dF
= −aF + bF R, a, b > 0
dt
dR
= −cR − dF R, c, d > 0.
dt
(March-June 2023) MA 102-ODE Lecture-14 5 / 49
Competition Model
Suppose two different species of animals occupy the same ecosystem, not as
predator and prey, but rather as competitors for the same resources (such as food
and living space) in the system. In absence of the other, let us assume that the
rate at which each population grows is given by
dx dy
= ax and = cy where a, c > 0.
dt dt
Since the two species compete, another assumption might be that each of these
rates is diminished simply by the influence or existence of other population.
dx
= ax − by,
dt
dy
= −dx + cy,
dt
dx1
= a11 (t)x1 + · · · + a1n (t)xn + b1 (t)
dt
dx2
= a21 (t)x1 + · · · + a2n (t)xn + b2 (t)
dt
... ...............
dxn
= an1 (t)x1 + · · · + ann (t)xn + bn (t)
dt
It can be written in the vector differential equation (VDE) form as
dx dX dX~
= Ax + f or =AX + F or ~X
=A ~ + F~ ,
dt dt dt
where x(t) = [x1 (t), . . . , xn (t)]T , f (t) = [b1 (t), . . . , bn (t)]T , and A(t) = [aij (t)]
is a n × n matrix.
Example 2: Here x(t) = (x(t), y(t)). That is, t is independent variable and x, y
are dependent variables.
dx
= t2 x + (t + 1) y + 3 t3
dt
dy
= 2 x + t3 y − e t
dt
dx dX dX~
= Ax or =AX or =A~X ~ ,
dt dt dt
a11 (t) · · · a1n (t) x1 (t)
.. .. .. .
A = A(t) = . . . & x = x(t) = .. .
an1 (t) · · · ann (t) xn (t)
(March-June 2023) MA 102-ODE Lecture-14 12 / 49
Example
Consider the homogeneous linear VDE x0 (t) = A(t) x(t) where
7 −1 6 x1 (t)
A(t) = −10 4 −12 and x(t) = x2 (t) .
−2 1 −1 x3 (t)
is to find a vector function x(t) ∈ C 1 that satisfies the system (1) on an interval I
and the initial conditions x(t0 ) = x0 = (x0,1 , . . . , x0,n )T , where t0 ∈ I and
x0 ∈ Rn .
Step 1:
The solution to the IVP (2) is equivalent to the solution of the vector integral
equation Z t
x(t) = x0 + A(s) x(s) ds (3)
t0
x0 (t) = x0 ,
Z t
xn+1 (t) = x0 (t) + A(s) xn (s) ds t ≥ t0 , t ∈ I (4)
t0
Since xn (t) → x(t) uniformly on I, we can take the limit under the integral sign
of (4) to get
Z t
x(t) = x0 (t) + A(s) x(s) ds
t0
which proves that x(t) is the solution of the integral equation (3) and hence x(t)
is a solution to the IVP (2).
Step 3: Uniqueness of the Solution
If x(t) and y(t) are the solutions of (2) then
Rt Rt
x(t)−y(t) = t0 A(s)(x(s)−y(s)) ds =⇒ |x(t)−y(t)| ≤ M t0 |x(s)−y(s)| ds.
Thus for any given >R0, we get from the above inequality,
t
|x(t) − y(t)| < + M t0 |x(s) − y(s)| ds. By Grownwall integral inequality, one
can get
|x(t) − y(t)| < exp(M (t − t0 )) for t ≥ t0 .
Since the above inequality is true for each > 0, we can get |x(t) − y(t)| = 0 and
hence x(t) = y(t) for t ≥ t0 .
Proof.
By the Theorem on Existence and Uniqueness of solution to the IVP involving
homogeneous linear VDE, the solution Φ(t) is a unique solution to the IVP
satisfying Φ(t0 ) = 0. Then as given in the proof of that theorem, one can
compute the Picard’s iterates (successive approximation) xn (t) and observe that
xn (t) = 0 for all t ∈ [a, b] and for each n ∈ N. Therefore, {xn } → Φ and
Φ(t) = 0 for all t ∈ [a, b].
x0 = A x
Proof.
Step 1: Let V = {x(t) : x0 (t) = A(t)x(t) for all t ∈ I}.
To show: V is a vector space over the real field.
Let x1 and x2 be any two elements in V .
To show: c1 x1 + c2 x2 ∈ V for any scalars c1 and c2 .
d dx1 dx2
(c1 x1 + c2 x2 ) = c1 + c2 = c1 Ax1 + c2 Ax2 = A(c1 x1 + c2 x2 ) .
dt dt dt
Therefore, c1 x1 + c2 x2 ∈ V . Hence V is a vector space.
c1 x1 + c2 x2 + · · · + cn xn = 0 .
To show: x1 , . . ., xn spans V .
Let Φ be the solution to the IVP x0 = Ax with x(t0 ) = x0 on I, where x0 ∈ Rn .
That is, Φ ∈ V .
Since x0 ∈ Rn , there exist unique scalars c1 , . . ., cn such that
c1 e1 + c2 e2 + · · · + cn en = x0 .
X n
Since xi (t0 ) = ei , we get x0 = ci xi (t0 ).
Pn i=1
Thus, the function Ψ(t) = i=1 ci xi (t) for t ∈ I is a solution to the IVP
x0 = Ax satisfying the initial condition Ψ(t0 ) = x0 .
Since the solution
Pn to the IVP is unique, we can conclude that
Φ ≡ Ψ = i=1 ci xi on I.
Therefore, every solution Φ in V is a linear combination of x1 , . . ., xn .
This completes the proof of the theorem.
is to find a vector function x(t) ∈ C 1 that satisfies the system (4) on an interval I
and the initial conditions x(t0 ) = x0 = (x1,0 , . . . , xn,0 )T , where t0 ∈ I and
x0 ∈ Rn .
where aii (t), i = 1, ..., n are the main diagonal elements of A(t).
c1 x1 (t) + · · · + cn xn (t)
Therefore
x(t) = Φ(t) c for t ∈ [a, b]
is the general solution of x0 = A(t) x.
is the unique solution to the IVP. Here Φ−1 (t0 ) denotes the inverse matrix of
Φ(t0 ).
(March-June 2023) MA 102-ODE Lecture-14 31 / 49
Example: The set {x1 , x2 , x3 }, where
2t
e −e−t 0
x1 = e2t , x2 = 0 , x3 = e−t ,
e2t e−t −e−t
0
is a fundamental
solution
set for the system x (t) = A(t)x(t) on R,
0 1 1
where A = 1 0 1 .
1 1 0
Note that Axi (t) = x0i (t), i = 1, 2, 3. Further,
e2t −e−t 0
2t
W (0) = e 0 e−t = −3 6= 0.
2t −t
e e −e−t t=0
Thus, the GS is
2t
e −e−t 0 c1
x(t) = c1 e2t +c2 0 +c3 e−t = Φ(t) c2 = Φ(t)c
e2t e−t −e−t c3
c1
where c = c2 .
c3