The Laplace transform converts a function of a real variable $t \ge 0$ — typically representing time — into a function of a complex variable $p$, replacing [differentiation](/page/Derivative) with multiplication and initial-value problems with algebraic equations. It is the principal tool of **operational calculus**: the systematic reduction of linear differential equations with constant coefficients to polynomial algebra in the transform variable. Unlike the [Fourier transform](/page/Fourier%20Transform), which requires global [integrability](/page/Integral) or square-integrability over the entire real line, the Laplace transform accommodates exponentially growing signals by restricting integration to the half-line $[0, \infty)$ and introducing an exponential damping factor $e^{-pt}$. This makes it the natural transform for causal systems — those whose behaviour depends only on the present and the past.

The theory rests on three pillars. First, the **forward transform** maps a function $f$ to its image $\hat{f}(p)$, which is holomorphic in a half-plane determined by the growth rate of $f$. Second, the **operational rules** express derivatives, integrals, convolutions, and shifts of $f$ as algebraic operations on $\hat{f}$. Third, the **inversion formula** recovers $f$ from $\hat{f}$ via a contour integral in the complex plane, connecting the Laplace transform to the theory of residues and [complex analysis](/page/Cambridge%20IB%20Complex%20Analysis).

[motivation] ## Motivation ### Direct Methods for Constant-Coefficient ODEs Consider a linear ordinary differential equation with constant coefficients, such as the second-order initial-value problem \begin{align*} y'' + 3y' + 2y &= e^{-3t}, \quad t > 0, \\ y(0) &= 1, \quad y'(0) = 0. \end{align*} The classical method proceeds in two stages: find the general solution to the homogeneous equation $y'' + 3y' + 2y = 0$ via the characteristic polynomial $\lambda^2 + 3\lambda + 2 = (\lambda + 1)(\lambda + 2) = 0$, then find a particular solution to the non-homogeneous equation by undetermined coefficients or [variation of parameters](/theorems/835), and finally impose initial conditions to determine the constants. This works, but it has two weaknesses. First, the initial conditions are imposed at the very end, after the general solution has been constructed — meaning the work of finding the full general solution is performed even though only one specific solution is needed. Second, the method scales poorly: for higher-order equations, systems of equations, or equations with discontinuous or impulsive forcing terms (such as the Heaviside step function or the Dirac delta), the machinery of undetermined coefficients becomes cumbersome or inapplicable entirely. What we want is a method that incorporates the initial conditions from the outset and reduces the entire problem — differential equation plus initial data — to algebra in one step. ### The Fourier Transform Attempt A natural candidate is the [Fourier transform](/page/Fourier%20Transform). If $f: \mathbb{R} \to \mathbb{C}$ is integrable, its Fourier transform $\hat{f}(\xi) = \int_{\mathbb{R}} f(t) e^{-i\xi t} \, d\mathcal{L}^1(t)$ converts differentiation into multiplication by $i\xi$, so a differential equation becomes a polynomial equation in frequency space. However, the Fourier transform requires $f \in L^1(\mathbb{R}, \mathcal{L}^1)$ (or $L^2(\mathbb{R}, \mathcal{L}^1)$ via the Plancherel extension), and the solutions to initial-value problems on $[0, \infty)$ typically grow — or at least do not decay — as $t \to \infty$. For example, the solution $y(t) = e^{-t}$ to the homogeneous problem $y' + y = 0$, $y(0) = 1$ does lie in $L^1([0,\infty))$, but the solution $y(t) = e^t$ to $y' - y = 0$, $y(0) = 1$ does not, and neither does the constant solution $y(t) = 1$ to $y' = 0$, $y(0) = 1$. The Fourier integral simply diverges for these functions. There is a second, subtler difficulty. Even when the integral converges, the Fourier transform is defined over $(-\infty, \infty)$ and treats the function as a two-sided signal. It has no mechanism for encoding initial conditions at $t = 0$ — the very data that distinguishes one solution of a linear ODE from another. Any transform that hopes to solve initial-value problems must be sensitive to the starting point $t = 0$. ### Exponential Damping and the Half-Line Both difficulties are resolved by a single modification: restrict the domain of integration to $[0, \infty)$ and introduce a decaying exponential weight $e^{-pt}$ with $\operatorname{Re} p$ large enough to kill the growth of $f$. If $f$ grows at most as fast as $e^{\sigma_0 t}$ for some $\sigma_0 \in \mathbb{R}$, then the product $f(t) e^{-pt}$ decays exponentially whenever $\operatorname{Re} p > \sigma_0$, and the integral \begin{align*} \hat{f}(p) = \int_0^\infty f(t) e^{-pt} \, d\mathcal{L}^1(t) \end{align*} converges absolutely. This is the Laplace transform. The half-line integration is not merely a convenience — it is the mechanism by which initial conditions enter the calculus. When we integrate $f'(t) e^{-pt}$ by parts over $[0, \infty)$, the boundary term at $t = 0$ produces exactly the initial value $f(0)$, so the derivative rule becomes $\mathcal{L}\{f'\}(p) = p\hat{f}(p) - f(0)$. The initial data is encoded automatically in the algebraic relation. [/motivation]

## The Laplace Transform

### Exponential Order

Before defining the transform, we must identify the class of [functions](/page/Function) for which the defining integral converges. The key condition is that $f$ does not grow faster than some exponential — without this constraint, no exponential damping factor $e^{-pt}$ can restore integrability.

[definition:Exponential Order] A function $f: [0, \infty) \to \mathbb{C}$ is of **exponential order** if there exist constants $M > 0$ and $\sigma_0 \in \mathbb{R}$ such that \begin{align*} |f(t)| \le M e^{\sigma_0 t} \quad \text{for all } t \ge 0. \end{align*} The infimum of all such $\sigma_0$ is called the **abscissa of exponential growth** of $f$. [/definition]

Functions of exponential order include all polynomials (with $\sigma_0 > 0$ arbitrarily small), all exponentials $e^{at}$ (with $\sigma_0 = \operatorname{Re} a$), all bounded functions (with $\sigma_0 = 0$), and products of these.

[example:Failure Of Exponential Order] The function $f: [0, \infty) \to \mathbb{R}$ defined by $f(t) = e^{t^2}$ is not of exponential order. To see this, suppose for contradiction that $|f(t)| \le M e^{\sigma_0 t}$ for all $t \ge 0$. Then $e^{t^2} \le M e^{\sigma_0 t}$, which gives $e^{t^2 - \sigma_0 t} \le M$ for all $t$. But \begin{align*} t^2 - \sigma_0 t = \left(t - \frac{\sigma_0}{2}\right)^2 - \frac{\sigma_0^2}{4} \to \infty \quad \text{as } t \to \infty, \end{align*} so $e^{t^2 - \sigma_0 t} \to \infty$, contradicting the bound. The super-exponential growth of $e^{t^2}$ cannot be tamed by any linear exponential $e^{\sigma_0 t}$. Consequently, the integral $\int_0^\infty e^{t^2} e^{-pt} \, d\mathcal{L}^1(t)$ diverges for every $p \in \mathbb{C}$: for any fixed $\operatorname{Re} p$, the integrand $e^{t^2 - (\operatorname{Re} p) t}$ grows without bound, and the Laplace transform of $e^{t^2}$ does not exist anywhere in the complex plane. [/example]

### Definition

With exponential order in hand, we can define the transform itself. The construction follows directly from the motivation: integrate $f$ against the kernel $e^{-pt}$ over the half-line, with $\operatorname{Re} p$ large enough to ensure convergence.

[definition:Laplace Transform] Let $f: [0, \infty) \to \mathbb{C}$ be a locally integrable function of exponential order with abscissa of exponential growth $\sigma_0$. The **Laplace transform** of $f$ is the function \begin{align*} \mathcal{L}\{f\}: \{p \in \mathbb{C} : \operatorname{Re} p > \sigma_0\} &\to \mathbb{C} \\ p &\mapsto \int_0^\infty f(t) e^{-pt} \, d\mathcal{L}^1(t). \end{align*} We write $\hat{f}(p) = \mathcal{L}\{f\}(p)$ for the transform of $f$ evaluated at $p$. [/definition]

The transform is defined on a half-plane $\operatorname{Re} p > \sigma_0$, not on a line as in the Fourier case. The larger $\operatorname{Re} p$, the stronger the damping — so the integral converges more easily. The [boundary](/page/Boundary) $\operatorname{Re} p = \sigma_0$ is the critical line: on it the integral may or may not converge, depending on the finer behaviour of $f$.

[example:Elementary Transforms] The simplest transforms can be computed directly from the definition. **The exponential.** Let $a \in \mathbb{C}$ and define $f: [0, \infty) \to \mathbb{C}$ by $f(t) = e^{at}$. Then $f$ is of exponential order with $\sigma_0 = \operatorname{Re} a$. For $\operatorname{Re} p > \operatorname{Re} a$: \begin{align*} \hat{f}(p) &= \int_0^\infty e^{at} e^{-pt} \, d\mathcal{L}^1(t) = \int_0^\infty e^{-(p-a)t} \, d\mathcal{L}^1(t) = \left[ \frac{e^{-(p-a)t}}{-(p-a)} \right]_0^\infty = \frac{1}{p - a}. \end{align*} The boundary term at $t = \infty$ vanishes because $\operatorname{Re}(p - a) > 0$ ensures $|e^{-(p-a)t}| = e^{-\operatorname{Re}(p-a) \cdot t} \to 0$. Setting $a = 0$ gives $\mathcal{L}\{1\}(p) = 1/p$ for $\operatorname{Re} p > 0$. **Powers of $t$.** Let $n \in \mathbb{N}$ and define $f(t) = t^n$. Repeated integration by parts yields \begin{align*} \mathcal{L}\{t^n\}(p) &= \int_0^\infty t^n e^{-pt} \, d\mathcal{L}^1(t) = \frac{n}{p} \int_0^\infty t^{n-1} e^{-pt} \, d\mathcal{L}^1(t) = \cdots = \frac{n!}{p^{n+1}}, \end{align*} valid for $\operatorname{Re} p > 0$. At each stage, the integration by parts uses [Integration By Parts](/theorems/210) with $u = t^k$ and $dv = e^{-pt} \, dt$; the boundary term at infinity vanishes because the exponential decay dominates any polynomial growth. **The sine function.** Let $b \in \mathbb{R} \setminus \{0\}$ and define $f(t) = \sin(bt)$. Since $|\sin(bt)| \le 1$, this function is of exponential order with $\sigma_0 = 0$. Using the identity $\sin(bt) = (e^{ibt} - e^{-ibt})/(2i)$ and the exponential transform: \begin{align*} \mathcal{L}\{\sin(bt)\}(p) &= \frac{1}{2i}\left(\frac{1}{p - ib} - \frac{1}{p + ib}\right) = \frac{1}{2i} \cdot \frac{2ib}{p^2 + b^2} = \frac{b}{p^2 + b^2}, \end{align*} valid for $\operatorname{Re} p > 0$. An analogous computation gives $\mathcal{L}\{\cos(bt)\}(p) = p/(p^2 + b^2)$. **The Heaviside step.** Define $H_a: [0, \infty) \to \mathbb{R}$ by $H_a(t) = \mathbb{1}_{[a, \infty)}(t)$ for $a > 0$. Then $H_a$ is bounded, so $\sigma_0 = 0$. For $\operatorname{Re} p > 0$: \begin{align*} \mathcal{L}\{H_a\}(p) &= \int_a^\infty e^{-pt} \, d\mathcal{L}^1(t) = \left[\frac{e^{-pt}}{-p}\right]_a^\infty = \frac{e^{-ap}}{p}. \end{align*} The factor $e^{-ap}$ encodes the delay: shifting a signal by $a$ units in the time domain multiplies its transform by $e^{-ap}$ in the frequency domain. [/example]

## Analytic Properties

The Laplace transform is not merely a formal device — it produces a holomorphic function in the half-plane of convergence, and this holomorphicity is the source of its power.

### Holomorphicity

A function defined by an integral depending on a complex parameter is often holomorphic in that parameter, provided the integrand is sufficiently well-behaved. The Laplace integral is a paradigmatic example: the integrand $f(t) e^{-pt}$ is entire in $p$ for each fixed $t$, and the exponential order condition provides the uniform domination needed to differentiate under the integral sign.

[theorem:Holomorphicity Of The Laplace Transform] Let $f: [0, \infty) \to \mathbb{C}$ be locally integrable and of exponential order with abscissa $\sigma_0$. Then $\hat{f}$ is holomorphic on the half-plane $\{p \in \mathbb{C} : \operatorname{Re} p > \sigma_0\}$, and its derivative is given by \begin{align*} \hat{f}'(p) = -\int_0^\infty t f(t) e^{-pt} \, d\mathcal{L}^1(t) = -\mathcal{L}\{t f(t)\}(p). \end{align*} More generally, for each $n \in \mathbb{N}$: \begin{align*} \hat{f}^{(n)}(p) = (-1)^n \mathcal{L}\{t^n f(t)\}(p). \end{align*} [/theorem]

The proof proceeds by verifying the hypotheses of the [Dominated Convergence Theorem](/theorems/4): for any compact subset $K$ of the half-plane $\operatorname{Re} p > \sigma_0$, there exists $c > \sigma_0$ such that $\operatorname{Re} p \ge c$ on $K$, and the integrand $|f(t) e^{-pt}| \le M e^{(\sigma_0 - c)t}$ provides an integrable dominating function. Differentiating under the integral sign is then justified, and each differentiation pulls down a factor of $-t$.