Skip to main content

Section 4.1 Homogeneous Linear Equations

A differential equation of the form
\begin{equation*} a(t) x'' + b(t) x' + c(t) x = g(t) \end{equation*}
is called a second-order linear differential equation. We will first consider the case
\begin{equation*} a x'' + b x' + cx = 0, \end{equation*}
where \(a\text{,}\) \(b\text{,}\) and \(c\) are constants and \(a \neq 0\text{.}\) An equation of this form is said to be homogeneous with constant coefficients. We already know how to solve such equations since we can rewrite them as a system of first-order linear equations. Thus, we can find the general solution of a homogeneous second-order linear differential equation with constant coefficients by computing the eigenvalues and eigenvectors of the matrix of the corresponding system.

Subsection 4.1.1 RLC Circuits

Recall the RC circuits that we studied earlier (see Section 1.3). Such circuits contained a voltage source, a capacitor, and a resistor. A battery or generator is an example of a voltage source, and a toaster or an electric stove is an example of something that might provide a resistance in a circuit. Capacitors store an electrical charge and are used in electronic flashes for cameras. We will now add an inductor such as a solenoid, a coil that generates a magnetic field. Inductor applications include transformers, power supplies, televisions, and radios. Our new circuit is called an RLC circuit (Figure 4.1.1).
an electrical circuit containing a voltage source, and resistor, a capacitor, and an inductor
Figure 4.1.1. An RLC Circuit
Current, \(I(t)\text{,}\) is the rate at which a charge flows through this circuit and is measured in amperes or amps. We assign a direction to the current, and a current flowing in the opposite direction will be given negative values. The impressed voltage, \(E(t)\text{,}\) is measured in volts, the resistance \(R\) is measured in ohms, and the capacitance \(C\) is measured in farads. The charge on the capacitor \(Q(t)\) at time \(t\) is measured in coulombs. Inductance on the coil, \(L\text{,}\) is measured in henrys.
The following laws from physics govern how our circuit behaves.
  • \(I = \dfrac{dQ}{dt}\text{.}\)
  • The voltage drop across a resistor is \(IR\) (Ohm’s Law).
  • The voltage drop across a capacitor is \(Q/C\text{.}\)
  • The voltage drop across an inductor is \(L (dI/dt)\text{.}\)
  • In a closed circuit the impressed voltage is equal to the sum of the voltage drops in the rest of the circuit (Kirchhoff’s Second Law).
Applying Kirchhoff’s Second Law to our circuit, we have the differential equation
\begin{equation} L \frac{dI}{dt} + RI + \frac{1}{C} Q = E(t)\tag{4.1.1} \end{equation}
or
\begin{equation*} LQ'' + RQ' + \frac{1}{C} Q = E(t). \end{equation*}
Differentiating both sides of (4.1.1), we have
\begin{equation*} L I'' + RI' + \frac{1}{C} I = E'(t). \end{equation*}
For example, we might consider an RLC circuit with \(R = 1\text{,}\) \(L = 1\text{,}\) and \(C = 1\text{.}\) At \(t = 0\) when both \(I(0) =0\) and \(I'(0) = Q(0) = 0\text{,}\) the impressed voltage on the circuit is given by \(E(t) = \sin(t)\text{.}\) Our equation becomes
\begin{equation*} I'' + I' + I = E'(t) = \cos t. \end{equation*}
This is an example of a second-order linear differential equation.

Subsection 4.1.2 Second-Order Linear Equations

Suppose that we have a homogeneous second-order linear differential equation with constant coefficients,
\begin{equation} a x'' + b x' + cx = 0.\tag{4.1.2} \end{equation}
The goal of this section is to be able to solve all such equations. However, we did a great deal of work finding unique solutions to systems of first-order linear systems equations in Chapter 3. Our efforts are now rewarded. Since each second-order homogeneous system with constant coefficients can be rewritten as a first-order linear system, we are guaranteed the existence and uniqueness of solutions. Indeed, we can rewrite (4.1.2) as a system of first-order linear equations,
\begin{align*} x' & = y\\ y' & = -\frac{c}{a} x - \frac{b}{a} y, \end{align*}
and then find the general solution by computing the eigenvalues and eigenvectors of the matrix of the corresponding system.

Example 4.1.2.

Solutions of a linear system \({\mathbf x}' = A {\mathbf x}\) often include terms of the form \(e^{r t}\text{.}\) It makes sense that solutions to equation (4.1.2) take the same form. Consider the equation
\begin{equation} x'' + 3 x' - 10 x = 0.\tag{4.1.3} \end{equation}
If we assume that a solution is of the form \(e^{rt}\text{,}\) we can substitute this expression into the left-hand side of (4.1.3) to obtain
\begin{align*} \frac{d^2}{dt^2} e^{rt} + \frac{d}{dt} 3 e^{rt} - 10 e^{rt} & = r^2 e^{rt} + 3 r e^{rt} - 10 e^{rt}\\ & = (r^2 + 3r - 10) e^{rt}\\ & = (r + 5)(r - 2) e^{rt}. \end{align*}
Since \(e^{rt}\) is never zero, we find that \((r + 5)(r - 2) = 0\) or \(r = -5\) or \(2\text{.}\) Thus, we have two solutions
\begin{equation*} x_1(t) = e^{-5t} \mbox{ and } x_2(t) = e^{2t}. \end{equation*}
By the Principle of Superposition,
\begin{equation} x(t) = c_1 x_1(t) + c_2 x_2(t) = c_1 e^{-5t} + c_2 e^{2t}\tag{4.1.4} \end{equation}
is a solution to \(x'' + 3x' - 10 x = 0\text{.}\)
Indeed, this is the general solution of our second-order equation since we have a one-to-one correspondence between the solutions of
\begin{equation*} x'' + 3x' - 10 x = 0 \end{equation*}
and the system
\begin{align*} x' & = y\\ y' & = 10x - 3 y. \end{align*}
The matrix associated with this system
\begin{equation*} A = \begin{pmatrix} 0 \amp 1 \\ 10 \amp -3 \end{pmatrix} \end{equation*}
has characteristic polynomial \(\lambda^2 + 3 \lambda - 10\text{.}\) The eigenvalues of \(A\) are \(\lambda_1 = -5\) and \(\lambda_2 = 2\) with eigenvectors \(\mathbf v_1 = (1, -5)\) and \(\mathbf v_2 = (1, 2)\text{,}\) respectively. Consequently, the solution to our system is
\begin{equation*} \begin{pmatrix} x(t) \\ y(t) \end{pmatrix} = c_1 e^{-5t} \begin{pmatrix} 1 \\ -5 \end{pmatrix} + c_2 e^{2t} \begin{pmatrix} 1 \\ 2 \end{pmatrix}, \end{equation*}
which agrees with (4.1.4).
In general, suppose that
\begin{equation*} a x'' + b x' + c x = 0, \end{equation*}
where \(a \neq 0\text{.}\) Applying the strategy in Example 4.1.2, we can find the general solution for this equation by finding the roots of the quadratic polynomial \(a \lambda^2 + b \lambda + c\text{,}\)
\begin{equation*} \lambda = - \frac{-b \pm \sqrt{b^2 - 4ac}}{2a}. \end{equation*}
If \(b^2 - 4ac \gt 0\text{,}\) we have real roots
\begin{equation*} \lambda_1 = \frac{-b + \sqrt{b^2 - 4ac}}{2a} \quad \text{and} \quad \lambda_2 = \frac{-b - \sqrt{b^2 - 4ac}}{2a}, \end{equation*}
and the solution to our second-order differential equation is
\begin{equation} x(t) = c_1 e^{\lambda_1 t} + c_2 e^{\lambda_2 t},\tag{4.1.5} \end{equation}
where \(c_1\) and \(c_2\) are arbitrary constants.
To prove that equation (4.1.5) is indeed the general solution to the second-order equation \(a x'' + b x' + c x = 0\text{,}\) we can study the equivalent system of linear equations. If we let \(y = x'\text{,}\) the corresponding linear system is \(\mathbf x' = A \mathbf x\text{,}\) where
\begin{equation*} A = \begin{pmatrix} 0 \amp 1 \\ -c/a \amp -b/a \end{pmatrix}. \end{equation*}
The characteristic polynomial of \(A\) is
\begin{equation*} p(\lambda) = \det(A - \lambda I) = \lambda^2 + \frac{b}{a} \lambda + \frac{c}{a}. \end{equation*}
The roots of \(p(\lambda)\) are the same as the roots of \(a\lambda^2 + b\lambda + c\text{.}\)
If \(b^2 - 4ac \gt 0\text{,}\) we have real roots
\begin{equation*} \lambda_1 = \frac{-b + \sqrt{b^2 - 4ac}}{2a} \quad \text{and} \quad \lambda_2 = \frac{-b - \sqrt{b^2 - 4ac}}{2a}. \end{equation*}
We can find eigenvectors
\begin{align*} \mathbf v_1 \amp = \begin{pmatrix} 1 \\ (-b + \sqrt{b^2 - 4ac}\,)/2a \end{pmatrix} = \begin{pmatrix} 1 \\ \lambda_1 \end{pmatrix}\\ \mathbf v_1 \amp = \begin{pmatrix} 1 \\ (-b - \sqrt{b^2 - 4ac}\,)/2a \end{pmatrix} = \begin{pmatrix} 1 \\ \lambda_2 \end{pmatrix} \end{align*}
for \(\lambda_1\) and \(\lambda_2\text{,}\) respectively. Thus, the general solution to the system of differential equations \(\mathbf x' = A \mathbf x\) is
\begin{equation*} \mathbf x(t) = \begin{pmatrix} x(t) \\ y(t) \end{pmatrix} = \begin{pmatrix} x(t) \\ x'(t) \end{pmatrix} = c_1 e^{\lambda_1 t} \begin{pmatrix} 1 \\ \lambda_1 \end{pmatrix} + c_2 e^{\lambda_2 t} \begin{pmatrix} 1 \\ \lambda_2 \end{pmatrix}, \end{equation*}
which agrees with (4.1.5).

Example 4.1.3.

Now let us solve the initial value problem
\begin{align*} x'' + 4 x' + 5 x & = 0\\ x(0) & = 1\\ x'(0) & = 1. \end{align*}
Again, we will assume that our solution has the form \(x(t) = e^{rt}\text{.}\) Substituting this function into our differential equation, we find that
\begin{equation*} 0 = x'' + 4 x' + 5 x = r^2 e^{rt} + 4r e^{rt} + 5 e^{rt} = (r^2 + 4r + 5) e^{rt}. \end{equation*}
As in Example 4.1.2, \(r^2 + 4r + 5 = 0\text{;}\) however, the roots of this polynomial are complex,
\begin{equation*} r = \frac{-4 \pm \sqrt{-4}}{2} = - 2 \pm i. \end{equation*}
Using Euler’s formula, we can find a complex solution
\begin{equation*} x(t) = e^{(-2 + i)t} = e^{-2t} e^{it} = e^{-2t} (\cos t + i \sin t). \end{equation*}
The real and imaginary parts of our solution are
\begin{align*} x_1(t) \amp = e^{-2t} \cos t\\ x_2(t) \amp = e^{-2t} \sin t, \end{align*}
respectively. We claim that both \(x_1(t)\) and \(x_2(t)\) are solutions to our differential equation. Indeed, since \(x(t) = x_1(t) + i x_2(t)\) is a solution,
\begin{align*} 0 \amp = ax'' + bx' + cx\\ \amp = a(x_1 + i x_2)'' + b(x_1 + i x_2)' + c(x_1 + i x_2)\\ \amp = (ax_1'' + bx_1' + cx_1) + i(ax_2'' + bx_2' + cx_2). \end{align*}
Since the real part and the imaginary part of \(x(t)\) must both be zero, we can conclude that \(ax_1'' + bx_1' + cx_1 = 0\) and \(ax_2'' + bx_2' + cx_2 = 0\text{.}\) Therefore, the general solution to our equation is
\begin{equation*} x(t) = c_1 e^{-2t} \cos t + c_2 e^{-2t} \sin t. \end{equation*}
To apply our initial conditions \(x(0) = 1\) and \(x'(0) = 1\text{,}\) we first calculate
\begin{equation*} x'(t) = -2 e^{-2t} (c_1 \cos t + c_2 \sin t) + e^{-2t} (- c_1 \sin t + c_2 \cos t). \end{equation*}
Thus,
\begin{align*} 1 \amp = x(0) = c_1\\ 1 \amp = x'(0) = -2c_1 + c_2, \end{align*}
and \(c_1 =1\) and \(c_2= 3\text{.}\) Hence, the solution to our initial value problem is
\begin{equation*} x(t) = e^{-2t} \cos t + 3 e^{-2t} \sin t. \end{equation*}
As before, the corresponding linear sytem is \(\mathbf x' = A \mathbf x\text{,}\) where
\begin{equation*} A = \begin{pmatrix} 0 \amp 1 \\ -c/a \amp -b/a \end{pmatrix}. \end{equation*}
If \(b^2 - 4ac \lt 0\text{,}\) the eigenvalue of \(A\) are \(\lambda = \alpha + i \beta\) and \(\overline{\lambda} = \alpha - i \beta\text{,}\) where
\begin{equation*} \alpha = - \frac{b}{2a} \quad \text{and} \quad \beta = \frac{\sqrt{4ac - b^2}}{2a}. \end{equation*}
The vector \(\mathbf v = (1, \alpha + i \beta)\) is an eigenvector for \(\lambda\text{.}\) Thus, a solution to our system of differential equations is
\begin{align*} \mathbf x(t) \amp = \begin{pmatrix} x(t) \\ y(t) \end{pmatrix}\\ \amp = \begin{pmatrix} x(t) \\ x'(t) \end{pmatrix}\\ \amp = e^{(\alpha + i \beta)t} \begin{pmatrix} 1 \\ \alpha + i \beta \end{pmatrix}\\ \amp = e^{\alpha t} (\cos \beta t + i \sin \beta t) \begin{pmatrix} 1 \\ \alpha + i \beta \end{pmatrix}\\ \amp = e^{\alpha t} \begin{pmatrix} \cos \beta t \\ \alpha \cos \beta t - \beta \sin \beta t \end{pmatrix} + i e^{\alpha t} \begin{pmatrix} \sin \beta t \\ \alpha \sin \beta t + \beta \cos \beta t \end{pmatrix}. \end{align*}
Taking the real and imaginary parts of \(x(t)\text{,}\) we obtain two real solutions to the system, \(x_1(t) = e^{\alpha t} \cos \beta t\) and \(x_2(t) = e^{\alpha t} \sin \beta t\text{.}\) Therefore, the general solution to \(ax_1'' + bx_1' + cx_1 = 0\) is
\begin{equation*} x(t) = c_1 e^{\alpha t} \cos \beta t + c_2 e^{\alpha t} \sin \beta t. \end{equation*}
Given a second-order linear differential equation with constant coefficients, \(ax'' + bx' + cx = 0\text{,}\) our strategy has been to solve the characteristic equation \(a \lambda^2 + b \lambda + c = 0\) to obtain two linearly independent solutions. We have covered the case where this equation has two distinct real solutions as well as when there are complex solutions, but what if there is only a single real solution \(\lambda = -b/2a\text{?}\)

Example 4.1.4.

Consider the equation
\begin{equation*} x'' + 2x' + x = 0. \end{equation*}
If we choose \(e^{\lambda t}\) as our guess, we find
\begin{align*} x'' + 2 x' + x & = \lambda^2 e^{\lambda t} + 2\lambda e^{\lambda t} + e^{\lambda t}\\ & = e^{\lambda t} (\lambda + 1)^2\\ & = 0. \end{align*}
Thus, \(\lambda = -1\) and we have a solution \(x_1(t) = e^{-t}\text{.}\)
In order to find a general solution to \(x'' + 2x' + x = 0\text{,}\) we must find a second solution that is not a multiple of \(x_1(t) = e^{-t}\text{.}\) Since we already know that \(c x_1(t)\) is a solution to our differential equation, we will try to generalize this observation by replacing \(c\) with a nonconstant function \(v(t)\) and then try to determine \(v(t)\) so that \(v(t) x_1(t)\) is a solution to \(x'' + 2x' + x = 0\text{.}\) Indeed, if
\begin{equation*} x(t) = v(t) x_1(t) = v(t) e^{-t}, \end{equation*}
then
\begin{equation*} x'(t) = v(t) x_1'(t) + v'(t) x_1(t) = -v(t) e^{-t} + v'(t) e^{-t} \end{equation*}
and
\begin{align*} x''(t) & = v''(t) x_1(t) + 2v'(t) x_1'(t) + v(t) x_1''(t)\\ & = v''(t) e^{-t} - 2 v'(t) e^{-t} + v(t) e^{-t}. \end{align*}
Consequently,
\begin{align*} x'' + 2 x' + x & = [v'' e^{-t} - 2 v' e^{-t} + v e^{-t}] + 2 [-v e^{-t} + v' e^{-t}] + [v e^{-t}]\\ & = e^{-t} v''\\ & = 0, \end{align*}
and \(v'' = 0\text{.}\) Therefore, \(v = c_1 t + c_2\text{.}\) Letting \(c_1 = 1\) and \(c_2 = 0\text{,}\) we can assume that \(v(t) = t\text{,}\) and the second solution to our equation is \(x = t e^{-t}\text{.}\) Hence, the general solution to \(x'' + 2x' + x = 0\) is
\begin{equation*} x(t) = c_1 e^{-t} + c_2 t e^{-t}. \end{equation*}
We leave it as an exercise to show that our solution agrees with the solution that we would obtain from solving the equivalent first-order linear system.
The technique that we have used in Example 4.1.4 is called reduction of order. We leave it as an exercise to show that this technique works in general. That is, given a second-order linear differential equation
\begin{equation*} a x'' + bx' + cx = 0 \end{equation*}
such that \(b^2 - 4ac = 0\text{,}\) then the general solution is given by
\begin{equation*} x(t) = c_1 e^{\lambda t} + c_2 t e^{\lambda t}, \end{equation*}
where \(\lambda = -b/2a\text{.}\)

Activity 4.1.1. Solving Second-Order Homogeneous Linear Differential Equations.

Solve each of the following initial value problems.
(a)
\begin{align*} x'' - 5 x' -14 x & = 0\\ x(0) & = 0\\ x'(0) & = 2. \end{align*}
(b)
\begin{align*} 6x'' - 11 x' + 3 x & = 0\\ x(0) & = 0\\ x'(0) & = 2. \end{align*}
(c)
\begin{align*} x'' - 14 x' + 53 x & = 0\\ x(0) & = 1\\ x'(0) & = 1. \end{align*}
(d)
\begin{align*} 4x'' - 20 x' + 41 x & = 0\\ x(0) & = 0\\ x'(0) & = 1. \end{align*}
(e)
\begin{align*} x'' -1 4 x' +49 x & = 0\\ x(0) & = 1\\ x'(0) & = 0. \end{align*}

Subsection 4.1.3 Classifying Harmonic Oscillators

Recall from Subsection 1.1.3 that we can model harmonic motion using the equation
\begin{equation*} m \frac{d^2 x}{dt^2} + b \frac{dx}{dt} + k x = 0. \end{equation*}
In the case of a spring-mass system, \(m\) is the oscillating mass, \(b\) is the damping coefficient, and \(k\) is the spring constant. It is important to remember that both \(m\) and \(k\) are positive constants and \(b \geq 0\text{.}\)
If \(b = 0\text{,}\) then the oscillator is undamped In this case,
\begin{equation} m \frac{d^2 x}{dt^2} + k x = 0.\tag{4.1.6} \end{equation}
The characteristic equation is
\begin{equation*} m \lambda^2 + k = 0, \end{equation*}
and we have eigenvalues \(\lambda = \pm i \sqrt{k/m}\text{.}\) Hence, the complex solution to our undamped oscillator is
\begin{equation*} e^{i \omega t} = \cos \omega t + i \sin \omega t, \end{equation*}
where \(\omega = \sqrt{k/m}\text{.}\) Since both the real and imaginary parts of the complex solution are also solutions to (4.1.6), the general solution to the undamped harmonic oscillator is
\begin{equation*} x(t) = c_1 \cos \omega t + c_2 \sin \omega t, \end{equation*}
giving us the position of the mass at time \(t\text{.}\) Now, of course, it is easy to determine the velocity of the mass at time \(t\) to be
\begin{equation*} v(t) = x'(t) = - c_1 \omega \sin \omega t + c_2 \omega \cos \omega t. \end{equation*}

Example 4.1.5.

Suppose that an undamped harmonic oscillator is modeled by the initial value problem
\begin{align*} \frac{d^2x}{dt^2} + 9 x \amp = 0\\ x(0) \amp = 2\\ x'(0) \amp = 1. \end{align*}
We can quickly determine the solution of this initial value problem to be
\begin{align*} x(t) \amp = 2 \cos 3t +\frac{1}{3} \sin 3t\\ v(t) \amp = \cos 3t - 6 \sin 3t, \end{align*}
where \(v(t) = x'(t)\) is the velocity of the oscillator (Figure 4.1.6). Examining the phase plane of the undamped oscillator, we find that the period of the oscillations is given by \(2 \pi / \omega = 2 \pi / 3 \approx 2.094\) (Figure 4.1.7).
position and velocity of an undamped harmonic oscillator plotted agains time
Figure 4.1.6. Position and velocity of an undamped harmonic oscillator
the velocity of a harmonic oscillator plotted against the position forming an ellipse in the position-velocity plane
Figure 4.1.7. The phase plane of an undamped harmonic oscillator
If we add damping to the oscillator, the equation becomes
\begin{equation} m \frac{d^2 x}{dt^2} + b \frac{dx}{dt} + k x = 0.\tag{4.1.7} \end{equation}
where \(b \gt 0\text{.}\) The charactersitic equation of (4.1.7) is
\begin{equation*} m\lambda^2 + b\lambda + k = 0\text{,} \end{equation*}
which has roots
\begin{equation} \lambda = \frac{-b \pm \sqrt{b^2 - 4mk}}{2m}.\tag{4.1.8} \end{equation}
There are three possible types of types of motion for the oscillator depending on the nature of the roots of (4.1.8).
  • If the damping value of \(b\) is small when compared to \(4mk\text{,}\) then \(b^2 - 4mk \lt 0\) and the roots of (4.1.8) will be complex. Furthermore, the real part of each root, \(-b/ 2m\text{,}\) is always negative. In such a situation, we say that the oscillator is under-damped.
  • If the damping value of \(b\) is large , then \(b^2 - 4mk \gt 0\text{,}\) and we obtain distinct real negative roots for (4.1.8). The oscillator is over-damped.
  • Finally, we say that the oscillator is critically-damped if \(b^2 - 4mk = 0\text{.}\)

Example 4.1.8. An Under-Damped Oscillator.

Suppose that an oscillator is modeled by the initial value problem
\begin{align*} \frac{d^2 x}{dt^2} + 0.4 \frac{dx}{dt} + 1.04 x \amp = 0\\ x(0) \amp = 0\\ x'(0) \amp = 1. \end{align*}
Notice that the damping \(b = 0.4\) is very small compared with the spring constant \(k = 1.04\text{.}\) The characteristic equation of the differential equation is \(\lambda^2 + 0.4 \lambda + 1.04 = 0\text{,}\) which has roots \(\lambda = -0.2 \pm i\text{.}\) Therefore, the complex solution must be
\begin{equation*} x(t) = e^{(-0.2 \pm i)t} = e^{-0.2t} ( \cos t + i \sin t), \end{equation*}
and the general solution must be
\begin{equation*} x(t) = c_1 e^{-0.2t} \cos t + c_2 e^{-0.2t} \sin t. \end{equation*}
Applying the initial conditions, our solution becomes
\begin{align*} x(t) \amp = e^{-0.2t} \sin t\\ v(t) \amp = x'(t) = e^{-0.2t} (\cos t - 0.2 \sin t). \end{align*}
Notice that the period of the oscillations, \(2 \pi / \omega = 2 \pi \approx 6.283\text{,}\) does not change; however, the amplitude slowly decreases (Figure 4.1.9 and Figure 4.1.10).
velocity and position plotted against time with oscillations becoming smaller as time increases
Figure 4.1.9. Position and velocity of an under-damped harmonic oscillator
the velocity of a harmonic oscillator plotted against the position spiraling towards the origin on the position-velocity plane
Figure 4.1.10. The phase plane of an under-damped harmonic oscillator

Example 4.1.11. An Over-Damped Oscillator.

We can expect a different type of behavior in the case of an over-damped oscillator. For example,
\begin{align*} \frac{d^2 x}{dt^2} + 4 \frac{dx}{dt} + 3 x \amp = 0\\ x(0) \amp = 0\\ x'(0) \amp = 1. \end{align*}
The characteristic equation of this initial value problem is
\begin{equation*} \lambda^2 + 4 \lambda + 3 = (\lambda + 3)(\lambda + 1) = 0; \end{equation*}
hence, we have the general solution
\begin{equation*} x(t) = c_1 e^{-t} + c_2 e^{-3t}. \end{equation*}
Applying the initial conditions, our solution is
\begin{align*} x(t) \amp = x(t) = \frac{1}{2}e^{-t} - \frac{1}{2} e^{-3t}\\ v(t) \amp = x'(t) = -\frac{1}{2}e^{-t} + \frac{3}{2} e^{-3t}. \end{align*}
Notice that the damping is too strong for any oscillations to occur (Figure 4.1.12 and Figure 4.1.13).
velocity and position plotted against time with the two curves approaching zero for large times
Figure 4.1.12. Position and velocity of an over-damped harmonic oscillator
the velocity of a harmonic oscillator plotted against the position with the solution curve approaching the origin on the position-velocity plane
Figure 4.1.13. The phase plane of an over-damped harmonic oscillator

Example 4.1.14. A Critically-Damped Oscillator.

As we increase the damping, the oscillations will cease to occur for some value of \(b\text{.}\) This will happen when \(b^2 - 4mk = 0\text{.}\) the At this point we have critical damping. Consider the system
\begin{align*} \frac{d^2 x}{dt^2} + 4 \frac{dx}{dt} + 4 x \amp = 0\\ x(0) \amp = 0\\ x'(0) \amp = 1. \end{align*}
The general solution to this initial value problem is
\begin{equation*} x(t) = c_1 e^{-2t} + c_2 t e^{-2t}. \end{equation*}
The solution to the initial value problem is
\begin{align*} x(t) \amp = x(t) = t e^{-2t}\\ v(t) \amp = x'(t) = (1 - 2t)e^{-2t}. \end{align*}
Although we see that no oscillations for this oscillator (Figure 4.1.15 and Figure 4.1.16), oscillations will commence as soon as we start to reduce the damping constant \(b = 4\text{.}\)
velocity and position plotted against time with the two curves approaching zero for large times
Figure 4.1.15. Position and velocity of a critically damped harmonic oscillator
the phase plane of a critically damped harmonic oscillator with the solution curve approaching the origin
Figure 4.1.16. The phase plane of a critically damped harmonic oscillator

Subsection 4.1.4 Important Lessons

  • A second-order linear differential equation with constant coefficients is an equation of the form
    \begin{equation*} a x'' + bx' + cx = 0. \end{equation*}
    We can guess the solution to this equation. Since we can rewrite this equation as a system of first-order linear differential equations, we can determine the general solution to \(a x'' + bx' + cx = 0\text{.}\)
  • Suppose that
    \begin{equation*} a x'' + b x' + c x = 0, \end{equation*}
    where \(a \neq 0\) and \(b^2 - 4ac \gt 0\text{.}\) If the roots of \(ar^2 + br + c\) are \(r_1\) and \(r_2\text{,}\) the general solution to this differential equation is
    \begin{equation*} x(t) = c_1 e^{r_1 t} + c_2 e^{r_2 t}. \end{equation*}
  • If \(b^2 - 4ac \lt 0\text{,}\) the differential equation
    \begin{equation*} a x'' + b x' + c x = 0 \end{equation*}
    has a general solution
    \begin{equation*} x(t) = c_1 e^{\alpha t} \cos \beta t + c_2 e^{\alpha t} \sin \beta t, \end{equation*}
    where \(\alpha \pm i \beta\) are the roots of \(a r^2 + br + c = 0\text{.}\)
  • If \(b^2 - 4ac = 0\text{,}\) the differential equation
    \begin{equation*} a x'' + b x' + c x = 0 \end{equation*}
    has a general solution
    \begin{equation*} x(t) = c_1 e^{- bt/2a} + c_2 t e^{-bt/2a}. \end{equation*}
  • A simple harmonic oscillator can be modeled by the equation
    \begin{equation*} m \frac{d^2 x}{dt^2} + b \frac{dx}{dt} + k x = 0, \end{equation*}
    where \(m \gt 0\text{,}\) \(k \gt 0\text{,}\) and \(b \geq 0\text{.}\) There are three possible types of motion for the oscillator depending on the sign of \(b^2 - 4mk\text{.}\)
    • If \(b^2 - 4mk \lt 0\text{,}\) the oscillator is under-damped.
    • If \(b^2 - 4mk \gt 0\text{,}\) the oscillator is over-damped.
    • If \(b^2 - 4mk = 0\text{,}\) the oscillator is critically damped.

Reading Questions 4.1.5 Reading Questions

1.

What is the characteristic equation of \(ax'' + bx' + cx = 0\text{?}\)

2.

Describe the possible types of damping of a harmonic oscillator?

Exercises 4.1.6 Exercises

Finding General Solutions.

Find the general solution for each equation in Exercise Group 4.1.6.1–10.
1.
\(\dfrac{d^2 y}{dx^2} - y = 0\)
2.
\(x'' - 2x' - 8x = 0\)
3.
\(y'' + 5y = 0\)
4.
\(\dfrac{d^2 x}{dt^2} + 6 \dfrac{dx}{dt} + 5x = 0\)
5.
\(x'' - 10x' + 25x = 0\)
6.
\(\dfrac{d^2 y}{dx^2} - 2\dfrac{dy}{dx} + 4y = 0\)
7.
\(y'' - 8y' + 4y = 0\)
8.
\(\dfrac{d^2 x}{dx^2} + 3\dfrac{dx}{dt} - 10x = 0\)
9.
\(\dfrac{d^2 Q}{dt^2} - 4 \dfrac{dQ}{dt} + 9Q = 0\)
10.
\(\dfrac{d^2 y}{dt^2} + 6 \dfrac{dy}{dt} + 9y = 0\)

Solving Initial Value Problems.

Solve the initial value problems in Exercise Group 4.1.6.11–20.
11.
\(\dfrac{d^2 y}{dx^2} - y = 0\text{,}\) \(y(0) = 1\text{,}\) \(y'(0) = 0\)
12.
\(x'' - 2x' - 8x = 0\text{,}\) \(x(0) = 1\text{,}\) \(x'(0) = 2\)
13.
\(y'' + 5y = 0\text{,}\) \(y(0) = 1\text{,}\) \(y'(0) = 1\)
14.
\(\dfrac{d^2 x}{dt^2} + 6 \dfrac{dx}{dt} + 5x = 0\text{,}\) \(x(0) = 2\text{,}\) \(x'(0) = -1\)
15.
\(x'' - 10x' + 25x = 0\text{,}\) \(x(0) = 1\text{,}\) \(x'(0) = 0\)
16.
\(\dfrac{d^2 y}{dx^2} - 2\dfrac{dy}{dx} + 4y = 0\text{,}\) \(x(0) = -1\text{,}\) \(x'(0) = 1\)
17.
\(y'' - 8y' + 4y = 0\text{,}\) \(y(0) = 1\text{,}\) \(y'(0) = -2\)
18.
\(\dfrac{d^2 x}{dx^2} + 3\dfrac{dx}{dt} - 10x = 0\text{,}\) \(x(0) = 1\text{,}\) \(x'(0) = 2\)
19.
\(\dfrac{d^2 Q}{dt^2} - 4 \dfrac{dQ}{dt} + 9Q = 0\text{,}\) \(Q(0) = -1\text{,}\) \(Q'(0) = 2\)
20.
\(\dfrac{d^2 y}{dt^2} + 6 \dfrac{dy}{dt} + 9y = 0\text{,}\) \(y(0) = 0\text{,}\) \(y'(0) = 0\)

Harmonic Oscillators.

Consider the harmonic oscillators with mass \(m\text{,}\) damping coeeficient \(b\text{,}\) and spring constant \(k\) in Exercise Group 4.1.6.21–28.
  1. Write the second-order initial value problem corresponding for the harmonic oscillator.
  2. Classify the oscillator as undamped, under-damped, over-damped, or critically damped.
  3. Solve the initial value problem.
  4. Sketch the \(x(t)\) and \(v(t)\)-graphs of the solution of the initial value problem.
  5. Sketch the phase portrait of the initial value problem.
21.
\(m = 1\text{,}\) \(b = 1\) \(k = 1\text{,}\) \(x(0) = 1\text{,}\) \(v(0) = 0\)
22.
\(m = 1\text{,}\) \(b = 2\) \(k = 3\text{,}\) \(x(0) = -3\text{,}\) \(v(0) = 4\)
23.
\(m = 1\text{,}\) \(b = 5\) \(k = 3\text{,}\) \(x(0) = 2\text{,}\) \(v(0) = -3\)
24.
\(m = 1\text{,}\) \(b = 0\) \(k = 25\text{,}\) \(x(0) = 2\text{,}\) \(v(0) = 0\)
25.
\(m = 2\text{,}\) \(b = 3\) \(k = 5\text{,}\) \(x(0) = 2\text{,}\) \(v(0) = -1\)
26.
\(m = 4\text{,}\) \(b = 4\) \(k = 1\text{,}\) \(x(0) = 2\text{,}\) \(v(0) = 1\)
27.
\(m = 3\text{,}\) \(b = 4\) \(k = 1\text{,}\) \(x(0) = 2\text{,}\) \(v(0) = 1\)
28.
\(m = 8\text{,}\) \(b = 4\) \(k = 1\text{,}\) \(x(0) = 2\text{,}\) \(v(0) = 1\)

Oscillations of a Hanging Mass.

In Exercise Group 4.1.6.29–31, we will consider the motion of a mass \(m\) hanging at the end of a vertical spring as in Figure 4.1.17. The mass stretches the spring in a downward (positive) direction by length \(L\text{.}\) There are two forces acting on the point where the mass is attached to the spring—the force exerted by the spring and gravity. The force of gravity, the weight of the mass, acts downward with a magnitude of \(mg\text{,}\) where \(g\) is the acceleration due to gravity. On the other hand, the force of the spring acts upward and is given by \(-kL\text{,}\) where \(L\) is the length of the spring.
 1 
This is commonly known as Hooke’s law. See Subsection 1.1.3.
When the mass is hanging in equilibrium the force of gravity and the force of the spring balance each other out; that is,
\begin{equation*} mg - kL = 0. \end{equation*}
three springs hanging from a beam with masses attacheds the right two springs and the rightmost spring extended
Figure 4.1.17. A spring-mass system
We would like to investigate the motion of the mass if is initially displaced or acted on by an external force. Let \(x(t)\) be the displacement of the mass from its equilibrium position, where a downward displacement is positive. The force acting on the mass are the weight of the mass \(mg\) and the force exerted by the spring, which is the total elongation of the spring, or
\begin{equation*} -k(L + x). \end{equation*}
By Newton’s second law of motion,
\begin{equation*} mx'' = mg - k(L + x). \end{equation*}
Furthermore, it is possible to add a damping term, \(bx'\text{,}\) or even an external force \(F(t)\) to obtain
\begin{equation*} mx'' = mg - k(L + x) + bx' + F(t). \end{equation*}
Since \(mg - kL = 0\text{,}\) we obtain the familiar equation
\begin{equation*} mx'' + bx' + kx = F(t). \end{equation*}
29.
Suppose that a mass of \(100\) grams stretches a spring \(2\) centimeters.
  1. Determine the spring constant \(k\text{.}\)
  2. If the mass is displaced an additional \(4\) centimeters and released, write an initial value problem that will model the motion of the oscillating mass.
  3. Solve the initial value problem.
30.
Suppose that a mass of \(1\) kilogram stretches a spring \(5\) centimeters.
  1. Determine the spring constant \(k\text{.}\)
  2. If the mass is displaced an additional \(5\) centimeters and released, write an initial value problem that will model the motion of the oscillating mass.
  3. Suppose the the spring-mass system is suspended in a fluid that exerts a resistance of \(0.25\) kilograms when the mass has a velocity of \(2\) centimeters per second. Modify the intial-value problem that you wrote in (b) to take this fact into account.
  4. Solve the initial value problem.
31.
Suppose that a mass weighing 4 lbs stretches a spring 3 inches.
  1. If \(g = 32 \text{ft/sec}^2\text{,}\) determine \(m\text{.}\)
  2. Determine the spring constant \(k\text{.}\)
  3. If the mass is displaced an additional 6 inches and released, write an initial value problem that will model the motion of the oscillating mass.
  4. Solve the initial value problem.
Hint.
Pay careful attention to units.

32.

Let \(a x'' + b x' + cx = 0\text{,}\) where \(a \neq 0\) and \(b^2 - 4ac = 0\text{.}\)
  1. Show that \(x_1(t) = e^{-bt/2a}\) is a solution to \(a x'' + b x' + cx = 0\text{.}\)
  2. Assume that
    \begin{equation*} y = v(t) x_1(t) = v(t) e^{-bt/2a} \end{equation*}
    is a solution to \(a x'' + b x' + cx = 0\) and show that \(v(t) = c_1 + c_2 t\text{.}\) Thus,
    \begin{equation*} x(t) = c_1 e^{-bt/2a} + c_2 t e^{-bt/2a} \end{equation*}
    is a general solution for \(a x'' + b x' + cx = 0\text{.}\)
Hint.
  1. Observe that
    \begin{align*} a x_1'' + b_1' + cx_1 & = a \left(\frac{-b}{2a}\right)^2e^{-bt/2a} + b \left( \frac{-b}{2a} \right) e^{-bt/2a} + c e^{-bt/2a}\\ & = e^{-bt/2a} \left( \frac{b^2}{4a} - \frac{b^2}{2a} + c \right)\\ & = e^{-bt/2a} \left( \frac{-b^2 + 4ac}{4a} \right)\\ & = 0. \end{align*}
  2. If \(y = v(t) x_1(t) = v(t) e^{-bt/2a}\) is a solution to our differential equation, then
    \begin{align*} a y'' + b y' + cy & = a (v''x_1 + 2 v'x_1' + vx_1'' ) + b(v' x_1 + v x_1') + cv x_1\\ & = a v''x_1 + 2a v'x_1' + bv' x_1 + v(a x_1'' +b x_1' + c x_1)\\ & = a v'' e^{-bt/2a} + \left[2a \left( \frac{-b}{2a} \right) e^{-bt/2a} + b e^{-bt/2a} \right] v'\\ & = a v'' e^{-bt/2a}\\ & = 0. \end{align*}
    Since \(a \neq 0\text{,}\) we know that \(v'' = 0\text{.}\) Hence, \(v(t) = c_1 + c_2 t\text{.}\)

33.

Suppose that \(x_1(t)\) is a solution (not identically zero) to the equation
\begin{equation*} x'' + p(t) x' + q(t) x = 0. \end{equation*}
  1. Assume that \(x(t) = v(t) x_1(t)\) is a solution to \(x'' + p(t) x' + q(t) x = 0\) and derive the equation
    \begin{equation} x_1 v'' +(2x_1' + px_1)v' = 0.\tag{4.1.9} \end{equation}
  2. Let \(u = v'\) and show that (4.1.9) is a first-order linear differential equation in \(u\text{.}\)
  3. Show that \(x_1(t) = 1/t\) is a solution to
    \begin{equation} 2 t^2 x'' + 3t x' - x = 0 \tag{4.1.10} \end{equation}
    for \(t \gt 0\) and find a second linearly independent solution to (4.1.10).
Hint.
  1. \begin{align*} x'' + px' + qx & = (v''x_1 + 2 v' x_1' + vx_1'') + p(v'x_1 + vx_1') + q(vx_1)\\ & = x_1 v'' + 2v' x_1' + p x_1 v' + v (x_1'' + p x_1' + q x_1)\\ & = x_1 v'' +(2x_1' + px_1)v' \\ & = 0. \end{align*}
  2. If \(u = v'\text{,}\) then \(x_1 u' +(2x_1' + px_1)u= 0\text{.}\)
  3. If \(x_1(t) = 1/t\text{,}\) then
    \begin{equation*} 2 t^2 x_1'' + 3t x_1' - x_1 = 2 t^2 \left(\frac{2}{t^3}\right) + 3t \left(\frac{-1}{t^2}\right) - \frac{1}{t} = 0. \end{equation*}
    If we assume that \(x = v/t\) is a second solution, then
    \begin{equation*} 2 t^2 x'' + 3t x' - x = 2tv'' - v' = 0. \end{equation*}
    If we let \(u = v'\text{,}\) then a solution of \(2tu' - u = 0\) is \(u = \sqrt{t}\) and \(v = \int \sqrt{t} \, dt = 2 t^{3/2} / 3\text{.}\) Therefore, the second solution to our equation is
    \begin{equation*} x = \frac{v}{t} = \frac{2}{3} \sqrt{t}. \end{equation*}

34.

Let \(a x'' + b x' + cx = 0\text{,}\) where \(a \neq 0\) and \(b^2 - 4ac = 0\text{.}\)
  1. Show that \(x_1(t) = e^{-bt/2a}\) is a solution to \(a x'' + b x' + cx = 0\text{.}\)
  2. Assume that
    \begin{equation*} y = v(t) x_1(t) = v(t) e^{-bt/2a} \end{equation*}
    is a solution to \(a x'' + b x' + cx = 0\) and show that \(v(t) = c_1 + c_2 t\text{.}\) Thus,
    \begin{equation*} x(t) = c_1 e^{-bt/2a} + c_2 t e^{-bt/2a} \end{equation*}
    is a general solution for \(a x'' + b x' + cx = 0\text{.}\)

35.

Consider the equation
\begin{equation*} y'' - (2 \alpha -1) y' + \alpha (\alpha - 1) y =0. \end{equation*}
Determine all values of \(\alpha\text{,}\) if any, for which all solutions tend toward zero as \(t \to \infty\text{.}\) Also, determine the values of \(\alpha\text{,}\) if any, for which all nonzero solutions become unbounded as \(t \to \infty\text{.}\)

36.

Solve each of the following initial value problems.
  1. \begin{align*} y'' - 2y' + 5y & = 0\\ y(\pi/2) & = 0\\ y'(\pi/2) & = 2 \end{align*}
  2. \begin{align*} 9y'' - 12y' + 4y & = 0\\ y(0) & = 2\\ y'(0) & = -1 \end{align*}
  3. \begin{align*} y'' + 8y' -9y & = 0\\ y(1) & = 1\\ y'(1) & = 0 \end{align*}
  4. \begin{align*} y'' + 2ay' + (a^2 + 1)y & = 0\\ y(0) & = 1\\ y'(0) & = 0 \end{align*}

37. Reduction of Order.

Suppose that \(x_1(t)\) is a solution (not identically zero) to the equation
\begin{equation*} x'' + p(t) x' + q(t) x = 0. \end{equation*}
  1. Assume that \(x(t) = v(t) x_1(t)\) is a solution to \(x'' + p(t) x' + q(t) x = 0\) and derive the equation
    \begin{equation} x_1 v'' +(2x_1' + px_1)v' = 0.\tag{4.1.11} \end{equation}
  2. Let \(u = v'\) and show that ((4.1.11)) is a first-order linear differential equation in \(u\text{.}\)
  3. Show that \(x_1(t) = 1/t\) is a solution to
    \begin{equation} 2 t^2 x'' + 3t x' - x = 0\tag{4.1.12} \end{equation}
    for \(t \gt 0\) and find a second linearly independent solution to ((4.1.12)).

38. Euler Equations.

An important class of second-order linear differential equations are equations of the form
\begin{equation*} t^2 y'' + \alpha t y' + \beta y = 0, \end{equation*}
where \(t \gt 0\) and \(\alpha\) and \(\beta\) are real constants. An equation of this form is called an Euler equation.
  1. Show that the substitution \(x = \ln t\) transforms an Euler equation into an equation of constant coefficients.
  2. Solve the equation
    \begin{equation*} t^2 y'' + 4ty' + 2y = 0. \end{equation*}
Hint.
Show that
\begin{align*} \frac{dy}{dt} & = \frac{dx}{dt} \frac{dy}{dx}\\ \frac{d^2y}{dt^2} & = \left( \frac{dx}{dt} \right)^2 \frac{d^2y}{dx^2} + \frac{d^2x}{dt^2} \frac{dy}{dx}. \end{align*}

39. Higher Order Linear Equqtions with Constant Coefficients.

The characteristic equation of the \(n\)th order linear differential equation
\begin{equation*} y^{(n)} + a_{n - 1} y^{(n-1)} + \cdots a_1 y' + a_0 y = 0 \end{equation*}
is
\begin{equation*} \lambda^n + a_{n - 1} \lambda^{n-1} + \cdots a_1 \lambda + a_0 = 0. \end{equation*}
For example, the characteristic equation of
\begin{equation*} \frac{d^3 y }{dt^3} - 6 \frac{d^2 y}{dt^2} + 11 \frac{dy}{dt} - 6y = 0 \end{equation*}
is
\begin{equation*} \lambda^3 - 6 \lambda^2 + 11 \lambda - 6 = 0.\text{.} \end{equation*}
This equation factors as \((\lambda - 1)(\lambda - 2)(\lambda - 3) = 0\text{.}\) Since the roots of the characteristic equation are \(\lambda_1 = 1\text{,}\) \(\lambda_2 = 2\text{,}\) and \(\lambda_3 = 3\text{.}\) As with second-order linear differntial equations, we can immediately write down our solution as
\begin{equation*} y = c_1 e^t + c_2 e^{2t} + c_3 e^{3t}. \end{equation*}
  1. Find the general solution of the equation \(y^{(4)} - 9 y'' + 20y = 0\text{.}\)
  2. Find the general solution of the equation \(y''' - 6 y'' + 2y' + 36y = 0\text{.}\)
  3. Find the general solution of the equation \(y''' - y'' + y' - y = 0\text{.}\)
In general, we can always factor the characteristic polynomial over the complex numbers; however, this may be very difficult to do in practice. We are often forced to numerical techniques, especially for systems of very high order.
Hint.
Show that
\begin{align*} \frac{dy}{dt} & = \frac{dx}{dt} \frac{dy}{dx}\\ \frac{d^2y}{dt^2} & = \left( \frac{dx}{dt} \right)^2 \frac{d^2y}{dx^2} + \frac{d^2x}{dt^2} \frac{dy}{dx}. \end{align*}

Subsection 4.1.7 Solving Second Order Linear Equations with Sage

Second order homogeneous linear differential equations with constant coefficeints can be solved sybolically using Sage. For example,
\begin{equation*} x'' + 2 x' + x = 0 \end{equation*}
can be solved using the following Sage commands.
We can even solve initial value problems such as
\begin{align*} x'' + 2 x' + x \amp = 0\\ x(0) \amp = 1\\ x'(0) \amp = 2 \end{align*}
Of course, we can add a forcing term to our initial value problem,
\begin{align*} x'' + 2 x' + x \amp = t \cos 2t\\ x(0) \amp = 1\\ x'(0) \amp = 2 \end{align*}

Exercises 4.1.8 Exercises

Solve each of the following differential equations using Sage.

1.

\(\dfrac{d^2 y}{dx^2} - y = 0\)

2.

\(\dfrac{d^2 y}{dx^2} - y = 0\text{,}\) \(y(0) = 1\text{,}\) \(y'(0) = 0\)

3.

\(\dfrac{d^2 y}{dx^2} - y = \cos 2t\)

4.

\(\dfrac{d^2 y}{dx^2} - y = \cos 2t\text{,}\) \(y(0) = 1\text{,}\) \(y'(0) = 0\)

5.

\(\dfrac{d^2 y}{dx^2} - y = \cos t\)

6.

\(\dfrac{d^2 y}{dx^2} - y = \cos t\text{,}\) \(y(0) = 1\text{,}\) \(y'(0) = 0\)

7.

\(x'' - 2x' - 8x = 0\)

8.

\(x'' - 2x' - 8x = 0\text{,}\) \(x(0) = 1\text{,}\) \(x'(0) = 2\)

9.

\(x'' - 2x' - 8x = e^{-2t}\)

10.

\(x'' - 2x' - 8x = e^{-2t} + e^{4t}\)