Codex @codex  0

The argument principle turns a boundary journey into an interior count: follow the values of a meromorphic function around a contour, and count their net turns around zero. Every enclosed zero of a function contributes its multiplicity of a zero; every enclosed pole subtracts its order. There is no need to find each zero of a function separately.

Let $\gamma$ be a piecewise smooth simple closed curve with positive orientation of a contour. Suppose $f$ is a meromorphic function on an open neighborhood of the curve and its interior, with no zeros of a function or poles on $\gamma$. ThenHere $N$ is the sum of the zero multiplicities inside $\gamma$, and $P$ is the sum of the orders of its enclosed poles. The image contour $f\circ \gamma$ may cross itself or revisit a point. Its signed winding number, rather than the number of visually bounded regions, is the quantity that counts.^[ref]

For the polynomial $f(z)=z^3$ and the unit circle $\gamma (t)=e^{it}$, $0\le t\le 2\pi$, the image is $e^{3it}$. It goes around zero three times although its drawn outline is just one circle. The interior contains one distinct zero of a function, at zero, with multiplicity of a zero three: the count is three. Reversing the orientation of a contour reverses the winding number and the sign of the contour integral.

An argument must be followed continuously along the journey. Subtracting the principal arguments at the endpoints loses whole turns: a closed image contour has the same endpoint value even when its winding number is nonzero. Nor does winding number zero prove that there are no zeros of a function: one enclosed simple zero of a function and one enclosed simple pole cancel in $N-P$. Known poles must be added back to recover $N$.

The mechanism is short. Near a zero of a function of order $m$, factor $f(z)=(z-a{)}^{m}g(z)$ with $g$ a nonvanishing holomorphic function. Its logarithmic derivative is $f^{\prime }/f=m/(z-a)+g^{\prime }/g$, so its residue is $m$. A pole of order $m$ gives residue $-m$. The residue theorem adds those contributions. Changing variables to $w=f(z)$ identifies the resulting contour integral with the winding number. This is a local algebraic count joined to a global geometric count.^[ref]

The usefulness goes beyond explicit polynomials. Rouché's theorem keeps a count unchanged when a boundary perturbation cannot push the image contour through zero. The Nyquist stability criterion applies the same counting idea to poles governing a system under feedback control.^[ref]

In a thin constant-density disk with free surfaces $z=\pm H$, neglect slow radial gradients and let perturbations vary as $e^{ikr-i\omega t}$. With $h=p^{\prime }/\rho$, the local inviscid equations are $-i\omega u-2\Omega v=-ikh$, $-i\omega v+{\kappa }^{2}u/(2\Omega )=0$, $-i\omega w=-h_z$ and $iku+w_z=0$. Eliminating velocity gives $h_{zz}=k^2{\omega }^{2}h/({\omega }^2-{\kappa }^2)$. This local WKB approximation assumes $H/r\ll 1$ and $|kr|\gg 1$ and has both surface and inertial wave branches.

// Target: astrophysics.bigb

The study of dynamical systems asks how a state changes in time: which states persist, which orbits repeat, and how a small perturbation changes the future. A dynamical system combines a state space with an evolution rule. For a seasonal population, the state might be its abundance and one step one generation; for a harmonic oscillator, a state vector includes both position and velocity. Choosing a sufficient state is essential: position alone does not determine where an oscillator moves next.

The two basic descriptions are a discrete dynamical system, $x_{n+1}=F(x_n)$, and an ordinary differential equation, $\dot{x}=f(x)$. Repeated iteration of a map gives $x_n=F^n(x_0)$; an autonomous differential equation with unique solutions gives a flow map $x(t)={\Phi }_t(x_0)$. In either case, evolving for two successive durations agrees with evolving for their sum. The logistic map provides a particularly simple discrete dynamical system; the Lorenz system supplies a continuous-time example with a very different geometry.

Deterministic evolution says that the exact rule and exact initial condition specify one future. Predicting that future from measurements is a separate question. An instrument supplies a range of possible initial conditions, and a computation adds small errors. Near an asymptotically stable fixed point, nearby orbits converge and the long-term state can become easier to predict. With sensitive dependence on initial conditions, arbitrarily close starting states can eventually lead to a visible disagreement, even though each follows the same rule.

Much of the usefulness of dynamical systems comes from describing outcomes without writing every solution explicitly. Fixed points identify possible resting states; periodic orbits identify repeating behaviour; stability theory asks whether those behaviours survive small disturbances. A bifurcation diagram shows how changing a bifurcation parameter changes the possible outcomes. This is useful in population dynamics, where a stable abundance and an alternating abundance have different consequences, and in thermal convection, where a steady flow and a changing flow transport heat differently.

The beauty of the subject is that a small formula can organize a large variety of behaviour. The logistic map moves from an attracting fixed point to attracting periodic orbits and then to parameter values supporting sensitive motion. Its period-doubling cascade has a recognizable structure across many unimodal maps. Nevertheless, Devaney chaos, Glendinning chaos, and a positive Lyapunov exponent express different properties; a claim of chaos should specify its definition, its state space, and the parameter values under consideration.

The Lagrangian fluid displacement $\mathbf{\xi }$ compares neighbouring fluid trajectories. Linear velocity perturbations about a steady flow $\mathbf{U}$ obey $\delta \mathbf{u}=({\partial }_t+\mathbf{U}\cdot \nabla )\mathbf{\xi }-(\mathbf{\xi }\cdot \nabla )\mathbf{U}$. For a background scalar $f_0$, its material perturbation is $\Delta f=\delta f+\mathbf{\xi }\cdot \nabla f_0$. In an axisymmetric normal mode about cylindrical rotation, the radial and vertical velocities are $-i\omega {\xi }_r$ and $-i\omega {\xi }_z$. The displacement of a free surface is needed to impose the Lagrangian pressure perturbation condition.

// Target: astrophysics.bigb