System of First Order Equations

Anil Kumar
CC-205

System of 1st order ODE

 General form of the system of simultaneous first order
ordinary differential equations
x1  f1 (t , x1 , x2 , xn )
x2  f 2 (t , x1 , x2 , xn )

xn  f n (t , x1 , x2 , xn )
where each xk is a function of t.
 If each fk is a linear function of x1, x2, …, xn, then the system
of equations is said to be linear, otherwise it is nonlinear.
 Systems of higher order differential equations can similarly
be defined.
•2 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 1

 Example 1. System of 1st order ODE
 Consider the differential equation
u (t )  16 u (t )  192u (t )  0
 This second order equation can be converted into a system
of first order equations by letting
x1 = u and x2 = u'.
Thus
x1  x2
x2  16 x2  192 x1  0

or x1  x2
x2  16 x2  192 x1

•3 •Anil Kumar, Dept. of Mathematics

nth Order ODEs and Linear 1st Order Systems

 The method illustrated in previous example can be used to
transform an arbitrary nth order equation

y ( n )  f  t , y , y, y, , y ( n 1)  (1)

into a system of n first order equations, first by defining

x1  y, x2  y, x3  y,  , xn  y ( n 1) (2)
Then
x1  x2
x2  x3
 (3)
xn 1  xn
xn  f (t , x1 , x2 , xn )
•4 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 2

 System of 1st order ODE

 (1) and (3) are equivalent in the sense

if y(t) is a solution of (1), then the functions x1(t), x2(t),

…, xn(t) defined by (2) satisfy (3), and conversely, if

x1(t), x2(t), …, xn(t) satisfy (3), then y(t) is a solution of

(1).

•5 •Anil Kumar, Dept. of Mathematics

Solutions of First Order Systems

 A system of simultaneous first order ordinary differential
equations has the general form
x1  f1 (t , x1 , x2 , xn )
x2'  f 2 (t , x1 , x2 , xn )
 (1)
xn  f n (t , x1 , x2 , xn ).
It has a solution on I: a ≤ t ≤ b if there exists n functions
x1 (t ), x2 (t ), , xn (t )
which are differentiable on I and satisfy the system of
equations at all points t in I.
 Initial conditions may also be prescribed to give an IVP:
x1 (t0 )  a1 , x2 (t0 )  a2 , , xn (t0 )  an (2)
•6 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 3

 Theorem 1:
(Existence and uniqueness of the solution)

Suppose f1,…, fn and the partial derivatives

 f1 f f f
, ..., 1 , ... , n , ..., n
 x1 xn  x1 xn
are continuous in a region R of (t, x1, x2, …,xn) space. If
 t0 , a1 , a2 , , an  is an interior point of R. Then the system (1) has
a unique solution x1 (t ), x2 (t ), , xn (t ) that satisfies the initial
conditions (2).

•7 •Anil Kumar, Dept. of Mathematics

System of 1st order ODE

 Theorem 1 is a special case of the following theorem for the

equation (n)
y  f  t , y , y, y, , y
( n 1)
 (*)

 Theorem 2: Let the function f and the partial derivatives

f/y, f/y',…, f/y(n-1) be continuous in a region R of (t, y,
y',…, y(n-1)) space. If  t0 , a1 , a2 , , an  is an interior point of R,
then the equation (*) has a unique solution y(t) that satisfies
the initial conditions y(t0) = a1, y'(t0) = a2, …, y(n-1)(t0) = an.

•8 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 4

 Examples
 Transform the following equations into a system of first
order equations.
1
(i ) y  y   2 y  0,
2
1
(ii) y  y  2 y  3sin t ,
2
 1
(iii ) t 2 y  ty   t 2   y  0,
 4
(iv) y  y  x ( y) ,
2 2

(v ) y  p (t ) y  q (t ) y  g (t ), y (0)   0 and y(0)  1.

•9 •Anil Kumar, Dept. of Mathematics

Exercises

1. Transform the given system into a single equation of second

order

(i ) x1  x2 , x2  2 x1  2 x2 .
(ii ) x1  2 x2 , x2  2 x1.

•10 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 5

 Linear systems
 The general linear system of first order differential
equations is of the form
 dx
  a1 (t ) x  b1 (t ) y  f1 (t )
dt

 dy
 a2 (t ) x  b2 (t ) y  f 2 (t )
 dt

where a1(t), b1(t), a2(t), b2(t), f1(t), f2(t) are continuous on

a closed interval I = [a, b] of the t-axis. If f1(t) = f2(t) = 0
then the system is called homogeneous, otherwise non-
homogeneous.
•Anil Kumar, Dept. of Mathematics •11

Solutions of the linear systems

 A solution of the above system on I is a pair of functions x(t),
y(t) satisfying both the equations throughout the interval.

 x  x (t )
We shall write such a solution in the form  .
 y  y (t )

 x  e4t x  e 2 t
Example:  and  2 t
are the solutions of the
y  e  y  e
4t
 dx
 dt  x  3 y
homogeneous system (on any interval I)  .
 dy  3 x  y
 dt

•Anil Kumar, Dept. of Mathematics •12

Dr. Anil Kumar, BITS Pilani Goa Campus 6

 Example
 The equation y   y  0, 0  t  2

can be written as system of first order equations by letting y1 = y

and y2 = y'. Thus '
y1  y2
y2'   y1
 A solution to this system is
y1  sin t , y2  cos t (0  t  2 )
which is a parametric description for the unit
circle.

•13 •Anil Kumar, Dept. of Mathematics

Linear Systems
Existence and uniqueness theorem:

Theorem 1: If t0 is any point of the interval I and if x0 and y0 are

any real numbers, then the system
 dx
 dt  a1 (t ) x  b1 (t ) y  f1 (t )

 dy  a (t ) x  b (t ) y  f (t )
 dt 2 2 2

 x  x(t )
has one and only one solution  valid throughout I such
that x(t0) = x0 and y(t0) = y0.  y  y (t )

•Anil Kumar, Dept. of Mathematics •14

Dr. Anil Kumar, BITS Pilani Goa Campus 7

 Linear Systems
 Task is to study the structure of the solutions of the
homogeneous system

 dx
 dt  a1 (t ) x  b1 (t ) y
 (1)
 dy  a (t ) x  b (t ) y
 dt 2 2

It is obvious that (x = 0, y = 0) is the trivial solution.

Hence, we seek nontrivial solutions.

•15 •Anil Kumar, Dept. of Mathematics

Linear Systems
 Theorem 2: If the above homogeneous system (1) has two
solutions
 x  x1 (t )  x  x2 (t )
 and  (2)
 y  y1 (t )  y  y2 (t )
 x  Ax1 (t )  Bx2 (t )
on I = [a, b], then  (3)
 y  Ay1 (t )  By2 (t )
is also a solution of (1) on I for any arbitrary constants A and B.
(Or any linear combination of two solutions of the homogeneous
system (1) is also a solution.
 Can we say (3) contains all the solutions of (1) or it is a general
solution of (1)?
 By theorem (3) will be general solution if the constants A and B
can be chosen so as to satisfy ICs.
•16 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 8

 Linear Systems

 Definition: Wronskian of the two solutions (2) of the homo-

geneous system (1) is defined as

x1 (t ) x 2 (t )
W(t) =
y1 (t ) y 2 (t )

 Theorem 3: If the Wronskian of the two solutions (2) of the

homogeneous system (1) doesn’t vanish on I then (3) is the
general solution of the given homogeneous system on this
interval.

•17 •Anil Kumar, Dept. of Mathematics

Linear Systems

 Remark: Vanishing or non-vanishing of the Wronskian of

two solutions does not depend on the choice of t.

 Therefore we have the following theorem:

 Theorem 4: The Wronskian of the two solutions (2) of the

homogeneous system (1) is either identically zero or never
zero on I.

•18 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 9

 Linear dependence/independence
 x  x1 (t )  x  x 2 (t )
 The solutions  and 
 y  y1 (t )  y  y 2 (t )
of the homogeneous system (1) are called linearly dependent on
I if x1(t) = k x2(t) or x2(t) = k x1(t)
y1(t) = k y2(t) y2(t) = k y1(t)
for some constant k and all t in I. Otherwise linearly
independent.
 Linear dependence is equivalent to the condition that there exist
two constants A and B, at least one of them is not zero, such that
Ax1(t) + B x2(t) = 0
Ay1(t) + B y2(t) = 0 for all t in I.
•Anil Kumar, Dept. of Mathematics •19

Linear dependence/independence

Theorem 5: If the two solutions (2) of the homogeneous

system (1) are linearly independent on I, then (3) is the
general solution of (1) on I.

Proof: To prove the theorem, it suffices to show that the solutions

(2) are linearly dependent iff their Wronskian W(t) is zero.

•20 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 10

 Homogeneous linear systems with constant
coefficients
 Find a general solution of the following homogeneous system of
differential equations
 dx
 dt  a1 x  b1 y
 (1)
 dy  a x  b y
 dt 2 2

where a1, b1, a2, b2 are given constants.

 One method of finding the general solution of above system is
based on differentiating one of the two equations given,
eliminating one of the dependent variables and solving the
resulting second order linear equation.

•Anil Kumar, Dept. of Mathematics •21

EXERCISES
 By differentiating an equation of the system w.r.t. t and
eliminating any one of the dependent variable find the
general solution of the following
 dx  dx
 dt  x  2 y  dt  x  y
(1)  (2) 
 dy  3 x  2 y  dy  y
 dt  dt

Any 2nd order equation obtained from the system is

 dx not equivalent to the system in the sense that it has
 dt  x solutions that are not part of any solution of the
(3)  system.
 dy  y
 dt Although, higher order equations are equivalent to the
systems, the reverse is not true and systems are more
general.
•Anil Kumar, Dept. of Mathematics •22

Dr. Anil Kumar, BITS Pilani Goa Campus 11

 Homogeneous linear systems
Method 2:
 We know that the exponential function has the property that its
derivatives are constant multiples of the function itself. Hence it is
natural to seek solution of (1) of the type x  Ae 
mt

 (2)
y  Be mt 
Substituting in (1) we get
Ame mt  a1 Ae mt  b1 Be mt and Bme mt  a2 Ae mt  b2 Bemt
and therefore, we get
(a1  m) A  b1 B  0 
 (3)
a2 A  (b2  m) B  0 

•Anil Kumar, Dept. of Mathematics •23

Homogeneous linear systems

 The above equations have nontrivial solutions whenever

a1  m b1
0  m 2 - (a1  b2 )m  (a1b2 - a2b1 )  0
a2 b2  m

 This equation is called auxiliary equation of the system (1).

 Let m1 and m2 be the roots of the auxiliary equation.
 When m = m1, system (3) gives solutions, say A = A1 and B = B1.
 And then the corresponding nontrivial solution of system (1) is
 x  A1e m1t

 y  B1e
m1t

•Anil Kumar, Dept. of Mathematics •24

Dr. Anil Kumar, BITS Pilani Goa Campus 12

 Homogeneous linear systems
 When m = m2 system (3) gives solutions, say A = A2 and B = B2.
 Then the corresponding nontrivial solution of system (1) is
x  A2e m2t 

y  B2 e m2t 
 Thus two solutions are x  A1e m1t  x  A2 em2t 
 and 
y  B1em1t  y  B2 e m2t 
 The two solutions, got above, in general may not be linearly
independent.
 Therefore, while constructing the general solution, three cases on
the nature of the roots m1 and m2 of auxiliary equation must be
dealt carefully. (Three cases)
•Anil Kumar, Dept. of Mathematics •25

Case 1: m1 ≠ m2 and real

 When m1 ≠ m2 clearly the above set of solutions are linearly
independent.
 Hence the general solution is given by

x  c1 A1e m1t  c2 A2 em2t  In this case the two solutions are

 linearly independent?
y  c1 B1e m1t  c2 B2 e m2t 
 Example: Solve the system of linear first order equations
dx  dx 
(i )  3 x  4 y  (ii )  4x  3 y
dt  dt 
 
dy dy
 2 x  3 y   8x  6 y 
dt  dt 
•Anil Kumar, Dept. of Mathematics •26

Dr. Anil Kumar, BITS Pilani Goa Campus 13

 Case 2: m1 = m2
 In this case, we get a repeated solution as

 x  Ae mt

 y  Be
mt

 Get the second linearly independent solution as

x  ( A1  A2t )e mt and y  ( B1  B2t )emt
and hence the general solution is
x  c1 ( A1  A2t )e mt  c2 Ae mt
y  c1 ( B1  B2t )e mt  c2 Be mt (*)

where c1 and c2 are arbitrary constants.

•Anil Kumar, Dept. of Mathematics •27

Examples

dx dx
 4 x  y  5x  4 y
dt dt
(i ) (ii )
dy dy
 x  2y  x  y
dt dt

•Anil Kumar, Dept. of Mathematics •28

Dr. Anil Kumar, BITS Pilani Goa Campus 14

 Remark
 In formula (*), the constants A2 and B2 satisfy the same
linear algebraic system as the constants A and B, and that
consequently we may put
A2 = A and B2 = B
without loss of generality.

Therefore, the general solution of the system is

x  c1 ( A1  At )e mt  c2 Ae mt
y  c1 ( B1  Bt )e mt  c2 Be mt

•29 •Anil Kumar, Dept. of Mathematics

Case 3: m1, m2 are complex

 Let m = a ± ib, (where a, b are real numbers and b ≠ 0)
 For m1 = a + ib, the expected solution is of the type
 x  ( A1  iA2 )e ( a  ib ) t x  ( A1  iA2 )e at (cos bt  i sin bt )
 ( a  ib ) t

 y  ( B1  iB2 )e y  ( B1  iB2 )e at (cos bt  i sin bt )

where A1, A2, B1, B2 are real constants such that neither both
of A1, A2, nor both of B1, B2 are zero. Therefore, we have
 x  e at ( A1 cos bt  A2 sin bt )  i ( A1 sin bt  A2 cos bt )

 y  e ( B1 cos bt  B2 sin bt )  i ( B1 sin bt  B2 cos bt )
at

•30 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 15

 Case 3: m1, m2 are complex…
 Easy to see that if a pair of complex-valued functions is a
solution of the system, in which the coefficients are real
constants, then their two real parts and their two
imaginary parts are real-valued solutions.
 From this fact, two real-valued solutions are
 x  e at ( A1 cos bt  A2 sin bt )  x  e at ( A1 sin bt  A2 cos bt )
 and 
 y  e ( B1 cos bt  B2 sin bt )  y  e ( B1 sin bt  B2 cos bt )
at at

 Therefore, the general solution of the system is

x  c1e at ( A1 cos bt  A2 sin bt )  c2 e at ( A1 sin bt  A2 cos bt )
y  c1e at ( B1 cos bt  B2 sin bt )  c2e at ( B1 sin bt  B2 cos bt )

•31 •Anil Kumar, Dept. of Mathematics

Case 3: m1, m2 are complex…

 Two solutions give the following pair of real solutions.

 x  e at ( A1 cos bt  A2 sin bt )  x  e at ( A1 sin bt  A2 cos bt )

 and 
 y  e ( B1 cos bt  B2 sin bt )  y  e ( B1 sin bt  B2 cos bt )
at at

 Wronskian of the two solutions is given by

W (t )  ( A1 B2  A2 B1 )e2 at  0 as A1 B2  A2 B1  0

and hence these two solutions are linearly independent.

•Anil Kumar, Dept. of Mathematics •32

Dr. Anil Kumar, BITS Pilani Goa Campus 16

 Example

dx dx
 x  2y  4x  2 y
dt dt
(i ) (ii )
dy dy
 4x  2 y  5x  2 y
dt dt

•Anil Kumar, Dept. of Mathematics •33

General solution of non homogeneous system of

equations.
 x  Ax1 (t )  Bx 2 (t )
Theorem 6: If  is the general solution of the
 y  Ay1 (t )  By 2 (t )
homogeneous part of the non- homogeneous system
dx
 a1 (t ) x  b1 (t ) y  f1 (t )
dt
(1)
dy
 a2 (t ) x  b2 (t ) y  f 2 (t )
dt

on I and if x  x p (t ) and y  y p (t ) is any particular solution of

(1) on I, then  x  Ax1 (t )  Bx 2 (t )  x p (t ) is the general solution of the
 y  Ay1 (t )  By 2 (t )  y p (t )
system (1).
•34 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 17

 Exercise: Solution to non-homogeneous
system
 dx
 dt  a1 (t ) x  b1 (t ) y  f1 (t )
 Consider a non-homogeneous linear system 
 dy  a (t ) x  b (t ) y  f (t )
 dt 2 2 2

 dx
 dt  a1 (t ) x  b1 (t ) y
and the corresponding homogeneous system 
 dy  a (t ) x  b (t ) y
 dt 2 2

 x  x1 (t )  x  x2 (t )
If  y  y (t ) and  are linearly independent solution
 1  y  y2 (t )
of homogeneous system so that

•Anil Kumar, Dept. of Mathematics •35

Exercise: Solution to non-homogeneous

system
 x  Ax1 (t )  Bx 2 (t )

 y  Ay1 (t )  By 2 (t ) is its general solution, show that
 x  v1 (t ) x1 (t )  v2 (t ) x2 (t )

 y  v1 (t ) y1 (t )  v2 (t ) y2 (t )
will be a particular solution of non-homogeneous system if the
functions v1(t) and v2(t) satisfy the system

v1/ x1  v 2/ x 2  f 1
v1/ y1  v 2/ y 2  f 2
(Proof is Just the verification)
•Anil Kumar, Dept. of Mathematics •36

Dr. Anil Kumar, BITS Pilani Goa Campus 18

 Just the verification
 Note: This method of getting the particular solution of a non-
homogeneous system is called method of variation of
parameters

 Example: Use method of variation of parameters to find the

particular solution of the non homogeneous system
dx
 x  y  5t  2
dt
dy
 4 x  2 y  8t  8
dt
Also find the general solution.

•Anil Kumar, Dept. of Mathematics •37

Example

 Use method of undetermined coefficients to find the

particular solution of the non homogeneous system
dx
 x  y  5t  2
dt
dy
 4 x  2 y  8t  8
dt

Also find the general solution.

•38 •Anil Kumar, Dept. of Mathematics

Dr. Anil Kumar, BITS Pilani Goa Campus 19