Second-Order Differential Equation

Interactive DE Solver & Plotter

Solve $ay'' + by' + cy = 0$

Initial Conditions (Optional)

Understanding (Second-Order) 2nd Order Differential Equations

A second-order linear homogeneous differential equation with constant coefficients is a cornerstone in the study of differential equations. It takes the general form:

$$ a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = 0 $$

Alternatively, using prime notation ($y' = \frac{dy}{dx}$, $y'' = \frac{d^2y}{dx^2}$):

$$ ay'' + by' + cy = 0 $$

Key characteristics of this equation type:

  • 2nd Order (Second Order): The highest derivative present is the second derivative ($y''$).
  • Linear: The dependent variable $y$ and its derivatives ($y', y''$) appear only to the first power and are not multiplied together.
  • Homogeneous: The right-hand side of the equation is zero. If it were a non-zero function of $x$, say $f(x)$, the equation would be non-homogeneous.
  • Constant Coefficients: The coefficients $a, b,$ and $c$ are real numbers (constants), with the crucial condition that $a \neq 0$ (otherwise, it wouldn't be a second-order equation).

These equations are incredibly powerful for modeling a wide array of physical phenomena. For instance, they describe the motion of a mass attached to a spring (simple harmonic motion, damped oscillations), the flow of charge in an RLC electrical circuit, and many other systems that exhibit oscillatory or exponential decay/growth behavior.

The Auxiliary (Characteristic) Equation Method

The standard approach to solving $ay'' + by' + cy = 0$ involves a clever substitution. We propose a solution of the exponential form $y = e^{mx}$, where $m$ is a constant to be determined.

Step 1: Assume a Solution Form

Let $y = e^{mx}$. Then, the derivatives are:

  • $y' = \frac{dy}{dx} = me^{mx}$
  • $y'' = \frac{d^2y}{dx^2} = m^2e^{mx}$

Step 2: Substitute into the DE

Plugging these into the original differential equation $ay'' + by' + cy = 0$ gives:

$$ a(m^2e^{mx}) + b(me^{mx}) + c(e^{mx}) = 0 $$

We can factor out $e^{mx}$ (which is never zero for any real $m$ and $x$):

$$ e^{mx}(am^2 + bm + c) = 0 $$

Step 3: Form the Auxiliary Equation

Since $e^{mx} \neq 0$, for the equation to hold, the term in the parenthesis must be zero. This gives us the auxiliary equation (also known as the characteristic equation):

$$ am^2 + bm + c = 0 $$

This is a simple quadratic equation in $m$. The roots of this auxiliary equation dictate the form of the general solution to the differential equation.

Step 4: Find the Roots of the Auxiliary Equation

We use the quadratic formula to find the roots $m_1$ and $m_2$:

$$ m = \frac{-b \pm \sqrt{b^2 - 4ac}}{2a} $$

The nature of these roots depends on the discriminant, $\Delta = b^2 - 4ac$. There are three distinct cases:

Case 1: Real and Distinct Roots ($\Delta > 0$)

If $b^2 - 4ac > 0$, the auxiliary equation has two distinct real roots, $m_1$ and $m_2$. In this scenario, $e^{m_1x}$ and $e^{m_2x}$ are two linearly independent solutions. The general solution is a linear combination of these:

$$ y(x) = Ae^{m_1x} + Be^{m_2x} $$

Where $A$ and $B$ are arbitrary constants. These constants are typically determined if initial conditions (e.g., $y(x_0) = y_0$ and $y'(x_0) = y'_0$) are provided.

Case 2: Real and Equal Roots ($\Delta = 0$)

If $b^2 - 4ac = 0$, the auxiliary equation has one real root (a repeated root), $m_1 = m_2 = m = -\frac{b}{2a}$. One solution is $y_1 = e^{mx}$. A second, linearly independent solution can be found (e.g., using the method of reduction of order) to be $y_2 = xe^{mx}$. The general solution is then:

$$ y(x) = (A + Bx)e^{mx} \quad \text{or} \quad y(x) = Ae^{mx} + Bxe^{mx} $$

$A$ and $B$ are arbitrary constants.

Case 3: Complex Conjugate Roots ($\Delta < 0$)

If $b^2 - 4ac < 0$, the auxiliary equation has two complex conjugate roots. These roots can be written as $m_1 = \alpha + i\beta$ and $m_2 = \alpha - i\beta$, where:

  • $\alpha = -\frac{b}{2a}$ (the real part)
  • $\beta = \frac{\sqrt{4ac - b^2}}{2a}$ (the imaginary part, $\beta > 0$)

The two linearly independent solutions would initially appear as $e^{(\alpha+i\beta)x}$ and $e^{(\alpha-i\beta)x}$. Using Euler's formula ($e^{i\theta} = \cos\theta + i\sin\theta$), these can be rewritten in terms of real-valued functions. The general solution becomes:

$$ y(x) = e^{\alpha x}(A\cos(\beta x) + B\sin(\beta x)) $$

$A$ and $B$ are arbitrary constants. This form of the solution often represents damped oscillations in physical systems.

Step 5: Apply Initial Conditions (if given)

If initial conditions such as $y(x_0) = y_0$ and $y'(x_0) = y'_0$ are provided, they are used to solve for the constants $A$ and $B$ in the general solution, yielding a unique particular solution.