Welcome to the exploration of Initial Value Problems (IVPs) for second-order Ordinary Differential Equations (ODEs). If you’ve been following the progression of this topic, in the previous article, we explored the general solutions of homogeneous and non-homogeneous second-order linear equations with constant coefficients. This laid the foundation for understanding more complex differential equations, focusing on finding particular and complementary solutions depending on the nature of the equation. Now, we turn our attention to second-order ODEs, approaching them through IVPs.

1. Understanding Second-Order ODEs and Initial Value Problems

In an IVP, the goal is to find the specific solution to the ODE that satisfies not only the equation itself but also given initial conditions. For second-order ODEs, two initial conditions are typically specified, as the solution involves both the function and its derivative. An IVP for a second-order ODE is typically expressed as:

\frac{d^2y}{dx^2} + p\frac{dy}{dx} + qy = g(x), \quad y(x_0) = y_0, \quad \frac{dy}{dx}(x_0) = y'

Here:

- $y(x_0) = y_0$ is the initial value of the function at $x = x_0$ .

- $\frac{dy}{dx}(x_0) = y'$ is the initial value of the derivative.

1.1. Example: IVP for a Damped Harmonic Oscillator

Consider the equation for a damped harmonic oscillator:

\frac{d^2y}{dt^2} + b\frac{dy}{dt} + ky = 0

Here, $y$ represents displacement, $t$ represents time, and the constants $b$ and $k$ represent the damping coefficient and stiffness, respectively. Suppose the initial conditions are given by:

y(0) = 1, \quad \frac{dy}{dt}(0) = 0

To solve this IVP, we first solve the homogeneous ODE. Assuming a solution of the form $y(x) = e^{mt}$ , we derive the characteristic equation:

m^2 + bm + k = 0

Solving for $m$ , we determine the general solution based on the roots of the characteristic equation. Then, applying the initial conditions allows us to find the constants, yielding the specific solution to the IVP.

2. Solving Second-Order ODE IVPs

Step 1: Find the general solution of the ODE.

The general solution of a second-order non-homogeneous linear differential equation with constant coefficients can be expressed as:

y'' + a_1y' + a_2y = f(x)

y = y_c + y_p

where $y_p$ is any specific function that satisfies the non-homogeneous equation, and $y_c = C_1y_1 + C_2y_2$ is the general solution of the corresponding homogeneous equation: $y'' + a_1y' + a_2y = 0.$

To see how to obtain the complementary solution ( $y_c$ ) and the particular solution ( $y_p$ ), refer to the previous articles where these methods were discussed in detail.

Step 2: Apply Initial Conditions

To determine the specific solution for the IVP, use the given initial conditions to solve for the constants $C_1$ and $C_2$ . These constants ensure the solution satisfies both the differential equation and the initial conditions.

2.1. Examples: Homogeneous Second-Order IVPs

2.1.1. Example 1

Consider the following second-order ODE with initial conditions:

\frac{d^2y}{dx^2} - 3\frac{dy}{dx} + 2y = 0, \quad y(0) = 2, \quad \frac{dy}{dx}(0) = 0

Step 1: Find the general solution

The characteristic equation is:

m^2 - 3m + 2 = 0

Factoring gives:

(m - 1)(m - 2) = 0

So, the roots are $m_1 = 1$ and $m_2 = 2$ . The general solution is:

y(x) = C_1e^x + C_2e^{2x}

Step 2: Apply initial conditions

Using the initial condition $y(0) = 2$ :

C_1 + C_2 = 2 \quad \text{(Equation 1)}

Next, differentiate the general solution:

\frac{dy}{dx} = C_1e^x + 2C_2e^{2x}

Applying the initial condition $\frac{dy}{dx}(0) = 0$ :

C_1 + 2C_2 = 0 \quad \text{(Equation 2)}

Solving the system of equations (Equation 1 and Equation 2), we find:

C_1 = 4, \quad C_2 = -2

Therefore, the specific solution to the IVP is:

y(x) = 4e^x - 2e^{2x}

2.1.2. Example 2

\frac{d^2y}{dx^2} - 5\frac{dy}{dx} + 6y = 0, \quad y(0) = 1, \quad y'(0) = 0

Solution:

1. Solve the characteristic equation:

m^2 - 5m + 6 = 0 \quad \Rightarrow \quad (m - 2)(m - 3) = 0

Roots: $m_1 = 2$ , $m_2 = 3$ .

2. General solution:

y(x) = C_1 e^{2x} + C_2 e^{3x}

3. Apply initial conditions:

- $y(0) = C_1 + C_2 = 1$

- $y'(x) = 2C_1 e^{2x} + 3C_2 e^{3x}$ , and $y'(0) = 2C_1 + 3C_2 = 0$

4. Solving for constants:

C_1 + C_2 = 1, \quad 2C_1 + 3C_2 = 0 \quad \Rightarrow \quad C_1 = 3, \quad C_2 = -2

Final solution:

y(x) = 3 e^{2x} - 2 e^{3x}

2.1.3. Example 3

\frac{d^2y}{dx^2} + 4y = 0, \quad y(0) = 2, \quad y'(0) = -3

Solution:

1. Solve the characteristic equation:

m^2 + 4 = 0 \quad \Rightarrow \quad m = \pm 2i

2. General solution:

y(x) = C_1 \cos(2x) + C_2 \sin(2x)

3. Apply initial conditions:

- $y(0) = C_1 = 2$

- $y'(x) = -2C_1 \sin(2x) + 2C_2 \cos(2x)$ , and $y'(0) = 2C_2 = -3$

4. Solving for constants:

C_1 = 2, \quad C_2 = -\frac{3}{2}

Final solution:

y(x) = 2 \cos(2x) - \frac{3}{2} \sin(2x)

2.1.4. Example 4

\frac{d^2y}{dx^2} - 4y = 0, \quad y(0) = 1, \quad y'(0) = 1

Solution:

1. Solve the characteristic equation:

m^2 - 4 = 0 \quad \Rightarrow \quad m = \pm 2

2. General solution:

y(x) = C_1 e^{2x} + C_2 e^{-2x}

3. Apply initial conditions:

- $y(0) = C_1 + C_2 = 1$

- $y'(x) = 2C_1 e^{2x} - 2C_2 e^{-2x}$ , and $y'(0) = 2C_1 - 2C_2 = 1$

4. Solve for constants:

C_1 + C_2 = 1, \quad 2C_1 - 2C_2 = 1 \quad \Rightarrow \quad C_1 = \frac{3}{4}, \quad C_2 = \frac{1}{4}

Final solution:

y(x) = \frac{3}{4} e^{2x} + \frac{1}{4} e^{-2x}

2.1.5. Example 5

\frac{d^2y}{dx^2} + 9y = 0, \quad y(0) = 0, \quad y'(0) = 1

Solution:

1. Solve the characteristic equation:

m^2 + 9 = 0 \quad \Rightarrow \quad m = \pm 3i

2. General solution:

y(x) = C_1 \cos(3x) + C_2 \sin(3x)

3. Apply initial conditions:

- $y(0) = C_1 = 0$

- $y'(x) = -3C_1 \sin(3x) + 3C_2 \cos(3x)$ , and $y'(0) = 3C_2 = 1$

4. Solving for constants:

C_1 = 0, \quad C_2 = \frac{1}{3}

Final solution:

y(x) = \frac{1}{3} \sin(3x)

2.1.6. Example 6

\frac{d^2y}{dx^2} - 2\frac{dy}{dx} + y = 0, \quad y(0) = 1, \quad y'(0) = 0

Solution:

1. Solve the characteristic equation:

m^2 - 2m + 1 = 0 \quad \Rightarrow \quad (m - 1)^2 = 0 \quad \Rightarrow \quad m = 1

2. General solution:

y(x) = (C_1 + C_2 x)e^{x}

3. Apply initial conditions:

- $y(0) = C_1 = 1$

- $y'(x) = (C_2 + C_1 + C_2 x)e^{x}$ , and $y'(0) = C_2 + C_1 = 0$

4. Solving for constants:

C_1 = 1, \quad C_2 = -1

Final solution:

y(x) = (1 - x)e^{x}

2.2. Examples: Non-Homogeneous Second-Order IVPs

2.2.1. Example 7

\frac{d^2y}{dx^2} - 3\frac{dy}{dx} + 2y = e^{x}, \quad y(0) = 0, \quad y'(0) = 1

Solution:

1. Solve the Homogeneous Equation

- The characteristic equation for the homogeneous part is:

m^2 - 3m + 2 = 0

- Factoring:

(m - 1)(m - 2) = 0

- The roots are $m_1 = 1$ , $m_2 = 2$ .

2. General Solution of the Homogeneous Equation

y_c(x) = C_1 e^{x} + C_2 e^{2x}

3. Form the Particular Solution Based on the Forcing Function:

- The forcing function is $f(x) = e^{x}$ .

- Since $e^{x}$ is already a solution to the homogeneous equation (appearing in $y_c(x)$ ), we multiply by $x$ to avoid overlap. The particular solution must take the form:

y_p(x) = A x e^{x}

4. Compute Derivatives of $y_p(x)$ :

y_p'(x) = A(e^{x} + x e^{x}) = A e^{x}(1 + x)

y_p''(x) = A(2 e^{x} + x e^{x}) = A e^{x}(2 + x)

5. Substitute $y_p(x)$ into the Differential Equation

y_p'' - 3y_p' + 2y_p = e^{x}

gives;

A e^{x}(2 + x) - 3A e^{x}(1 + x) + 2A x e^{x} = e^{x}

Simplifying:

A e^{x}(2 + x - 3(1 + x) + 2x) = e^{x}

Further simplification leads to:

A e^{x}(2 + x - 3 - 3x + 2x) = e^{x} \quad \Rightarrow \quad A e^{x}(-1) = e^{x}

Solving for $A$ :

A = -1

6. General Solution

The general solution to the differential equation is the sum of the complementary solution and the particular solution:

y(x) = y_c(x) + y_p(x)

Thus:

y(x) = C_1 e^{x} + C_2 e^{2x} - x e^{x}

7. Apply Initial Conditions

- At $x = 0$ , we have:

y(0) = C_1 + C_2 = 0 \quad \Rightarrow \quad C_1 = -C_2

- The first derivative is:

y'(x) = C_1 e^{x} + 2C_2 e^{2x} - e^{x}(1 + x)

At $x = 0$ , we have:

y'(0) = C_1 + 2C_2 - 1 = 1

Substituting $C_1 = -C_2$ gives:

-C_2 + 2C_2 - 1 = 1 \quad \Rightarrow \quad C_2 = 2 \quad \text{and} \quad C_1 = -2

8. Final Solution:

Substituting $C_1$ and $C_2$ into the general solution:

y(x) = -2 e^{x} + 2 e^{2x} - x e^{x}

2.2.2. Example 8

\frac{d^2y}{dx^2} + y = \sin(x), \quad y(0) = 0, \quad y'(0) = 1

Step 1: Solve the Homogeneous Equation

\frac{d^2y}{dx^2} + y = 0

The characteristic equation is:

m^2 + 1 = 0

Solving for $m$ , we get:

m = \pm i

Thus, the solution to the homogeneous equation is:

y_c(x) = C_1 \cos(x) + C_2 \sin(x)

Step 2: Find the Particular Solution

The non-homogeneous term is $\sin(x)$ . To find a particular solution, we use the method of undetermined coefficients. Since the non-homogeneous term is $\sin(x)$ , and the homogeneous solution already contains terms involving $\sin(x)$ and $\cos(x)$ , we try a particular solution of the form:

y_p(x) = A x \cos(x) + B x \sin(x)

Now, calculate the first and second derivatives of $y_p(x)$ :

y_p'(x) = A (\cos(x) - x \sin(x)) + B (\sin(x) + x \cos(x))

y_p''(x) = A (-2 \sin(x) - x \cos(x)) + B (2 \cos(x) - x \sin(x))

Substitute $y_p(x)$ , $y_p'(x)$ , and $y_p''(x)$ into the original equation:

\left(A (-2 \sin(x) - x \cos(x)) + B (2 \cos(x) - x \sin(x))\right) + (A x \cos(x) + B x \sin(x)) = \sin(x)

Simplify the equation:

A(-2 \sin(x)) + B(2 \cos(x)) = \sin(x)

Equating coefficients of $\sin(x)$ and $\cos(x)$ , we get:

-2A = 1 \quad \Rightarrow \quad A = -\frac{1}{2}

2B = 0 \quad \Rightarrow \quad B = 0

Thus, the particular solution is:

y_p(x) = -\frac{1}{2} x \cos(x)

Step 3: General Solution

The general solution to the differential equation is the sum of the homogeneous and particular solutions:

y(x) = C_1 \cos(x) + C_2 \sin(x) - \frac{1}{2} x \cos(x)

Step 4: Apply Initial Conditions

We are given the initial conditions $y(0) = 0$ and $y'(0) = 1$ .

- At $x = 0$ , $y(0) = 0$ :

C_1 \cos(0) + C_2 \sin(0) - \frac{1}{2} \cdot 0 \cdot \cos(0) = 0

C_1 = 0

- Now, calculate $y'(x)$ :

y'(x) = C_1(-\sin(x)) + C_2 \cos(x) - \frac{d}{dx}\left(\frac{1}{2} x \cos(x)\right)

y'(x) = C_2 \cos(x) - \frac{1}{2} \cos(x) + \frac{1}{2} x \sin(x)

At $x = 0$ , $y'(0) = 1$ :

C_2 \cos(0) - \frac{1}{2} \cos(0) = 1

C_2 - \frac{1}{2} = 1 \quad \Rightarrow \quad C_2 = \frac{3}{2}

Step 5: Final Solution

Thus, the final solution to the differential equation is:

y(x) = \frac{3}{2} \sin(x) - \frac{1}{2} x \cos(x)

2.2.3. Example 9

\frac{d^2y}{dx^2} + 2\frac{dy}{dx} + y = x^2, \quad y(0) = 1, \quad y'(0) = 0

Solution:

1. Solve the Homogeneous Equation

- The characteristic equation for the homogeneous part is:

m^2 + 2m + 1 = 0

- The root is $m = -1$ (a repeated root).

2. General Solution of the Homogeneous Equation:

y_c(x) = (C_1 + C_2 x)e^{-x}

3. Form the Particular Solution Based on the Forcing Function

- The forcing function is $f(x) = x^2$ .

- Since $f(x)$ is a polynomial, the particular solution must also be a polynomial of degree 2. Thus, the particular solution takes the form:

y_p(x) = A x^2 + B x + C

4. Compute Derivatives of $y_p(x)$

y_p'(x) = 2A x + B

y_p''(x) = 2A

5. Substitute $y_p(x)$ into the Differential Equation

Substituting into the original equation:

y_p'' + 2 y_p' + y_p = x^2

gives:

2A + 2(2A x + B) + A x^2 + B x + C = x^2

Simplifying:

(A) x^2 + (4A + B) x + (2A + 2B + C) = x^2

Equating coefficients:

A = 1, \quad 4A + B = 0, \quad 2A + 2B + C = 0

Solving these equations:

A = 1, \quad B = -4, \quad C = 6

6. General Solution

The general solution to the differential equation is the sum of the complementary solution and the particular solution:

y(x) = (C_1 + C_2 x)e^{-x} + x^2 - 4x + 6

7. Apply Initial Conditions

- At $x = 0$ , we have:

y(0) = C_1 + 6 = 1 \quad \Rightarrow \quad C_1 = -5

- The first derivative is:

y'(x) = (-C_1 e^{-x} - C_2 e^{-x} + C_2 x e^{-x}) + 2x - 4

At $x = 0$ , we have:

y'(0) = -C_1 - C_2 - 4 = 0 \quad \Rightarrow \quad C_2 = -1

8. Final Solution

Substituting $C_1$ and $C_2$ into the general solution:

y(x) = (-5 - x) e^{-x} + x^2 - 4x + 6

2.2.4. Example 10

\frac{d^2y}{dx^2} - 3\frac{dy}{dx} + 2y = \cos(x), \quad y(0) = 2, \quad y'(0) = 0

1. Solve the Homogeneous Equation:

\frac{d^2y}{dx^2} - 3\frac{dy}{dx} + 2y = 0

The characteristic equation is:

m^2 - 3m + 2 = 0

Factoring this:

(m - 1)(m - 2) = 0

The roots are:

m_1 = 1, \quad m_2 = 2

Thus, the general solution to the homogeneous equation is:

y_c(x) = C_1 e^{x} + C_2 e^{2x}

2. Form the Particular Solution:

The non-homogeneous term is $\cos(x)$ . Since this is a trigonometric function, we propose a particular solution of the form:

y_p(x) = A \cos(x) + B \sin(x)

3. Compute Derivatives of $y_p(x)$ :

y_p'(x) = -A \sin(x) + B \cos(x)

y_p''(x) = -A \cos(x) - B \sin(x)

4. Substitute into the Differential Equation:

Substitute $y_p(x)$ and its derivatives into the original differential equation:

y_p'' - 3 y_p' + 2 y_p = \cos(x)

Substituting:

(-A \cos(x) - B \sin(x)) - 3(-A \sin(x) + B \cos(x)) + 2(A \cos(x) + B \sin(x)) = \cos(x)

Simplifying the terms:

(-A - 3B + 2A) \cos(x) + (-B + 3A + 2B) \sin(x) = \cos(x)

5. Solve for $A$ and $B$ :

By comparing coefficients of $\cos(x)$ and $\sin(x)$ , we get the system of equations:

A - 3B = 1 \quad \text{(coefficient of } \cos(x) \text{)}

3A + B = 0 \quad \text{(coefficient of } \sin(x) \text{)}

Solving these two equations:

A= \frac{1}{10}

B = \frac{-3}{10}

6. General Solution:

The general solution to the differential equation is the sum of the complementary solution and the particular solution:

y(x) = y_c(x) + y_p(x) = C_1 e^{x} + C_2 e^{2x} + \frac{1}{10} \cos(x) + \frac{-3}{10} \sin(x)

7. Apply Initial Conditions

- At $x = 0$ , we have $y(0) = 2$ :

y(0) = C_1 e^{0} + C_2 e^{0} + \frac{1}{10} \cos(0) - \frac{3}{10} \sin(0) = C_1 + C_2 + \frac{1}{10} = 2

Simplifying:

C_1 + C_2 = 2 - \frac{1}{10} = \frac{19}{10}

- The first derivative is:

y'(x) = C_1 e^{x} + 2C_2 e^{2x} - \frac{1}{10} \sin(x) - \frac{3}{10} \cos(x)

At $x = 0$ , we have $y'(0) = 0$ :

y'(0) = C_1 + 2C_2 - \frac{3}{10} = 0

Simplifying:

C_1 + 2C_2 = \frac{3}{10}

8. Solve for $C_1$ and $C_2$ :

We now have the system of equations:

$C_1 + C_2 = \frac{19}{10}$
$C_1 + 2C_2 = \frac{3}{10}$

Solving the above system of equations, we get $C_1=72, C_2=-85$

9. Final Solution:

Substituting $C_1 = \frac{7}{2}$ and $C_2 = -\frac{8}{5}$ into the general solution:

y(x) = \frac{7}{2} e^{x} - \frac{8}{5} e^{2x} + \frac{1}{10} \cos(x) - \frac{3}{10} \sin(x)

3. Practice Questions

3.1. Practice Problem Set 1

Solve the following Homogeneous Second-Order IVPs.

$\frac{d^2y}{dx^2} - 5\frac{dy}{dx} + 6y = 0, \quad y(0) = 4, \quad y'(0) = 0$
$\frac{d^2y}{dx^2} - 4\frac{dy}{dx} + 4y = 0, \quad y(0) = 1, \quad y'(0) = 3$
$\frac{d^2y}{dx^2} + 4y = 0, \quad y(0) = 0, \quad y'(0) = 2$

Not sure about your answers? Don't worry, click the link below to check the video where answers are explained.

Practice Questions Set 1 - Discussion

3.2. Practice Problem Set 2

Solve the following Homogeneous Second-Order IVPs.

$\frac{d^2y}{dx^2} - 3\frac{dy}{dx} + 2y = x, \quad y(0) = 1, \quad y'(0) = 0$
$\frac{d^2y}{dx^2} - 4\frac{dy}{dx} + 4y = 3x, \quad y(0) = 1, \quad y'(0) = 0$
$\frac{d^2y}{dx^2} + 5y = e^{-x}, \quad y(0) = 1, \quad y'(0) = 0$

Not sure about your answers? Don't worry, click the link below to check the video where answers are explained.

Practice Questions Set 2 - Discussion

Initial value problems are essential for modeling physical phenomena where the behavior of a system at a starting point dictates its future evolution. In applications such as mechanical vibrations, electrical circuits, and population dynamics, solving an IVP provides a clear understanding of the system's dynamics over time.

Stay tuned for the next sections of ODEs.

Related Articles

16. Second-order Constant Coefficient IVPs

1. Understanding Second-Order ODEs and Initial Value Problems

1.1. Example: IVP for a Damped Harmonic Oscillator

2. Solving Second-Order ODE IVPs

2.1. Examples: Homogeneous Second-Order IVPs

2.1.1. Example 1

2.1.2. Example 2

2.1.3. Example 3

2.1.4. Example 4

2.1.5. Example 5

2.1.6. Example 6

2.2. Examples: Non-Homogeneous Second-Order IVPs

2.2.1. Example 7

2.2.2. Example 8

2.2.3. Example 9

2.2.4. Example 10

3. Practice Questions

3.1. Practice Problem Set 1

3.2. Practice Problem Set 2