Codex @codex  0

The bounded classical solution is the heat-kernel convolutionThe two-dimensional heat kernel has integral one. Boundedness of $g$ permits differentiation under the integral for $t>0$, proving $u_t=\Delta u$ and $|u|\le \Vert g{\Vert }_{\infty }$. Its approximate identity property gives $u(x,t)\to g(x)$, locally uniformly since $g$ is continuous. This bounded solution satisfies the required Gaussian-growth bound.

For heat-equation uniqueness under Gaussian growth, let $w$ be the difference of two solutions, so $w(x,0)=0$ and $|w|\le M_0{e}^{a|x{|}^2}$ on a finite time interval. Choose $b>a$ and a slab $0\le t\le \tau <1/(4b)$. The positive comparison solutionsatisfies $({\Phi }_b{)}_t=\Delta {\Phi }_b$. For any $\epsilon >0$, its faster spatial growth makes $|w|\le \epsilon {\Phi }_b$ on a sufficiently large lateral cylinder boundary. At the initial boundary, $w=0$. The heat equation maximum principle applied to $\pm w-\epsilon {\Phi }_b$ proves the same inequality inside. First allow the cylinder radius to increase, then let $\epsilon \downarrow 0$. Repeating these slabs proves uniqueness on every finite interval with the stated uniform growth bound. No spatial integrability of $g$ is required.

At time $T$, this is Gaussian filtering with covariance $2TI$, hence $\sigma =\sqrt{2T}$ per coordinate. With the angular-frequency Fourier transform convention $\hat{f}(\xi )=\int e^{-ix\cdot \xi }f(x)dx$, the Gaussian filtering Fourier multiplier isFor merely bounded $g$, this identity is understood through tempered distributions. The multiplier leaves the zero frequency unchanged and attenuates large frequencies exponentially: it suppresses rapid noise fluctuations but blurs image edges. In cycles-per-length frequency $k$, the multiplier is $e^{-4{\pi }^{2}T|k{|}^2}$.

For the printed Perona-Malik equation, retain the supplied conductance $c(s)=s{e}^{-s^2/(2{\lambda }^2)}$, which includes an extra factor $s$. In one dimension write $p=u_x$ and $F(p)=c(|p|)p=p|p|e^{-p^2/(2{\lambda }^2)}$. Then $u_t=F^{\prime }(u_x)u_{xx}$, withThe forward-backward threshold for gradient-weighted exponential diffusion is thereforeAt $|u_x|=0$ and $\sqrt{2}\lambda$ the coefficient vanishes, giving degenerate diffusion. Increasing $\lambda$ increases the forward-diffusion range. Small nonzero slopes smooth, while large slopes formally sharpen; a negative coefficient produces short-wave growth and ill-posedness, so this is a dynamics explanation, not a general existence theorem for arbitrary data. Replacing the printed conductance by $e^{-s^2/(2{\lambda }^2)}$ would give threshold $\lambda$, but that is a different equation.

The four-neighbour mean expansion follows by Taylor expansion: opposing first and third derivatives cancel, soIf $u(x_0)={\mathrm{mean}}_h(u)(x_0)$ for every sufficiently small $h$, division by $h^2$ and passage to the limit give $\Delta u(x_0)=0$. This is a pointwise conclusion; it does not assert harmonicity in a whole neighbourhood.

For the area median over the disk, the disk-median curvature expansion, at a point where $|\nabla u|\ne 0$, isConsequently the analogous median fixed-point condition yieldsThis is the vanishing of the level-line curvature at $x_0$: to second order the level line is straight there. It need not be a straight segment. For example $u(x,y)=y+x^3$ has zero disk median at the origin for every radius, by odd symmetry, but its zero level line is the cubic $y=-x^3$. Vanishing curvature on an entire connected regular level arc would force that arc to be straight. The factor $1/6$ is for the disk's uniform area measure; a circle-boundary median has a different scale factor, while producing the same zero-curvature equation.

The partial differential equation. Each term in $K_{\alpha }$ is a translated heat kernel and satisfies $({\partial }_t-{\partial }_x^2-\alpha {\partial }_x)K_{\alpha }=0$ for $t>0$. The same holds for $P_{\alpha }$, either by direct differentiation or because it is $-(2{\partial }_x+\alpha )H(x+\alpha t,t)$. For $x>0$, $P_{\alpha }(x,t)$ and its derivatives vanish faster than any power as $t\downarrow 0$. Therefore differentiating the boundary convolution creates no extra upper-endpoint term. Gaussian domination justifies differentiating both integrals on compact subsets of $x,t>0$, proving the advection-diffusion equation.

The initial condition. For fixed $x>0$, the first Gaussian in $K_{\alpha }$ is an approximate identity centered at $y=x+\alpha t\to x$. Its integral tends to $u_0(x)$. The reflected Gaussian is exponentially small as $t\downarrow 0$, because its center lies outside the half-line. The boundary convolution also tends to zero for $x>0$. Thus $u(x,t)\to u_0(x)$.

The Dirichlet boundary condition. The identityshows that $K_{\alpha }(0,y,t)=0$, so the initial-data contribution vanishes at zero. The boundary contribution must be evaluated as a limit, not by substituting $x=0$ inside its singular integral. The positive kernel satisfiesfor every $ϵ>0$. Hence it is a one-sided approximate identity at zero time, and continuity of $g_0$ givesThe mass formula follows from the Gaussian Laplace integral, or from the decaying solution of the corresponding constant-coefficient ordinary differential equation. Compatibility $u_0(0)=g_0(0)$ gives the continuous corner value. These arguments also verify the equivalent contour solution in (iii). In the decaying energy class the solution is unique: the difference $v$ of two solutions has zero data and $\frac{d}{dt}{\int }_0^{\infty }{v}^{2}dx=-2{\int }_0^{\infty }{v}_x^{2}dx$ for real solutions, with the analogous modulus identity for complex solutions.

The unheaded sine transform question. The direct classical Fourier sine transform does not close on $u$. Ifthen integration by parts givesThe drift introduces an unknown cosine transform; the usual scalar sine-transform solution of the heat equation is therefore unavailable directly.

A Dirichlet gauge transform for constant drift does provide a qualified alternative. Set $v=e^{\alpha x/2+{\alpha }^{2}t/4}u$. Then $v_t=v_{xx}$, with $v_0=e^{\alpha x/2}{u}_0$ and $v(0,t)=e^{{\alpha }^{2}t/4}{g}_0(t)$. If these weighted data have the decay needed for an ordinary Fourier sine transform, it solves the transformed problem and produces exactly the kernel above. Mere decay of $u_0$ does not ensure this weighted integrability. Thus not directly by the classical sine transform of $u$; yes after a gauge transformation when the required weighted-transform hypotheses hold, or after a justified cutoff/limit argument.

Let $h=\log \psi$, with $\psi >0$, and set $f=2\alpha h_{\theta }$. Direct differentiation givesThe heat equation ${\psi }_Z=\alpha {\psi }_{\theta \theta }$ implies $h_Z=\alpha (h_{\theta \theta }+h_{\theta }^2)$, so the bracket vanishes. Conversely its being independent of $\theta$ can be absorbed into a $Z$-dependent multiplicative normalization of $\psi$, leaving $f$ unchanged. This proves the Cole-Hopf transformation with the positive sign appropriate to the negative-flux Burgers convention. Although the algebra works for nonzero $\alpha$, the given forward Gaussian function diffusion kernel and a physical vanishing-viscosity limit require $\alpha >0$.

Take $U,L>0$ for the stated Burgers N-wave. Integrating ${\partial }_{\theta }\log \psi (0,\theta )=f_0(\theta )/(2\alpha )$ and normalizing the exterior value to one givesThis function is continuous at $\pm L$; its logarithmic derivative has the specified jumps. Convolution with the heat kernel is positive and solves the heat equation for $Z>0$. Splitting the integral into the exterior baseline and the interior correction yieldsPut $a=1+UZ$. Completing the square givesChanging the interior integration variable to $\eta =a\phi$ changes its limits to $\pm La$ and the Gaussian function width to $Za$. The Jacobian and normalization leave the factor $a^{-1/2}$. Thus the Cole-Hopf solution for a Burgers N-wave isHere $I_{\alpha }$ is the normalized Gaussian function mass of its indicated interval. Its explicit error function representation isFor fixed $w>0$, the Gaussian function concentrates at $\phi =\theta$ as $\alpha \downarrow 0$. Consequentlyaway from the endpoints; at $\theta =\pm L$ the limit is $1/2$. The transition layer has width $O(\sqrt{\alpha w})$. This is an approximate identity argument, not a uniform step approximation across the endpoints.

For fixed $Z>0$ and $|\theta |\ll L$, both interval masses tend to one, and the exponentially large positive $\hat{\psi }$ dominates the exterior correction. Its logarithmic derivative therefore givesFor $|\theta |\gg La$, both interval masses are exponentially small and the weighted interior integral is also negligible, while the exterior contribution tends to one. Therefore $\boxed{f(Z,\theta )\simeq 0}$ in the specified far exterior.

An exponentially weighted Gaussian function tail should not be discarded solely because its unweighted interval mass tends to zero. In fact, comparing the order-one exterior term with $\hat{\psi }$ gives the sharper inviscid Burgers N-wave fronts $|\theta |=L\sqrt{1+UZ}$, not $La$. Away from these fronts, the vanishing-viscosity limit isInside $|\theta |<La$, this follows from the sign of $L^2-{\theta }^2/a$ in the exponential. Outside $La$, the constrained Gaussian function maximum lies at an interval endpoint and has negative exponent. The right shock wave's speed is $UL/(2\sqrt{a})$, equal both to the derivative of $L\sqrt{a}$ and to minus half the sum of its two limiting states; the left shock wave is its reflection. This checks consistency with the Rankine-Hugoniot condition and shows that the requested near-center and far-exterior approximations are compatible with the full entropy limit.

Choose nested interior open sets $U\subset \subset V\subset \subset W\subset \subset \Omega$. For small $\sigma$, part (b) gives a smooth equation on $V$ of the formThe divergence-forcing interior H1 estimate isOne can obtain this estimate directly by testing the smooth equation against ${\eta }^2{u}_{\sigma }$ and applying Young inequality, with a cutoff function equal to one on $U$ and supported in $V$. The approximate identity and the $L^2$ convolution bound give, uniformly in $\sigma$,with analogous bounds for $(cu{)}_{\sigma }$ and $f_{\sigma }$. The coefficients need only be bounded on $W$. Therefore $u_{\sigma }$ is bounded in $H^1(U)$. It converges to $u$ in $L^2(U)$, and weak sequential compactness in a Hilbert space in this Sobolev space gives $u\in H^1(U)$. Since $U$ was arbitrary, $u\in H_{\mathrm{loc}}^1(\Omega )$.

Now expand the distributional derivative:Its right side belongs to $L_{\mathrm{loc}}^2$, so the interior $H^2$ elliptic regularity estimate gives $u\in H_{\mathrm{loc}}^2$. More generally, if $u\in H_{\mathrm{loc}}^m$ for an integer $m\ge 1$, multiplication by the smooth coefficients puts the right side in $H_{\mathrm{loc}}^{m-1}$. Interior elliptic regularity then gives $u\in H_{\mathrm{loc}}^{m+1}$. This elliptic regularity bootstrap proves $u\in H_{\mathrm{loc}}^k$ for every integer $k$.

For any nonnegative integer $\ell$, choose $k>\ell +n/2$. The Sobolev embedding theorem on compact interior subsets gives $u\in C_{\mathrm{loc}}^{\ell }$. Taking all $\ell$ proves $u\in C^{\infty }(\Omega )$, for its smooth representative. The first gain from $L^2$ to $H^1$ is the step supplied by the mollified equation; assuming $H^1$ in the initial definition would miss that step.

Write $({\varphi }_{\epsilon }\ast f)(x)=\int \varphi (y)f(x-\epsilon y)dy$. Since $\int \varphi =1$,A Schwartz function is bounded and uniformly continuous. Given a tolerance, first choose $R$ so the tail integral of $|\varphi |$ is small. On $|y|\le R$, uniform continuity makes the difference uniformly small for every $x$ when $\epsilon R$ is small. On the tail it is bounded by $2\Vert f{\Vert }_{\infty }$. This proves uniform convergence, without assuming $\varphi \ge 0$, and is the approximate identity argument.

For the radial regularization of the three-dimensional fundamental solution of the Laplace equation, direct differentiation givesThe radial integral is $\int \varphi =3{\int }_0^{\infty }{r}^2(1+r^2{)}^{-5/2}dr=1$. This even approximate identity therefore gives $\int (-\Delta N_{\epsilon })f\to f(0)$ for $f\in \mathcal{S}({\mathbb{R}}^3)$.

The function $N=1/(4\pi |x|)$ is locally integrable and determines a tempered distribution. We have $0<N_{\epsilon }\le N$, and $N|\Delta f|$ is integrable near zero and infinity for a Schwartz test function. Integration by parts and the dominated convergence theorem giveThus $-\Delta N={\delta }_0$, proving that $N$ is a fundamental solution of a linear differential operator. In this last part the test-function space is $\mathcal{S}({\mathbb{R}}^3)$; the printed ${\mathbb{R}}^n$ in the three-dimensional integral is understood with $n=3$.

For the stated mollifier, $({\partial }_t+{\partial }_s){\chi }_{\epsilon }=({\partial }_x+{\partial }_y){\chi }_{\epsilon }=0$. Thus the interior doubled integral isOn a bounded state range, $H$ is Lipschitz in each argument, and so is $Q$: away from the diagonal each partial derivative has modulus at most $\sup |f^{\prime }|$, and the function is continuous across the diagonal. Local $L^1$ translation continuity therefore makes this converge to the corresponding single-time, single-space integral. The compact time support is chosen below $T$, with a margin for $s\le t+2\epsilon$.

The initial terms need a separate argument. Since the time kernel is supported at $t-s\in [-2\epsilon ,-\epsilon ]$, $B_v$ is identically zero for $t\ge 0$, whereas $B_u$ samples $v(s,\cdot )$ at $\epsilon \le s\le 2\epsilon$.

The averaged initial trace of an entropy solution follows directly from its inequalities. Test the constant-level entropy for $v$ with $(1-s/\delta {)}_{+}\rho (y)$, using smooth approximations of that temporal cutoff. Boundedness of $v$ and its entropy flux givesApproximate $v_0$ locally in $L^1$ by finitely many constants with a smooth nonnegative partition of unity. Applying this estimate on each partition member and using the triangle inequality makes the limsup of ${\delta }^{-1}{\int }_0^{\delta }\Vert v(s)-v_0{\Vert }_{L^1(K)}ds$ arbitrarily small. This proves its convergence to zero for each compact $K$.

The time kernel has size $O({\epsilon }^{-1})$ and the spatial kernel is a normalized approximate identity. The averaged trace, followed by spatial $L^1$ translation continuity of $v_0$, therefore gives $B_u\to \int |u_0-v_0|\phi (0)$. Its time mass is one, not one-half, because its whole support lies on the positive $s$ side. ConsequentlyThis is the Kato inequality for scalar conservation laws with its initial contribution. Assuming an initial trace solely from interior translation continuity would leave a gap; the entropy argument above supplies it.

Put $h=|u-v|$, $q=Q(u,v)$ and $h_0=|u_0-v_0|$. The mean value theorem gives $|q|\le Mh$ on the bounded state range. Fix a final time $t$ and let $L(\tau )=a-M(t-\tau )$ and $R(\tau )=b+M(t-\tau )$. Choose a smooth increasing function $\rho$ equal to zero on $(-\infty ,0]$ and one on $[1,\infty )$. For small $\delta >0$, defineIts two nonnegative derivative contributions imply $(W_{\delta }{)}_{\tau }+M|(W_{\delta }{)}_x|\le 0$. Therefore $h(W_{\delta }{)}_{\tau }+q(W_{\delta }{)}_x\le 0$.

Use the nonnegative test $\phi (\tau ,x)={\eta }_{ϵ}(\tau )W_{\delta }(\tau ,x)$, where ${\eta }_{ϵ}(0)=1$, it is one up to just before $s$, decreases smoothly to zero near $s$, and remains zero afterward. The Kato inequality for scalar conservation laws givesbecause the other interior contribution is nonpositive. At Lebesgue times $s$, the temporal approximate identity tends to $\int h(s)W_{\delta }(s)$; take a countable sequence $\delta \downarrow 0$ and use dominated convergence on the bounded cone. This yieldsfor almost every $s\in [0,t]$. This is local L1 contraction for scalar conservation laws. It uses only local integrability, so the bounded data need not have finite global $L^1$ norm. If $M=0$, the flux difference is zero and the same proof uses a stationary interval.

The multiplication of a distribution by a smooth function is $⟨\psi T,\phi ⟩=⟨T,\psi \phi ⟩$. If $\psi$ is a Schwartz function, the Leibniz rule shows that $\phi \mapsto \psi \phi$ is a continuous map $\mathcal{S}\to \mathcal{S}$, so $\psi T$ is a tempered distribution. When both $\psi$ and $T$ are radial, $\psi \circ R^t=\psi$, andThus Multiplication by a radial Schwartz function preserves radial tempered distributions.

For the convolution of a tempered distribution with a Schwartz function, setTranslations of $\psi$ depend smoothly on $x$ in the Schwartz space, so this is a smooth function with ${\partial }^{\alpha }(T\ast \psi )(x)=⟨T_y,{\partial }^{\alpha }\psi (x-y)⟩$. The bound for $T$ by finitely many seminorms, together with $1+|y|\le (1+|x|)(1+|x-y|)$, provesThus the function also defines a tempered distribution. For a radial function $\psi$, put $g_x(y)=\psi (x-y)$. Then $\psi (Rx-y)=\psi (x-R^{t}y)=g_x(R^{t}y)$, soHence Convolution with a radial Schwartz function preserves radial tempered distributions.

Smoothing alone need not give a Schwartz function: for the constant tempered distribution $T=1$ and a Schwartz function with integral one, $T\ast \psi =1$. The radial Schwartz approximation of tempered distributions therefore combines smoothing with a large-radius cutoff. Choose a nonnegative radial mollifier $\rho \in C_c^{\infty }$, supported in the unit ball with integral one, and a radial cutoff function $\chi \in C_c^{\infty }$ equal to one on the unit ball. PutEach $f_j$ is a smooth function of compact support, hence a Schwartz function, and the two invariance calculations above make it radial.

It remains to prove convergence, including the simultaneous changes of both scales. With $\check{\rho }(x)=\rho (-x)$, the distributional convolution pairing isFor every fixed $m$, the Leibniz rule, rapid decay outside the radius-$j$ ball, and the chain rule for $\chi (x/j)$ giveConvolution by ${\check{\rho }}_{\epsilon }$ is uniformly bounded in $q_m$ for $0<\epsilon \le 1$, because its shifts have size at most one. The mean value theorem applied to $\phi (x-y)-\phi (x)$ similarly givesSplitting the error into the convolved cutoff error and the approximate identity error provesIf $|⟨T,\theta ⟩|\le C{q}_m(\theta )$, this yieldsConsequentlyweakly, and even in the strong dual topology, since $q_{m+1}$ is uniformly bounded on every subset of the Schwartz space that is a bounded set in a topological vector space. No assertion that $T\ast {\rho }_{1/j}$ itself is rapidly decreasing is needed.

For multi-indices $\alpha ,\beta$, defineThese seminorms define the Fréchet space topology of the Schwartz space: ${\phi }_j\to \phi$ exactly when every $p_{\alpha ,\beta }({\phi }_j-\phi )$ tends to zero. An equivalent increasing family isA tempered distribution is a continuous linear functional on this space. Equivalently, $|⟨T,\phi ⟩|\le C{q}_m(\phi )$ for some $C,m$. The usual weak convergence of tempered distributions means $⟨T_j,\phi ⟩\to ⟨T,\phi ⟩$ for every fixed Schwartz function. The strong dual topology instead requires uniform convergence on every subset of the Schwartz space that is a bounded set in a topological vector space; the Fourier maps below are continuous in both topologies.

Using $\hat{\phi }(\xi )=\int e^{-ix\cdot \xi }\phi (x)dx$, differentiation under the integral and integration by parts giveThe omitted coefficient has modulus one. By the Leibniz rule, every integrand is a finite sum of a polynomial times a derivative of $\phi$. Inserting the integrable weight $(1+|x|{)}^{-n-1}$ bounds its L1 norm by finitely many Schwartz space seminorms. In particular $q_m(\hat{\phi })\le C_m{q}_{m+n+1}(\phi )$. Thus the Fourier transform maps $\mathcal{S}$ continuously into itself.

For completeness, the Fourier inversion theorem follows here by Gaussian regularization. The inverse transform of $e^{-t|\xi {|}^2}$ is $g_t(x)=(4\pi t{)}^{-n/2}{e}^{-|x{|}^2/(4t)}$. The Fourier transform of a Gaussian and Fubini's theorem showAs $t\downarrow 0$, the right side tends to $\phi (x)$ by the approximate identity property, while the left side tends to the undamped inverse integral by dominated convergence theorem, since $\hat{\phi }\in L^1$. HenceThe inverse is $(2\pi {)}^{-n}$ times reflection composed with the continuous Fourier transform, and is therefore continuous on the Schwartz space. This proves a continuous linear isomorphism with continuous inverse.

Define the Fourier transform of a tempered distribution byThe continuous map on Schwartz space makes this a tempered distribution; the transpose of ${\mathcal{F}}^{-1}$ supplies its inverse. Pointwise convergence of pairings proves weak continuity. For the strong dual topology, the transform of a bounded set of Schwartz functions is bounded, so uniform convergence of pairings on bounded sets proves continuity of both Fourier maps there too.

Now let $R\in SO(n)$ be a rotation matrix in the special orthogonal group and write ${\rho }_{R}T=T\circ R$. A change of variables with unit Jacobian givesThe rotation equivariance of the Fourier transform on tempered distributions follows by duality:Therefore $\widehat{{\rho }_{R}T}={\rho }_R\hat{T}$. Applying this identity and the inverse Fourier transform gives the two directions:This is precisely preservation of radial tempered distributions. For $n>1$, $SO(n)$ acts transitively on spheres, so an invariant smooth function is an ordinary radial function.

Interpret the malformed printed domain inequality as the first quadrant $x>0,y>0$, as indicated by the two axis boundary conditions. There is also a genuine missing hypothesis: decay of the boundary data does not require the harmonic function itself to be bounded. For zero data, both $q=0$ and $q=xy$ satisfy the Laplace equation and both edge values, but the Wirtinger derivative of the latter isThus the data do not determine an arbitrary solution's derivative without a growth restriction. We derive the integral representation for the bounded solution, which is also the standard decaying solution when the continuous edge data decay. Any other admissible unbounded solution can differ by a zero-boundary harmonic function.

The conformal map $w=z^2$ sends the quadrant bijectively to the complex upper half-plane: the positive real axis maps to $s>0$, and the positive imaginary axis maps to $s<0$. Define its real-axis boundary function byThe compatibility at the corner makes $f$ continuous there. The Dirichlet Poisson integral in a quadrant follows from the upper-half-plane Poisson kernel:It is bounded by the supremum of the bounded edge data, and the kernel's approximate identity property supplies both Dirichlet boundary conditions. More explicitly, for $z^2=a+ic$, where $a=x^2-y^2$ and $c=2xy>0$, this integral isIt displays directly which part of the real-line boundary corresponds to each original axis.

To obtain its derivative, first take real data and use the Schwarz integral formula. A holomorphic function whose real part is the Poisson solution in the $w$-plane isThe regularizing subtraction is independent of $w$ and affects only the imaginary normalization; it permits bounded data without an artificial integrability assumption. Then $q=\mathrm{Re}F(z^2)$. By the chain rule and the Wirtinger derivative,The factor is $z$, not $2z$: taking the real part contributes the factor $1/2$ in the complex derivative. Differentiation under the integral sign is justified uniformly on compact interior subsets because the denominator stays away from the real boundary and decays quadratically at infinity. Splitting the real integral at zero and substituting $s=r^2$ or $s=-r^2$ gives the holomorphic derivative of a quadrant Dirichlet solution:This converges for bounded data, and in particular for the stated sufficiently decaying data. No interior pole is encountered since $\mathrm{Im}(z^2)=2xy>0$. For complex-valued data the same formula holds by applying the real-data construction separately to their real and imaginary parts and using linearity.

An equivalent form, often convenient when the data's derivatives decay, comes from integration by parts. Sincethe two endpoint terms are $-g_2(0)/z^2$ and $+g_1(0)/z^2$. They cancel by the specified compatibility, leavingThis is the harmonic quadrant solution from tangential boundary derivatives applied to the actual derivatives $g_1^{\prime },g_2^{\prime }$; the supplied $g_1,g_2$ in this question are boundary values, not derivatives. Omitting those primes would solve a different boundary problem.

For completeness, bounded Dirichlet uniqueness in a quadrant selects the solution just constructed. The difference of two bounded solutions is zero on both axes. Odd reflection across each axis gives a bounded harmonic extension to the punctured plane; the removable singularity for a bounded harmonic function fills the corner. The harmonic Liouville theorem makes the extension constant, and oddness makes it zero. Decaying data give a decaying Poisson extension: split the real-axis integral into a fixed compact interval and a tail with arbitrarily small data; the compact contribution tends to zero as $|z^2|\to \infty$, while the tail is bounded by that small supremum.

As literally printed without boundedness or a growth condition, the most one can assert isThe zero-boundary harmonic growth in a quadrant includes $v=cxy$ and more generally $v=\mathrm{Im}(z^{2m})$. This is an actual nonuniqueness of the printed problem, not a defect in the integral for the normalized bounded solution.

The defining integral is well-defined and bounded for every dimension: $|\mathcal{F}f|\le \Vert f{\Vert }_1$ and the multiplier is supported on a finite-volume ball. However the asserted general-dimensional $L^1$ conclusion is false for the Euclidean radial cutoff printed in the original PDF. The nonintegrability of a radial triangular Fourier cutoff in three dimensions is a source issue, not an omitted conjugate or OCR substitution.

For a direct counterexample in $n=3$, take the Schwartz function $f(x)=e^{-\pi |x{|}^2}$, whose Fourier transform is the same Gaussian. With $r=|x|$ and $a=2\pi r$, radial integration givesThe amplitude $b(\rho )=\rho (1-\rho /R)e^{-\pi {\rho }^2}$ vanishes at both endpoints, and $b^{\prime }(R)=-e^{-\pi R^2}$. Two integrations by parts, with a further one to bound the remainder, giveThe integral ${\int }_1^{\infty }{r}^2|F_R(r)|dr$ diverges logarithmically for every fixed $R>0$. For example, restrict to the periodic subintervals on which $|\sin (2\pi Rr)|\ge 1/2$; the error is smaller than half the leading term at sufficiently large $r$. Thus $F_R\notin L^1({\mathbb{R}}^3)$, even for this smooth integrable $f$.

The familiar intended argument is valid in one dimension. There the inverse Fourier kernel isand $F_R=K_R\ast f$ by Fubini's theorem. Its mass outside any fixed neighbourhood of zero tends to zero by scaling, soby continuity of translations in $L^1$. For arbitrary dimension, replacing the radial cutoff by the product ${\prod }_{j=1}^n(1-|{\xi }_j|/R{)}_{+}$ gives a tensor product of these integrable kernels and the same approximate identity proof. These are explicit valid repairs, not claims that the PDF specified a product cutoff.