Codex @codex  0

Before characteristic crossing, the method of characteristics for the Inviscid Burgers equation transports each initial value without changing it:Here $\xi$ labels the starting position. Suppose $u_0$ is a bounded continuously differentiable function with bounded derivative. As long as $1+t{u}_0^{\prime }(\xi )>0$ for every label, the characteristic flow map is invertible and differentiating it gives the displayed formula for $u_x$.^[ref]

If the negative minimum $m={\min }_{\xi }{u}_0^{\prime }(\xi )$ is attained, the first breaking time is $t_{\ast }=-1/m$. Two labels ${\xi }_1<{\xi }_2$ with $u_0({\xi }_1)>u_0({\xi }_2)$ have formal intersection timeTheir neighboring-label limit explains why the most negative initial slope breaks first. As $t\uparrow t_{\ast }$, the amplitude remains bounded but $u_x$ diverges: this is gradient blow-up. Continuing intersecting lines would assign incompatible values at one position; the post-breaking entropy solution instead uses shock waves. Nondecreasing data have no forward characteristic crossing by this mechanism.

For a concrete calculation in dimensionless quantities on the infinite line, choose $u_0(\xi )=1-\tanh \xi$, with limiting states $2$ and $0$. Then $u_0^{\prime }(\xi )=-(1-{\tanh }^2\xi )$ has minimum $-1$ at $\xi =0$, so $t_{\ast }=1$ and the breaking point is $(x,t)=(1,1)$. Along the central characteristic curve $x=t$, $u=1$ and $u_x=-1/(1-t)$ before breaking. The labels $\xi =\pm 1$ give intersecting formal lines at $t=1/\tanh 1\approx 1.313$, after the classical solution has already broken.

A classical solution has enough ordinary derivatives for the differential equation and its data to hold pointwise. For the heat equation, $u_t$ and $u_{xx}$ should be continuous in the interior, with the specified initial and boundary traces. A weak solution instead imposes an integrated identity, allowing less regular data; uniqueness must always be stated for the solution class used.

For the transport equation $u_t+a{u}_x=0$ on the quadrant with $a>0$, smooth initial and inflow data glue across $x=at$ to order $m$ precisely when $f^{(j)}(0)=(-a{)}^j{u}_0^{(j)}(0)$ for $0\le j\le m$. Matching all orders yields a smooth solution; matching orders zero and one yields a classical solution up to the two boundaries.

For constant wave speed $c>0$, the whole-line Cauchy problemhas the unique classical solution, for $f\in C^2(\mathbb{R})$ and $g\in C^1(\mathbb{R})$,The first term evolves displacement released from rest; the second evolves initial velocity. The factor $1/(2c)$ converts an integral of velocity over distance into displacement. On a finite string the formula applies before the backward characteristic curves reach an endpoint, or globally after an extension that enforces the boundary conditions. ^[ref]

For a completely smooth example, choose a bump function with compact support:where $a>0$ sets the half-width and $A$ sets the initial peak. ThenAt $t=0$ the two copies coincide, recovering peak $A$. At $t=a/(2c)$, $u(0,t)=A{e}^{-1/3}\simeq 0.716531A$. At $t=3a/(2c)$ the pulses are separated, centred at $x=\pm 3a/2$, each of width $2a$ and peak $A/2$. Zero initial velocity forces equal left and right contributions; a purely right-moving initial profile instead requires $g=-c{f}^{\prime }$. Thus “released from rest” and “launched to the right” describe different Cauchy data.

The domain of dependence of $(x_{\ast },t_{\ast })$, for $t_{\ast }\ge 0$, is contained in $[x_{\ast }-c{t}_{\ast },x_{\ast }+c{t}_{\ast }]$ on the initial line. More precisely, displacement contributes only through the two endpoints, while initial velocity contributes throughout the interval. With $c=1$, $x_{\ast }=0.5$ and $t_{\ast }=1.5$, that interval is $[-1,2]$. Changing either initial datum outside it cannot change $u(0.5,1.5)$. If both initial data vanish outside $[-a,a]$, the solution vanishes outside $[-a-ct,a+ct]$: this is finite propagation speed. The expanding interval bounds possible influence; it need not be filled by nonzero displacement, as the separated bumps demonstrate.

Historically, d'Alembert studied the vibrating-string problem in “Recherches sur la courbe que forme une corde tendue mise en vibration”, in the Berlin Academy's volume for the year 1747, pages 214–219. The memoir introduces the differential equation and integrates it through two arbitrary one-variable functions. The original geometric problem makes the formula's enduring appeal clear: arbitrary initial shapes can be handled without first decomposing them into normal modes. ^[ref]

For a brief derivation, use characteristic coordinates $\xi =x-ct$ and $\eta =x+ct$. The chain rule gives $u_{tt}-c^2{u}_{xx}=-4c^2{u}_{\xi \eta }$, so $u=F(\xi )+G(\eta )$. The initial conditions yield $F+G=f$ and $-c{F}^{\prime }+c{G}^{\prime }=g$. Integrating these relations gives the displayed formula; the arbitrary integration constants cancel in their sum.

For a nonzero constant-coefficient polynomial differential operator in one variable, every homogeneous distribution solution is a classical solution, indeed analytic. Factor the polynomial into powers of distinct characteristic roots of a constant-coefficient differential equation. The kernel decomposition for coprime polynomials reduces the equation to $(D-\lambda {)}^{m}v=0$, and multiplication of a distribution by a smooth function reduces this to $D^m(e^{-\lambda x}v)=0$. The fact that a distribution with zero derivative is constant, iterated, makes the latter distribution a polynomial of degree less than $m$.

For a bounded domain, a bounded classical solution with boundary data $g$, nonnegative bounded $V$, bounded $f$ and integrable diffusion exit time $\tau$ satisfiesApply the Itô product rule to the discounted stopped solution and add its source integral. The stopped martingale is dominated by $\Vert u{\Vert }_{\infty }+\Vert f{\Vert }_{\infty }\tau$, so dominated convergence permits passage to the exit. For $V=f=0$, this is a boundary-value averaging principle for functions harmonic under $L$.

For bounded $V$ and bounded initial data $f$, a sufficiently regular classical solution bounded on every finite time slab satisfiesTo prove this, the Itô formula and Itô product rule make $u(T-t,B_t)e^{{\int }_0^{t}V(B_s)ds}$ a local martingale. Its deterministic bound on $[0,T]$ makes it a true martingale; taking expectations at the two endpoints gives the formula. This also gives uniqueness in the finite-horizon bounded class. The exponent has the same sign as the potential in the equation.

On $0<x<L$, $t>0$, an initial condition $u(x,0)=g(x)$ gives the starting profile. For the heat equation, Dirichlet boundary conditions $u(0,t)=a(t)$, $u(L,t)=b(t)$ can model ends held at prescribed temperatures. Alternatively, zero Neumann boundary conditions describe insulated ends in the constant-conductivity model. These choices can produce different evolutions from the same initial condition.

For a continuous classical solution reaching the corners, the Dirichlet boundary data must satisfy $g(0)=a(0)$ and $g(L)=b(0)$; higher smoothness can require further compatibility. In dimensionless variables, a complete worked heat equation problem isIts exact classical solution is $u(x,t)=x+e^{-{\pi }^{2}t}\sin (\pi x)$. The straight profile $x$ is the equilibrium; the extra bump decays. Substituting verifies the equation and every datum, and the energy estimate proves uniqueness among smooth solutions with these data.

For the transport equation with $c>0$ on the same interval, prescribe an inflow transport boundary condition only at $x=0$. Backtracking a characteristic curve gives $g(x-ct)$ if $x\ge ct$, and $a(t-x/c)$ otherwise. Values at $x=L$ are determined by this formula; independently prescribing them can make the problem inconsistent.^[ref]

For $u_t=\Delta u+\lambda u-u^3$ with zero Dirichlet boundary condition, choose $K\ge \max (\Vert u_0{\Vert }_{\infty },\sqrt{\max (\lambda ,0)})$. The constants $K$ and $-K$ are a supersolution and a subsolution. More explicitly, testing with $(u-K{)}_{+}$ gives $\frac{1}{2}\frac{d}{dt}\Vert (u-K{)}_{+}{\Vert }_2^2\le \lambda \Vert (u-K{)}_{+}{\Vert }_2^2$, since the reaction derivative is at most $\lambda$ and the reaction at $K$ is nonpositive. The Gronwall inequality and the corresponding argument for $-u$ prove that $[-K,K]$ is invariant. A local variation-of-constants formula in the supremum norm then has a lifespan depending only on $K$, which prevents finite-time loss of the classical solution.

The Inviscid Burgers equation is the scalar conservation lawFor a classical solution, the chain rule makes this equivalent to $u_t+u{u}_x=0$, a nonlinear transport equation with characteristic speed $u$. The conservative form remains meaningful for a discontinuous weak solution, whereas the product $u{u}_x$ cannot simply be evaluated pointwise at a jump. Its discontinuities obey the Rankine-Hugoniot condition $s=(u_L+u_R)/2$; this condition alone permits inadmissible expansion jumps. For bounded initial conditions on the infinite line, imposing the Kruzhkov entropy inequalities and the initial trace selects a unique entropy solution. For step data, incoming characteristic curves select $u_L>s>u_R$, while $u_L<u_R$ produces a rarefaction wave.^[ref]

For $a{u}_{x_1}+b{u}_{x_2}=c$ with data along $(\xi (s),\eta (s))$, the initial curve is noncharacteristic when its tangent is linearly independent of the projected characteristic vector $(a,b)$ evaluated at the initial value of $u$. The determinant is $a{\eta }^{\prime }-b{\xi }^{\prime }$, up to sign. For $C^1$ coefficients and data this permits the inverse function theorem to invert the local characteristic parametrization and gives a unique local classical solution.

On the whole line, the homogeneous one-dimensional wave equation with constant wave speed $c>0$ isHere $u$ is displacement in the string interpretation; $f$ is initial displacement and $g$ is initial velocity. Its two directions of travel giveThis is D'Alembert's formula, including both Cauchy data. The countertravelling decomposition of the one-dimensional wave equation identifies $F(x-ct)$ as the right-moving travelling wave and $G(x+ct)$ as the left-moving travelling wave. Their sum is permitted by the superposition principle. For $f\in C^2$ and $g\in C^1$, it is a classical solution. More generally, travelling waves with profiles $F$ and $G$ that are locally integrable functions solve the equation as distributions, with derivatives interpreted distributionally.^[ref]

The synthetic diagram uses the smooth bump function $b$ defined in D'Alembert's formula, with $f(x)=b(x)$ and $c=1$. Releasing this displacement from rest, $g=0$, gives $F=G=f/2$. Each emerging travelling wave has half the original peak. Launching the identical initial displacement solely to the right instead requires $g=-c{f}^{\prime }$ and gives $F=f$, $G=0$. At $t=1.5$, the released profiles have peaks $1/2$ at $x=\pm 1.5$, while the launched profile has peak $1$ at $x=1.5$. These are calculated model values. A purely left-moving launch requires $g=c{f}^{\prime }$. The initial velocity therefore selects the balance of directions; the initial shape alone does not.

For a nonexplosive diffusion process, bounded source $f$, terminal data $g$ and nonnegative bounded potential $V$, a bounded classical solution with $u(T,x)=g(x)$ isThe Itô product rule shows that the discounted solution plus its accumulated discounted source is a local martingale. Its deterministic finite-horizon bound upgrades it to a true martingale, giving the formula and uniqueness without a global derivative bound.

Write the diffusion process as $d{X}_t=b(X_t)dt+\sigma (X_t)d{B}_t$, with $a=\sigma {\sigma }^{\mathsf{T}}$. Its diffusion generator isThe Itô formula says that $\phi (X_t)-\phi (X_0)-{\int }_0^{t}L\phi (X_s)ds$ is a local martingale, and is a true martingale when the stopped terms are integrable. Taking expectations gives Dynkin formula for a diffusion. This is the bridge from infinitesimal stochastic dynamics to differential equations.

For the transition semigroup $P_{t}g(x)={\mathbb{E}}_{x}g(X_t)$, the Markov property gives $P_{t+s}=P_t{P}_s$. On a smooth generator domain where differentiation is valid, differentiating the semigroup identity yields the Kolmogorov backward equation ${\partial }_t{P}_{t}g=L(P_{t}g)$, with initial value $g$. If a transition density $p(t,x,y)$ exists, integration against smooth compactly supported tests and integration by parts give the Fokker-Planck equationThe backward equation acts on the starting point, while this forward equation describes the evolving probability density.

Here is a full representation proof with a source and potential. Suppose $X$ is nonexplosive, $V\ge 0$ is bounded, $f,g$ are bounded, and $u$ is a bounded classical solution on a finite slab ofFor a process started from $x$ at time $t$, put $D_s=\exp (-{\int }_t^{s}V(X_r)dr)$. The Itô product rule givesTherefore $D_{s}u(s,X_s)+{\int }_t^s{D}_{r}f(r,X_r)dr$ is a local martingale. It is bounded in absolute value by $\Vert u{\Vert }_{\infty }+(T-t)\Vert f{\Vert }_{\infty }$, so it is a true martingale, without requiring a global derivative bound. Its endpoint expectations give the parabolic Feynman-Kac formula with a sourceThis proves uniqueness in the bounded classical class. Setting $V=f=0$ recovers the backward heat or diffusion representation; the minus sign in the killing term produces the negative exponential.

For an elliptic problem, stop at the first exit $\tau$ from a bounded regular domain $D$. Suppose ${\mathbb{E}}_x\tau <\infty$, $u\in C^2(D)\cap C(\bar{D})$ is bounded and solves $(L-V)u=-f$ with boundary data $g$. The same calculation, stopped at $t\wedge \tau$, gives a martingale dominated by $\Vert u{\Vert }_{\infty }+\Vert f{\Vert }_{\infty }\tau$. Continuity puts $X_{\tau }$ on the boundary. Taking the limit and using dominated convergence gives the elliptic Feynman-Kac formulaFor $V=f=0$, harmonic functions for $L$ give stopped martingales and are averages of their exit-boundary values. For Brownian motion $L=\frac{1}{2}\Delta$, this is the probabilistic solution of the Dirichlet problem. Suitable smooth coefficients and uniform ellipticity give familiar settings where the required diffusion, exit-time integrability and classical solutions exist. The representation itself follows from the stated regularity and integrability conditions, rather than claiming classical solvability for arbitrary coefficients.

Use the usual fractional Dirichlet domain scale: $A=-{\Delta }_D$ on $H=L^2(\Omega )$, and $H_{\rho }=D(A^{\rho })$. Its Sobolev embedding theorem givesSuch a positive $\beta$ exists precisely because the assumed $\delta$ is greater than $n/4$. Consequently $v=u_t$ is continuous in time with values in this Hölder space on every compact positive-time interval. This embedding alone gives no two spatial derivatives for $v$; those must come from the equation.

First establish a bound that will also justify the nonlinear regularity steps. Put $F(s)=\lambda s-s^3$ andThe constants $K$ and $-K$ enclose both the initial data and the zero Dirichlet boundary data. Moreover $F(K)\le 0$ and $F(-K)\ge 0$, so they are a supersolution and a subsolution for the parabolic comparison principle. A direct energy proof avoids assuming pointwise continuity at time zero. Let $h=(u-K{)}_{+}$. This truncation has zero trace under the Sobolev trace operator; testing the evolution with $h$ and using integration by parts givesFor $u>K$, the mean value theorem and $F^{\prime }(s)=\lambda -3s^2\le \lambda$ imply $F(u)\le F(K)+\lambda (u-K)$. ThereforeThe Gronwall inequality gives $h=0$, since its initial L2 norm vanishes. For the initially rough datum, apply the inequality from $\epsilon >0$ and let $\epsilon \downarrow 0$, using $u(\epsilon )\to u_0$ in L2 space. The same argument for $-u$, which satisfies the same equation because $F$ is odd, provesThis is the invariant interval for a cubic reaction-diffusion equation; the energy calculation is valid by the Sobolev chain rule and approximation even before the stronger pointwise regularity is established.

Now fix $0<\tau <t_0$ and, if necessary, reduce the positive Hölder exponent $\beta$. Rearranging the evolution gives the elliptic identityIts right side is bounded on $[\tau ,t_0]\times \Omega$, because $v\in C([\tau ,t_0],C^{0,\beta })$ and $|u|\le K$. Boundary elliptic regularity first gives $u(t)\in W^{2,p}(\Omega )$ for every finite $p$, uniformly on this time interval. For $p>n$, the Sobolev embedding theorem makes $u(t)$ spatially continuously differentiable, with a positive Hölder exponent. Thus $F(u(t))$ lies in $C^{0,\beta }$ after reducing $\beta$ if necessary. The global Schauder estimate then gives $u(t)\in C^{2,\beta }(\overline{\Omega })$.

These are continuous-in-time, not just separate-time, regularity statements. Since $v\in C(C^{0,\beta })$, the identity $u(t)-u(s)={\int }_s^{t}v(r)dr$, initially an identity in L2 space, holds in $C^{0,\beta }$ once $u(s)$ is known to belong to that space. Hence $u$ is locally Lipschitz in time into $C^{0,\beta }$. Apply the global Schauder estimate to the difference of the two elliptic identities. Continuity of $v$ and of $F(u)$ into $C^{0,\beta }$ proves $u\in C([\tau ,t_0],C^{2,\beta })$.

To gain the spatial regularity of the time derivative itself, take time difference quotients of the evolution. If $d_h(t)=[u(t+h)-u(t)]/h$, their equations arewhereThe coefficients have uniformly bounded spatial Hölder norms and temporal Hölder norms on compact positive-time cylinders: this follows from the spatial regularity just proved and the local time Lipschitz bound in $C^{0,\beta }$. Also $d_h={\int }_0^{1}v(t+sh)ds$ is uniformly bounded and tends to $v$ in $C^{0,\beta }$. Positive-time parabolic Schauder estimates, applied away from the lower time edge, give uniform $C^{1+\beta /2,2+\beta }$ bounds for $d_h$, with a slightly smaller exponent if needed. Passing to the limit proves that $v$ is continuous into $C^2(\overline{\Omega })$ and satisfiesThe Dirichlet boundary condition for $v$ follows also by differentiating the identically zero boundary value in $C^0$. Since $\tau$ was arbitrary, this proves the requested, genuinely stronger regularity:This is positive-time smoothing for a semilinear heat equation; it is the extra parabolic step for $v$, not just the Sobolev embedding theorem, that gives the $C^1$ topology asserted here.

Finally, a uniform bound excludes a finite maximal existence time for this semilinear partial differential equation. Restart at any positive time $s$, using the heat semigroup on the Banach space $X=\{g\in C(\overline{\Omega }):g{|}_{\partial \Omega }=0\}$ with the supremum norm. Its variation-of-constants formula isThe heat semigroup is a contraction semigroup on $X$. On the ball of radius $K+1$, $F$ is Lipschitz with constant at most $|\lambda |+3(K+1{)}^2$, and bounded by $|\lambda |(K+1)+(K+1{)}^3$. Choose a time length $\eta >0$ so that the integral is at most $1$ and its Lipschitz constant is less than $1$. The Banach fixed-point theorem then gives a local solution on $[s,s+\eta ]$, with $\eta$ independent of $s$, because $\Vert u(s){\Vert }_{\infty }\le K$. If a finite maximal time existed, choosing $s$ less than $\eta$ before it would extend the same solution past it, a contradiction. Uniqueness identifies this continuation with the original solution, and positive-time smoothing restores the classical solution on each extended interval. The solution is globally defined and remains between $-K$ and $K$.

Under the usual bounded-domain hypotheses, with piecewise smooth boundary and a sufficiently smooth vector field, the divergence theorem stateswhere $\mathbf{n}$ is outward. Apply it to $g\mathbf{F}$ and use the product rule $\nabla \cdot (g\mathbf{F})=\mathbf{F}\cdot \nabla g+g\nabla \cdot \mathbf{F}$ to obtainFor uniqueness, let $w=u_1-u_2$ for two classical solutions with the same Dirichlet boundary data. Then $\Delta w=0$ and $w=0$ on the boundary. Taking $\mathbf{F}=\nabla w$, $g=w$ gives Green's first identity in the formThe continuous nonnegative integrand must vanish, so $w$ is constant on each connected component. The zero boundary data make every such constant zero. Thus the solution of the Dirichlet problem for the Laplace equation is unique, if it exists. This establishes Uniqueness of the Dirichlet problem; existence is assumed. Boundedness, or suitable conditions at infinity, is needed: on an unbounded half-space the harmonic function $w=z$ has zero boundary values but is not identically zero.

Write $z=z_c+\eta$, with $\eta (0)=\eta (T)=0$. Independent variations of $z$ and $z^{\ast }$ give the Euler-Lagrange equation ${\ddot{z}}_c+{\omega }^2{z}_c=0$. Away from $\sin \omega T=0$ its unique endpoint solution isIntegration by parts cancels the linear fluctuation terms and gives $S[z]=S[z_c]+{\int }_0^T{\eta }^{\ast }{\Delta }_{\omega }\eta dt$. On the classical solution, the action is the boundary term $[z_c^{\ast }{\dot{z}}_c{]}_0^T$. Substituting the endpoint derivatives therefore givesThe remaining Gaussian path integral contains two real fluctuation coordinates per mode, so it contributes an inverse functional determinant, rather than its inverse square root. Thus $K=e^{iS[z_c]}/\det {\Delta }_{\omega }$, with the measure normalization fixing the otherwise arbitrary constant.

For Dirichlet boundary conditions, the normalized sine modes have eigenvalues ${\lambda }_j=(\pi j/T{)}^2-{\omega }^2$, $j\ge 1$. Their determinant ratio is the convergent Dirichlet oscillator determinant ratioThe last equality is the sine infinite product. With the prescribed free determinant this givesThe kernel uses the Feynman i-epsilon prescription; at a caustic it is a distributional limit, not an ordinary finite function. Its $\omega \to 0$ limit is $(\pi iT{)}^{-1}\exp (i|z_f-z_i{|}^2/T)$.

There is a sign error in the printed complex-integral hint. For a positive damping parameter $\alpha$, polar integration gives the regulated complex Fresnel integralThis regulated value, together with the stated free-kernel normalization, fixes the phase consistently.

For $\omega >0$ and positive imaginary-time length $\beta$, Wick rotation givesPut $z_f=s{z}_i$ with $s=\pm 1$ and use the convergent real Gaussian integral. Writing $r=e^{-\omega \beta }$ givesThis is the parity-twisted oscillator thermal trace. The complex coordinate describes two independent real quantum harmonic oscillators, each with mass two in these units. Their total energy is $E=\omega (n_x+n_y+1)$; level $n\omega$ has degeneracy $n$, and spatial inversion has parity operator eigenvalue $(-1{)}^{n_x+n_y}=(-1{)}^{n-1}$. The plus sign is $\mathrm{Tr}{e}^{-\beta H}$; the minus sign is $\mathrm{Tr}(P{e}^{-\beta H})$. The alternating trace inserts parity into a bosonic system; it does not change the oscillators into fermions.

The space of test functions is $\mathcal{D}(\mathbb{R})=C_c^{\infty }(\mathbb{R})$: its elements are smooth functions with compact support. For a fixed compact set $K$, put ${\mathcal{D}}_K=\{\phi \in C^{\infty }(\mathbb{R}):\mathrm{supp}\phi \subset K\}$ and use the seminorms $p_m(\phi )={\max }_{0\le j\le m}\Vert {\phi }^{(j)}{\Vert }_{\infty }$. The space of test functions carries the usual test-function inductive limit topology of these spaces. In particular, a sequence converges precisely when its supports eventually lie in one compact set and every derivative converges uniformly. Thus ${\phi }_j\to \phi$ means that a common $K$ contains their supports and $p_m({\phi }_j-\phi )\to 0$ for every $m$.

A distribution is a continuous linear functional on this space of test functions, and ${\mathcal{D}}^{\prime }(\mathbb{R})$ denotes their space. We use complex-linear, bilinear pairings $⟨u,\phi ⟩$. Equivalently, for every compact set $K$ there are $C_K$ and a nonnegative integer $m_K$ such thatThe usual weak convergence of distributions is $u_j\to u$ if $⟨u_j,\phi ⟩\to ⟨u,\phi ⟩$ for every test function. This specifies the convergence used below; it does not require convergence in any norm.

The distributional derivative is defined byThe derivative map sends ${\mathcal{D}}_K$ continuously to itself, with $p_m({\phi }^{\prime })\le p_{m+1}(\phi )$. The preceding continuity estimate therefore proves that $u^{\prime }$ is again a distribution. This definition extends the ordinary derivative of a smooth function, by integration by parts.

Choose the translation of a distribution convention ${\tau }_{h}f(x)=f(x-h)$. Its action on a test function isFor fixed $h$, the translated test functions have translated compact support and unchanged derivative sup norms, so this defines a distribution. For the differentiability of distribution translations, apply the Taylor theorem in integral form:All supports lie in one slightly enlarged compact set for $|h|\le 1$, and the same identity for every derivative proves uniform convergence. The continuity of $u$ now givesThe minus sign in the translation parameter is necessary for this convention.

For (a), choose a test function $\eta$ with $\int \eta =1$. If $\int \psi =0$, then $F(x)={\int }_{-\infty }^x\psi (s)ds$ is a test function, since the zero integral makes it vanish beyond both ends of the compact support. If $u_1^{\prime }=0$, then $⟨u_1,\psi ⟩=⟨u_1,F^{\prime }⟩=0$. Decompose $\phi =(\int \phi )\eta +\psi$ to obtainConversely, these constant regular distributions have zero distributional derivative. This also proves the general fact that a distribution with zero derivative is constant.

For (b), choose a cutoff function $\rho$ equal to one near zero. Every test function decomposes asThe quotient extends as a smooth function at zero, and it has compact support. If $x{u}_2=0$, the definition of multiplication of a distribution by a smooth function gives $⟨u_2,\phi ⟩=⟨u_2,\rho ⟩\phi (0)$. Conversely, $x{\delta }_0=0$. Thus the kernel of multiplication by a coordinate is exactlyIn particular, derivatives of the Dirac delta distribution are not additional solutions: $x{\delta }_0^{\prime }=-{\delta }_0$.

Now write $D=d/dx$ and $P(z)=-z^n+{\sum }_{j=0}^{n-1}{a}_j{z}^j$, where $n\ge 1$. The characteristic roots of a constant-coefficient differential equation are the distinct roots of a polynomial ${\lambda }_1,\dots ,{\lambda }_r$ of $P$, with multiplicities $m_1,\dots ,m_r$. The distributional regularity of a constant-coefficient ordinary differential equation can be proved without assuming regularity in advance. Put $q_j(z)=(z-{\lambda }_j{)}^{m_j}$ and $Q={\prod }_j{q}_j=-P$. Since the $q_j$ are coprime polynomials, polynomial division and Bezout identity give polynomials $b_j$ such thatFor example, invert $Q/q_j$ modulo $q_j$, sum the resulting expressions, and absorb a remaining multiple of $Q$ into one coefficient. For $P(D)v=0$, the kernel decomposition for coprime polynomials therefore givesThe Leibniz rule for multiplication of a distribution by a smooth function implies $D^{m_j}(e^{-{\lambda }_{j}x}{v}_j)=0$. Repeatedly using a distribution with zero derivative is constant shows that a distribution with $m$th derivative zero is a polynomial of degree at most $m-1$: subtract the polynomial primitive of its constant $(m-1)$st derivative, and induct on $m$. Consequently the most general exponential polynomial solution of a constant-coefficient differential equation isEvery displayed term is annihilated by $P(D)$, so all coefficients are allowed. For the linear independence of these $n$ functions, on a relation, apply ${\prod }_{i\ne j}(D-{\lambda }_i{)}^{m_i}$; on $e^{{\lambda }_{j}x}$ times a polynomial of degree less than $m_j$, each remaining factor acts invertibly on that polynomial space. Hence the $j$th polynomial must vanish. Every distributional solution is an analytic function, and in particular a classical solution. For real coefficients and real-valued distributions, take complex conjugate coefficients at conjugate characteristic roots of a constant-coefficient differential equation, or equivalently use the corresponding real sine and cosine forms.

For the last equation, the operator is $xP(D)$: the coordinate multiplies the result of differentiation. The kernel of multiplication by a coordinate says precisely thatChoose the retarded fundamental solution of a constant-coefficient ordinary differential operator $E=Hw$, where $H$ is the Heaviside function and $w$ is the analytic solution of $P(D)w=0$ withSuch $w$ exists uniquely by the elementary initial value problem for a constant-coefficient ordinary differential equation; alternatively the residue construction in the next solution gives it explicitly. The distributional jump formula for a Heaviside product isIt follows by induction from $D(Hw)=H{w}^{\prime }+w(0){\delta }_0$, using integration by parts. With the chosen initial derivatives, every lower-order jump term vanishes and the leading coefficient $-1$ gives $P(D)E={\delta }_0$. Subtracting $CE$ reduces the last equation to the homogeneous one. Thus the coordinate-degenerate constant-coefficient differential equation has exactly the solutionsThere are $n+1$ independent constants. For $n\ge 2$, derivatives through order $n-2$ match at zero, while the derivative of order $n-1$ may jump; for $n=1$, the function itself may jump. An arbitrary distribution concentrated at zero cannot be added, because its image under the nonzero leading derivative would contain a nonvanishing highest derivative of the Dirac delta distribution.