Skip to main content

Handout Lessons 22 and 23, Systems of Differential Equations and Applications

Textbook Section(s).

This lesson is based on Sections 4.1 and 4.2 of your textbook by Edwards, Penney, and Calvis.

Systems of Differential Equations.

Definition 143.

A system of differential equations is a set of differential equations.
A solution of a system of differential equations is a set of functions that satisfy all of the differential equations in the system over some interval of the independent variable.

Example 144. Verifying a solution of a system of differential equations.

(Based on Number 23 from Section 4.1 of your textbook by Edwards, et.al.)
Consider the system of differential equations
\begin{align} \begin{cases} x' \amp = y \\ y' \amp = 6x-y \end{cases}\tag{#} \end{align}
The independent variable for this system is . Verify that \(x(t)=e^{2t}\text{,}\) \(y(t)=2e^{2t}\) is a solution of the system of differential equations in (#).
Systems of differential equations arise from:

Example 145. Multiple Springs.

(See Example 1 from Section 4.1 of your textbook by Edwards, et.al.)
Consider a horizontal dual spring system without damping or friction. The first spring is attached to a vertical wall and a mass \(m_1\) and has spring constant \(k_1\text{.}\) The second spring is attached to the mass \(m_1\) and a second mass \(m_2\text{.}\) The second spring has spring constant \(k_2\text{.}\) An external force \(f(t)\) is being applied to the system at the second mass. Figure 4.1.1 from your textbook by Edwards, et.al. represents this dual spring system.
described in detail following the image
Spring 1 attached to wall and mass \(m_1\text{.}\) Spring 2 attached to mass \(m_1\) and mass \(m_2\text{.}\) Mass \(m_2\) acted on by external force \(f(t)\text{.}\)
Figure 146.
In this diagram:
  • \(x(t) =\) the displacement to the right of \(m_1\) from its static equilibrium at time \(t\text{.}\)
  • \(y(t) =\) the displacement to the right of \(m_2\) from its static equilibrium at time \(t\text{.}\)
The system of differential equations representing this dual spring system is
\begin{align*} \begin{cases} m_1x'' \amp = -k_1x+k_2(y-x) \\ m_2y'' \amp = -k_2(y-x)+f(t) \end{cases} \end{align*}

Converting a differential equations into a system of first-order differential equations.

There are computer algorithms for solving systems of first-order differential equations. Therefore, if we can convert a single differential equaiton in to a system of first-order differential equations, then we can use the computer to solve higher order differential equations.

Example 147. Converting a higher order DE into a system of first order DE’s.

Find a system of first-order differential equatons that is equivalent to
\begin{gather} x^{(4)}+2x^{(3)}-x=e^{5t}\tag{†} \end{gather}
General Procedure: Convert
\begin{gather} x^{(n)}=f(t,x,x', \dots, x^{(n-1)})\tag{✢} \end{gather}
to
\begin{align} \begin{cases} x_1' \amp = x_2 \\ x_2' \amp = x_3 \\ \amp \vdots \\ x_{n-1}' \amp = x_n \\ x_n' \amp = f(t,x_1,x_2, \dots, x_n) \end{cases}\tag{✢✢} \end{align}
Solve (✢✢). Then \(x(t)=x_1(t)\) is a solution of (✢).

Example 148. Converting a system of higher order DE’s into a systemm of first order DE’s.

(Number 10 from Section 4.1 of your textbook by Edwards, et.al.)
Transform the system of differential equations into an equivalent system of first-order differential equations.
\begin{align*} \begin{cases} x''-5x+4y \amp = 0 \\ y''+4x-5y \amp = 0 \end{cases} \end{align*}

Solving two-dimensional systems of differential equations.

We now turn our attention to solving some systems of differential equations. In this course, we will focus on solutions of two-dimensional systems of differential equations. These systems have two unknown functions. Normally we will denote these unknown fuctions by \(x(t)\) and \(y(t)\text{,}\) or just \(x\) and \(y\) for short. In this section of the notes, we will use a substitution approach and direction fields to help us understand the solutions of these systems of equations. In the last section of these notes, we will use the method of elimination to solve a special category of these systems of equations.

Substitution.

In the next example, we will use a substitution approach to solve a two-dimensional system of differential equations. As we are working through this problem, notice that the technique does not produce a clear-cut algorithm for solving these systems. It requires a bit of creativity and is very dependent on the equations in the system.

Example 149. Solving a system using substitution.

(Number 22 from Section 4.1 of your textbook by Edwards, et.al.)
Solve the system of differential equations.
\begin{align} \begin{cases} x' \amp = 8y \\ y' \amp = -2x \end{cases}\tag{‑} \end{align}

Using direction fields to understand solutions of two-dimensional systems of differential equations.

There are various programs for plotting direction fields of two-dimensional systems of differential equations. I will describe two of them in these notes.

Darryl Nester’s Slope and Direction Field Plotter.

  • Then click on the direction fields several times to show some solutions. (The points where you click are providing initial values.) This isn’t perfect. You will notice that near the center of the graph we have some "fat ellipses." Away from the center, we appear to get some jagged lines that stop abruptly. Even with the imperfections, it does give us a sense that possibly the solution curves are elliptical.
    described in detail following the image
    The phase plane and some solution curves.

pplane by John C. Polking in the Department of Mathematics at Rice University.

  • This program is available on the department website. I will also make a link to it in Brightspace with this lecture.
  • There are two versions of the program, a MATLAB version and a Java version. I will discuss the Java version in these notes.
  • Once you download the java version, you will have a pplane.jar file on your computer. The directions on the course website say to double-click this file in your computer directory. I am running Linux. Before I could run the file, I had to change the permissions so that it was executable. I do not know if this will be a problem on Windows or Macs, but I wanted to warn you about this possibility.
  • Once the file could be run, I double-clicked and 4 windows opened up. The first window contains a copyright statement. You can safely close this window.
  • The PPLANE Equation Window allows you to enter the system and to change the viewing window. I entered my system and changed the viewing window to \([-4,4] \times [-4,4]\text{.}\)
    described in detail following the image
    The PPLANE Equation Window.
    I then clicked on the Graph Phase Plane button.
  • The PPLANE Phase Plane window then plots the phase plane for the system that I entered. Clicking on points in the PPLANE Phase Plane window produced specific solutions. These solutions are nicer than those that appeared on Nester’s website, but Nester’s website is easy to access when you are on a public computer.
    described in detail following the image
    A phase plane for the system.
    described in detail following the image
    A phase plane for the system with some solution curves.

Example 150. Verifying that the solution curves are ellipses.

From the direction field plotters, we have an idea that our solution curves could be ellipses:
\begin{align*} x(t) \amp = A\cos(4t)+B\sin(4t)\\ y(t) \amp = -\frac{1}{2}A\sin(4t)+\frac{1}{2}B\cos(4t) \end{align*}
Can we verify this algebraically?

Solving systems of linear differential equations with constant coefficients by elimination.

In general, solving systems of differential equations can be difficult and may require some creativity. When the system contains only linear differential equations with constant coefficients, then we can apply the Elimination Algorithm to solve the system. This is similar to the elimination method you used in high school to solve systems of linear equations.
  • Combine equations to eliminate unknown functions until you have a single differential equation containing only one unknown function.
  • Solve that equation for the unknown function.
  • Use that solution to help you find the other unknown functions.
This approach is more than substitution or the direction field tools that we used previously. For those of you who have already done your supplemental work on section 5.1, you may notice that we are gearing up for some matrix notation soon. The key to the elimination method is the behavior of .
In general, differential operators do commute. The fact that polynomial differential operators do commute is a very special case.
In this set of notes, \(L_i\) will always refer to a polynomial differential operator.
The Elimination Algorithm can be used to solve a system of differential equations of the form
\begin{align} \begin{cases} L_1x+L_2y \amp = f_1(t) \\ L_3x+L_4y \amp = f_2(t) \end{cases}\tag{✠} \end{align}
We could spend a lot of time describing the algorithm, but I think that the algorithm becomes clear after you see a few examples.

Example 152. Solving systems of DE’s by elimination.

(Number 8 from Section 4.2 of your textbook by Edwards, et.al.)
Solve the system of differential equations.
\begin{align*} \begin{cases} x' \amp = 2x+y \\ y' \amp = x+2y-e^{2t} \end{cases} \end{align*}

Example 153. Solving systems of DE’s by elimination.

(Number 12 from Section 4.2 of your textbook by Edwards, et.al.)
Solve the system of differential equations.
\begin{align*} \begin{cases} x'' \amp = 6x-2y \\ y'' \amp = -3x+7y \end{cases} \end{align*}
I would encourage you to take another look at number 22 from Section 4.1 of your textbook by Edwards, et. al. We did this problem in ExampleΒ 149. Can you use elimination to find the solution we found before?