Chapter 2
Chapter 2
2023
Index:
• Definitions and basic concepts.
1/30
Definitions and basic concepts
The initial condition of an ODE specifies the time and the phase space point
(t0 , x0 ) of a solution. Then ϕ(t0 ) = x0 .
ODE may involve higher derivatives and definitions of Eqs.(1-2) are easily
generalised. Explicit high order ODE with autonomous vector field can be
always reduced to a system of 1st order ODE.
2
4 ẍ = −kx, k>0
Vector Length
1
3
0 is equivalent to the system
y
−1 2
ẋ = y
−2
−3
1
ẏ = −kx.
−4 0
−4 −3 −2 −1 0
x
1 2 3 4 with vector field f (x, y) = (y, −kx).
Note that scaling time we can generically assume that k = 1. The integral
curves are circles centered at the origin x2 + y 2 = R2 . If at t0 = 0 the
position and the velocity of the particle are (x0 , ẋ0 ) = (x0 , y0 ) = (1, 0) the
solution is ϕ(t) = (cos t, sin t) (inner circle of the figure).
4/30
Separation of variables
Example:
Solve y ′ − 2xy = x with y(0) = 5.
11 x2
Solution: y(x) = 2 e − 12 .
5/30
Example: Population models in biology
2.5
ṫ = 1
2
1 ẋ = x.
1.5
0 1
with vector field f (t, x) = (1, x).
0 1 2 3 4 6/30
t
Example: Population models in biology
Verhulst (1836) improved Malthus’ model, taking into account finite available
food:
ẋ = x(a − bx), with a, b > 0.
By scaling t and x the general logistic equation ẋ = x(1 − x) is obtained.
Note that for x = 1, the population remains constant since ẋ = 0. Moreover,
since x > 0 then the population increases (ẋ > 0) if x < 1. Conversely, the
population decreases.
Considering x0 = x(0) (the initial population) and by separation of variables
the population at each time is
t
ke
1+ket if x0 < 1
x(t) = 1 if x0 = 1
−ket
if x0 > 1
1−ket
x0
where k = |1−x0 | .
7/30
Example: Population models in biology
f (t, x) = (1, x(1 − x)), (t, et /(1 + et )), (t, 1), (t, −5et /2(1 − 5et /2))
2 2.2
Vector Length
1.8
1
x
1.6
1.4
1.2
0 1
0 1 2 3
t
0.4
Phase portrait (x, ẋ) displaying the
fixed points (ẋ = f (t, x) = 0). The
ẋ
0.2
point x = 1 is stable and x = 0 is
0•
0 0.5
•1 1.5
unstable.
x 8/30
Example: Population models in biology
Exercice:
Consider now that part of the individuals are emigrating at a constant rate
c > 0. By scaling t and x the general migration equation (also called
harvesting)
ẋ = x(1 − x) − c
is obtained.
Determine the fixed points and plot the corresponding vector field with the
integral curves as well as the phase portrait. Discuss the results as function
of the migration parameter c.
Perform the same analysis assuming now the migration rate is relative to
the population:
ẋ = x(1 − x) − cx
9/30
Homogeneous differential equations
dy y
y′ = =f ,
dx x
P (1,y/x)
with f (y/x) = − Q(1,y/x) , and the solution is given by separation of variables:
du dx
=−
u+ P (1,u) x
Q(1,u)
10/30
Example: Aircraft with side-wind
Q •S
11/30
Example: Aircraft with side-wind
x − yy′ w
−−→ = v .
|P O|
w −−→ p
Considering α = v and with |P O| = x2 + y 2 an homogeneous ODE is
obtained:
y p p
2 + y 2 =⇒ ydx + α x2 + y 2 − x dy = 0
x− = α x
y′
12/30
Example: Aircraft with side-wind
dy du
=− √
y α 1 + u2
√
which by integration gives (cy)−α = u + 1 + u2 .
Finally, since
p p p
u + 1 + u2 u − 1 + u2 = −1 =⇒ u − 1 + u2 = −(cy)α ,
13/30
Linear ODE
Given a(x) and b(x) an ODE of the form y ′ + a(x)y = b(x) is called linear.
The solution with initial condition y0 = y(x0 ) is:
Z x
−A(x) A(t)
y(x) = e y0 + e b(t)dt , (3)
x0
where Z x
A(x) = a(s)ds. (4)
x0
Example: Solve 2x2 y ′ + 4xy = 7x4 with y0 = y(x0 ). This ODE is equivalent
to y ′ + x2 y = 72 x2 . We first find A(x) = 2 ln xx0 to solve the integral
Z x Z x
2 ln xt 7 2 7
A(t)
e b(t)dt = e 0 t dt = (x5 − x50 ).
x0 x0 2 10x20
The solution is then
7 7 3 c 7 5
− 2 ln xx
y(x) = e 0 y0 + (x5 − x50 ) = x + 2, c = y0 x20 − x
10x20 10 x 10 0
14/30
Example: Chemical dilutions
15/30
Example: Chemical dilutions
Using Eqs. (3-4) we find the solute quantity at each instant t assuming
ve ̸= vs . If ve = vs (V (t) = V0 , ∀t) the ODE for the solute quantity can be
solved by separation of variables. If s0 is the initial (t = 0) quantity of the
solute, the general solution is
vs
V0 ve −vs (s − C V ) + C V (t) if v ̸= v
V (t) 0 e 0 e e s
s(t) =
(s − C V ) e −v V0
et
+C V if v = v
0 e 0 e 0 e s
• If ve > vs and assuming a tank with very large volume the solute
quantity s(t) tends to Ce V (t) for large t, i.e, C(t) → Ce .
V0
• If ve < vs then the tank will be empty after t = vs −ve minutes.
• If ve = vs the solute quantity s(t) tends to Ce V0 for large t, i.e,
C(t) → Ce .
16/30
Exact differential equations
∂P ∂Q
The ODE P (x, y)dx + Q(x, y)dy = 0 is exact if and only if ∂y = ∂x .
origin.
17/30
Example: Building a radio telescope
18/30
Example: Building a radio telescope
Solving Eq. (7) for y ′ (2nd degree equation) the ODE is rewritten as
p
′ −2x ± 4x2 + 4y 2 p
y = =⇒ ydy + xdx ∓ x2 + y 2 dx = 0 =⇒
2y
!
x y
=⇒ ∓ p + 1 dx ∓ p dy = 0,
x2 + y 2 x2 + y 2
which is an exact ODE with
!
x y
P (x, y) = ∓p +1 and Q(x, y) = ∓ p ,
x2 + y2 x2 + y2
verifying
∂P ∂Q 3
= = ±xy(x2 + y 2 )− 2
∂y ∂x
19/30
Example: Building a radio telescope
which is equivalent to
y2 k
x= − .
2k 2
This family of curves correspond to parabolas with their focus in the origin.
20/30
Bibliography recommended
21/30
Qualitative Theory of ODE. Basics
ϕ(T, x∗ , p) = x∗
22/30
Qualitative Theory of ODE. Basics
x − ϕ(T, x, p) = 0
g(x, p) = 0,
23/30
Qualitative Theory of ODE. Basics
asymptotically stable. 0
max max
0.3
The eigenvalues of
max max
M (T ) of largest 0.2
max max
-0.5
modulus have to be
max max
computed to asses the 0.1
stability of periodic (a) (b)
orbits. -1 0
0 0.5 1 0.7 0.8 0.9 1
26/30
ℜ(λ) ℜ(λ)
Qualitative Theory of ODE. Bifurcations
27/30
ments
Qualitative
10
(a)
Theory of ODE. Bifurcation Diagrams
1
P2
P1
variables of the solution).
Example: Bifurcation
k bran h 2
diagrams in a problem of
Rayleigh-Bénard rotating
1
10
-3
10
-2
10
-1
10
0 convection in a sphere.
f = Ra/Rac − 1
Ra 28/30
Qualitative Theory of ODE. Fold and Hopf Bifurcations
29/30
Bibliography recommended
30/30