Solving IVPs via Laplace — Solving Linear ODE IVPs Using Laplace Transforms - Differential Equations

SECTION 1

Historical Context & Motivation

Differential equations describe how quantities change over time — from the swing of a pendulum to the charge in an electrical circuit. For centuries, mathematicians solved these equations by guessing solution forms and checking them, a process that could be tedious and error-prone. The Laplace transform changed all of that by converting calculus problems into algebra problems, making solutions far more systematic.

An initial value problem (IVP) is a differential equation paired with specific starting conditions — for instance, knowing where an object is and how fast it is moving at time t = 0. The Laplace transform method provides a recipe-like procedure: transform the equation, solve for the unknown in the new domain, and then transform back. This approach is especially powerful for engineers and scientists who need reliable, repeatable methods, and it is a standard technique in college-level differential equations courses, building directly on the calculus and algebra you have already studied.

💡

Prerequisites

This lesson assumes you have completed a calculus sequence (differentiation and integration, including improper integrals) and are comfortable with algebraic techniques such as factoring polynomials and manipulating rational expressions. These are the standard prerequisites for a first course in differential equations, and the worked examples build directly on that foundation.

1737 • Euler's Integral Ideas

Leonhard Euler explored integral transforms in his work on differential equations, laying groundwork for later transform methods.

1782 • Laplace Introduces the Transform

Pierre-Simon Laplace developed the integral transform that bears his name while studying probability theory and celestial mechanics.

1816 • Poisson and Cauchy Refine the Method

French mathematicians extended Laplace's ideas, developing the inverse transform and making the method more rigorous for solving differential equations.

1910s • Heaviside Popularizes Operational Calculus

Oliver Heaviside used similar transform techniques to solve circuit equations in electrical engineering, bringing the method into wide practical use.

1940s–Today • Standard Tool in College Engineering & Science

The Laplace transform method became a core technique taught in college-level engineering and mathematics programs worldwide, essential for courses in differential equations, circuits, and applied mathematics.

The key question this lesson addresses is: How can we systematically solve a linear differential equation with initial conditions by transforming it into an algebraic equation? By the end of this lesson, you will have a clear, step-by-step procedure for doing exactly that.

SECTION 2

Core Principles & Definitions

Before diving into the method, you need to understand a few foundational ideas. The Laplace transform converts a function of time, f(t), into a function of a new variable s, written as F(s). Think of s simply as an algebraic variable — like x in a polynomial — that helps us work in a different mathematical setting where calculus operations become easier. This conversion turns derivatives into algebraic expressions involving s, which is the entire secret of the method. Once the equation is algebraic, you solve for F(s), and then use the inverse Laplace transform to convert back to the time domain.

Laplace Transform

The operation ℒ{f(t)} = F(s) converts a time-domain function into the s-domain. It replaces calculus operations (derivatives) with algebra operations (multiplication by s).

Initial Value Problem (IVP)

A differential equation together with starting conditions, such as y(0) = 2 and y′(0) = −1. These conditions specify the exact solution among infinitely many possibilities.

s-Domain

The algebraic world where the transformed equation lives. Here, derivatives become polynomials in s, making the equation much simpler to manipulate. You can think of s as just an algebraic variable, similar to x in ordinary polynomial algebra.

Partial Fraction Decomposition

An algebra technique that breaks a complicated fraction into simpler pieces. For example, (s + 3)/((s+1)(s+2)) can be split into A/(s+1) + B/(s+2). Each simple piece can then be matched to a known inverse Laplace transform in a table.

Inverse Laplace Transform

The reverse operation ℒ⁻¹{F(s)} = f(t) that converts the algebraic answer back to a function of time — the actual solution to the original problem.

✦

KEY TAKEAWAY

Think of the Laplace transform like translating a difficult sentence from English into a simpler language, solving the puzzle there, and then translating the answer back. Differential equations are hard to solve directly in the "time language," but once you translate them into the "s-language," they become straightforward algebra problems. After solving, you translate back to get your answer in terms of time.

SECTION 3

Visual Explanation — The Laplace Transform Pipeline

The diagram below illustrates the complete pipeline for solving an IVP using Laplace transforms. Instead of trying to solve the differential equation directly (the hard way), you follow the three-step path around the diagram: transform, solve algebra, and inverse-transform.

The Laplace transform pipeline: instead of solving the differential equation directly (the dashed arrow, the hard way), follow the three solid arrows. Step 1 transforms to the s-domain, Step 2 solves algebra, and Step 3 converts back to the time-domain answer.

Notice how the three-step path goes around the outside of the diagram in a clockwise direction. The dashed arrow down the left side represents trying to solve the differential equation directly — which is often much harder. The beauty of this method is that each of the three steps uses techniques you already know or can learn quickly: applying transform formulas, doing algebra, and looking up inverse transforms in a table.

SECTION 4

Mathematical Framework

The mathematical machinery behind this method rests on a few key formulas. The most important ones are the transform rules for derivatives, because these are what convert a differential equation into an algebraic equation. Let's look at each formula you'll need.

LAPLACE TRANSFORM DEFINITION

ℒ{f(t)} = F(s) = ∫₀∞ e⁻ˢᵗ f(t) dt

This integral defines the transform: it multiplies f(t) by a decaying exponential e⁻ˢᵗ and integrates from 0 to infinity. In practice, you rarely compute this integral by hand — instead, you use a table of known transforms. The variable s is treated here as a positive real number that is large enough to make the integral converge (produce a finite value). The condition written as "s > a" in the transform table simply means s must be greater than some minimum value for the formula to be valid.

FIRST DERIVATIVE RULE

ℒ{y′(t)} = sY(s) − y(0)

Taking the Laplace transform of a first derivative produces sY(s) minus the initial value y(0). Notice how the initial condition is automatically built in — this is a major advantage of the method.

SECOND DERIVATIVE RULE

ℒ{y″(t)} = s²Y(s) − s·y(0) − y′(0)

For a second derivative, you get s²Y(s) minus terms involving both initial conditions y(0) and y′(0). Higher derivatives follow a similar pattern, picking up more initial condition terms.

COMMON INVERSE TRANSFORMS

ℒ⁻¹{1/(s − a)} = eᵃᵗ ℒ⁻¹{n!/(s − a)ⁿ⁺¹} = tⁿeᵃᵗ

These are the formulas you use in Step 3 to convert back from the s-domain. The first says that a simple fraction 1/(s − a) transforms back to the exponential eᵃᵗ. You look these up in a Laplace transform table rather than memorizing all of them.

💡

Why Initial Conditions Matter

When you solve a differential equation using other methods, you typically find a general solution with arbitrary constants (like C₁ and C₂) and then use the initial conditions to find those constants in a separate step. With the Laplace method, the initial conditions are incorporated from the very beginning through the derivative rules. This means you get the specific solution directly — no extra step needed.

SECTION 5

Essential Laplace Transform Pairs

A Laplace transform table is your most important reference tool when using this method. It lists functions f(t) alongside their transforms F(s), so you can quickly convert in either direction. Below is a table of the most commonly needed pairs for solving IVPs, followed by a diagram showing how partial fractions connect the s-domain solution to these table entries.

f(t) — Time Domain	F(s) — s-Domain	Conditions
1	1/s	s > 0
t	1/s²	s > 0
tⁿ	n!/sⁿ⁺¹	s > 0, n = 1, 2, 3, …
eᵃᵗ	1/(s − a)	s > a
sin(bt)	b/(s² + b²)	s > 0
cos(bt)	s/(s² + b²)	s > 0
eᵃᵗ sin(bt)	b/((s − a)² + b²)	s > a
eᵃᵗ cos(bt)	(s − a)/((s − a)² + b²)	s > a

Essential Laplace transform pairs for solving linear ODE IVPs

💡

Partial Fractions: A Key Algebra Skill

Partial fraction decomposition is a technique from algebra (often first seen in a precalculus or calculus course) that breaks a complex rational expression into a sum of simpler fractions. For example: (s + 3)/((s+1)(s+2)) = A/(s+1) + B/(s+2). To find A and B, multiply both sides by the full denominator, then substitute convenient values of s to solve. This skill is central to the Laplace method — practice it on simple examples first if it is new to you.

This diagram traces how a complicated fraction Y(s) is broken down step by step via partial fraction decomposition until each piece matches a known entry in the Laplace transform table, allowing you to apply the inverse transform.

The key insight is that partial fraction decomposition is the bridge between the algebraic solution Y(s) and the table of inverse transforms. If you can factor the denominator and split the fraction, you can always convert back to the time domain. This is why strong algebra skills — especially factoring polynomials and working with fractions — are essential for this method.

SECTION 6

Worked Example

Let's solve a complete initial value problem from start to finish. We'll work through every step so you can see exactly how the method operates.

ℹ

Problem Statement

Solve the IVP: y″ + 3y′ + 2y = 0 with initial conditions y(0) = 1 and y′(0) = 0.

Solving y″ + 3y′ + 2y = 0 via Laplace Transforms

1
Step 1 — Apply the Laplace Transform to Both Sides

Take ℒ of each term: ℒ{y″} + 3·ℒ{y′} + 2·ℒ{y} = ℒ{0}. Using the derivative rules with y(0) = 1 and y′(0) = 0, we substitute: [s²Y(s) − s·y(0) − y′(0)] + 3[sY(s) − y(0)] + 2Y(s) = 0. Plugging in initial values gives: [s²Y(s) − s·1 − 0] + 3[sY(s) − 1] + 2Y(s) = 0.

s²Y(s) − s + 3sY(s) − 3 + 2Y(s) = 0
2
Step 2 — Collect Terms and Solve for Y(s)

Group all terms with Y(s) on one side: (s² + 3s + 2)Y(s) = s + 3. Now solve for Y(s) by dividing both sides by the polynomial in front of Y(s).

Y(s) = (s + 3) / (s² + 3s + 2)
3
Step 3 — Factor the Denominator

Factor s² + 3s + 2. We need two numbers that multiply to 2 and add to 3. Those numbers are 1 and 2, so: s² + 3s + 2 = (s + 1)(s + 2).

Y(s) = (s + 3) / ((s + 1)(s + 2))
4
Step 4 — Perform Partial Fraction Decomposition

Write (s + 3)/((s + 1)(s + 2)) = A/(s + 1) + B/(s + 2). Multiply both sides by (s + 1)(s + 2) to get: s + 3 = A(s + 2) + B(s + 1). Set s = −1: (−1 + 3) = A(−1 + 2) + B(0), so 2 = A. Set s = −2: (−2 + 3) = A(0) + B(−2 + 1), so 1 = −B, meaning B = −1.

Y(s) = 2/(s + 1) + (−1)/(s + 2)
5
Step 5 — Apply the Inverse Laplace Transform

Use the table: ℒ⁻¹{1/(s + 1)} = e⁻ᵗ and ℒ⁻¹{1/(s + 2)} = e⁻²ᵗ. By linearity of the inverse transform: y(t) = 2·e⁻ᵗ + (−1)·e⁻²ᵗ.

y(t) = 2e⁻ᵗ − e⁻²ᵗ
6
Step 6 — Verify the Initial Conditions

Check: y(0) = 2e⁰ − e⁰ = 2 − 1 = 1 ✓. For the derivative: y′(t) = −2e⁻ᵗ + 2e⁻²ᵗ, so y′(0) = −2 + 2 = 0 ✓. Both initial conditions are satisfied, confirming our solution is correct.

Solution verified: y(t) = 2e⁻ᵗ − e⁻²ᵗ ✓

SECTION 7

Strengths & Limitations of the Laplace Method

The Laplace transform method is powerful, but it is not always the best choice. Understanding when to use it — and when another method might be simpler — is an important part of building your problem-solving toolkit.

Aspect	Strengths	Limitations
Initial Conditions	Automatically incorporates y(0), y′(0), etc., from the start — no extra step to find constants.	Requires initial conditions at t = 0. Boundary value problems (conditions at two different points) need other methods.
Equation Types	Handles linear ODEs with constant coefficients very efficiently, including those with discontinuous or impulsive forcing functions.	Does not directly handle nonlinear equations or variable-coefficient equations well.
Procedure	Systematic and recipe-like — follow the same steps every time. Reduces the need for clever guesses.	Partial fraction decomposition can become algebraically messy for higher-order equations.
Discontinuous Inputs	Excels at problems with step functions, impulses, or piecewise forcing — areas where other methods struggle.	Requires learning additional tools (Heaviside function, Dirac delta) for these advanced inputs.

Comparing strengths and limitations of the Laplace transform method

✦

KEY TAKEAWAY

Think of the Laplace method as a GPS for solving differential equations. A GPS doesn't find the most creative scenic route — it gives you reliable, turn-by-turn directions that always get you to the destination. Similarly, the Laplace method gives you a consistent procedure that works every time for linear constant-coefficient problems, even if there might be a quicker shortcut for simple cases. It's especially valuable when the road gets complicated — like when the forcing function has sudden jumps or spikes.

SECTION 8

Connection to Advanced Topics

The Laplace transform method you've learned here is the gateway to a rich set of advanced techniques used throughout engineering and applied mathematics. Understanding where this topic leads can help you appreciate its importance and prepare for later coursework.

What You've Learned	Where It Leads
Solving homogeneous IVPs (right side = 0)	Non-homogeneous IVPs with forcing functions like sin(t), step functions, or impulses on the right side
Partial fraction decomposition with distinct real roots	Repeated roots and complex conjugate roots, requiring modified partial fraction techniques
Single ODE problems	Systems of coupled differential equations, where multiple unknowns interact
Laplace transforms of basic functions	Transfer functions and frequency response — key concepts in electrical engineering and control systems courses
Transform tables for known functions	The convolution theorem, which handles products of transforms and integral equations

How today's skills connect to advanced topics

In later engineering courses, you'll encounter the concept of a transfer function, which is essentially the Laplace transform of a system's response — a way of describing how a system reacts to different inputs using algebra in the s-domain. Engineers use transfer functions to design everything from autopilot systems in aircraft to noise-canceling headphones. This algebraic perspective — treating systems as input-output relationships — becomes the foundation of advanced coursework in control theory and signal processing, typically studied in electrical and mechanical engineering programs.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain in your own words why the Laplace transform method automatically incorporates initial conditions, whereas the method of undetermined coefficients requires a separate step to find constants C₁ and C₂. What specific formula makes this possible?

Show answer

The Laplace transform derivative rules — ℒ{y′} = sY(s) − y(0) and ℒ{y″} = s²Y(s) − s·y(0) − y′(0) — contain the initial values y(0) and y′(0) explicitly. When you apply the transform to the ODE, these values are plugged in immediately, so the resulting algebraic equation already reflects the specific initial conditions. In contrast, methods like undetermined coefficients first find a general solution with arbitrary constants, and then substitute the initial conditions in a separate step to determine those constants. The Laplace method skips straight to the particular solution.

Question 2: BASIC CALCULATION

Solve the IVP: y′ + 2y = 0, with y(0) = 5. Use the Laplace transform method to find y(t).

Show answer

Apply ℒ to both sides: sY(s) − y(0) + 2Y(s) = 0. Substitute y(0) = 5: sY(s) − 5 + 2Y(s) = 0. Factor: (s + 2)Y(s) = 5, so Y(s) = 5/(s + 2). Using the table entry ℒ⁻¹{1/(s − a)} = eᵃᵗ with a = −2, the inverse transform gives y(t) = 5e⁻²ᵗ. You can verify: y(0) = 5e⁰ = 5 ✓, and y′ + 2y = −10e⁻²ᵗ + 10e⁻²ᵗ = 0 ✓.

Question 3: INTERMEDIATE

Solve the IVP: y″ + 5y′ + 6y = 0, with y(0) = 2 and y′(0) = −3. Show all steps including the partial fraction decomposition.

Show answer

Applying ℒ: [s²Y − 2s − (−3)] + 5[sY − 2] + 6Y = 0, which simplifies to (s² + 5s + 6)Y = 2s + 7. So Y(s) = (2s + 7)/(s² + 5s + 6) = (2s + 7)/((s + 2)(s + 3)). For partial fractions: (2s + 7)/((s + 2)(s + 3)) = A/(s + 2) + B/(s + 3). Setting s = −2: (−4 + 7) = A(1), so A = 3. Setting s = −3: (−6 + 7) = B(−1), so B = −1. Therefore Y(s) = 3/(s + 2) − 1/(s + 3). Taking ℒ⁻¹: y(t) = 3e⁻²ᵗ − e⁻³ᵗ. Verification: y(0) = 3 − 1 = 2 ✓, y′(0) = −6 + 3 = −3 ✓.

Question 4: APPLIED

A simple RC circuit has the equation for charge: q′ + (1/RC)q = 0, where R = 500 Ω, C = 0.002 F, and the initial charge is q(0) = 10 C. Use the Laplace transform to find q(t), the charge as a function of time. After how many seconds does the charge drop to about 37% of its initial value?

Show answer

First compute 1/(RC) = 1/(500 × 0.002) = 1. The equation becomes q′ + q = 0 with q(0) = 10. Applying ℒ: sQ(s) − 10 + Q(s) = 0, so (s + 1)Q(s) = 10, giving Q(s) = 10/(s + 1). The inverse transform yields q(t) = 10e⁻ᵗ. The charge drops to 37% (which is approximately 1/e ≈ 0.368) of its initial value when e⁻ᵗ = 1/e, meaning t = 1 second. This time constant τ = RC = 1 second is a fundamental concept in circuit analysis, and the Laplace method finds it naturally.

Question 5: CRITICAL THINKING

Consider the IVP: y″ + 4y = 0, with y(0) = 0 and y′(0) = 6. Solve using Laplace transforms. Then explain why the solution involves trigonometric functions instead of exponentials, and relate this to the nature of the roots of the characteristic equation s² + 4 = 0.

Show answer

Applying ℒ: s²Y(s) − 0 − 6 + 4Y(s) = 0, so (s² + 4)Y(s) = 6, giving Y(s) = 6/(s² + 4). Recognizing this matches the table entry b/(s² + b²) with b = 2, we get y(t) = 3 sin(2t). The solution is trigonometric rather than exponential because the characteristic equation s² + 4 = 0 has roots s = ±2i — that is, the solutions for s are the imaginary numbers 2i and −2i, where i² = −1. When the roots of the characteristic equation are real numbers, you get exponential solutions; when they are purely imaginary (±bi), you get oscillatory solutions involving sin(bt) and cos(bt). Physically, this corresponds to undamped oscillation — like a frictionless spring — where energy is conserved and the system oscillates forever. The Laplace transform seamlessly handles both cases through the same algebraic procedure.

SUMMARY

Lesson Summary

Summary

The Laplace transform method provides a systematic, three-step procedure for solving linear ODE initial value problems. First, apply the Laplace transform to convert the differential equation into an algebraic equation in the s-domain, using the derivative rules ℒ{y′} = sY(s) − y(0) and ℒ{y″} = s²Y(s) − s·y(0) − y′(0), which automatically incorporate the initial conditions. Second, solve the resulting algebra for Y(s). Third, use partial fraction decomposition and a transform table to apply the inverse Laplace transform and obtain the solution y(t).

This method is especially powerful because it handles initial conditions from the start (no separate step to find constants), follows a consistent recipe every time, and extends naturally to problems with discontinuous forcing functions. The core algebraic skills needed are factoring polynomials and performing partial fractions. This topic is a standard part of college-level differential equations courses, and mastering it prepares you for advanced study in control theory, signal processing, and engineering design.