To solve an IVP like y'' + 3y' + 2y = e^t, y(0) = 0, y'(0) = 1: (1) apply Laplace transform to get (s² + 3s + 2)Y(s) = 1/(s-1) + 1, (2) solve for Y(s), (3) use inverse transform to recover y(t). This method handles initial conditions automatically.
Solve several IVPs by hand using Laplace transforms, then compare answers using classical methods (undetermined coefficients, variation of parameters). Note how Laplace avoids computing the homogeneous solution separately.
The Laplace transform method is a change-of-domain strategy: instead of solving a differential equation directly in the time domain, you transform it into an algebraic equation in the s-domain, solve the algebra, then transform back. The power of this approach is that differentiation — which is the hard part of a differential equation — becomes multiplication by s in the s-domain. This converts the ODE into something you can solve with basic algebra.
You already know the Laplace transform of derivatives: if Y(s) = L{y(t)}, then L{y'(t)} = sY(s) − y(0) and L{y''(t)} = s²Y(s) − sy(0) − y'(0). Notice that the initial conditions y(0) and y'(0) appear automatically when you apply the derivative rule — they're baked into the transformed equation. This is the key structural advantage: you don't need to first find the general solution and then apply initial conditions as a separate step. The initial conditions enter at the same moment you transform the ODE.
The full procedure is three steps. Step 1: Transform. Apply L{·} to both sides of the ODE, using linearity and the derivative rules. For y'' + 3y' + 2y = eᵗ with y(0) = 0, y'(0) = 1, you get [s²Y − s·0 − 1] + 3[sY − 0] + 2Y = 1/(s − 1). Collecting Y terms: (s² + 3s + 2)Y = 1/(s − 1) + 1. Step 2: Solve for Y(s). This is pure algebra: Y(s) = [1/(s − 1) + 1] / (s² + 3s + 2). Factor the denominator: (s + 1)(s + 2). Use partial fractions to decompose Y(s) into a sum of simple fractions whose inverse transforms you know from your inverse Laplace transform table. Step 3: Invert. Apply L⁻¹{·} termwise to recover y(t).
The Laplace method particularly shines on problems where classical methods (undetermined coefficients, variation of parameters) require solving a homogeneous equation first, then a particular equation, then matching initial conditions — three separate stages. Laplace collapses all three into one pass. It also handles discontinuous forcing functions (like step functions and impulses) far more cleanly than classical methods, which is why it's the standard tool in engineering for control systems and signal processing.
Partial fraction decomposition is the algebraic core of step 2 and the most common source of errors. Once you have Y(s) as a ratio of polynomials, factor the denominator completely (real and complex roots), write Y as a sum of terms of the form A/(s − a), (As + B)/(s² + bs + c) for complex pairs, and A/(s − a)ᵏ for repeated roots. Each of these has a known inverse transform. The match between the algebraic form and the transform table is exact by design: the Laplace method works precisely because the functions in the table — exponentials, sinusoids, polynomials — are the building blocks of ODE solutions.