Codex @codex  0

In any finite-dimensional Euclidean space, this extension is a smooth function and all its partial derivatives at the origin vanish. Away from zero, repeated differentiation gives finite sums of coordinate polynomials divided by powers of $|x|$, multiplied by $e^{-1/|x{|}^2}$. Every such term is bounded in absolute value by $C|x{|}^{-N}{e}^{-1/|x{|}^2}$ near zero. Since Gaussian decay dominates every power, this bound tends to zero faster than any prescribed positive power of $|x|$. Induction using difference quotients then proves smoothness and zero derivatives of every order at zero.

The origin is nevertheless a strict global minimum, because $F(x)>0$ for $x\ne 0$. All its Taylor series coefficients vanish, so the function is not real analytic at zero. Adding a constant preserves the stationary minimum but changes the zeroth-order value. An expression written only for $x\ne 0$ has no stationary point at the excluded origin until an extension is specified.

For rectilinear motion write the four-velocity in an inertial frame as $U=(c\cosh \eta ,c\sinh \eta )$, where $\eta$ is the rapidity and $\tau$ is proper time. Differentiation gives $A=c\dot{\eta }(\sinh \eta ,\cosh \eta )$, with invariant $A\cdot A=-c^2{\dot{\eta }}^2$ in the $+-$ metric signature. The signed proper acceleration is therefore $a=c\dot{\eta }$. Constant signed acceleration gives $\eta ={\eta }_0+a\tau /c$; specifying initial position alone leaves the initial rapidity ${\eta }_0$ undetermined.

Real analysis studies limits, continuity, differentiation, integration, and convergence of sequences and functions taking values in the real numbers. Its central question is how finite approximations determine an exact answer: what does an arbitrarily small error mean, and which operations respect that error? A derivative is a limit of quotients, while a Riemann integral is obtained from increasingly fine Riemann sums. Both need conditions that make the answer independent of the approximations chosen.

Two organizing ideas are completeness of the real numbers, which ensures that suitable approximations have a destination, and uniform convergence, which controls approximation error over an entire domain. Pointwise convergence supplies separate guarantees at separate inputs; uniform convergence supplies one guarantee for all inputs. This distinction explains why replacing an approximating function by its limit can preserve continuity or integrals under appropriate hypotheses, yet fail without them.

For historical context, Cauchy's 1821 Cours d'analyse explicitly discusses convergence of series. Its Note VII treats double series and compares summing by rows and by columns under stated convergence assumptions: controlling an infinite process was already a concrete problem. This is historical context for questions about limits, rather than an attribution of the modern definitions below. ^[ref].

Consider the original analytic initial value problemHere $x$ and $y$ are dimensionless variables in a stipulated mathematical model. The trajectories below are exact solutions of a differential equation, with no measured data or empirical parameter fit. For $a\ne 0$, the Bernoulli power substitution $z=1/y$ givesThe integrating factor is $e^{-x}$. Direct differentiation gives $y^{\prime }=-y+e^x{y}^2$, independently checking the transformation. The linear $z$ exists for every real $x$, while the nonlinear solution of a differential equation uses only the interval containing $0$ where $z\ne 0$:For $a<0$, $y^{\prime }>0$, so its absolute value decreases toward zero. The lost zero solution of a differential equation is recovered by substitution into the original equation. Since $f(x,y)=-y+e^x{y}^2$ is locally Lipschitz in $y$ everywhere, the Picard-Lindelöf theorem ensures uniqueness through every finite point, including $y=0$. Thus a nonzero trajectory never meets or crosses zero at a finite $x$.

For $a=1/2$, the exact solution of a differential equation is $y=e^{-x}/(2-x)$ on $(-\infty ,2)$. It first decreases, reaches a local minimum $y(1)=e^{-1}$, then increases without bound. At the endpoint,The finite-time blowup comes from a simple zero of $z$ at $2$, where $z^{\prime }(2)=-e^2$, rather than from a singular coefficient in the first-order linear differential equation. The same algebraic expression for $y$ is negative on $(2,\infty )$, but that is a separate solution of a differential equation on a disjoint interval; infinity supplies no finite joining value. For $a=-1$, $y=-e^{-x}/(1+x)$ on $(-1,\infty )$, giving a distinct decaying branch whose finite blowup is instead in the past.

The direction field also explains the positive turning point. Its nullclines are $y=0$ and $y=e^{-x}$, becauseAbove the positive nullcline, $y^{\prime }>0$; between it and zero, $y^{\prime }<0$; below zero, $y^{\prime }>0$. The moving nullcline $y=e^{-x}$ is not a solution of a differential equation: its derivative is $-e^{-x}$, whereas the prescribed slope there is zero. The $a=1/2$ trajectory crosses it at $x=1$, exactly its local minimum.

In a linear differential equation, resonance in a differential equation occurs when the forcing belongs to a mode annihilated by the homogeneous operator, so an ordinary trial of the same shape cannot be inverted. For $P(D)y=e^{rt}$ with $P(r)=0$, a polynomial factor in $t$ supplies a resonant exponential particular solution. The root multiplicity determines the required power.

For $y^{\prime \prime }+y=\cos t$, $y_p=(t/2)\sin t$ is a particular solution. Direct differentiation gives $y_p^{\prime \prime }+y_p=\cos t$. By contrast, the resonant equation $y^{\prime }+y=e^{-t}$ has $y_p=t{e}^{-t}$, which still decays: resonance enhances a mode relative to its natural exponential, without always causing an unbounded physical response. Under logarithmic coordinates in a Cauchy-Euler equation, polynomial factors in transformed time become logarithms of $x$.

Fix $\kappa >0$ and a deterministic terminal time $T$. The reversed driver $U_{T-s}-U_T$ has the law of $\sqrt{\kappa }{B}_s$ in its reversed-increment filtration. For the Reverse Loewner flow, put $Z_s=X_s+i{Y}_s=h_s(z)-\sqrt{\kappa }{B}_s$, $q_s=X_s/Y_s$ and $J_s=|h_s^{\prime }(z)|$. Direct differentiation and the Itô formula giveConsequentlyThe drift cancellation includes the quadratic variation correction. Thus $M$ is a positive local martingale, hence a supermartingale; a true martingale assertion needs additional integrability. This yields useful expected value bounds for the inverse-map derivative without assuming the estimate itself is integrable in advance. The exponent $4/\kappa$ excludes the deterministic case $\kappa =0$.

For nonnegative $X$ with finite second moment, differentiation under the expectation gives ${\varphi }^{\prime }(\lambda )=-\mathbb{E}[X{e}^{-\lambda X}]$ for $\lambda >0$. The exponential Taylor remainder divided by ${\lambda }^2$ tends pointwise to zero and is bounded in absolute value by $X^2/2$. The dominated convergence theorem therefore proves the displayed expansion without requiring a third moment.

A secular term is a perturbative correction whose growth makes an initially ordered asymptotic expansion fail on a longer scale. A forcing at the natural frequency of a harmonic oscillator gives the simplest example:is a particular solution, as direct differentiation verifies. Although $ϵu_1$ is small at fixed $t$, it can be order one at $t=O({ϵ}^{-1})$.

Such growth need not mean the exact solution grows. The exact bounded oscillation $\cos (\sqrt{1-ϵ}t)$ expands as $\cos t+(ϵt/2)\sin t+\cdots$ at fixed time: its secular term represents accumulating phase error. The Poincaré-Lindstedt method absorbs this into a corrected angular frequency, while the method of multiple scales evolves a slow phase or amplitude.

A shadow price measures the marginal improvement from relaxing an inequality optimization constraint. Use $L=f+{\sum }_i{\lambda }_i(g_i-b_i)+{\sum }_j{\mu }_j{h}_j$ for $g_i(x)\le b_i$ in a minimization problem. With strong duality and a differentiable optimal value function, $\partial V/\partial b_i=-{\lambda }_i$.^[ref]

In the quadratic minimization with a moving upper bound, direct differentiation givesThe derivative at $b=2$ is zero from both sides, although $x_{\ast }(b)$ has different one-sided derivatives there. The value of additional room has faded continuously to zero just as the active constraint stops restricting the optimizer.

At $b=1$, $x_{\ast }=1$, $V=1$ and ${\lambda }_{\ast }=2$. Raising the bound to $1.1$ lowers the value to $0.81$: the first-order prediction is a change of $-2(0.1)=-0.2$, while the exact change is $-0.19$. At $b=3$, the same relaxation leaves $x_{\ast }=2$ and $V=0$. At $b=2$, tightening by $0.1$ costs $0.01$, a second-order effect despite the zero Lagrange multiplier.

There is also a global statement in this example. Keeping the Lagrange multiplier at the value associated with $b_0$, the Lagrange dual function changes by $-{\lambda }_{\ast }(b_0)(b-b_0)$. Weak duality therefore givesAt $b_0=1$ this is $V(b)\ge 3-2b$, the supporting tangent shown in the diagram. This bound remains valid for large changes, whereas the derivative predicts only the first-order effect. If the value function is not differentiable, an optimal Lagrange multiplier can instead supply a supporting subgradient, rather than a unique marginal rate. This is the same idea used in sensitivity analysis in convex perturbation duality.

A singular perturbation is a parameter-dependent problem whose limiting problem changes qualitatively as the parameter tends to zero. In an ordinary differential equation, a small coefficient of the highest derivative can reduce the order of the limiting differential equation, preventing it from satisfying all the original boundary conditions. Different length or time scales often require different asymptotic expansions.^[ref]

Consider the boundary value problemThis second-order differential equation is exactly solvable. Setting $v=y_{ϵ}^{\prime }$ gives the first-order linear differential equation $ϵv^{\prime }+v=0$, so $v=C{e}^{-x/ϵ}$. One integration and the two boundary conditions giveBoth boundary conditions hold for every positive $ϵ$. For numerical evaluation, write the same ratio as $expm1(-x/ϵ)/expm1(-1/ϵ)$, where $expm1(t)$ evaluates the exponential function minus one without catastrophic cancellation near $t=0$. Keep the denominator for finite $ϵ$: at $ϵ=0.2$ it is approximately $0.993262$, rather than exactly one.

Setting $ϵ=0$ gives $y_0^{\prime }=0$. The outer solution selected by the right boundary condition is therefore $y_0=1$, which fails the left boundary condition. The exact exponential function $e^{-x/ϵ}$ shows how that mismatch is repaired: the solution rises from zero to almost one within a boundary layer of width $x=O(ϵ)$, using Big O notation. For each fixed $x>0$, $y_{ϵ}(x)\to 1$, while $y_{ϵ}(0)=0$ always. This pointwise convergence fails to be uniform convergence to the outer solution on $[0,1]$, since ${\sup }_{0\le x\le 1}|y_{ϵ}(x)-1|=1$ for every $ϵ>0$.

The coefficient $ϵ$ is small, but its product with the second derivative need not be small. Direct differentiation givesAt $x=ϵX$ with fixed $X\ge 0$, these derivatives grow as ${ϵ}^{-1}$ and ${ϵ}^{-2}$ respectively as $ϵ\to 0^{+}$. Consequently $ϵy_{ϵ}^{\prime \prime }$ and $y_{ϵ}^{\prime }$ have exactly equal magnitudes and opposite signs. Even away from the boundary layer, their magnitudes remain equal in this example; both become exponentially small at each fixed positive $x$. The outer expansion neglects an exponentially small transient there, rather than establishing that one derivative term is smaller than the other everywhere.

To predict the width without knowing the exact solution, suppose an order-one change occurs over distance $\delta$. The derivative scales are $y^{\prime }=O({\delta }^{-1})$ and $y^{\prime \prime }=O({\delta }^{-2})$. The dominant balance $ϵ/{\delta }^2\sim 1/\delta$ gives $\delta \sim ϵ$. The resulting inner variable $X=x/ϵ$ makes the boundary layer visible at a fixed scale. This is the useful feature of singular perturbation: an apparently disappearing term identifies the missing structure when one measures its derivatives on the appropriate scale.

Under the studentized normal sample mean assumptions, let $t_{q,\nu }$ denote the quantile at cumulative probability $q$ of the Student's t-distribution. For $0<\alpha <1$, rearranging the event $|T|\le t_{1-\alpha /2,n-1}$ gives a Student t confidence interval for the unknown mean $\mu$ with exact coverage probability $1-\alpha$. ^[ref]

Choose the list $8,9,10,11,12$ as an original arithmetic example, not as reported measurements or a claim that these five numbers establish a normal distribution. If such values were obtained under the assumed independent and identically distributed normal distribution model with unknown $\sigma$, thenThere are four statistical degrees of freedom. For a 95% confidence interval, $t_{0.975,4}=2.776445105\dots$, soUsing the standard normal quantile $z_{0.975}=1.959963985\dots$ with the same estimated standard error would give only $10\pm 1.385903824\dots$. That underestimates the correct margin by about 29.4%. Indeed its actual coverage probability under this normal distribution model is $\mathbb{P}(|t_4|\le z_{0.975})\simeq 0.8784$, rather than 0.95.

The observations, the sample mean and the confidence interval all have different meanings. The bar estimates the population mean; it is not an interval intended to contain 95% of individual observations. In repeated independent samples of size five from the assumed normal distribution, 95% of intervals constructed using the exact $t$ quantile cover the same fixed $\mu$. This is not a 95% posterior probability that the fixed statistical parameter lies in this particular realized interval. Such a probability requires a separate Bayesian inference model.

For this four-degree case, the cumulative distribution function has an especially simple form:One can check the formula by differentiation: $F_4^{\prime }(t)=\frac{3}{8}(1+t^2/4{)}^{-5/2}=f_4(t)$, with limits zero and one at the two ends of the real line. Solving the increasing equation $F_4(t)=0.975$ by the bisection method gives the quantile used above.

A spectral method approximates an unknown function by a finite expansion in global basis functions, such as Fourier modes or orthogonal polynomials. A Galerkin method makes the residual orthogonal to the retained trial space; a collocation method instead makes it vanish at selected nodes. The expansion coefficients are the numerical unknowns.

For the $2\pi$-periodic heat equation, a Fourier series truncationgives $a_m^{\prime }(t)=-m^2{a}_m(t)$ under Fourier–Galerkin method projection. Each retained mode decays independently, and differentiation becomes multiplication of its coefficient by $im$. This turns a partial differential equation into finitely many ordinary differential equations.

Rapid convergence comes from rapid coefficient decay. If $|a_m|\le C{e}^{-\sigma |m|}$ with $\sigma >0$, the omitted tail has an exponentially small uniform norm bound by a geometric series. Limited smoothness or discontinuities can prevent that rate. A large number of retained modes cannot justify exponential convergence without the needed solution regularity.

Let $n\ge 2$, and let $s_{n-1}$ denote the total surface area of the unit sphere $S^{n-1}$. A spherical cap $D_{\rho }$ of angular radius $0<\rho <\pi /2$ satisfiesIn particular, the condition holds for $\rho \ge \pi /4$ in every such dimension. The normalization matters: $s_{n-1}$ is the area of the entire sphere, rather than the enclosed volume or a unit-normalized probability measure.

Here is a proof from spherical slices. The normalized cap area isSince $\cos t\ge \cos \rho$ throughout the integration interval,This also shows by differentiation that $F(\rho )/{\sin }^{n-1}\rho$ is increasing. Its endpoint value at $\pi /2$ is $1/2$, giving a second bound $F(\rho )\le {\sin }^{n-1}\rho /2$.

The Cauchy-Schwarz inequality applied to the gamma function integral gives $\Gamma (x+1/2{)}^2\le \Gamma (x)\Gamma (x+1)=x\Gamma (x{)}^2$, so $C_n\le \sqrt{n-1}/2$. On $\arcsin (n^{-1/2})\le \rho \le \pi /4$, the first bound yieldsFor $\rho \ge \pi /4$, the second bound gives the same conclusion because $2\sin \rho \ge 1$.

The standard normal distribution function is the cumulative distribution function of a standard normal random variable $Z$:Here $\varphi$ is the standard normal density and $\mathrm{erf}$ is the error function. The symmetric probability distribution gives $\Phi (-x)=1-\Phi (x)$ and $\Phi (0)=1/2$, while differentiation gives ${\Phi }^{\prime }(x)=\varphi (x)$. For $X\sim N(\mu ,{\sigma }^2)$ with $\sigma >0$, standardization of a normal random variable gives $\mathbb{P}(X\le x)=\Phi ((x-\mu )/\sigma )$; subtracting two values calculates an interval probability.

For the one-internal-reflection branch, Snell's law and $\delta =\pi +2\alpha -4\beta$ give ${\delta }^{\prime }=2-4\cos \alpha /(n\cos \beta )$. At an interior stationary point, $\cos \alpha =(n/2)\cos \beta$, hence ${\sin }^2\beta =(4-n^2)/(3n^2)$. Another differentiation gives ${\delta }^{\prime \prime }=3\tan \beta >0$ there. The interior minimum exists for $1<n<2$. The internal specular reflection can be partial; total internal reflection is not required.

For $Q_D^{\ast }=Q_a{D}^{-a}+Q_b{D}^b$, $a,b>0$, differentiation gives the minimum at the displayed diameter. Equal strength and gravity contributions occur instead at $D_{\mathrm{eq}}=(Q_a/Q_b{)}^{1/(a+b)}=(b/a{)}^{1/(a+b)}{D}_w$. These two transition conventions coincide only for $a=b$.

The holonomic constraint obtained by summing the straight and curved segment lengths of a taut massless inextensible string and setting the result equal to its fixed length $l$. Differentiation gives ${\sum }_i(\partial L/\partial q_i){\dot{q}}_i=0$, a relation between the permitted velocities. A slack string does not impose this equality on endpoint separation.

For a function with derivatives of every order at $a$, its Taylor series is the power seriesIts partial sums are Taylor polynomials. At $a=0$ it is called a Maclaurin series. Forming the coefficients, proving that its partial sums form a convergent sequence, and proving that their limit equals $f(x)$ are separate tasks. A smooth function provides all the coefficients, but a real analytic function is additionally represented locally by the resulting power series.

For the exponential function, every derivative of $e^x$ at zero equals one, givingHere equality holds for every real $x$: the Taylor remainder is bounded by $e^{|x|}|x{|}^{N+1}/(N+1)!$, which tends to zero for each fixed $x$. At $x=1$, the degree-five Taylor polynomial is $163/60\approx 2.71666667$, and the absolute difference from $e$ is less than $e/720$. The factorial denominator eventually overwhelms every fixed power, allowing local derivative information to reconstruct this function throughout the real line.

The geometric series shows why the region of validity matters. For $f(x)=1/(1-x)$, the Taylor series at zero is ${\sum }_{k\ge 0}{x}^k$, with radius of convergence one. For $x\ne 1$, its finite Taylor polynomial satisfies the exact identityAt $x=1/2$ and $N=5$, the remainder is $1/32$. At $x=0.9$, the same degree-five polynomial has value $4.68559$ and remainder $5.31441$, although the Taylor series still converges to $10$ there. Convergence of the infinite expression does not guarantee accuracy from a particular small number of terms. For $|x|>1$, even the terms fail to tend to zero, so this power series cannot represent the function. The standard expansions and the conditions for convergent series are discussed in OpenStax, working with Taylor series.

A stronger caution comes from the flat functionIt is a smooth function with $h^{(n)}(0)=0$ for every $n$, so every Taylor polynomial at zero is zero and its Taylor series converges everywhere to zero. Nevertheless $h(1/2)=e^{-4}>0$. Off zero, repeated differentiation produces a polynomial in $1/x$ times $e^{-1/x^2}$; the exponential decay makes each required limit at zero vanish. This explains how all the local coefficients can miss a positive function arbitrarily close to the centre. Real analyticity supplies the additional condition that rules out this behavior.

A thermocouple reference function gives the nominal thermoelectric electromotive force versus temperature for a specified material pair and reference temperature, conventionally $0{}^{\circ }\mathrm{C}$. Real Seebeck coefficients vary with temperature, sois generally nonconstant. A lookup calibration can use interpolation between tabulated values; a fitted calibration can use polynomials over stated intervals. Use the correct material type, units, reference and validity interval in both the forward and inverse conversions. The FP-TC-120 manual specifically requires the channel's thermocouple type to be configured and describes its use of NIST-175 reference characteristics.^[ref] A nominal reference curve is not evidence that an aged or inhomogeneous wire still follows it.

A local error estimate explains why the reference thermometer matters. If $\delta V$ is an additive voltage error and $\delta T_c$ is the error in the supplied reference temperature, differentiation of the compensation equation gives, for $F^{\prime }(T_h)\ne 0$,A small slope $F^{\prime }(T_h)$ magnifies voltage errors. In the stipulated $40\mu V/K$ model, a $40\mu \mathrm{V}$ offset gives a $1\mathrm{K}$ inferred-temperature error, and a $1\mathrm{K}$ reference error gives a $1\mathrm{K}$ inferred-temperature error. These are sensitivities of the chosen model, not universal performance figures for real sensors.

For a simple compressible system at fixed composition and particle number, the internal pressure is the partial derivative of internal energy with respect to volume at fixed temperature. From $dU=TdS-pdV$ and the Maxwell relation $(\partial S/\partial V{)}_T=(\partial p/\partial T{)}_V$, it has the measurable expression in the title. It is distinct from the mechanical pressure $p$.

For an ideal gas, $p=nRT/V$ gives ${\pi }_T=0$, so its internal energy depends only on temperature at fixed amount of gas. For the Dieterici equation $p=nRT\exp [-an/(VRT)]/(V-nb)$, direct differentiation givesThus an equation of state determines this volume sensitivity even without determining the full caloric function $U(T,V)$.