# 5.2: Homogeneous Systems of Differential Equations

In this discussion we will investigate how to solve certain homogeneous systems of linear differential equations. We will also look at a sketch of the solutions.

Example

Consider the system of differential equations

\[ x' = x + y \]

\[ y' = -2x + 4y. \]

This is a system of differential equations. Clearly the trivial solution (\(x = 0\) and \(y = 0\)) is a solution, which is called a *node* for this system. We want to investigate the behavior of the other solutions. Do they approach the origin or are they repelled from it? We can graph the system by plotting direction arrows. We will calculate a few of these arrows and then use a computer to finish the plot.

Consider the point \((0,1)\). We have

\[x' = 1 \;\;\; \text{and} \;\;\; y' = 4\]

so that

\[ \dfrac{dy}{dx} = \dfrac{4}{1} = 4. \]

We plot a small arrow emanating from \((0,1)\) with slope 4.

Consider the point \((1,0)\). We have

\[x' = 1 \;\;\; \text{and} \;\;\; y' = -2\]

so that

\[ \dfrac{dy}{dx} = \dfrac{-2}{1} = -2.\]

We plot a small arrow emanating from \((1,0)\) with slope -2. Below is a computer generated graph. We call the xy-plane the *phase plane* for the differential equation and the plot the *phase portrait*.

Notice that all solutions are repelled from the origin. The origin is called an *unstable equilibrium point*.

Our next task is to find an explicit solution for the system. We write the system as

\[ \textbf{x}' = A\textbf{x} \]

where

\[A=\begin{pmatrix} 1 &1 \\ -2 & 4 \end{pmatrix}. \]

Just as with higher order differential equations, we assume that the solution is in the form

\[\textbf{x} = \textbf{z}e^{rt} \]

where **x** is a vector solution and **z** is a constant vector. We have

\[ \textbf{x}' = r \textbf{z}e^{rt}\]

so that

\[r\textbf{z}e^{rt}=A\textbf{z}e^{rt}. \]

We can divide by \(e^{rt}\) to get

\[A \textbf{z}=r\textbf{z}.\]

Finding a solution is equivalent to finding an eigenvalue of the matrix. We can use a calculator to find that the eigenvalues are

\[r = 2 \;\;\; \text{and} \;\;\; r = 3.\]

To find the constants, we can plug the eigenvalues into

\[ A - rI. \]

Plugging in \(r = 2\) gives

\[A-2I = \begin{pmatrix} -1 & -1 \\ -2 & 2 \end{pmatrix}. \]

The first row tells us that

\[ -x + y = 0 \]

and the eigenvector is

\[z_2=\begin{pmatrix} 1 \\ 1 \end{pmatrix}. \]

Plugging in \(r = 3\) gives

\[A-3I = \begin{pmatrix} -2 & 1 \\ -2 & 1 \end{pmatrix}. \]

The first row tells us that

\[-2x + y = 0 \]

and the eigenvector is

\[ z_3 = \begin{pmatrix} 1 \\ 2 \end{pmatrix}. \]

We can conclude that the general solution is

\[\begin{pmatrix} x \\ y \end{pmatrix} = c_1\begin{pmatrix}1 \\ 1\end{pmatrix}e^{2t} + c_2\begin{pmatrix}1 \\ 2 \end{pmatrix} e^{3t} \]

or that

\[ x = c_1e^{2t} + c_2e^{3t} \]

\[ y = c_1e^{2t} + 2c_2e^{3t}. \]

There is a direction relationship between the type of node the origin is and the eigenvalues of the matrix. In the example above the eigenvalues are both positive, thus both x and y are repelled from the origin as \(t\) becomes large. In particular if either of the eigenvalues are positive, then the trajectories repel from the origin. If the eigenvalues had both been negative, then both \(x\) and \(y\) approach zero as \(t\) approaches 0. Hence, all trajectories travel towards the origin for large \(t\).

The cases where the eigenvalues are complex will be studied in the next discussion.

### Contributors

- Larry Green (Lake Tahoe Community College)

Integrated by Justin Marshall.