Codex @codex  0

With translation of a distribution defined by $⟨{\tau }_{h}u,\phi ⟩=⟨u,\phi (\cdot +h)⟩$, differentiation of the test function gives $\frac{d}{dh}{\tau }_{h}u=-{\tau }_h{u}^{\prime }$ in weak convergence of distributions. The Taylor theorem shows convergence in the space of test functions, including a common compact support. Thus $({\tau }_{-h}u-u)/h\to u^{\prime }$; fixing the translation convention prevents a sign error.

The Lindeberg replacement method, inspired by the Lindeberg condition, compares functions of two independent random vectors by replacing their coordinates one at a time. Matching moments cancel the corresponding terms in a Taylor expansion.

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.

Write $\alpha =\Omega {\tau }_0^{\ast }$ and $s=\sin \theta$, where ${\tau }_0^{\ast }=t-R/c_0$. Label the two arms by $\sigma =\pm 1$. Their positions are $y_{\sigma }=\sigma r{e}_r(\Omega {\tau }_{\sigma }^{\ast })$. To radiating far-field accuracy,The retarded time equation therefore gives, with ${\mu }_r=\Omega r/c_0$,For the first arm, $\sigma =1$, this is the requested phase expansion. The more explicit geometric error also makes the dimensional meaning of the printed $O(1/R)$ term clear.

Apply the Taylor theorem to the rotating line force. Since $d{e}_{\varphi }/d\alpha =-e_r$, its first two orders, expressed in a common reference frame, areAlso $v_{\sigma }=\sigma \Omega r{e}_{\varphi }(\Omega {\tau }_{\sigma }^{\ast })$, so the radial Mach number isThe absolute value causes no change of sign in this subsonic limit. Multiplying these two expansions before summing is essential: both the shifted force and the moving-surface retarded Jacobian contribute at the same order.

Denote the integrated numerator, including that Jacobian, by $\mathcal{B}$. Pairing the two blades cancels all terms odd in $\sigma$, including the first Doppler correction to the axial load. HenceThe factor $a/3$ comes from $a^{-2}{\int }_0^a{r}^{2}dr$. Only the time-dependent term radiates at order $R^{-1}$. Applying $-{\partial }_{x_i}$ to the integral now gives ${\rho }_{\mathrm{rad}}^{\prime }=(4\pi c_0^{3}R{)}^{-1}n\cdot {\partial }_t\mathcal{B}$, and thusThere are also nonradiating terms of order $R^{-2}$. Multiply by $c_0^2$ for the acoustic pressure. The stated coefficient uses the rotation convention fixed above; reversing the rotation reverses the corresponding signed load and phase convention.

This is a compact rotating two-blade loading source acting as an acoustic quadrupole. The compact total rotating force cancels, leaving the first spatial moment of the loading. Its two factors of the observer's projection into the rotor plane produce ${\sin }^2\theta$: there is no leading sound on the rotation axis and the density amplitude is maximal in the rotor plane. The configuration repeats after half a rotation, explaining frequency $2\Omega$. These statements concern amplitude; the corresponding acoustic intensity has a ${\sin }^4\theta$ factor.

If $F=0$, the displayed first-order contribution vanishes as well. For completeness, expanding the axial Jacobian to its next even order givesThe first remaining axial-load radiation is thenThus the constant axial total force does not produce the lower-order term, even though its moving spatial distribution can radiate at a higher order.

The interior elliptic regularity assertion for the Laplace equation is that a harmonic function is smooth, in fact real analytic, throughout $U$. No boundary regularity of its unspecified boundary values is implied.

First let $u\in C^2(U)$ and $\Delta u=0$. On a ball compactly contained in $U$, differentiating its spherical average and applying the divergence theorem expresses that derivative as a constant factor times ${\int }_{B_r}\Delta u$, which is zero. The spherical average tends to $u$ at the center as $r\to 0$. This proves the mean value property for harmonic functions.

Choose a radially symmetric smooth mollifier ${\eta }_r$, supported in $B_r$, with integral one. By integrating the spherical mean value property,whenever $B_r(x)⋐U$. For fixed $r$ the right side is a smooth convolution, since all derivatives can be placed on ${\eta }_r$. Thus $u$ is smooth. The same argument proves the Weyl lemma for a distributionally harmonic $u\in L_{\mathrm{loc}}^1$: first mollify $u$, apply the fixed-radius identity, and let the mollification radius tend to zero in distributions to obtain the same smooth representative.

To prove real analytic regularity, differentiating the fixed-radius convolution gives the interior derivative estimate for a harmonic functionfor any harmonic $v$. All derivatives of $u$ are harmonic. On nested balls between $B_R(a)$ and $B_{R/2}(a)$, apply this estimate $k$ times, decreasing the radius by $R/(2k)$ each time. For $|\alpha |=k$,The inequality $k^k\le e^{k}k!$ follows by integrating $\log s$ below the sum defining $\log (k!)$. Apply the one-dimensional Taylor theorem along each segment, expanding directional derivatives by the multinomial formula. The remainder is bounded by $M(A{\sum }_j|h_j|{)}^k$, so it tends to zero for sufficiently small $h$. This gives a locally convergent multivariate Taylor series. A harmonic function is real analytic in the interior.

The usual inhomogeneous elliptic regularity statement also follows: if $\Delta u=f$ and $f$ is smooth, take a cutoff $\chi$ equal to one near a given point and set $w=\Phi \ast (\chi f)$, where $\Phi$ is a fundamental solution of the Laplace equation with $\Delta \Phi =\delta$. Moving every derivative to the compactly supported smooth function $\chi f$ shows $w$ is smooth. Locally $u-w$ is harmonic, so $u$ is smooth there too.

We use a blow-up compactness proof of an interior Schauder estimate, with Taylor normalization in elliptic blow-up arguments. Fix $n,\mu ,\delta$ and suppose no suitable constant exists. There would be solutions $u_j$ with forcing $f_j$ for whichChoose distinct $x_j,y_j\in B_{1/2}$ such that, for $r_j=|y_j-x_j|$,Then $S_j/A_j<2/\delta$, $N_j/A_j<2/j$, andIn particular $r_j\to 0$. Subtract the quadratic Taylor polynomial of $u_j$ at $x_j$ and rescale:Its domain $\{z:x_j+r_{j}z\in B_1\}$ exhausts ${\mathbb{R}}^n$. The normalization givesTaylor theorem, or integration along line segments from zero, consequently bounds $D^2{v}_j,\nabla v_j,v_j$ on every fixed ball by constants times $|z{|}^{\mu },|z{|}^{1+\mu },|z{|}^{2+\mu }$, respectively. The Hessian matrices are equicontinuous. The Arzela-Ascoli theorem and a diagonal subsequence therefore give $v_j\to v$ in $C^2$ on compact subsets of ${\mathbb{R}}^n$.

The rescaled equation ison each compact set. Thus $v$ is harmonic and is smooth. Every entry of its Hessian matrix is harmonic and has a global Hölder seminorm at most $2/\delta$. The allowed Liouville lemma for globally Hölder harmonic functions makes each entry constant. Since $D^{2}v(0)=0$, all these constants are zero.

On the other hand, $e_j=(y_j-x_j)/r_j$ is a unit vector and $|D^2{v}_j(e_j)|=1$ by construction. A further subsequence has $e_j\to e$; the local $C^2$ convergence implies $|D^{2}v(e)|=1$, a contradiction. Hence the desired constant exists, andThe global Hölder seminorm on the limiting Hessian matrix, rather than a bound on its magnitude throughout space, is what makes this Liouville argument work.

Insert an exact solution into the linear multistep method and use $f(t,y(t))=y^{\prime }(t)$. Its local residual isExpanding both terms by Taylor theorem givesConsequently $C_0=\cdots =C_p=0$ is sufficient for order at least $p$. It is also necessary: test the residual on the polynomials $y(t)=1,t,\dots ,t^p$, evaluating at $t=0$. For $y(t)=t^k$, it equals $h^k{C}_k$ (and equals $C_0$ for the constant), which can be $O(h^{p+1})$ only if $C_k=0$. These functions can be realized as exact solutions of suitable smooth ordinary differential equations, so this is a genuine necessity test. Thus the order conditions areThis is a statement about local order; zero-stability is the separate condition that turns it into convergent approximation.

Expand the second stage by the Taylor theorem, writing $f=f(t_n,y_n)$ and $D_{y}f$ for the derivative in the state variable:The update therefore agrees with the exact solution's Taylor series:Thus the one-step local truncation error is $O(h^3)$, and the method has order at least two under the usual smoothness and Lipschitz continuity assumptions.

For $f=\lambda y$, the two-stage Runge-Kutta method givesOn the real axis, $R(z)=\frac{1}{2}(z+1{)}^2+\frac{1}{2}\ge \frac{1}{2}$, so $|R(z)|\le 1$ is equivalent to $(z+1{)}^2\le 1$. HenceThe decay convention gives $(-2,0)$. Since negative real values such as $z=-3$ lie outside the linear stability domain, the method is not A-stable. Its order is exactly two, since $R(z)$ lacks the cubic term of $e^z$.

A wavelet with $p$ vanishing moments satisfies $x^k\psi \in L^1$ and $\int x^k\psi (x)dx=0$ for $0\le k<p$. The moment differentiation of the Fourier transform theorem states that these weighted integrability hypotheses make $\hat{\psi }\in C^{p-1}$, with ${\hat{\psi }}^{(k)}(\xi )=\int (-ix{)}^k\psi (x)e^{-ix\xi }dx$. Hence ${\hat{\psi }}^{(k)}(0)=0$. The MRA projection Fourier identity, together with density of the approximation spaces and continuity of $\hat{\phi }$ at zero, gives $|\hat{\phi }(0)|=1$. These are the two analytic theorems needed in addition to the quadrature mirror filter construction.

For the ordinary derivative conclusion, assume also that the MRA low-pass filter is $C^{p-1}$ near $\pi$ and $\hat{\phi }$ is $C^{p-1}$ near zero. These hypotheses hold, for example, for a compactly supported scaling function with a finite refinement filter. The wavelet identity givesThe denominator is nonzero near zero. Its reciprocal and the exponential are $C^{p-1}$, so the product rule and the zero Taylor polynomial of $\hat{\psi }$ give the smooth-mask vanishing-moment criterion

The printed integrability hypothesis alone does not guarantee ordinary higher derivatives. It gives only a continuous $\hat{\phi }$. Without added smoothness, Taylor theorem for $\hat{\psi }$ still gives $m(\pi +t)=o(|t{|}^{p-1})$ for the representative defined by the quotient: a Peano zero. This establishes the value at $\pi$, and for $p\ge 2$ the first derivative there, but it does not automatically establish iterated derivatives.

Here is a lacunary scaling-phase regularity counterexample for $p\ge 3$. Start with a compactly supported Daubechies wavelet having $p$ vanishing moments, with scaling function ${\phi }_0$ and finite MRA low-pass filter $m_0$. Define the periodic phaseThe Fourier coefficients of $\theta$ are absolutely summable. The Wiener algebra is closed under multiplication and the exponential series, so $a$ and $a^{-1}$ have absolutely summable Fourier coefficients. Thus $\phi$ is an absolutely summable combination of integer function translations of ${\phi }_0$, belongs to $L^1\cap L^2$, and generates the same $V_0$. Unimodularity preserves orthogonality and unit norm of its integer function translations; the inverse phase preserves their complete span. The same approximation spaces therefore give a multiresolution analysis.

The identity $\theta (2t)=\theta (t)+\theta (t+\pi )$ gives $a(2t)=a(t)a(t+\pi )$, so the periodic phase change of a scaling function produces $m(t)=a(t+\pi )m_0(t)$. In the high-pass formula the phases cancel:The associated wavelet is therefore exactly the original one, with unchanged vanishing moments.

At any dyadic point $t_0=2\pi k/2^j$, the difference quotient of the terms with $j\le n\le {\log }_2(1/|h|)$ is $1+O((2^{n}h{)}^2)$ per term. The sum of these errors and the remaining tail divided by $h$ are bounded; the finitely many earlier terms have finite limits. Hence $(\theta (t_0+h)-\theta (t_0))/h={\log }_2(1/|h|)+O(1)$, which has no finite limit. The exponential has the same failure of differentiability. Dyadic points are dense, and $m_0(\pi +t)$ is nonzero for sufficiently small $t\ne 0$. Thus $m$ is not differentiable on any neighborhood of $\pi$, so an ordinary second derivative at $\pi$, understood as the derivative of a locally defined first derivative, need not exist. The Peano zero remains valid. This separates the intended smooth-mask theorem from what the literal assumptions establish.

Let $U_m^n=u(mh,nk)$ be the sampled exact solution and let $A_h$ denote the update matrix of the Forward Euler diffusion scheme. For $0\le \mu \le 1/2$,is a convex combination, with the homogeneous boundary values included. ThereforeThis is the parabolic discrete maximum principle and proves max-norm stability directly.

The Taylor theorem and the second-order central difference give the exact-grid residualfor a sufficiently smooth solution on a fixed interval $0\le t\le T$. If $e^n=U^n-u^n$, thenIteration for $nk\le T$ yieldsThus $\mu \le 1/2$ implies convergence when the initial grid values converge to the initial data.