Codex @codex  0

Charge conservation says that the total signed electric charge of an isolated collection stays constant. More generally, a collection closed to electric charge exchange has constant total electric charge, even if energy crosses its boundary. This is the charge bookkeeping relevant to a closed system or an isolated system; isolation from energy exchange is not required. Positive and negative electric charge count with opposite signs. Creating or removing an Electron and a Positron together is compatible with charge conservation, because their combined electric charge is zero. ^[ref].

Imagine a fixed box initially containing electric charge $Q(0)=10\mu \mathrm{C}$. Stipulate incoming conventional current $I_{\mathrm{in}}=3\mu \mathrm{A}$ and outgoing conventional current $I_{\mathrm{out}}=1\mu \mathrm{A}$ throughout $0\le t\le 2\mathrm{s}$, with no other boundary transport. Since one ampere is one coulomb per second, the stored electric charge grows at $2\mu C/s$:If instead the outgoing electric current also equals $3\mu \mathrm{A}$, the stored electric charge remains $10\mu \mathrm{C}$ while electric charge passes through. These are stipulated transport examples; constant terminal electric currents alone do not guarantee a time-independent interior charge density.

Locally, charge conservation is the charge continuity equation, an instance of the general continuity equation:Here $\rho (\mathbf{r},t)$ is signed volume charge density and $\mathbf{J}(\mathbf{r},t)$ is current density. A positive divergence of current density means net outward transport from a small neighborhood, so its charge density decreases. For a fixed control volume $V$ with closed surface $\partial V$ and outward unit normal $\mathbf{n}$, the equivalent integrated balance isThe surface integral is signed: incoming positive electric charge contributes negatively to outward electric current, and outgoing positive electric charge contributes positively. Integrate the charge continuity equation over fixed $V$ and apply the divergence theorem to obtain this balance. Zero net boundary transport gives constant $Q_V$; the global and local statements describe the same conservation law at different scales. ^[ref].

In a steady transport state, charge density and current density are time-independent, so the charge continuity equation reduces to $\nabla \cdot \mathbf{J}=0$. At an electrical circuit junction, choose each terminal electric current $I_k$ positive outward and assume the terminals account for all boundary transport. ThenKirchhoff's current law, ${\sum }_k{I}_k=0$, follows when accumulation of electric charge at the junction is negligible. During storage, retain the same balance with its nonzero derivative; charge conservation continues to hold. This makes charge conservation useful for checking electrical circuit models and for tracking the charging of a capacitor. ^[ref].

The beautiful consistency with Maxwell equations is that the electric field and magnetic field already enforce this transport balance. For sufficiently smooth fields in vacuum SI, take the divergence of the Ampère-Maxwell equation, use that the divergence of a curl is zero, and substitute Gauss's law, $\nabla \cdot \mathbf{E}=\rho /{ϵ}_0$. Constant vacuum permittivity ${ϵ}_0$ givesThe displacement current term supplies the part needed when electric charge accumulates. Historically, Maxwell added that term to extend the steady-current equation; see the existing charge conservation from Maxwell equations for this connection. ^[ref].

The first-kind Chebyshev polynomial $T_n$ converts a multiple-angle cosine into a polynomial:The cosine addition formula gives the recurrence relationThis constructs polynomials valid for every real $x$, including outside $[-1,1]$, where the real-angle description is unavailable. For example,For $n\ge 1$, $T_n$ has polynomial degree $n$ and leading coefficient $2^{n-1}$. On $[-1,1]$ its uniform norm is one, even as its oscillations become more numerous.^[ref][ref]

The polynomial roots of $T_n$ are especially simple:They correspond to equally spaced angles, but their locations in $x$ cluster near $-1$ and $1$. For $T_3$, the polynomial roots are $\sqrt{3}/2,0,-\sqrt{3}/2$. For $T_8$, the positive polynomial roots are approximately $0.980785,0.831470,0.555570,0.195090$, with their negatives completing the set. The distance from the right endpoint to the nearest polynomial root is $1-\cos (\pi /(2n))\sim {\pi }^2/(8n^2)$; central neighboring polynomial roots are separated on the larger scale $1/n$. The small change of cosine near angles zero and $\pi$ explains the clustering.

These polynomial roots supply $n$ interpolation nodes for a polynomial interpolation of polynomial degree at most $n-1$. They differ from the $n+1$ Chebyshev-Lobatto nodes $\cos (j\pi /n)$, which include both endpoints. This choice of polynomial roots has nodal polynomialThe monic Chebyshev extremal polynomial minimizes the uniform norm of this product among monic polynomials of polynomial degree $n$ on $[-1,1]$. Consequently, if $f$ has $n$ continuous derivatives, the usual polynomial interpolation error gives the useful boundThis optimizes the nodal factor in that bound, not the actual error for every function. Endpoint clustering helps control the large endpoint oscillations of equally spaced interpolation associated with Runge's phenomenon. The two node conventions and their geometric construction are discussed in Trefethen's authored text.^[ref]

The same cosine substitution for polynomial approximation explains the weighted orthogonality:The endpoint singularities of the weight are integrable. Thus the Chebyshev polynomials are orthogonal polynomials for this weight, making coefficients in a least-squares polynomial in an orthogonal-polynomial basis separate cleanly. Their polynomial roots also give Chebyshev–Gauss quadrature, connecting polynomial interpolation to Gaussian quadrature.^[ref]

Accuracy still depends on the target function. The absolute value function $|x|$ has a corner at zero, so the hypothesis requiring continuous derivatives in the displayed polynomial interpolation error bound fails. Endpoint clustering cannot remove that corner. A jump poses a stronger obstruction: polynomials are continuous functions, so they cannot converge in the uniform norm to a discontinuous function on the whole interval. Chebfun's documented examples show persistent overshoot for such jumps; these are a limitation of global polynomial interpolation, not a universal promise of rapid accuracy from Chebyshev polynomials.^[ref]

At chemical equilibrium, a reversible chemical reaction has equal forward and reverse reaction rates, so its macroscopic composition stays constant. Individual molecules can keep changing identity: the reaction has not stopped. This dynamic balance is the definition described by the ^[ref]. Equal reaction rates also need not mean equal concentrations. A small population converting rapidly can balance a large population converting slowly.

This distinction gives two useful questions. Which composition is favored at the specified temperature and other constraints? Thermodynamics answers through chemical potentials and an equilibrium constant. How quickly does a prepared mixture approach that composition? Chemical kinetics answers through reaction rates. Knowing the destination does not determine the journey's duration. A composition that barely changes during observation may reflect extremely slow chemical kinetics, so constancy alone does not establish chemical equilibrium.

Consider the idealized elementary chemical reactionAssume a well-mixed closed system at fixed temperature and volume, with first-order law of mass action rates in both directions and no competing chemical reactions. Write $a=[\mathrm{A}]$ and $b=[\mathrm{B}]$ for concentrations in $\mathrm{M}=\mathrm{mol}/\mathrm{L}$. The total $a+b=C=1\mathrm{M}$ is conserved because each conversion replaces one molecule by one molecule. This conservation is special to the chosen stoichiometry.

The forward and reverse reaction rates, measured as conversions per unit volume per unit time, are $J_f=k_{f}a$ and $J_r=k_{r}b$. HenceSetting this derivative to zero givesWith $a(0)=1\mathrm{M}$ and $b(0)=0$, solving the ordinary differential equation yieldsAt $t=5\mathrm{s}$, these are approximately $0.3515\mathrm{M}$ and $0.6485\mathrm{M}$. Initially $J_f=0.3M/s$ and $J_r=0$. At chemical equilibrium, both are $0.075M/s$, while the net reaction rate $J_f-J_r$ is zero. The persistent exchange is visible in the lower panel even after the upper panel has nearly flattened.

Why does chemical equilibrium also involve a free energy minimum? For a closed system at fixed temperature and pressure, without imposed non-expansion work, stable equilibrium minimizes Gibbs free energy over compositions allowed by conserved elemental inventories. See W. C. Carter, MIT thermodynamics notes, sections 2-1.14 and 2-1.19. This is constrained Gibbs minimization for chemical equilibrium: freely minimizing over unrelated species amounts would change the inventory. At fixed temperature and volume, the appropriate potential is Helmholtz free energy. In either setting, varying an allowed extent of reaction gives the corresponding free energy slope ${\sum }_i{\nu }_i{\mu }_i$, with molar chemical potentials and signed stoichiometric vector entries. An interior equilibrium has zero slope; stability requires a minimum. For A converting to B, this says ${\mu }_{\mathrm{B}}={\mu }_{\mathrm{A}}$.

The beauty of this connection is that molecular exchange, a macroscopic composition, and a geometric minimum describe the same balance. In the ideal concentration model, let $x=b/C$ and choose the standard molar chemical potential difference consistent with $K=3$. Up to a composition-independent constant, the fixed-volume Helmholtz free energy then obeyswhere $n$ is the conserved total amount of substance and $R$ the molar gas constant. Its derivative is $\ln [x/(1-x)]-\ln 3$, vanishing at $x=3/4$. Its second derivative is $1/x+1/(1-x)>0$ inside $0<x<1$, so this balance is the unique minimum. The mixing entropy terms explain why the favored composition still contains A.

Rigid compartments impose a different constraint: fixed volume rather than fixed pressure. With Helmholtz free energy $F=U-TS$,This follows immediately by differentiating $F=U-TS$ in the chemical potential differential above. It expresses the same chemical potential through the thermodynamic potential appropriate to fixed temperature and volume.

Consider a closed system consisting of two rigid compartments of equal volume $V$, maintained at a common fixed temperature $T$. A valve permits exchange of the same ideal gas; the partition stays in place. Initially $n_A=9\mathrm{mol}$ and $n_B=3\mathrm{mol}$, with conserved total amount of substance $n_{\mathrm{tot}}=12\mathrm{mol}$. These are macroscopic moles, modeled continuously, not twelve individual particles. By the ideal gas law, $p_B/p_A=n_B/n_A=1/3$, so initiallyAt $300\mathrm{K}$ with the stipulated $R$, this is approximately $-2740.311\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}$, favoring transfer from A to B.

At fixed $T$ and equal $V$, integrating the ideal-gas chemical potential gives the total Helmholtz free energy, up to a term independent of the allocation,Here $n_{\mathrm{ref}}$ is a positive reference amount of substance; any omitted common term linear in $n_A+n_B$ is constant along this exchange. Let $x=n_A/n_{\mathrm{tot}}$ be the dimensionless fraction in A. The relative Helmholtz free energy then takes the particularly simple formIts derivative along the allowed transfer and its second derivative areThus the Helmholtz free energy is a strictly convex function of the allocation for positive amounts. Its unique minimum has $x=1/2$, or $n_A=n_B=6\mathrm{mol}$, and equal chemical potentials. Initially its excess above this minimum is $n_{\mathrm{tot}}RT[\frac{3}{4}\ln (3/2)+\frac{1}{4}\ln (1/2)]\simeq 3915.474\mathrm{J}$ at $300\mathrm{K}$. Relaxation decreases the total Helmholtz free energy by this amount in the model; the final pressures are also equal because both compartments have the same volume.

The beautiful connection is that the direction of matter exchange and the slope of a free energy curve are the same calculation. The entropy gained by redistributing the ideal gas makes equal allocation favorable here. As with chemical equilibrium, a free energy minimum identifies the favored state, while transport properties determine how quickly it is reached. The fixed-temperature, fixed-pressure Gibbs free energy criterion is summarized in 16.4 Free Energy; the rigid experiment instead uses Helmholtz free energy.

A chemostat steady state has constant nutrient concentration and biomass concentration even while nutrient and cells flow through the vessel. It is an equilibrium of the chemostat ordinary differential equations, rather than a culture in which all growth has stopped. Assume $D>0$ and the Monod equation. Setting both derivatives to zero gives the always-existing chemostat washout equilibrium $(s_{\ast },x_{\ast })=(s_{\mathrm{in}},0)$. A positive biomass concentration instead requires growth to balance removal, $\mu (s_{\ast })=D$, whenceThe positive chemostat steady state exists exactly whenThe weaker condition $D<{\mu }_{\max }$ only makes the denominator positive; it does not ensure $s_{\ast }<s_{\mathrm{in}}$ or $x_{\ast }>0$. Because the Monod equation is increasing, the feed must permit a specific growth rate greater than $D$. At a positive chemostat steady state, increasing feed nutrient concentration at fixed $D$ increases the biomass concentration without changing $s_{\ast }$, within the one-nutrient assumptions.

For a concrete chemostat, let ${\mu }_{\max }=1{\mathrm{h}}^{-1}$, $K_s=1\mathrm{g}{\mathrm{L}}^{-1}$, $s_{\mathrm{in}}=9\mathrm{g}{\mathrm{L}}^{-1}$ and $Y=0.5\mathrm{g}\text{ biomass}/(\mathrm{g}\text{ nutrient})$. At $D=0.5{\mathrm{h}}^{-1}$,The nutrient supply minus outflow is $D(s_{\mathrm{in}}-s_{\ast })=4\mathrm{g}{\mathrm{L}}^{-1}{\mathrm{h}}^{-1}$, matching nutrient consumption $\mu (s_{\ast })x_{\ast }/Y=4\mathrm{g}{\mathrm{L}}^{-1}{\mathrm{h}}^{-1}$. Biomass production and removal both equal $2\mathrm{g}{\mathrm{L}}^{-1}{\mathrm{h}}^{-1}$. The chemostat washout threshold is $D_c=\mu (9)=0.9{\mathrm{h}}^{-1}$. For example, $D=0.95{\mathrm{h}}^{-1}<{\mu }_{\max }$ formally gives $s_{\ast }=19\mathrm{g}{\mathrm{L}}^{-1}$ and $x_{\ast }=-5\mathrm{g}{\mathrm{L}}^{-1}$, which cannot describe a physical culture.

There is a useful exact simplification: the nutrient-equivalent concentration $z=s+x/Y$ satisfies $\dot{z}=D(s_{\mathrm{in}}-z)$, so $z(t)=s_{\mathrm{in}}+[z(0)-s_{\mathrm{in}}]e^{-Dt}$. This combines dissolved nutrient and the nutrient equivalent of biomass; it does not assert that these two variables account for every chemical constituent. In coordinates $(z,x)$, the Jacobian matrix at the positive chemostat steady state has eigenvalues $-D$ and $-{\mu }^{\prime }(s_{\ast })x_{\ast }/Y$, both negative. Hence the positive chemostat steady state is locally asymptotically stable. In the worked example the eigenvalues are $-0.5{\mathrm{h}}^{-1}$ and $-2{\mathrm{h}}^{-1}$.

Use circulation to summarize how much a velocity field runs along an oriented closed curve $C$, such as a circle around a fluid vortex. Its definition is the line integralwhere $d\mathbf{l}$ is tangent to $C$ in the chosen direction. Reversing that direction reverses the sign of circulation. Its units are ${\mathrm{m}}^2{\mathrm{s}}^{-1}$, because velocity is multiplied by a length. Counterclockwise circulation in the horizontal plane is positive when the chosen unit normal vector is $+{\mathbf{e}}_z$, according to the right-hand rule.

For solid-body rotation $\mathbf{u}=(-\Omega y,\Omega x,0)$, a counterclockwise circle of radius $R$ has constant tangential velocity $\Omega R$. Multiplying by its circumference givesFor the potential line vortex $\mathbf{u}=K{\mathbf{e}}_{\theta }/r$, the same circle instead gives ${\Gamma }_C=(K/R)(2\pi R)=2\pi K$, independent of $R$. Stipulate $\Omega =1{\mathrm{s}}^{-1}$, $R=1\mathrm{m}$ and $K=1{\mathrm{m}}^2{\mathrm{s}}^{-1}$: both velocity fields have ${\Gamma }_C=2\pi {\mathrm{m}}^2{\mathrm{s}}^{-1}\approx 6.283185{\mathrm{m}}^2{\mathrm{s}}^{-1}$, although their vorticities in the fluid are $2{\mathbf{e}}_z{\mathrm{s}}^{-1}$ and zero respectively. On a circle of radius $2\mathrm{m}$, their circulations are $8\pi {\mathrm{m}}^2{\mathrm{s}}^{-1}$ and $2\pi {\mathrm{m}}^2{\mathrm{s}}^{-1}$. The comparison plot in vorticity shows this difference in radius dependence.

The bridge between circulation and vorticity is Stokes theorem:Here the velocity field must be continuously differentiable on a neighborhood of the entire spanning surface $S$ inside the fluid, and its boundary orientation must match the unit normal vector $\mathbf{n}$. For solid-body rotation, the disk is admissible and its surface integral is $2\Omega \pi R^2$. For the potential line vortex, the origin is singular and excluded from the fluid: a disk crossing it is not admissible. On an annulus, Stokes theorem instead gives ${\Gamma }_{R_{\mathrm{out}}}-{\Gamma }_{R_{\mathrm{in}}}=0$, with both circulations here measured counterclockwise; the induced inner boundary direction is clockwise. Nonzero circulation therefore does not imply nonzero local curl in the available fluid. ^[ref]

This is useful because Stokes theorem lets a local vorticity measurement determine a loop's circulation when a suitable surface exists, while Kelvin's circulation theorem tracks circulation along a moving material curve. The beauty is the precise link between a local derivative and a global line integral, with the domain deciding which conclusions are valid.

A classical solution has enough ordinary derivatives for the differential equation and its data to hold pointwise. For the heat equation, $u_t$ and $u_{xx}$ should be continuous in the interior, with the specified initial and boundary traces. A weak solution instead imposes an integrated identity, allowing less regular data; uniqueness must always be stated for the solution class used.

The Clapeyron equation connects the slope of a phase coexistence curve to the latent heat and the change in volume during a first-order phase transition. A steep boundary can indicate a large latent heat, a small change in volume, or both. This lets a measured boundary constrain caloric properties and predicts how changing pressure shifts a melting temperature.

For a pure substance in bulk phase coexistence, label the initial and final thermodynamic phases $\alpha$ and $\beta$. Both have the same absolute temperature $T$, pressure $p$, and chemical potential. With $s$, $h$, and $v$ denoting entropy, enthalpy, and volume on one common basis, define $\Delta s=s_{\beta }-s_{\alpha }$, $\Delta h=h_{\beta }-h_{\alpha }$, and $\Delta v=v_{\beta }-v_{\alpha }$. ThenThe derivative follows the phase coexistence curve; it is not a partial derivative within one thermodynamic phase. Use molar volume $v_m$ with molar latent heat $L_m$ in $\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}$, or specific volume $v$ with latent heat $L$ in $\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$. Both give $dp/dT$ in $\mathrm{Pa}{\mathrm{K}}^{-1}$. For molar mass $M$, $L_m=ML$ and $\Delta v_m=M\Delta v$, so the predicted slope is independent of the chosen basis. This formulation assumes a simple compressible system, fixed composition, and negligible interface contributions; a curved interface can instead require the Young–Laplace equation.

For positive latent heat in the stated conversion direction, the sign of the slope is the sign of $\Delta v$. Converting liquid to gas normally increases volume and gives a positive slope. The negative melting-point slope of water illustrates the opposite sign. At a thermodynamic critical point, the coexisting gas and liquid merge: both $L$ and $\Delta v$ vanish, so their ratio needs a limiting analysis rather than substitution of zero into the Clapeyron equation.^[ref]

A short derivation uses the molar Gibbs free energy $g_m$, equal to the chemical potential for a pure substance. Equality persists along the boundary, so $d{g}_{m,\alpha }=d{g}_{m,\beta }$ with $d{g}_m=-s_{m}dT+v_{m}dp$. Therefore $-\Delta s_{m}dT+\Delta v_{m}dp=0$. Equality of $g_m=h_m-T{s}_m$ also gives $\Delta h_m=T\Delta s_m$. Eliminating $\Delta s_m$ yields the Clapeyron equation.

For a scalar function in the BV space, derivative variation is the integral of the relative perimeter measures of its superlevel sets. The formula holds on every Borel subset $A$ of the domain, not just for total mass. It converts the continuous ROF denoising problem into linked minimization problems for finite-perimeter sets.

If a velocity operator satisfies $\mathrm{Re}⟨Lh,h⟩\ge \delta \Vert h{\Vert }_2^2$, then sufficiently regular solutions of ${\partial }_{t}f+v\cdot {\nabla }_{x}f+Lf=0$ obey the displayed estimate when spatial boundary fluxes vanish. The free-transport part contributes zero to the energy derivative. Integrating the resulting differential inequality avoids dividing by a possibly zero norm.

For $\dot{x}=\lambda /(1+y^m)-x$, $\dot{y}=\lambda /(1+x^n)-y$, the equilibrium slope product is $K=mn{x}^n{y}^m/[(1+x^n)(1+y^m)]$. A value $K>1$ gives a saddle separating attracting alternatives. Both exponents greater than one allow this at sufficiently large common $\lambda$. If both are at most one it is impossible. If $0<m\le 1<n$, define $Y(x)$ as the unique positive solution of $Y^{-1}+Y^{m-1}=x^{-1}+x^{n-1}$ and set $K_{m,n}(x)=mn{x}^{n}Y(x{)}^m/[(1+x^n)(1+Y(x{)}^m)]$. All equilibria are parameterized by $\lambda =x(1+Y^m)$, whose derivative is $(1+Y^m)(1-K)/(1-m{Y}^m/(1+Y^m))$. Thus multistability for some common strength occurs exactly when ${\max }_{x>0}{K}_{m,n}(x)>1$. For $m=1$ the maximum is $(n-1{)}^2[(n+1)/(n-1){]}^{(n+1)/n}/(4n)$. Consequently $mn>1$ alone is not sufficient when the two synthesis strengths are required to coincide.

Let the real spline $s={\sum }_{i=p}^q{c}_i{N}_i$ have distinct knots, order $k\ge 2$, nonzero endpoint coefficients, and no identically zero knot cell in its support. It and its first $k-2$ derivatives vanish at both support endpoints. If $Z$ is its number of interior distinct zeros, there are initially $Z+2$ zero components. Each differentiation through order $k-2$ increases this count by at least one: a nonzero gap between zero components has an interior extremum, and the two endpoint zeros remain separate. The final derivative is continuous piecewise linear on $q-p+k$ cells, and each cell can meet at most one zero component unless it is identically zero, in which case all its zeros belong to the same component. It therefore has at most $q-p+k$ zero components. Thus $Z+k\le q-p+k$, proving the bound.

For an order-$m$ operator with smooth coefficients on a bounded piecewise smooth domain, integration by parts has no boundary remainder when at each boundary point either the entire boundary jet of $u$ through order $m-1$ vanishes or the entire such jet of $v$ vanishes. Every boundary summand pairs derivatives of the two functions of order at most $m-1$. It is insufficient to require only that, for each same multi-index $\beta$, one of ${\partial }^{\beta }u$ and ${\partial }^{\beta }v$ vanish: for $P={\partial }_x^2$, $u=x(a-x)\chi (y)$ and $v=1$ on a rectangle, all those same-index products vanish on the boundary but $\int (Pu)v=-2a\int \chi \ne 0$.

Complex analysis studies functions of a complex number, especially holomorphic functions. Its distinctive feature is that a single complex derivative controls variation in every direction of the complex plane. The map $f(z)=z^2$ gives a first picture: straight coordinate lines become curved, yet small perpendicular arrows remain perpendicular wherever $f^{\prime }(z)\ne 0$.

The bottom panels show displacements rather than absolute positions. Multiplication by $f^{\prime }(1+i)=2+2i$ predicts the image of each small arrow; the exact image differs by $h^2$. The origin is exceptional because $f^{\prime }(0)=0$. This distinction between local conformal maps and a critical point is central: complex differentiability imposes strong geometric structure, without making every holomorphic function globally one-to-one.

The remarkable next step is that holomorphic functions are determined inside a circle by their values on its boundary. Precisely, let $R>0$, let $D=\{z:|z-a|<R\}$, and suppose $f$ is a holomorphic function on an open set containing the closed disc $\overline{D}$. Orient its boundary counterclockwise. The Cauchy integral formula states that, for every $w\in D$,The contour integral combines boundary values with a kernel depending on the requested interior point. It is a reconstruction formula, not merely a trick for evaluating an integral. The hypothesis about the closed disc ensures that all boundary values and the contour integral are well defined. ^[ref]

At the centre, substituting $z=a+R{e}^{i\theta }$ cancels the denominator against $dz=iR{e}^{i\theta }d\theta$:Thus the centre value is the angular average of the boundary values. For the polynomial $f(z)=z^2$, $a=1+i$ and $R=1/2$, expand the boundary data asBoth oscillating terms have zero average, leaving $f(a)=2i$. Away from the centre the weights change: on the unit circle,With $f(z)=z^2$ and $w=1/2$, expanding the denominator as a geometric series leaves only its second-power term after integration, giving $f(1/2)=1/4$. This also explains how the Cauchy integral formula recovers an entire Taylor series from boundary data.

When a holomorphic function has isolated singularities, a different boundary calculation detects their residues. Suppose $\Gamma$ is a counterclockwise simple closed contour and $f$ is holomorphic on a neighbourhood of it and its interior except at finitely many interior isolated singularities $a_j$. The residue theorem givesA residue is the coefficient of $(z-a_j{)}^{-1}$ in the Laurent series. A short proof idea is to deform the contour to small circles around the isolated singularities using the Cauchy integral theorem. Integrating the Laurent series on each circle kills every power except $-1$. ^[ref]

Here is how a complex pole determines a real improper integral. Set $g(z)=1/(1+z^2)$ and use the contour consisting of the segment from $-R$ to $R$, followed by the upper semicircle $A_R$, with $R>1$. Only the simple pole $i$ lies inside, andThe arc must actually disappear. On $A_R$, the reverse triangle inequality gives $|1+z^2|\ge R^2-1$, and its length is $\pi R$, soThe real integrand is bounded by $x^{-2}$ for $|x|\ge 1$, so both tails converge separately. Taking the limit therefore proves the ordinary improper integralFor example, at $R=10$ the real segment contributes $2\arctan (10)\approx 2.942255$, and the arc contributes $0.199337$, below the bound $10\pi /99\approx 0.317333$. The usefulness of the residue theorem is that the same strategy works for many rational functions whose real antiderivatives are cumbersome: find the enclosed poles, compute their residues, and control the connecting arc. ^[ref]

Historically, Cauchy's 1825 memoir on integrals with imaginary limits develops complex contour integration and its use in transforming and evaluating definite integrals. The opening sections define complex path integrals, discuss changing the path, and examine what happens near singular points. These are recognizable ancestors of the modern link between holomorphic functions, contours and residues. ^[ref]