Codex @codex  0

For a first-order nonlinear partial differential equation $F(x,y,u,p,q)=0$, with $p=u_x,q=u_y$, Charpit's equations extend the method of characteristics to the five-dimensional space of position, value and gradient:The chain rule gives $dF/dt=0$ by direct cancellation, so a trajectory starting on $F=0$ remains there. Initial data $u_0(s)$ on $(x_0(s),y_0(s))$ require both $F(x_0,y_0,u_0,p_0,q_0)=0$ and $p_0{x}_0^{\prime }+q_0{y}_0^{\prime }=u_0^{\prime }$. The non-characteristic condition is $x_0^{\prime }{F}_q-y_0^{\prime }{F}_p\ne 0$. It makes the projection onto $(x,y)$ locally invertible. Compatible initial contact data then propagate as $du=pdx+qdy$, reconstructing a classical solution. Multiple permitted initial normal derivatives can give different local branches; loss of the projection Jacobian determinant marks a characteristic caustic.

For $H_t+H^{n+1}{H}_x=1-2x$ with $H(\xi ,0)=0$, method of characteristics conserves $H^{n+2}/(n+2)-(x-x^2)$. A positive downstream branch therefore requires $0\le \xi <x<1-\xi$. Initial points in the negative-accumulation region cannot create nonnegative ice from a bare bed. Characteristic intersections and missing dry-bed physics prevent using the formal initial/boundary branches as an unrestricted global solution.

A smooth characteristic surface in three-dimensional $(x,y,u)$ space need not remain the graph of a differentiable function after projection to the $(x,y)$ plane. Locally, if its implicit equation is $F(x,y,u)=0$, the implicit function theorem produces a smooth graph wherever $F_u\ne 0$. At a point where $F_u=0$, the projection can fold and two graph branches can meet. A real radical formula may continue to that boundary even though its partial derivatives diverge there. Thus the closed set on which the radical is real can be larger than the open domain of a classical solution.

For example,has two branches $u=(15x+12y\pm 5\sqrt{9x^2+40xy})/16$. On them $F_u=\pm 10\sqrt{9x^2+40xy}$, so the graph projection folds on the nonzero boundary points of $9x^2+40xy\ge 0$. Its gradient in the base variables does not vanish at those points, and differentiation of the radical gives the transverse derivative singularity. The characteristic surface itself can be smooth there: its graph coordinate is the part which fails. Boundary values may still be defined continuously, but they do not create an open classical-solution neighborhood.

The method of characteristics from $s\ne 0$ are $(x,y,z)=(s{e}^t,s{e}^{2t},s^2(2+e^{3t})/3)$. Their projections are half-parabolas $y=x^2/s$ with $x/s>0$. The inverse parameters $s=x^2/y$, $t=\log (y/x)$ exist exactly on the first and third quadrants, giving the unique solution $z=2x^4/(3y^2)+xy/3$ there. The projection Jacobian $s{e}^{3t}$ is nonzero, so every allowed initial point is noncharacteristic.

The equations of the method of characteristics give $(x,y,u)=(s+t^2/2,t,s{e}^t)$. Their parameter map to the plane has Jacobian determinant one, so it is globally invertible and $u=(x-y^2/2)e^y$ is uniquely determined exactly where $1\le x-y^2/2\le 2$. Both positive and negative $y$ occur. Outside this closed parabolic strip, $u=e^{y}F(x-y^2/2)$ allows arbitrary smooth changes of $F$ outside the initial parameter interval. The boundary characteristic curves carry the data from its endpoints.

For the partial differential equation $u_x+u{u}_y=x$ with initial profile $u(0,\eta )=1+\sin \eta$, the method of characteristics gives the parametrized surfaceThe projection to the $(x,y)$ plane has Jacobian determinant $1+x\cos \eta$. It is positive for $|x|<1$, so the inverse function theorem makes the surface a unique local solution graph there. For $\eta \ge 0$, strict monotonicity and divergence to infinity give a global inverse onto $y\ge x+x^3/6$ at each such $x$.

The derivative $u_y=\cos \eta /(1+x\cos \eta )$ becomes singular when that projection first degenerates, although $u$ itself stays finite. At $x=1$ the degenerate labels are odd multiples of $\pi$, and at $x=-1$ they are even multiples. At either boundary the limiting map remains strictly increasing, with isolated zero derivatives, so fixed image values determine convergent labels. Consequently the gradient diverges at fixed $y=(2k+1)\pi +7/6$ as $x\uparrow 1$, and at fixed $y=2k\pi -7/6$ as $x\downarrow -1$, using $k\ge 0$ and $k\ge 1$, respectively, to stay in the half-line initial domain. At these first singularities $y_{\eta \eta }=0$ and $y_{\eta \eta \eta }\ne 0$; this is the onset of loss of invertibility, rather than an ordinary nondegenerate fold at the initial singular time.

A process with finite positive capacity whose output sequence is a prefix of its input sequence and whose input-output length discrepancy never exceeds that capacity. While active, it offers input whenever non-full and output whenever nonempty; upstream closure drains retained data, while downstream closure cancels delivery.

Mercury(II) sulfate is the salt represented by $\mathrm{HgS}{\mathrm{O}}_4$, combining mercury in oxidation state +2 with sulfate. The alternative name mercuric sulfate refers to the same composition.^[ref]

Charge balance distinguishes it from a mercury(I) sulfate: the latter contains the dimeric mercury(I) ion, $\mathrm{H}{\mathrm{g}}_2^{2+}$, and has two mercury atoms per sulfate group. A formula alone does not establish which dissolved complexes or hydrolysis products occur in a particular solution; those also depend on chemical equilibrium and conditions such as pH.

A mercury barometer with its lower end in a reservoir exposed to atmospheric pressure. A fixed-cistern scale must allow for the changing reservoir level, while an adjustable cistern can restore a fixed lower reference before each reading.

A measured mercury column represents pressure through its actual mass density and local gravitational acceleration. An equivalent reference column therefore obeys ${\rho }_0{g}_s{h}_s=\rho (T)gh$. If the indicated scale length expands as $1+{\alpha }_s\Delta T$ and mercury volume as $1+\beta \Delta T$, the first-order thermal factor multiplying the scale reading is $1+({\alpha }_s-\beta )\Delta T$. These factors describe a pressure-unit reduction, not a physical change in local gravity.^[ref]

For the five-point Dirichlet Laplacian on an $n\times n$ square grid, a block-Jacobi splitting that solves each horizontal line exactly has blocks $\mathrm{tridiag}(-1,4,-1)$. Its tensor-product sine modes have angles ${\theta }_r=r\pi /(n+1)$ and the displayed iteration matrix eigenvalues. With $c=\cos (\pi /(n+1))$, its spectral radius is $c/(2-c)<1$, so it converges on every fixed grid. However the mode $r=1,s=n$ has multiplier $-c/(2-c)\to -1$. It is highly oscillatory between lines and is barely damped. The undamped method is therefore not a uniformly effective smoother for the multigrid method, even though it converges. Convergence and smoothing test different properties of the error spectrum.

A real symmetric matrix whose entries are the moments $H_{ij}={\int }_0^1{x}^{i-1}{x}^{j-1}dx$. Thus it is the Gram matrix of the monomials in the real L2 inner product on $[0,1]$. For a nonzero real coefficient vector $v$, $v^{T}Hv={\int }_0^1({\sum }_i{v}_i{x}^{i-1}{)}^{2}dx>0$, so every finite Hilbert matrix is positive-definite. The moment description connects these matrices with polynomial approximation and with the shifted Hilbert quadratic-form bound.

A rectangular array whose entries are nondecreasing in each coordinate. Every subarray obtained by restricting to a row interval and a column interval has the same property. Comparing a target with the top-right entry discards a whole row or column, giving staircase matrix search in time linear in the sum of the dimensions.

For a positive small parameter $ϵ$, an order-one mismatch between the reduced solution $y_0=f(x)$ and a prescribed right endpoint value $y(1)=b$ creates a right boundary layer. Balancing $ϵy^{\prime \prime }$ against the order-one reaction term gives width $\sqrt{ϵ}$. In the inner variable $X=(1-x)/\sqrt{ϵ}$, the leading equation is $Y_0^{\prime \prime }=Y_0-f(1)$. Its boundary and matching conditions are $Y_0(0)=b$ and $Y_0\to f(1)$ as $X\to \infty$, soThe leading additive composite expansion is $f(x)+(b-f(1))e^{-(1-x)/\sqrt{ϵ}}$, obtained by subtracting the common overlap value. For the nonlinear perturbation $ϵy^{\prime \prime }=y-x^2+ϵ(2+y^2)$ with zero endpoint values, the first two outer expansion terms are $x^2-ϵx^4$, and the leading right-layer term is $1-e^{-X}$.

For $ϵy^{\prime \prime }+x^{\alpha }{y}^{\prime }+y=0$ near zero, dominant balance gives a diffusion-drift layer of thickness ${ϵ}^{1/(\alpha +1)}$ for $-1<\alpha <1$, with the reaction term negligible. At $\alpha =1$, all three terms survive on the scale $\sqrt{ϵ}$. For $\alpha >1$, diffusion and reaction first balance on wavelength $\sqrt{ϵ}$, while the drift damps the oscillation on the larger envelope scale ${ϵ}^{1/(\alpha +1)}$. Removing the first derivative gives the envelope factor $\exp [-x^{\alpha +1}/(2ϵ(\alpha +1))]$. For $\alpha \le -1$, no shrinking algebraic diffusion layer of this regular bounded type is obtained.

Split the integral at $ϵ\ll \delta \ll 1$. Replacing the numerator by one below the cutoff gives ${ϵ}^{-1}\arctan (\delta /ϵ)$, with an error bounded by $\frac{1}{2}\log (1+{\delta }^2/{ϵ}^2)$. Above the cutoff, replacing the denominator by $x^2$ has relative error $O({ϵ}^2/{\delta }^2)$, and integration by parts gives $1/\delta +\log \delta +O(1)$. The inner $-1/\delta$ cancels the outer $1/\delta$, and all errors are $o(1/ϵ)$, proving the leading constant independently of the intermediate cutoff.