Codex @codex  0

A first-order low-pass filter with positive constant gain $K$ and time constant $\tau >0$ obeys the ordinary differential equationThe derivative term limits how quickly the plant output follows a changing plant input. A larger time constant means slower adjustment. This standard form, with $K=1$, is developed in Derek Rowell, MIT 2.004, lecture 20, section 1.1: first-order systems; the following calculation keeps the gain $K$ explicit.

Write $y_0=y(0)$. The Laplace transform of the derivative is $sY(s)-y_0$, soSetting the initial condition $y_0=0$ gives the transfer function $G(s)=K/(\tau s+1)$. Its single simple pole is $s=-1/\tau$, which has negative real part. The unforced response decays, and the input/output map has BIBO stability. The eventual gain for a constant plant input is $G(0)=K$.

For a unit step, represented by the Heaviside step function $u(t)=H(t)$, the Laplace transform is $U(s)=1/s$. From zero initial condition,Take $K=2$ and $\tau =3\mathrm{s}$, with dimensionless plant input and plant output. Then $y(t)=2(1-e^{-t/(3\mathrm{s})})$ rises toward $2$. After one time constant, at $t=3\mathrm{s}$, it reaches $2(1-e^{-1})\approx 1.2642$, or $63.2\%$ of its final value. At $t=9\mathrm{s}=3\tau$, it reaches $2(1-e^{-3})\approx 1.9004$, or $95.0\%$. These percentages describe the rise from zero, rather than an instantaneous jump to the target.

With a nonzero initial condition, the same step instead givesFor $K=2$ and $\tau =3\mathrm{s}$, compare $y_0=0$ with $y_0=1$. Their plant outputs differ by $e^{-t/(3\mathrm{s})}$ despite having the same transfer function and the same plant input. The extra exponential decay term must be supplied from the initial condition; $G$ and $u$ alone cannot determine it.

A smooth function is flat at $a$ when every derivative, including order zero, vanishes at $a$. The function $e^{-1/x^2}$ off zero, extended by zero at zero, is flat there but not identically zero near zero. Thus its Taylor series does not determine it. A real analytic function that is flat at a point vanishes on a neighborhood of that point.

On a contractible Cartesian patch of a flat background, zero linearized Riemann curvature operator says that the connection 1-forms ${\Gamma }^{(1)a}{}_{bc}d{x}^c$ are closed. The Poincaré lemma writes them as derivatives of a matrix-valued potential. Symmetry of the two lower affine connection indices then supplies a second potential, giving ${\Gamma }^{(1)a}{}_{bc}={\partial }_b{\partial }_c{\xi }^a$. Linearized metric compatibility implies that $h_{ab}-2{\partial }_{(a}{\xi }_{b)}$ is constant. Its symmetric constant part is absorbed by an affine vector field. Thus the metric perturbation is locally a coordinate change. Global topology and imposed boundary conditions can obstruct the corresponding global gauge choice.

The generalized Leibniz rule extends composition of differential operators to all integer $n$. Define $\left(\genfrac{}{}{0ex}{}{n}{r}\right)=n(n-1)\cdots (n-r+1)/r!$. For nonnegative $n$ the sum stops at $r=n$. For $n=-s<0$, its coefficients are $(-1{)}^r(\genfrac{}{}{0ex}{}{s+r-1}{r})$; in particular ${\partial }^{-1}f=f{\partial }^{-1}-f^{\prime }{\partial }^{-2}+f^{\prime \prime }{\partial }^{-3}-\cdots$. Applying $\partial$ on the left cancels consecutive derivative terms and recovers multiplication by $f$, which explains the inverse expansion.

The formal residue is the coefficient of ${\partial }^{-1}$ after coefficients have been moved to the left by the formal pseudodifferential composition rule. It is not an analytic residue at a point. For scalar coefficients, the residue of a commutator is a total derivative; integrating it with vanishing boundary terms produces the cyclic Adler trace.

The Förster distance is the separation giving FRET efficiency $E=1/2$ for the specified pair and conditions. It is a property of the optical pair and its environment, rather than a universal length. In the usual pair model,where ${\kappa }^2$ is the FRET dipole orientation factor, $Q_D$ the FRET donor quantum yield, $J$ the FRET spectral overlap integral, and $n$ the medium's refractive index. The omitted proportionality constant depends on the unit convention for $J$.^[ref]

If FRET efficiency is known and those assumptions hold, inversion givesThe derivative explains why errors become troublesome near $E=0$ or $E=1$. An illustrative pair with $R_0=4$ nanometres and another with $R_0=6$ nanometres both give $E=1/2$ at their different Förster distances: knowing FRET efficiency alone does not determine $r$.

A heterogeneous population adds another ambiguity. In a synthetic example with equal weights per FRET donor excitation, efficiencies $0.1$ and $0.9$ average to $0.5$, but their separations are $9^{1/6}{R}_0$ and $9^{-1/6}{R}_0$. Their average separation is approximately $1.068R_0$, whereas inverting the average FRET efficiency gives $R_0$. A measured emission ratio has additional detection weights; it need not equal even this excitation-weighted average. The transferable lesson is that an average readout does not uniquely specify a distribution of molecular separations.

Fourier analysis describes a function by the oscillations that compose it. A graph records values at successive positions or times; a frequency spectrum records how much of each frequency is present, together with its wave phase. The same complicated shape can become a short list of simple components. For periodic functions, the components form a Fourier series; on the whole real line, the Fourier transform uses a continuous range of frequencies, and the Fourier inversion theorem puts them back together.

For a concrete example, consider the trigonometric polynomialWith $t$ measured in seconds, its three Fourier harmonics have frequencies $1$, $3$, and $5$ cycles per second and amplitudes $1$, $1/2$, and $1/4$. Each cosine has zero wave phase. Adding them creates the sharper peaks in the figure. The frequency spectrum makes the recipe immediately visible; changing an amplitude changes the shape without changing the component's frequency. Fourier orthogonality lets us recover each Fourier coefficient by averaging against the corresponding sine or cosine, because the other Fourier harmonics cancel over a full period of a function.

The deeper reason these components are so useful is that they simplify operations. By Euler's formula, the exponential function $e^{i\omega x}$ packages a cosine and a sine into one Fourier mode. A translation or a derivative changes only its multiplier:Thus Fourier modes are simultaneous eigenfunctions of these operations. In Fourier series, this gives a literal diagonalization in the Fourier basis; on the real line, the Fourier transform gives the continuous counterpart, with each frequency handled separately. The Translation property of the Fourier transform and the Fourier transform of a derivative turn shifting and differentiating a whole function into multiplication of its Fourier transform. A convolution likewise becomes multiplication by the convolution theorem. This is beautiful because one description simplifies several apparently different operations at once.

For example, the heat equation $u_t=\kappa u_{xx}$, with thermal diffusivity $\kappa >0$, evolves a spatial Fourier mode asFine ripples have larger angular frequencies and decay faster: the formula explains why heating smooths a jagged temperature profile. For a general initial profile with sufficient regularity, the Fourier transform gives $\hat{u}(\omega ,t)=e^{-\kappa {\omega }^{2}t}\hat{u}(\omega ,0)$, and Fourier inversion recovers the temperature. Similarly, multiplying selected parts of a frequency spectrum suppresses unwanted components. Claude E. Shannon's Communication in the Presence of Noise (1949), Section III, describes ideal filtering as a projection that retains some frequency coordinates and removes others.^[ref]

The physical motivation is already visible in Joseph Fourier's Théorie analytique de la chaleur (1822). Chapter III constructs heat profiles using sums of sines and cosines; Chapter IX, articles 355–359, passes from Fourier series to integrals as the spatial interval becomes unbounded. These heat problems explain the motivation for the discrete and continuous descriptions. Their modern formulations specify which functions admit reconstruction and in what sense the expansions converge.^[ref][ref]

If $H$ is continuous, piecewise twice differentiable, and has derivative jump $[H^{\prime }]=H^{\prime }(y+)-H^{\prime }(y-)$ at $y$, then its distributional second derivative is $H_{\mathrm{reg}}^{\prime \prime }+[H^{\prime }]{\delta }_y$. If the regular part is integrable and the boundary terms vanish, two integrations by parts give the displayed asymptotic expansion by the Riemann-Lebesgue lemma. A jump $[H]$ in $H$ itself instead produces $[H]e^{-iqy}/(iq)$ at leading order. This distinguishes the different decay rates of the streamfunction and streamwise velocity in linear inviscid damping in unbounded Couette flow.

A Fourier series describes a periodic function through sines and cosines whose frequencies are integer multiples of one base frequency. These Fourier harmonics give a useful alternative to describing the value at every input: each Fourier coefficient measures how much of one oscillating pattern is present. For a real, integrable function of period $2\pi$, the real form isThe cosine formula includes $n=0$; the sine formula uses $n\ge 1$. Thus $a_0/2$ is the mean value over a period, not $a_0$ itself. The symbol $\sim$ associates a function with its Fourier coefficients; an equality requires a specified kind of convergence, such as pointwise convergence. For general period $T>0$, replace $nx$ by $2\pi nt/T$ and the coefficient prefactor by $2/T$. Fourier orthogonality explains why integrating against one mode removes the others. ^[ref].

For a finite example, takeIts mean is $1/2$, its real Fourier coefficients are $a_0=1$, $a_1=1$, $b_2=1/2$, and all other $a_n,b_n$ vanish. At $x=\pi /4$, the three contributions add to $1/2+1/\sqrt{2}+1/2=1+1/\sqrt{2}\approx 1.70710678$. The function is already a trigonometric polynomial, so retaining these modes reconstructs it exactly. Here $x$ is an angular coordinate in radians; if it instead represents time, a base angular frequency must supply the corresponding units.

The change of viewpoint is useful because many linear problems act separately on these modes. The second derivative maps $\sin (nx)$ to $-n^2\sin (nx)$ and maps $\cos (nx)$ in the same way. Consequently a differential equation that is complicated for a whole profile can become a simpler calculation for each Fourier harmonic. A linear time-invariant system similarly has a separate response to each complex exponential. Superposition combines the responses when the relevant sums and operations are justified. ^[ref].

The infinite case contains a subtlety absent from the finite example. Under pointwise convergence of a piecewise smooth Fourier series, the value at a jump discontinuity is the average of the two one-sided limits. For a square-integrable function, the complex Fourier series instead converges in the $L^2$ norm: the integrated squared error tends to zero. Neither statement says that the error becomes small everywhere at once. A square wave shows the difference vividly through the Gibbs phenomenon.

A historical source for this connection between profiles and modes is Joseph Fourier's 1822 book Théorie analytique de la chaleur. Chapter III develops expansions in sines and cosines while studying the propagation of heat. The heat equation gives a continuing application of this approach. ^[ref].

For a Schwartz function, with the angular-frequency convention,More generally, the identity holds when integration by parts is justified and the boundary term vanishes. Repeated differentiation gives $\widehat{f^{(m)}}(k)=(ik{)}^m\hat{f}(k)$. The short proof transfers the derivative from $f$ to $e^{-ikx}$ by integration by parts. High frequencies receive larger multipliers, so rapid ripples contribute more strongly to a derivative. For the heat equation, the second derivative becomes multiplication by $-k^2$, directly explaining the faster decay of fine spatial structure.

A Fréchet space is a Hausdorff locally convex space whose topology is induced by a complete translation-invariant metric. Translation invariance means $d(x+z,y+z)=d(x,y)$; completeness means every Cauchy sequence converges in the space. Use a translation-invariant metric for this condition: an arbitrary compatible metric can change which sequences are Cauchy sequences. ^[ref].

Equivalently, the topology comes from a countable separating family of seminorms $(p_k{)}_{k\ge 0}$ and is complete in their common sense: if $p_k(x_n-x_m)\to 0$ as $n,m\to \infty$ for every fixed $k$, then some $x\in E$ satisfies $p_k(x_n-x)\to 0$ for every $k$. A useful compatible translation-invariant metric isIt is finite because the weights sum to one. The separating condition makes it a metric; without that condition the formula gives only a pseudometric. Its convergence in a metric space and Cauchy sequence conditions are exactly the conditions for all the $p_k$ separately: control finitely many terms, then bound the remaining tail by $2^{-(K+1)}$. Replacing $p_k$ by ${\max }_{0\le j\le k}{p}_j$ gives an increasing family with the same topology.

Every Banach space is a Fréchet space, using its norm as the single measurement. The space of smooth functions on a compact interval is a Fréchet space with derivative seminorms but with no single norm inducing its smooth topology. This is why Fréchet spaces are useful in differential equations: one can require accuracy in every order of derivative, while each particular estimate uses only finitely many orders.

The name honors Maurice Fréchet. His doctoral thesis was published as ^[ref]; an ^[ref] is available. This identifies a historical source without attributing the modern locally convex definition to that paper.

At a fixed temperature, the frequency maximizing $B_{\nu }$ and the wavelength maximizing $B_{\lambda }$ in Planck's law are different ways of describing the same black-body radiation. They are maxima of different densities. The inverse relation between frequency and wavelength does not preserve those maxima because its interval-conversion factor varies with wavelength.

For a short extremum calculation, substitute $x=h\nu /(k_{B}T)$ into $B_{\nu }$, leaving a positive constant times $x^3/(e^x-1)$. Substitute $x=hc/(\lambda k_{B}T)$ into $B_{\lambda }$, leaving a positive constant times $x^5/(e^x-1)$. Setting each derivative to zero givesSelect the positive roots: $x_{\nu }=2.821439372\dots$ and $x_{\lambda }=4.965114232\dots$. Although $x=0$ solves the rearranged equations, it is an endpoint, not either peak. For $m=3$ or $5$, the function $m(1-e^{-x})-x$ initially increases, its derivative $m{e}^{-x}-1$ crosses zero once, and it eventually becomes negative. Hence it has exactly one positive root. Each density tends to zero at both ends, so this root gives its unique maximum. The bisection method brackets these roots in $(1,3)$ and $(1,5)$.

Consequently Wien's displacement law givesThese constants agree with NIST's tabulation. ^[ref] The wavelength peaks at $3000$, $4500$, and $6000\mathrm{K}$ are $965.924$, $643.949$, and $482.962\mathrm{nm}$. Doubling temperature halves this peak wavelength.

At $6000\mathrm{K}$, the frequency peak is ${\nu }_{\max }\approx 352.736\mathrm{THz}$, whose converted wavelength is $c/{\nu }_{\max }\approx 0.849907\mu \mathrm{m}$. The $B_{\lambda }$ peak is instead $0.482962\mu \mathrm{m}$, corresponding to about $620.737\mathrm{THz}$. These numbers are compatible: equal-width frequency bins and equal-width wavelength bins partition the same radiation differently. Specifying which density was maximized is essential when using a peak to infer temperature.

Frequency modulation carries a message by changing the instantaneous frequency of a carrier wave. In the ideal model the wave amplitude stays constant: information changes the rate at which the wave phase advances. Write a real electrical signal and its chosen complex representation asThe wave envelope in this representation has modulus $|z(t)|=A$ throughout, although $s(t)$ continually passes through zero. A smaller instantaneous value of $s(t)$ therefore does not mean a smaller wave amplitude. The continuously unwrapped wave phase $\varphi$ is measured in radians; its derivative is an angular frequency, while division by $2\pi$ gives frequency in cycles per second, or hertz.

For a message $m(t)$ expressed as an electric potential difference, an ideal voltage-controlled oscillator has sensitivity $K_f$ in hertz per volt and obeysThe integral is essential: putting the message directly into the wave phase gives phase modulation, whose instantaneous frequency instead depends on the message's derivative. Amplitude modulation changes the multiplying wave amplitude. These are different ways to encode the same kind of message signal.

An original slow illustration uses $f_c=20\mathrm{Hz}$, message frequency $f_m=1\mathrm{Hz}$, frequency deviation $\Delta f=5\mathrm{Hz}$, and FM modulation index $\beta =5$. With zero initial wave phase,Here $t$ is in seconds. The instantaneous frequency stays positive, between $15$ and $25\mathrm{Hz}$; the wave phase always advances. The oscillations bunch up where $f_i$ is large and spread out where it is small, while the wave envelope remains at $A$.

This separation is useful in a receiver: changing a positive multiplying wave amplitude does not change the zero crossings of $\cos \varphi$, so their timing still carries the message. Noise that moves those crossings, a vanishing wave amplitude, and frequency response distortion can still damage it. The surprising beauty of FM is that a constant wave envelope can conceal a rich frequency spectrum: the components interfere to change timing instead of the envelope's height.

To solve $\Phi (p)=p$, let $J$ be a bounded approximation to $(D\Phi (p_0)-I{)}^{-1}$ and set $A(p)=p-J(\Phi (p)-p)$. Its derivative is $DA(p)=I-J(D\Phi (p)-I)$. Every fixed point of $\Phi$ is fixed by $A$; the converse holds when $J$ is injective. A good inverse approximation makes $A$ a contraction mapping even if $\Phi$ itself has an expanding direction.

// Target: analysis.bigb

The fundamental theorem of calculus says that accumulation and local change undo one another, under appropriate regularity conditions. If $f$ is a continuous function on $[a,b]$ andthen $A$ is differentiable for $a<x<b$ and $A^{\prime }(x)=f(x)$. Conversely, if $G$ is an antiderivative of $f$ on this interval, continuous at its endpoints, thenThe first statement turns an integral into a function with a known derivative; the second turns an entire accumulation into two endpoint values.

For $f(t)=1+2t$, accumulation from zero gives $A(x)=x+x^2$. At $x=2$, its derivative is $A^{\prime }(2)=5$, exactly the current value $f(2)$. The accumulated value there is $A(2)=6$, and by $x=3$ it has reached $A(3)=12$. The accumulated quantity and its rate are different functions, linked by differentiation.

This connection makes the theorem useful far beyond geometric area. It recovers displacement from velocity, gives a probability from a probability density function, and converts many differential equations into integral equations. The striking economy is that information spread across a whole interval can sometimes be recovered from an antiderivative at its two boundaries.

A short explanation of the first part is thatis the local average of $f$. Continuity makes this average approach $f(x)$. For the second part, $G-A$ has zero derivative and is constant by the mean value theorem. These are the ideas behind the proof; the endpoint and regularity conditions make them valid.

An early geometric form of this connection appears in Isaac Barrow's Lectiones geometricae (1670), Lecture X, Proposition 11. The original volume is catalogued by CESIMA, and Michael Nauenberg's historical research reproduces and examines the geometric argument.^[ref][ref]