The Laplace Equation

Introduction

The heat equation and the wave equation were each introduced with a promise that the present page now redeems: the classical trichotomy of linear second-order PDEs — parabolic, hyperbolic, elliptic — has a third member, the Laplace equation, and we committed to developing it in a subsequent page. To present that third member as another point on the same axis, however, would be misleading. Heat and wave are initial-value problems in time: at \(t = 0\) something is prescribed, and the equation propagates it forward, parabolically in the heat case and hyperbolically in the wave case. The Laplace equation has no time variable. The problem it solves is structurally different — given values of \(u\) on the boundary \(\partial\Omega\) of a spatial region, what is \(u\) inside? — and the time axis along which heat and wave were arrayed is replaced here by the spatial boundary.

A second characterization stands alongside the first. Setting the time derivative of the heat equation to zero, or both time derivatives of the wave equation to zero, reduces each equation to a single common form: \[ \Delta u \;=\; \partial_{xx} u + \partial_{yy} u \;=\; 0. \] The Laplace equation is the time-evolution endpoint of the other two classical PDEs — the steady state to which a heat distribution settles after long time, and the equation satisfied by a vibrating string at rest. The same equation arises in physics whenever an equilibrium of a diffusive or elastic medium is sought: the electrostatic potential in a region free of charge, the shape of a soap film stretched across a wire frame. The solutions of \(\Delta u = 0\) — called harmonic functions — are the static, equilibrium analogues of the dynamical objects studied in the heat and wave pages.

A third characterization, the deepest of the three, asks what controls the solution. For the heat equation, the maximum principle gives \(\inf f \leq u(x, t) \leq \sup f\): the solution is bounded pointwise by the initial data. For the wave equation, no such pointwise bound exists; the conserved quantity is the global mechanical energy, an integral invariant rather than a pointwise estimate. For the Laplace equation, as the present page will prove, a pointwise bound returns — \(\min_{\partial\Omega} u \leq u(x) \leq \max_{\partial\Omega} u\) — but with the controlling data shifted from the time-zero slice to the spatial boundary. The three classical PDEs realize three answers to a single question, "what controls the solution?": initial data, global energy, boundary data. That this three-way pattern is more than a coincidence of names is a claim we are not yet prepared to substantiate; the proof of the maximum principle in the section on harmonic function properties will articulate the three answers side by side as a concrete object.

Laplace Equation on the Disk

The first concrete setting in which we solve the Laplace equation is the unit disk \(D = \{(x, y) : x^2 + y^2 < 1\}\) with boundary \(\partial D\) the unit circle. The Dirichlet problem on \(D\) prescribes the values of \(u\) on \(\partial D\) and asks for the harmonic extension into the interior: \[ \begin{cases} \Delta u(x, y) \;=\; 0, & (x, y) \in D, \\ u(x, y) \;=\; f(x, y), & (x, y) \in \partial D, \end{cases} \] where \(f\) is a prescribed continuous function on the unit circle. The structural parallel with the heat and wave equations on the bounded interval is direct: in each case a second-order linear PDE is supplemented with data on the boundary of the spatial domain. The difference is the absence of a time variable and the absence of initial data — what was a one-dimensional spatial domain \((0, L)\) with two endpoint conditions is now a two-dimensional spatial domain \(D\) with a one-dimensional boundary \(\partial D\) carrying a full continuum of conditions. The natural choice of coordinates is polar.

Polar Coordinates and Separation of Variables

Setting \(x = r\cos\theta\) and \(y = r\sin\theta\) with \(r \in [0, 1]\) and \(\theta \in [-\pi, \pi]\), a direct chain-rule computation rewrites the Laplacian in the form \[ \Delta u \;=\; \partial_{rr} u + \frac{1}{r}\partial_r u + \frac{1}{r^2}\partial_{\theta\theta} u. \] The Dirichlet problem becomes \[ \begin{cases} \partial_{rr} u + \dfrac{1}{r}\partial_r u + \dfrac{1}{r^2}\partial_{\theta\theta} u \;=\; 0, & 0 < r < 1,\ -\pi < \theta \leq \pi, \\ u(1, \theta) \;=\; f(\theta), & -\pi < \theta \leq \pi, \end{cases} \] where, with a slight abuse of notation, \(f(\theta)\) denotes the boundary datum as a function on the circle and is assumed continuous and \(2\pi\)-periodic in \(\theta\). Following the technique used for the heat equation on the bounded interval, we seek separated solutions \[ u(r, \theta) \;=\; R(r)\, \Theta(\theta). \] Substituting into the polar Laplace equation and multiplying through by \(r^2 / \bigl(R(r)\,\Theta(\theta)\bigr)\) — valid wherever neither factor vanishes — separates the variables: \[ \frac{r^2 R''(r) + r R'(r)}{R(r)} \;=\; -\frac{\Theta''(\theta)}{\Theta(\theta)}. \] The left-hand side depends only on \(r\), the right only on \(\theta\); both must equal a common constant, which we denote \(\lambda\). The Dirichlet problem on the disk thus splits into a pair of ordinary differential equations coupled through the spectral parameter \(\lambda\): \[ \begin{align*} \Theta''(\theta) + \lambda\, \Theta(\theta) &\;=\; 0, \\\\ r^2 R''(r) + r R'(r) - \lambda\, R(r) &\;=\; 0. \end{align*} \]

The Angular Eigenvalue Problem

The angular equation \(\Theta'' + \lambda \Theta = 0\) is the same eigenvalue ODE that arose in the heat equation on the bounded interval, but the boundary conditions are different. There \(\Theta\) was required to vanish at two endpoints of an interval; here \(\theta\) parametrizes a closed curve and \(\Theta\) must instead be \(2\pi\)-periodic: \[ \Theta(\theta + 2\pi) \;=\; \Theta(\theta) \qquad \text{for all } \theta \in \mathbb{R}. \] This periodicity replaces the case analysis on the sign of \(\lambda\) that the Dirichlet endpoint conditions required. For \(\lambda < 0\), the general solution \(\Theta(\theta) = c_1 e^{\sqrt{-\lambda}\,\theta} + c_2 e^{-\sqrt{-\lambda}\,\theta}\) grows or decays exponentially and admits no nontrivial periodic representative; for \(\lambda = 0\), \(\Theta(\theta) = c_1 + c_2\theta\) is periodic only when \(c_2 = 0\), giving the constant solution \(\Theta_0 = 1\); for \(\lambda > 0\), the general solution \(\Theta(\theta) = c_1 \cos\sqrt{\lambda}\,\theta + c_2 \sin\sqrt{\lambda}\,\theta\) is \(2\pi\)-periodic if and only if \(\sqrt{\lambda}\) is a positive integer. The eigenvalues are therefore \[ \lambda_n \;=\; n^2, \qquad n \in \{0, 1, 2, \dots\}, \] with corresponding two-dimensional eigenspaces spanned by \(\{\cos n\theta, \sin n\theta\}\) for \(n \geq 1\) (or, in complex form, by \(e^{in\theta}\) and \(e^{-in\theta}\)) and one-dimensional for \(n = 0\) (spanned by the constant). Combining the positive and negative integer indices into a single signed parameter, we record the angular eigenfunctions in complex form as \[ \Theta_n(\theta) \;=\; e^{-in\theta}, \qquad n \in \mathbb{Z}. \]

The Radial Equation

With \(\lambda = n^2\) determined, the radial equation becomes \[ r^2 R''(r) + r R'(r) - n^2 R(r) \;=\; 0. \] This is the Cauchy–Euler equation, whose general solution is obtained by trying \(R(r) = r^\alpha\): substitution yields \(\alpha(\alpha - 1) + \alpha - n^2 = \alpha^2 - n^2 = 0\), so \(\alpha = \pm n\) for \(n \geq 1\). The general solution is \[ R_n(r) \;=\; A_n\, r^{n} + B_n\, r^{-n} \qquad (n \geq 1), \] with the degenerate case \(n = 0\) handled separately: the radial equation \(r^2 R'' + r R' = 0\) reduces to \((r R')' = 0\), giving \(R_0(r) = A_0 + B_0 \log r\). In both cases, one of the two basis solutions is singular as \(r \to 0^+\) — \(r^{-n}\) blows up for \(n \geq 1\), and \(\log r\) diverges to \(-\infty\) for \(n = 0\). We seek bounded harmonic functions on the entire open disk, including the origin; this regularity requirement forces \(B_n = 0\) for every \(n \geq 0\), leaving \[ R_n(r) \;=\; A_n\, r^{|n|}, \qquad n \in \mathbb{Z}, \] where the absolute value covers the negative-index case uniformly (the radial factor depends on \(|n|\) only, while the angular factor distinguishes \(n\) from \(-n\)).

The Series Solution

The product solutions \(R_n(r)\, \Theta_n(\theta) = r^{|n|} e^{-in\theta}\) are each harmonic on the disk. By superposition, any convergent series of the form \[ u(r, \theta) \;=\; \sum_{n \in \mathbb{Z}} c_n\, r^{|n|}\, e^{-in\theta} \] formally satisfies the Laplace equation. The remaining task is to choose the coefficients \(c_n\) so that the boundary condition \(u(1, \theta) = f(\theta)\) holds. Setting \(r = 1\) gives \[ f(\theta) \;=\; \sum_{n \in \mathbb{Z}} c_n\, e^{-in\theta}, \] which is precisely the complex Fourier series of the boundary datum \(f\). The coefficients are therefore the Fourier coefficients \[ c_n \;=\; \frac{1}{2\pi}\int_{-\pi}^{\pi} f(\phi)\, e^{+in\phi}\, d\phi. \] The boundary datum is decomposed into angular modes \(e^{in\phi}\), and each mode is propagated into the interior by the radial damping factor \(r^{|n|}\), which decays geometrically as \(r\) moves inward from the boundary at \(r = 1\) toward the origin. This is the structural mechanism of the disk Dirichlet problem: high-frequency boundary oscillations are damped most strongly toward the interior — the factor \(r^{|n|}\) shrinks faster for larger \(|n|\) — while low-frequency modes persist deeper into the disk.

The Poisson Integral Formula

Substituting the Fourier coefficients back into the series and interchanging the order of summation and integration — justified by the absolute convergence of \(\sum r^{|n|}\) for \(r < 1\) — yields \[ u(r, \theta) \;=\; \sum_{n \in \mathbb{Z}} \left(\frac{1}{2\pi}\int_{-\pi}^{\pi} f(\phi)\, e^{+in\phi}\, d\phi\right) r^{|n|}\, e^{-in\theta} \;=\; \int_{-\pi}^{\pi} f(\phi) \left(\frac{1}{2\pi}\sum_{n \in \mathbb{Z}} r^{|n|}\, e^{-in(\theta - \phi)}\right) d\phi. \] The inner sum is a function of \(r\) and \(\theta - \phi\) alone; writing \(\psi = \theta - \phi\) and splitting at \(n = 0\), \[ \sum_{n \in \mathbb{Z}} r^{|n|}\, e^{-in\psi} \;=\; 1 + \sum_{n=1}^{\infty} r^n\, e^{-in\psi} + \sum_{n=1}^{\infty} r^n\, e^{in\psi}, \] each tail is a geometric series in the ratio \(r e^{\mp i\psi}\), of modulus \(r < 1\). Summing, \[ \sum_{n=1}^{\infty} \bigl(r e^{-i\psi}\bigr)^n \;=\; \frac{r e^{-i\psi}}{1 - r e^{-i\psi}}, \qquad \sum_{n=1}^{\infty} \bigl(r e^{i\psi}\bigr)^n \;=\; \frac{r e^{i\psi}}{1 - r e^{i\psi}}. \] Adding these to \(1\), placing over a common denominator \((1 - r e^{-i\psi})(1 - r e^{i\psi}) = 1 - 2r\cos\psi + r^2\), and simplifying the numerator yields \[ \sum_{n \in \mathbb{Z}} r^{|n|}\, e^{-in\psi} \;=\; \frac{1 - r^2}{1 - 2r\cos\psi + r^2}. \] The right-hand side depends on \(\psi\) only through \(\cos\psi\), so the kernel is even in \(\psi\) and the sign in the exponent of the original sum did not affect the final closed form. Dividing by \(2\pi\), we obtain the kernel.

The Poisson kernel is positive on its domain (the denominator \(1 - 2r\cos\psi + r^2 = (1 - r)^2 + 2r(1 - \cos\psi) \geq (1 - r)^2 > 0\) is strictly positive, and the numerator \(1 - r^2 > 0\) for \(r < 1\)), and a direct integration of its series representation against \(d\psi\) gives \(\int_{-\pi}^{\pi} P_r(\psi)\, d\psi = 1\) (only the \(n = 0\) term survives after integration). With the kernel identified, the series solution takes integral form.

The Poisson kernel \(P_r(\theta - \phi)\) belongs to a broader family of integral kernels that arise in the theory of linear elliptic boundary-value problems. The fundamental object of that theory is the Green's function \(G(x, y)\) for the Laplace operator on a domain \(\Omega\), defined as the solution to \(\Delta_x G(x, y) = \delta_y(x)\) (the response to a point source at \(y\)) with \(G(x, y) = 0\) for \(x \in \partial\Omega\). The Green's function solves the inhomogeneous problem \(\Delta u = f\) with homogeneous boundary data; the Poisson kernel is its companion for the complementary problem treated on this page — the homogeneous equation \(\Delta u = 0\) with inhomogeneous boundary data \(u|_{\partial\Omega} = f\) — and is recovered from the Green's function by taking its normal derivative on the boundary, \(P(x, y) = -\partial G(x, y)/\partial n_y\) for \(y \in \partial\Omega\). The two kernels are therefore distinct objects: \(G\) satisfies an inhomogeneous equation involving the Dirac delta and is not a harmonic function, while \(P\) is itself harmonic on the open domain. The Green's function for the disk and its derivation via the method of images, along with the full distributional theory underlying \(\Delta G = \delta_y\), belong to a more advanced treatment than the present one and are not developed here; the name is introduced now because the same structural pair — a boundary kernel on a domain, arising as the normal derivative of an underlying Green's function — will reappear in the next section in a different geometric setting.

Laplace Equation on the Half-Plane

Replacing the unit disk by the upper half-plane \(\mathbb{H} = \{(x, y) : y > 0\}\) changes the analytic situation in the same structural way that passing from the bounded interval to the real line did for the heat equation: the boundary loses its compactness, the angular Fourier series of the disk becomes the continuous Fourier transform on \(\mathbb{R}\), and the discrete sum over modes is replaced by an integral over a continuous frequency variable. The Fourier transform developed earlier is the natural tool, and the parallel with the disk treatment will be exact: a Fourier-side computation that diagonalizes the spatial derivative will produce, after inversion, a Poisson integral representation with a kernel of different geometric form but the same structural role.

The Dirichlet problem on the upper half-plane is \[ \begin{cases} \Delta u(x, y) \;=\; 0, & (x, y) \in \mathbb{H}, \\ u(x, 0) \;=\; f(x), & x \in \mathbb{R}, \end{cases} \] with \(f\) a continuous and suitably decaying boundary datum on \(\mathbb{R}\), and \(u(x, y)\) required to decay as \(y \to \infty\) and to remain bounded as \(|x| \to \infty\). These decay requirements pick out the natural function class for the problem and rule out, for instance, the addition of a constant (which trivially satisfies \(\Delta u = 0\) but fails the decay condition unless the constant is zero) or polynomially growing harmonic backgrounds; the same role is played here by decay-at-infinity that was played on the disk by regularity at the origin.

Fourier Transform in the Horizontal Variable

Define the spatial Fourier transform of \(u(\cdot, y)\) for each fixed \(y > 0\): \[ \hat{u}(\xi, y) \;=\; \int_{-\infty}^{\infty} u(x, y)\, e^{ix\xi}\, dx. \] Applying the transform to the Laplace equation and using the differentiation property twice in \(x\) — which produces the factor \((-i\xi)^2 = -\xi^2\) — converts the partial differential equation in \((x, y)\) into an ordinary differential equation in \(y\) alone, with \(\xi\) entering as a parameter:

The structural parallel with the disk is exact. There the angular Fourier coefficients \(c_n\) of the boundary datum were propagated inward by the radial factor \(r^{|n|}\), which decays as \(r\) moves from the boundary \(r = 1\) toward the interior origin \(r = 0\). Here the continuous Fourier transform \(\hat{f}(\xi)\) of the boundary datum is propagated upward by the factor \(e^{-|\xi| y}\), which decays as \(y\) moves from the boundary \(y = 0\) into the interior \(y > 0\). In both cases, high-frequency boundary content is damped fastest — \(r^{|n|}\) shrinks rapidly for large \(|n|\), and \(e^{-|\xi| y}\) shrinks rapidly for large \(|\xi|\) — and only low-frequency content survives deep into the interior. The damping factor in each case is parametrized by the distance from the boundary, measured along the natural interior coordinate.

The Half-Plane Poisson Kernel

To recover \(u(x, y)\) we invert the Fourier transform of the right-hand side, applying the Fourier inversion formula: \[ u(x, y) \;=\; \frac{1}{2\pi}\int_{-\infty}^{\infty} \hat{f}(\xi)\, e^{-|\xi| y}\, e^{-ix\xi}\, d\xi. \] The identification of this expression as a convolution of \(f\) with a kernel requires the inverse transform of \(e^{-|\xi| y}\) itself, which is computable in closed form by a direct split-integral computation. For fixed \(y > 0\), \[ \frac{1}{2\pi}\int_{-\infty}^{\infty} e^{-|\xi| y}\, e^{-ix\xi}\, d\xi \;=\; \frac{1}{2\pi}\left[\int_0^{\infty} e^{-\xi y - ix\xi}\, d\xi + \int_{-\infty}^{0} e^{\xi y - ix\xi}\, d\xi\right]. \] Each integral is elementary: the first equals \(\frac{1}{y + ix}\), the second \(\frac{1}{y - ix}\). Their sum is \[ \frac{1}{y + ix} + \frac{1}{y - ix} \;=\; \frac{(y - ix) + (y + ix)}{(y + ix)(y - ix)} \;=\; \frac{2y}{x^2 + y^2}. \] Dividing by \(2\pi\) gives the kernel.

The half-plane Poisson kernel is positive on its domain (\(y > 0\) makes both numerator and denominator positive), and a direct integration gives \[ \int_{-\infty}^{\infty} P_y(x)\, dx \;=\; \frac{1}{\pi}\int_{-\infty}^{\infty} \frac{y}{x^2 + y^2}\, dx \;=\; \frac{1}{\pi}\left[\arctan\!\frac{x}{y}\right]_{-\infty}^{\infty} \;=\; 1. \] With the kernel identified, the inverse transform expresses the solution as a convolution.

Parallel Structure: Disk and Half-Plane

The two Dirichlet problems just solved exhibit a structural symmetry that goes beyond shared terminology. In both cases the boundary datum is decomposed into Fourier modes — discrete \(e^{-in\theta}\) on the circle, continuous \(e^{-ix\xi}\) on the line — and each mode is propagated into the interior by a damping factor parametrized by the distance from the boundary. On the disk, the interior coordinate is the radius \(r \in (0, 1)\) and the damping factor is \(r^{|n|}\); on the half-plane, the interior coordinate is the vertical distance \(y > 0\) and the damping factor is \(e^{-|\xi| y}\). The result in both cases is the same structural object: a Poisson kernel — a positive harmonic kernel on the interior that integrates to one over the boundary and concentrates onto a single boundary point as the interior coordinate approaches the boundary. The disk kernel \(P_r(\psi) = \frac{1}{2\pi}\frac{1 - r^2}{1 - 2r\cos\psi + r^2}\) and the half-plane kernel \(P_y(x) = \frac{1}{\pi}\frac{y}{x^2 + y^2}\) are two realizations of this structural pattern, differing in geometric form because the two domains differ in geometry, but identical in role; both arise, by the relation noted at the close of the previous section, as the boundary normal derivative of an underlying Green's function on their respective domains.

A further parallel reaches back to the heat and wave equations. The heat equation on the real line was solved by a kernel parametrized by the time variable \(t\) — the heat kernel \(K_t(x) = (4\pi k t)^{-1/2} e^{-x^2/(4kt)}\), which concentrates onto the initial datum as \(t \to 0^+\) and provides the solution by convolution. The wave equation admitted no such kernel representation: d'Alembert's formula is a different structural object, not a convolution with a single approximate identity. The Laplace equation on the disk and on the half-plane both admit kernel representations — Poisson kernels — parametrized by the boundary-distance coordinate, which plays in the elliptic setting the role that time played in the parabolic setting. Where heat has a one-parameter family of kernels concentrating onto the time-zero slice, Laplace has a one-parameter family of kernels concentrating onto the spatial boundary; where heat propagates information forward in time by convolution with a Gaussian, Laplace propagates information inward from the boundary by convolution with a Poisson kernel.

Harmonic Function Properties

The Poisson integral formulas just constructed are tied to two specific domains — the disk and the half-plane — and to the coordinate systems adapted to those domains. We now turn to two qualitative properties that hold for every harmonic function on any open subset of the plane, irrespective of domain shape: the mean value property and the maximum principle. Both will be derived from the Poisson integral on the disk, but their conclusions are domain-independent and supply the analytic structure that, in the heat-and-wave parallel, was carried by smoothing and energy conservation respectively. The maximum principle in particular is the signature property of the elliptic case — the pointwise bound on the solution by its boundary data — and is the result that discharges the marker placed in the heat-equation treatment, where the boundary-data form of the principle was promised but not derived.

The Mean Value Property

The Poisson integral formula on the unit disk has a particularly transparent consequence when evaluated at the centre. There the Poisson kernel \(P_r(\psi) = \frac{1}{2\pi}\frac{1 - r^2}{1 - 2r\cos\psi + r^2}\) reduces, at \(r = 0\), to the constant \(\frac{1}{2\pi}\), and the integral formula collapses to a pure arithmetic average. The same conclusion, transferred to a disk of arbitrary centre and radius by translation and dilation (both of which preserve harmonicity), yields the mean value property in its general form.

The Maximum Principle

The mean value property has an immediate qualitative consequence: a harmonic function cannot attain a strict interior maximum on a connected domain. Any point at which the function takes its maximum value would be forced, by the average being equal to the maximum, to be a value that every nearby point on a surrounding circle also takes; iterating this observation spreads the maximum value through the connected component. The precise statement is the elliptic counterpart of the heat-equation pointwise bound, with the initial-time data slice replaced by the spatial boundary.

The maximum principle yields uniqueness for the Dirichlet problem as an immediate corollary, paralleling the role energy conservation played for the wave equation. Suppose \(u_1\) and \(u_2\) are two solutions of the Dirichlet problem on a bounded \(\Omega\) with the same continuous boundary datum \(f\) on \(\partial\Omega\). Their difference \(w = u_1 - u_2\) is harmonic on \(\Omega\), continuous on \(\overline{\Omega}\), and satisfies \(w|_{\partial\Omega} = 0\). The maximum principle gives \(\max_{\overline{\Omega}} w = \max_{\partial\Omega} w = 0\) and \(\min_{\overline{\Omega}} w = \min_{\partial\Omega} w = 0\), so \(w \equiv 0\) on \(\overline{\Omega}\); hence \(u_1 = u_2\). The Poisson integral formulas constructed earlier therefore give the solution of the Dirichlet problem on the disk and on the half-plane respectively, not merely a solution. This is the elliptic analogue of the energy-based uniqueness established for the wave equation, with the controlling object replaced from a global quadratic invariant to a pointwise bound by the boundary data.

The maximum principle, in the boundary-data form just proved, completes a pattern initiated in the treatment of the heat equation. There the principle took the form \(\inf f \leq u(x, t) \leq \sup f\), bounding the solution pointwise by the initial data. The wave equation admits no such pointwise estimate; its quantitative invariant is the global mechanical energy. The Laplace equation here recovers a pointwise bound, but with the controlling data shifted from the time-zero slice to the spatial boundary. The three classical PDEs realize three answers to the question "what controls the solution?": initial data for the heat equation, global energy for the wave equation, boundary data for the Laplace equation. The structural reason for this three-way pattern — and its connection to the algebraic classification of second-order linear PDEs — is the subject of the closing section.

The Algebraic Origin of the Trichotomy

The three classical PDEs treated on this site share a common algebraic origin. A general second-order linear PDE in two variables has principal part \(A u_{xx} + B u_{xy} + C u_{yy}\), and its qualitative behaviour at each point is governed by the discriminant \(B^2 - 4AC\) — exactly the same expression that classifies plane conic sections \(A x^2 + B xy + C y^2 + \ldots = 0\). Negative discriminant gives ellipses and the elliptic Laplace equation \(\partial_{xx} u + \partial_{yy} u = 0\), whose discriminant equals \(-4\); zero discriminant gives parabolas and the parabolic heat equation \(\partial_t u - \partial_{xx} u = 0\), whose principal part \(-\partial_{xx} u\) yields discriminant \(0\); positive discriminant gives hyperbolas and the hyperbolic wave equation \(\partial_{tt} u - c^2\, \partial_{xx} u = 0\), whose discriminant equals \(4c^2\).

The discriminant is the surface manifestation of a deeper algebraic fact. The principal part \(A u_{xx} + B u_{xy} + C u_{yy}\) is the quadratic form associated with the symmetric coefficient matrix \[ M \;=\; \begin{pmatrix} A & B/2 \\ B/2 & C \end{pmatrix}, \] and \(B^2 - 4AC = -4 \det M\): the discriminant records the sign of the determinant of \(M\), which in turn records the sign pattern of its eigenvalues. Negative discriminant means \(\det M > 0\), so the two eigenvalues of \(M\) share the same sign — the matrix is definite, and the PDE is elliptic; positive discriminant means \(\det M < 0\), so the eigenvalues have opposite signs — the matrix is indefinite, and the PDE is hyperbolic; zero discriminant means \(\det M = 0\), so one eigenvalue vanishes — the matrix is degenerate, and the PDE is parabolic. The Laplace matrix \(\operatorname{diag}(1, 1)\) has eigenvalues \((+1, +1)\); the wave matrix \(\operatorname{diag}(-c^2, 1)\) (with the second variable identified with \(t\)) has eigenvalues \((-c^2, +1)\); the heat matrix \(\operatorname{diag}(-1, 0)\) has eigenvalues \((-1, 0)\). The shared names — elliptic, parabolic, hyperbolic — are not metaphor: they record the eigenvalue sign pattern of a \(2 \times 2\) symmetric matrix, and the same classification extends verbatim to higher-dimensional PDEs, where the principal symbol is a larger symmetric matrix and the discriminant \(B^2 - 4AC\) gives way to a richer signature.

The discriminant has a second meaning, geometric rather than algebraic. It governs the number of real solutions of the characteristic equation \(A(dy/dx)^2 - B(dy/dx) + C = 0\) for the level curves \(\phi(x, y) = \text{const}\) along which information propagates. For the wave equation, identifying the second variable \(y\) with time \(t\), the principal part \(\partial_{tt} u - c^2 \partial_{xx} u\) gives \(A = -c^2\), \(B = 0\), \(C = 1\), and the characteristic equation \(-c^2(dt/dx)^2 + 1 = 0\) has the two real roots \(dt/dx = \pm 1/c\), or equivalently \(dx/dt = \pm c\). These define two real characteristic families — the lines \(x \pm ct = \text{const}\) along which d'Alembert's formula carried the travelling waves of the wave page. For the heat equation, with principal part \(-\partial_{xx} u\) (the first-order term \(\partial_t u\) does not enter the classification) and degenerate quadratic in \(dt/dx\), the characteristic condition reduces to \(\phi_x = 0\), giving a single family of characteristics \(t = \text{const}\): information at any instant is coupled across the entire spatial axis at once — the infinite propagation speed of diffusion, recovered here as a geometric consequence of the parabolic signature. For the Laplace equation, no real characteristics exist at all; the equation \(\phi_x^2 + \phi_y^2 = 0\) admits only the trivial solution over the reals. There is no curve along which information flows, no time direction, no notion of forward evolution. The absence of characteristics has a second, qualitative consequence: singularities have no channel along which to propagate from the boundary into the interior. A discontinuous boundary datum on the wave equation persists as a discontinuity travelling along the characteristic line, but a merely continuous datum imposed on the Laplace equation produces an interior solution that is smooth — in fact analytic — at every interior point, no matter how irregular the boundary values are away from those points. Elliptic regularity is, in this sense, the structural shadow of no-characteristic ellipticity: with no curve to carry singularities inward, the interior cannot help but be smooth. The absence of a time axis announced at the start of this page — the orthogonality of Laplace to the heat–wave dichotomy — is not a modelling choice but the consequence of an absent characteristic structure.

The three control mechanisms catalogued in the maximum-principle section — initial data for heat, global energy for wave, boundary data for Laplace — are three faces of this single algebraic fact: the eigenvalue sign pattern of the principal-part matrix. Two-characteristic propagation makes the wave equation conserve mechanical energy along its trajectories; degenerate-characteristic diffusion makes heat bounded pointwise by initial data through the maximum principle on the real line; no-characteristic ellipticity makes Laplace bounded pointwise by boundary data through the boundary maximum principle. The three classical PDEs realize three paradigms of information flow — propagation, dissipation, harmonic interpolation — and each admits a discrete dual treated elsewhere on this site: the graph Laplacian as the discrete Dirichlet energy operator, and the factorization \(L = BB^\top\) in terms of the incidence matrix of a graph as the discrete counterpart of \(\Delta = \nabla \cdot \nabla\).

The parallel running through the trichotomy is summarized below.

One algebraic fact — the eigenvalue sign pattern of a \(2 \times 2\) symmetric matrix — produces, across its three possible patterns \((+,+)\), \((+,-)\), \((+,0)\), the analytic, geometric, and discrete structure of the entire trichotomy.