For positive endpoints $a,b$, halfway along a logarithmic scale meansHere $\beta >1$ is the base. Equivalently, using the natural logarithm, $2\ln (m/x_0)=\ln (a/x_0)+\ln (b/x_0)$. The product of numbers identity for logarithms gives $m^2=ab$, so the midpoint is the geometric mean $\sqrt{ab}$. The two endpoints lie at equal multiplicative distances because $m/a=b/m$.

For $a=10$ and $b=100$, this midpoint is $\sqrt{1000}\approx 31.6228$, giving a factor of $\sqrt{10}$ on each side. The arithmetic mean $55$ is halfway on a linear scale, but its two ratios are $55/10=5.5$ and $100/55\approx 1.81818$. Thus reading halfway between the marks $10$ and $100$ as $55$ on a logarithmic scale gives the wrong value.

For positive quantities, a difference of one order of magnitude means a factor of ten; a difference of $n$ orders of magnitude means a factor of $10^n$. The logarithmic difference ${\log }_{10}(x_2/x_1)$ records this comparison, without subtracting the quantities themselves. For example, $1000$ is three orders of magnitude above one, regardless of a common change of units. Roughly assigning a single value to an order of magnitude also requires a rounding convention; powers of ten are the reference levels.

For $t>0$, the natural logarithm satisfiesTo prove it, put $h(t)=\ln (1+t)-2t/(2+t)$. Then $h(0)=0$ and the derivative is $h^{\prime }(t)=t^2/[(1+t)(2+t{)}^2]>0$. More precisely,This positive integral controls errors in symmetric approximations to the logarithm. After writing $x=1+t$, the bound is $\ln x>2(x-1)/(x+1)$ for $x>1$; the inequality reverses for $0<x<1$ by replacing $x$ with $1/x$.

For an oriented parametric curve $\mathbf{r}(t)$ and scalar field $V$, define the componentwise line integralIts value is a vector and reversing orientation changes its sign. It differs from the scalar arc length integral ${\int }_{C}Vds=\int V(\mathbf{r}(t))|{\mathbf{r}}^{\prime }(t)|dt$, which is independent of orientation. A request to integrate a scalar field along a curve should specify which differential is intended.

For the counterclockwise circle $\mathbf{r}(\varphi )=(R\cos \varphi ,R\sin \varphi ,z_0)$, suppose the restricted field has Fourier series $V(\varphi )=A_0+{\sum }_{n\ge 1}(A_n\cos n\varphi +B_n\sin n\varphi )$. Orthogonality givesThe oriented vector integral selects the first angular harmonics, whereas the scalar arc-length integral selects the mean.

A line integral can be defined as an improper integral of its parametrized integrand when that integrand is undefined or unbounded at isolated parameter values. At an interior singular parameter, split the integral and require both one-sided limits to exist separately; their cancellation in a symmetric Cauchy principal value is a different convention. For a scalar arc length integral, use $f(\mathbf{r}(t))|{\mathbf{r}}^{\prime }(t)|$ in place of the vector-field dot product.

If the restricted integrand has a removable singularity, integrating its continuous extension gives the same value. This does not require the original vector field to extend continuously to a full neighbourhood of the curve. Likewise, a differential identity may hold only on the regular segments: a zero integrand there has improper integral zero, although its printed classical expression may not be defined at the omitted points.

An original toy model separates two ways to lose a central cathode ray signal. Let $x$ be forward depth into a uniform gas. Normalize the entering electron population to one. Write $U(x)$ for the fraction still in a narrow forward channel, and $D(x)$ for the fraction diverted into a wider collection channel. Assume constant removal rate $a$ per unit depth in both channels, and constant transfer rate $b$ from $U$ to $D$, with no return. Removal means that an electron stops or is trapped before reaching the collecting plane; it does not mean that its electric charge disappears. The wide ideal collector accepts all surviving members of both channels. This is a stipulated two-population differential equation model, not a general law of electron transport:The first equation gives exponential decay. Adding the equations gives $W^{\prime }=-aW$ for $W=U+D$, henceThese formulas also give a conservation check: $U+D+(1-W)=1$. The fraction $1-W$ has been removed before the plane, while $D$ remains present but misses the narrow collector. An exponential fit to the narrow signal therefore measures $a+b$, not the removal rate $a$ by itself. This is an identifiability issue: one declining signal cannot distinguish the two unknown rates.

For a worked example, choose $a=0.15\mathrm{c}{\mathrm{m}}^{-1}$ and $b=0.35\mathrm{c}{\mathrm{m}}^{-1}$. These are invented teaching inputs, not measured coefficients for air or Lenard's apparatus. At $x=2.00\mathrm{cm}$,Thus about $63.2\%$ misses the narrow collector, although only about $25.9\%$ has been removed in this model. Treating the whole narrow-signal deficit as stopping would confuse redistribution with removal.

Suppose a calibrated readout is proportional to arriving electron number with the same response for either channel. Choose an empty-path signal $R_0=1020$ arbitrary units, of which $B=20$ is background; the two collection geometries have the same calibrated reference. The calculated readings at $2\mathrm{cm}$ are $R_{\mathrm{n}}=20+1000e^{-1}\simeq 387.88$ and $R_{\mathrm{w}}=20+1000e^{-0.3}\simeq 760.82$. After background subtraction, define $T_j=(R_j-B)/(R_0-B)$ for each geometry. Then the original inverse problem has the solutionA single narrow reading would instead allow every nonnegative pair with $a+b=0.50\mathrm{c}{\mathrm{m}}^{-1}$. Using the raw ratios without subtracting $B$ gives an erroneous narrow rate of about $0.4834\mathrm{c}{\mathrm{m}}^{-1}$. In a real fluorescent screen, the assumed response must be tested: deposited energy and conversion to light need not remain constant as electrons slow. Repeating the comparison at several depths tests the constant slopes, although straight slopes alone do not establish the microscopic mechanism.

Actual electrons change kinetic energy and direction. Electron stopping power $S(K)=-dK/ds>0$ describes mean kinetic energy loss per travelled path length $s$; collisional losses include ionization, while radiative losses include bremsstrahlung. The continuous slowing-down approximation range integrates reciprocal total electron stopping power, neglecting fluctuations. It approximates travelled path length, not forward penetration depth.^[ref]As a separate original limiting example, stipulate a straight path, constant $S=2\mathrm{keV}\mathrm{c}{\mathrm{m}}^{-1}$, and entering kinetic energy $K_0=6\mathrm{keV}$. Then $K(s)=6\mathrm{keV}-Ss$ until stopping, and $R_{\mathrm{CSDA}}=3\mathrm{cm}$. In that deterministic model every entering electron reaches each depth below $3\mathrm{cm}$ and none crosses a plane beyond it. There is no exponential population tail. Real fluctuations and electron scattering make the situation more complicated, but the example already proves that continuous energy loss does not imply constant-rate removal.

The window also changes the entering population: source kinetic energy before the foil need not equal $K_0$ just outside it. For an electron initially at rest, an ideal static accelerating electric potential difference $V>0$ supplies electrical work $eV$, where $e$ is the elementary charge. If losses before and within the foil total $\Delta K\ge 0$, a transmitted electron enters the observation chamber with $K_0=eV-\Delta K\le eV$. A screen's brightness determines neither $V$ nor $\Delta K$. These formulas specify a modern ideal work-energy theorem argument; no numerical historical voltage or foil thickness is inferred from the diagram.

For nonnegative integer $n$, a second solution of the Legendre differential equation is obtained by reduction of order from $P_n$ as $P_n(x)\int dx/[(1-x^2)P_n(x{)}^2]$. Normalize it as $Q_0=\frac{1}{2}\log |(1+x)/(1-x)|$ and $Q_1=x{Q}_0-1$, with the usual Legendre polynomial recurrence relation for higher degrees. Inside $(-1,1)$, $Q_n\sim -\frac{1}{2}\log (1-x)$ near $1$ and $Q_n\sim (-1{)}^n\frac{1}{2}\log (1+x)$ near $-1$. These follow from $P_n(1)=1$, $P_n(-1)=(-1{)}^n$ in the integral. Poles in the integral at interior polynomial zeros cancel in the product, by reduction of order across a zero of the known solution. On $x>1$, the integral from infinity defines the decaying normalization $Q_n(x)\sim 2^n(n!{)}^2{x}^{-n-1}/(2n+1)!$. For $n\ge 1$ all bounded exterior solutions decay this way; for $n=0$ an additional constant solution remains bounded.

Choose a primitive cubic lattice with $a=0.400\mathrm{nm}$, one identical point scatterer per primitive unit cell, incident direction ${\mathbf{s}}_{\mathrm{in}}=(0,0,1)$, and the continuous wavelength band $0.080\le \lambda \le 0.160\mathrm{nm}$. This is a hypothetical model, not the zinc blende specimen of the Laue X-ray diffraction experiment of 1912. Write a candidate reciprocal lattice vector as $\mathbf{G}=(2\pi /a)(h,k,-m)$ with $m>0$, and set $S=h^2+k^2+m^2$. The Laue equations and elastic scattering giveA transmission photographic plate at $z=L>0$ intercepts only $s_{\mathrm{out},z}>0$. Its coordinates relative to the direct beam areFor $(h,k,-m)=(2,1,-1)$ and $L=50\mathrm{mm}$, $S=6$, $\lambda =0.13333\dots \mathrm{nm}$ and ${\mathbf{s}}_{\mathrm{out}}=(2/3,1/3,2/3)$, so $(X,Y)=(50,25)\mathrm{mm}$. Its length is one, and each component difference satisfies the corresponding Laue equations. This verifies a direction and a selected wavelength, without predicting an exposure or measured wave intensity.

The first panel has fourfold symmetry because its chosen Bravais lattice, beam axis and square detector share that symmetry. For the tilted panel, apply a rotation matrix $R$ to each index vector, set $v=R(h,k,l)$, and use $\lambda =-2a{v}_z/|v{|}^2$ and ${\mathbf{s}}_{\mathrm{out}}=(0,0,1)+(\lambda /a)v$. All indices with $|(h,k,l)|\le 2a/{\lambda }_{\min }=10$ suffice: elastic scattering bounds $|\mathbf{G}|\le 2|{\mathbf{k}}_{\mathrm{in}}|$. Thus the plotted spectral and angular selection is complete within this ideal point model and detector window. Multiples of an index vector select inversely scaled wavelengths and can land on the same spot; a white-beam spot need not identify a unique diffraction order.

For diffusion with a constant diffusion coefficient $D>0$ and no sources, the model flux is $J=-D\nabla u$. A continuity equation gives $u_t=-\nabla \cdot J=D\Delta u$. This proportionality is a constitutive assumption about transport; the continuity equation expresses conservation. If $D$ varies with position, it must remain inside the divergence, giving $u_t=\nabla \cdot (D\nabla u)$ instead.^[ref]

An exact constructed heat equation example on the whole plane uses dimensionless coordinates with $D=1$ and the smooth initial Gaussian function $u(x,y,0)=e^{-r^2}$. Set $s=1+4t$. For $t\ge 0$,Direct differentiation verifies both the initial condition andThe center decreases, while the region $r^2>s$ increases at that instant. The sign-change circle expands with radius $\sqrt{s}$. A fixed outer location can initially increase and later decrease as this circle passes it.

At time $t$, the peak is $1/s$, and the radius where the profile is $1/e$ times its peak is $\sqrt{s}$. The width grows as the height falls. Integration in polar coordinates givesThus a decreasing peak is compatible with a conserved spatial integral. Dividing $u$ by $\pi$ makes a probability density function with coordinate variances $s/2=1/2+2t$, connecting the same spreading to Brownian motion with diffusion coefficient one. Equivalently the profile is evolved by the Gaussian heat kernel with time parameter $2t$, since that article uses generator $\frac{1}{2}\Delta$.

The Laplacian detects whether a value lies below or above nearby values, to second order. For a scalar field whose partial derivatives through order four are continuous functions, take the arithmetic mean of four equally spaced neighbors:Applying Taylor's theorem in each coordinate cancels the first- and third-order terms between opposite neighbors. The two pairs leaveThis is an asymptotic expansion, rather than an equality for arbitrary $h$. Rearranging gives the familiar central finite difference approximationThis finite difference replaces partial derivatives by five sampled values. The arithmetic mean lies above the center for a sufficiently small $h$ if the Laplacian is strictly positive there, and below if it is strictly negative. A zero Laplacian cancels this second-order discrepancy without requiring equality of individual values.

At $(0,0)$, the bowl $x^2+y^2$ has center value zero and four neighbor values $h^2$, so the discrepancy is exactly $h^2$. The saddle $x^2-y^2$ has neighbor values $h^2,h^2,-h^2,-h^2$, whose arithmetic mean is exactly zero. For the Gaussian function $e^{-x^2-y^2}$ the discrepancy is $e^{-h^2}-1=-h^2+O(h^4)$, agreeing with $\Delta f(0,0)=-4$.

A harmonic function need not equal this four-neighbor arithmetic mean for a finite $h$. The polynomial $f(x,y)=x^4-6x^2{y}^2+y^4$ has $\Delta f=0$ everywhere, but $A_{h}f(0,0)=h^4$. Higher-order terms remain even when the second-order terms cancel.

At a local minimum, each coordinate second derivative is at least zero, hence the Laplacian is at least zero. At a local maximum the inequality reverses. The sum alone does not characterize an extremum: cancellation at a saddle point of a scalar function can also produce zero. The local-average interpretation explains why the heat equation smooths differences, and helps motivate the heat equation maximum principle.

For a constant orthogonal matrix $R$, the Laplacian is unchanged by the coordinate transformation $\xi =Rx$. If $u$ is twice continuously differentiable, the chain rule givesbecause ${\sum }_i{R}_{ji}{R}_{ki}={\delta }_{jk}$ by orthogonality of the rows. Here ${\delta }_{jk}$ is the Kronecker delta. In two dimensions, the rotation matrix $\left(\begin{array}{cc}a& -b\\ b& a\end{array}\right)$ with $a^2+b^2=1$ makes the mixed partial derivatives cancel explicitly. The same argument works for reflections and translations, so the Laplacian depends on the Euclidean geometry rather than the orientation or origin of the coordinates.

Choose the five numbers $1,2,2,5,10$ as an original illustration. Their sample median is 2 and their sample mean is 4. The sum of absolute losses at the sample median isIf these numbers were modeled as an independent and identically distributed Laplace distribution sample, the joint maximum-likelihood estimates would therefore be $\hat{\mu }=2$ and $\hat{b}=12/5=2.4$. Using the sample mean instead gives $S(4)=14$, so it does not maximize this likelihood function.

Now replace only 10 by 100. The sample median stays at 2, while the sample mean becomes $(1+2+2+5+100)/5=22$. The minimized sum of absolute losses rises to $S(2)=102$, and the fitted scale rises to $\hat{b}=102/5=20.4$. Stability of the fitted location therefore does not mean stability of every fitted parameter.

The left panel subtracts each objective's minimum so the two absolute loss curves can be compared despite the different totals. The right panel tracks the sample mean and sample median as the largest chosen value changes. This is an outlier sensitivity illustration, not actual observations or a check that a Laplace distribution fits any population. In applications, regression residuals and the data-generating assumptions need separate examination; a stable sample median alone does not validate a statistical model.

The Laguerre polynomial normalized by $L_n(0)=1$ has coefficients $L_n(x)={\sum }_{j=0}^n(-1{)}^j(\genfrac{}{}{0ex}{}{n}{j})x^j/j!$. By the Leibniz rule, this is equivalent to Rodrigues' formula, $e^{-x}{L}_n(x)=D^n(e^{-x}{x}^n)/n!$. For any polynomial $q$ of degree below $n$, $n$ applications of integration by parts give ${\int }_0^{\infty }{e}^{-x}{L}_{n}qdx=0$. All endpoint terms vanish: derivatives of $e^{-x}{x}^n$ of order below $n$ vanish at zero and decay at infinity. Since $L_n^{(n)}=(-1{)}^n$, the same calculation givesThus these orthogonal polynomials are already an orthonormal set for the weighted inner product with weight $e^{-x}$.

Set $v=x^n{e}^{-x}$. Direct differentiation gives $x{v}^{\prime }=(n-x)v$. Apply the Leibniz rule to differentiate $n+1$ times:Hence $x{v}^{(n+2)}+(x+1)v^{(n+1)}+(n+1)v^{(n)}=0$. By Rodrigues' formula, $v^{(n)}=n!e^{-x}{L}_n$, so the next two derivatives are $n!e^{-x}(L_n^{\prime }-L_n)$ and $n!e^{-x}(L_n^{\prime \prime }-2L_n^{\prime }+L_n)$. Substitution gives the displayed Laguerre differential operator on polynomials eigenvalue equation. This derivation also covers $n=0$.

A Frobenius solution $y={\sum }_{k\ge 0}{a}_k{x}^{k+\sigma }$ has the repeated indicial root $\sigma =0$. Equating coefficients gives the displayed recurrence relation. The resulting analytic solution is an entire function: its consecutive term ratio has magnitude asymptotic to $|x|/k$. For $\nu =n\ge 0$ an integer, the recurrence terminates and $L_n(x)={\sum }_{k=0}^n(-1{)}^k(\genfrac{}{}{0ex}{}{n}{k})x^k/k!$. The second local solution is independent and has a logarithmic term at zero; repeated indicial roots do not supply two independent analytic series.

Putting $Y=e^{z}V$ in the displayed Kummer differential equation gives $z{V}^{\prime \prime }+(z+b)V^{\prime }=0$. Thus $V^{\prime }=C{z}^{-b}{e}^{-z}$ on a chosen branch, yielding $Y=e^z[A+B\int z^{-b}{e}^{-z}dz]$. For $b=1/2$ and $z=x/2$, one solution is $e^{x/2}$ and a normalized second Frobenius solution is $\sqrt{x}{\sum }_{n\ge 0}{x}^n/(2n+1)!!$. The analytic factors in both solutions have infinite radius of convergence, since their coefficient ratios tend to zero. The second solution still has a square-root branch at zero; an entire series factor does not remove a nonintegral Frobenius exponent.

Choose a positive constant $\Sigma$ and the illustrative interface lawHere $T_h$ and $T_c$ are the hot and cold interface-side temperatures, not distant thermometer readings. This is a low-temperature phonon model under restrictive assumptions, not a universal fit to Kapitza's measured exponent or a material-property table. Its coefficient $\Sigma$ is not the Stefan-Boltzmann constant: the carriers in this model are phonons, rather than photons crossing a vacuum.

The fourth power has a simple counting explanation. In the three-dimensional low-temperature Debye model, a linear phonon dispersion relation gives a density of states proportional to ${\omega }^2$. A phonon has energy $\hslash \omega$ and Bose-Einstein distribution occupation $[\exp (\hslash \omega /k_{B}T)-1{]}^{-1}$. If propagation speeds and the angle-averaged transmission introduce no additional frequency scale over the occupied modes, the outgoing energy flux has the formExtending the upper integration limit assumes occupied modes lie well below the cutoff. Balancing the two outgoing fluxes, with coefficients satisfying zero net heat flux at equal temperatures, gives the chosen difference of fourth powers. Frequency-dependent transmission, different excitations or additional transport channels can change this conclusion.

The finite-drop area-specific thermal resistance follows by factoring a difference of powers:Near thermal equilibrium, write $T_h=T+\delta T$ and $T_c=T$. Expansion gives $q=4\Sigma T^3\delta T+O((\delta T{)}^2)$, so the differential area-specific Kapitza resistance isAt a finite bias, varying only the hot side gives $(\partial q/\partial T_h{)}^{-1}=1/(4\Sigma T_h^3)$; this is generally different from $r_{K,\sec }$. Stating what is held fixed makes a differential resistance reproducible.

For an original normalized example, choose reference scales $T_0>0$ and $A_0>0$, and setFor contact area $A=2A_0$, the total heat transfer rate is $\dot{Q}=2.1472A_0\Sigma T_0^4$. The finite-drop resistance satisfies $\Sigma T_0^3{r}_{K,\sec }=0.2/1.0736=125/671\approx 0.18629$. Using the equilibrium tangent at $T_0$ instead would predict $j=4(0.2)=0.8$, illustrating why a substantial temperature rise requires the nonlinear law. Doubling the common equilibrium temperature reduces $R_K$ by eight in this chosen model.

A series check can be done without guessing any material properties. Stipulate a hot-side slab with dimensionless areal bulk resistance $b_s=\Sigma T_0^3{L}_s/k_s=0.10$, and neglect the cold-side bulk drop. At the same $j=1.0736$, the hot reservoir must have $T_H/T_0=1.2+b_{s}j=1.30736$. The imposed reservoir drop splits into $0.10736$ in the bulk and $0.20$ at the interface. Applying the fourth-power law directly to $T_H$ would assign the bulk drop to the interface and overestimate the heat flux.

Finally, transfer heat at rate $Aq$ between maintained interface-side temperatures $T_h>T_c>0$. The hot side loses entropy at rate $Aq/T_h$, and the cold side gains it at rate $Aq/T_c$. With no interface entropy accumulation, the entropy production per unit area isFor the normalized example it is ${\dot{S}}_{\mathrm{prod}}/(A\Sigma T_0^3)=1.0736(1-1/1.2)=0.178933\dots$. The same energy crosses both sides in a steady state, while total entropy increases, consistently with the Second law of thermodynamics. A finite temperature jump during heat transfer therefore expresses dissipation, not a failure of energy conservation.

Pyotr Leonidovich Kapitza's 1941 study suspended cleaned, annealed copper blocks in liquid helium. Sealed, gas-filled cavities held a heater and a bronze resistance thermometer. Electrical readout tracked the heating-induced temperature rise; the earlier readout method assumed a stable bath without a separate bath thermometer. A nearby unheated block warmed much less, localizing the drop near the heated surface. Gold plating and coatings tested surface effects; tight cavities prevented liquid ingress. These were solid-body readings, not probes resolving the microscopic interfacial layer. ^[ref]

The measurement logic suggests a modern control sequence: calibrate both temperature readouts, establish their zero-load difference, increase heater power, wait for a steady state, and repeat while monitoring the bath. Determine lead losses and the bulk temperature gradients before using $r_B=(T_s-T_l)/q$. These are proposed controls derived from the energy balance, rather than additional reported steps in the 1941 apparatus. A nearby unheated solid also has its own interface: its reading need not equal the liquid's local temperature exactly. The experiment supports a drop concentrated near the surface, while a microscopic thickness and a universal interfacial law require further evidence.

For the Jacobi iteration for a three-variable tridiagonal system,An orthonormal eigenbasis iswith eigenvalues $q=\sqrt{2}/4$, $0$, and $-q$, respectively. These are error modes: their shapes persist under the iteration matrix, while their coefficients change. Direct multiplication verifies them, or expansion of the determinant gives $\det (tI-B)=t(t^2-1/8)$. Thus $\rho (B)=q<1$ and convergence holds for every starting vector.

With $x^{(0)}=0$, the initial forward error has coefficientsFor $k\ge 1$,The zero mode disappears after one sweep; the negative mode alternates in sign. The Euclidean norm drops by exactly $q$ per sweep from $k=1$ onward. The first norm drop is different because $\Vert e^{(0)}{\Vert }_2=\sqrt{14}$ includes the zero mode. A starting error proportional to $u_0$ would disappear entirely in one sweep, faster than the spectral radius predicts for generic errors.

The curves are computed iterates and exact formulas for this linear system. Since $A=4(I-B)$, the same eigenvectors have eigenvalues $4-\sqrt{2}$, $4$, and $4+\sqrt{2}$ for $A$. They prove that $A$ is a positive-definite matrix, and give the condition number used in residual and forward error in a stationary iteration.

Here $V_d$ is the volume of the unit ball in $d$-dimensional Euclidean space. Rotational symmetry makes the integral an isotropic second-rank tensor $C{\delta }_{ij}$. Taking its trace gives $dC={\int }_{|x|\le R}|x{|}^2{d}^{d}x=d{V}_d{\int }_0^R{r}^{d+1}dr$. This proves the formula. For a constant matrix $M$, contraction gives $\int x^{T}Mx{d}^{d}x=V_d{R}^{d+2}\mathrm{tr}M/(d+2)$. In three dimensions the coefficient is $4\pi R^5/15$.