*i.e.*you are probably on a mobile phone). Due to the nature of the mathematics on this site it is best views in landscape mode. If your device is not in landscape mode many of the equations will run off the side of your device (should be able to scroll to see them) and some of the menu items will be cut off due to the narrow screen width.

### Section 5.4 : Systems of Differential Equations

In the introduction to this section we briefly discussed how a system of differential equations can arise from a population problem in which we keep track of the population of both the prey and the predator. It makes sense that the number of prey present will affect the number of the predator present. Likewise, the number of predator present will affect the number of prey present. Therefore the differential equation that governs the population of either the prey or the predator should in some way depend on the population of the other. This will lead to two differential equations that must be solved simultaneously in order to determine the population of the prey and the predator.

The whole point of this is to notice that systems of differential equations can arise quite easily from naturally occurring situations. Developing an effective predator-prey system of differential equations is not the subject of this chapter. However, systems can arise from \(n^{\text{th}}\) order linear differential equations as well. Before we get into this however, let’s write down a system and get some terminology out of the way.

We are going to be looking at first order, linear systems of differential equations. These terms mean the same thing that they have meant up to this point. The largest derivative anywhere in the system will be a first derivative and all unknown functions and their derivatives will only occur to the first power and will not be multiplied by other unknown functions. Here is an example of a system of first order, linear differential equations.

\[\begin{align*}{{x'}_1} & = {x_1} + 2{x_2}\\ {{x'}_2} & = 3{x_1} + 2{x_2}\end{align*}\]We call this kind of system a **coupled** system since knowledge of \(x_{2}\) is required in order to find \(x_{1}\) and likewise knowledge of \(x_{1}\) is required to find \(x_{2}\). We will worry about how to go about solving these later. At this point we are only interested in becoming familiar with some of the basics of systems.

Now, as mentioned earlier, we can write an \(n^{\text{th}}\) order linear differential equation as a system. Let’s see how that can be done.

*order differential equation as a system of first order, linear differential equations. \[2y'' - 5y' + y = 0\hspace{0.25in}y\left( 3 \right) = 6\,\,\,\,\,\,\,\,\,y'\left( 3 \right) = - 1\]*

^{nd}We can write higher order differential equations as a system with a very simple change of variable. We’ll start by defining the following two new functions.

\[\begin{align*}{x_1}\left( t \right) & = y\left( t \right)\\ {x_2}\left( t \right) & = y'\left( t \right)\end{align*}\]Now notice that if we differentiate both sides of these we get,

\[\begin{align*}{{x'}_1} & = y' = {x_2}\\ {{x'}_2} & = y'' = - \frac{1}{2}y + \frac{5}{2}y' = - \frac{1}{2}{x_1} + \frac{5}{2}{x_2}\end{align*}\]Note the use of the differential equation in the second equation. We can also convert the initial conditions over to the new functions.

\[\begin{align*}{x_1}\left( 3 \right) & = y\left( 3 \right) = 6\\ {x_2}\left( 3 \right) & = y'\left( 3 \right) = - 1\end{align*}\]Putting all of this together gives the following system of differential equations.

\[\begin{align*}{{x'}_1} & = {x_2} & \hspace{0.25in}{x_1}\left( 3 \right) & = 6\\ {{x'}_2} & = - \frac{1}{2}{x_1} + \frac{5}{2}{x_2} & \hspace{0.25in}{x_2}\left( 3 \right) & = - 1\end{align*}\]We will call the system in the above example an **Initial Value Problem** just as we did for differential equations with initial conditions.

Let’s take a look at another example.

^{th}order differential equation as a system of first order, linear differential equations. \[{y^{\left( 4 \right)}} + 3y'' - \sin \left( t \right)y' + 8y = {t^2}\hspace{0.25in}y\left( 0 \right) = 1\,\,\,\,y'\left( 0 \right) = 2\,\,\,\,y''\left( 0 \right) = 3\,\,\,\,y'''\left( 0 \right) = 4\]

Just as we did in the last example we’ll need to define some new functions. This time we’ll need 4 new functions.

\[\begin{align*}{x_1} & = y & \Rightarrow \hspace{0.25in}{{x'}_1} & = y' = {x_2}\\ {x_2} & = y' & \Rightarrow \hspace{0.25in}{{x'}_2} & = y'' = {x_3}\\ {x_3} & = y'' & \Rightarrow \hspace{0.25in}{{x'}_3} & = y''' = {x_4}\\ {x_4} & = y''' & \Rightarrow \hspace{0.25in}{{x'}_4}& = {y^{\left( 4 \right)}} = - 8y + \sin \left( t \right)y' - 3y'' + {t^2} = - 8{x_1} + \sin \left( t \right){x_2} - 3{x_3} + {t^2}\end{align*}\]The system along with the initial conditions is then,

\[\begin{align*}{{x'}_1} & = {x_2} & \hspace{0.25in}{x_1}\left( 0 \right) & = 1\\ {{x'}_2} & = {x_3} & \hspace{0.25in}{x_2}\left( 0 \right) & = 2\\ {{x'}_3} & = {x_4} & \hspace{0.25in}{x_3}\left( 0 \right) & = 3\\ {{x'}_4} & = - 8{x_1} + \sin \left( t \right){x_2} - 3{x_3} + {t^2} & \hspace{0.25in}{x_4}\left( 0 \right) & = 4\end{align*}\]Now, when we finally get around to solving these we will see that we generally don’t solve systems in the form that we’ve given them in this section. Systems of differential equations can be converted to **matrix form** and this is the form that we usually use in solving systems.

First write the system so that each side is a vector.

\[\left( {\begin{array}{*{20}{c}}{{{x'}_1}}\\{{{x'}_2}}\end{array}} \right) = \left( {\begin{array}{*{20}{c}}{4{x_1} + 7{x_2}}\\{ - 2{x_1} - 5{x_2}}\end{array}} \right)\]Now the right side can be written as a matrix multiplication,

\[\left( {\begin{array}{*{20}{c}}{{{x'}_1}}\\{{{x'}_2}}\end{array}} \right) = \left( {\begin{array}{*{20}{c}}4&7\\{ - 2}&{ - 5}\end{array}} \right)\left( {\begin{array}{*{20}{c}}{{x_1}}\\{{x_2}}\end{array}} \right)\]Now, if we define,

\[\vec x = \left( {\begin{array}{*{20}{c}}{{x_1}}\\{{x_2}}\end{array}} \right)\]then,

\[\vec x' = \left( {\begin{array}{*{20}{c}}{{{x'}_1}}\\{{{x'}_2}}\end{array}} \right)\]The system can then be written in the matrix form,

\[\vec x' = \left( {\begin{array}{*{20}{c}}4&7\\{ - 2}&{ - 5}\end{array}} \right)\vec x\]We’ll start with the system from Example 1.

\[\begin{align*}{{x'}_1} & = {x_2}\hspace{0.25in}{x_1}\left( 3 \right) = 6\\ {{x'}_2} & = - \frac{1}{2}{x_1} + \frac{5}{2}{x_2}\hspace{0.25in}{x_2}\left( 3 \right) = - 1\end{align*}\]First define,

\[\vec x = \left( {\begin{array}{*{20}{c}}{{x_1}}\\{{x_2}}\end{array}} \right)\]The system is then,

\[\vec x' = \left( {\begin{array}{*{20}{c}}0&1\\{ - \frac{1}{2}}&{\frac{5}{2}}\end{array}} \right)\vec x\hspace{0.25in}\vec x\left( 3 \right) = \left( {\begin{array}{*{20}{c}}{{x_1}\left( 3 \right)}\\{{x_2}\left( 3 \right)}\end{array}} \right) = \left( {\begin{array}{*{20}{c}}6\\{ - 1}\end{array}} \right)\]Now, let’s do the system from Example 2.

\[\begin{align*}{{x'}_1} & = {x_2}\hspace{0.25in}{x_1}\left( 0 \right) = 1\\ {{x'}_2} & = {x_3}\hspace{0.25in}{x_2}\left( 0 \right) = 2\\ {{x'}_3} & = {x_4}\hspace{0.25in}{x_3}\left( 0 \right) = 3\\ {{x'}_4} & = - 8{x_1} + \sin \left( t \right){x_2} - 3{x_3} + {t^2}\hspace{0.25in}{x_4}\left( 0 \right) = 4\end{align*}\]In this case we need to be careful with the *t ^{2}* in the last equation. We’ll start by writing the system as a vector again and then break it up into two vectors, one vector that contains the unknown functions and the other that contains any known functions.

Now, the first vector can now be written as a matrix multiplication and we’ll leave the second vector alone.

\[\vec x' = \left( {\begin{array}{*{20}{c}}0&1&0&0\\0&0&1&0\\0&0&0&1\\{ - 8}&{\sin \left( t \right)}&{ - 3}&0\end{array}} \right)\vec x + \left( {\begin{array}{*{20}{c}}0\\0\\0\\{{t^2}}\end{array}} \right)\hspace{0.25in}\vec x\left( 0 \right) = \left( {\begin{array}{*{20}{c}}1\\2\\3\\4\end{array}} \right)\]where,

\[\vec x\left( t \right) = \left( {\begin{array}{*{20}{c}}{{x_1}\left( t \right)}\\{{x_2}\left( t \right)}\\{{x_3}\left( t \right)}\\{{x_4}\left( t \right)}\end{array}} \right)\]Note that occasionally for “large” systems such as this we will go one step farther and write the system as,

\[\vec x' = A\vec x + \vec g\left( t \right)\]The last thing that we need to do in this section is get a bit of terminology out of the way. Starting with

\[\vec x' = A\vec x + \vec g\left( t \right)\]we say that the system is **homogeneous** if \(\vec g\left( t \right) = \vec 0\) and we say the system is **nonhomogeneous** if \(\vec g\left( t \right) \ne \vec 0\).