Rack // Equations

Every Law on the Rack

616 lawsConstellation Gallery
616 / 616
  • The abelian propertyDhar 1990
    order-independent ⟹ group + fractal identity

    The final pile and the total number of topplings do not depend on the order the grains fall in. That commutativity turns the recurrent states into an abelian group — whose identity element is itself a striking fractal.

  • The adaptive belief marketfundamentalists against chartists
    rt+1 = μ[nfφ(F − pt) + ncδ tanh(g rt/δ)] + σεt

    A market maker moves the price by excess demand. Fundamentalists buy what is cheap against the fundamental F; chartists buy the trend they read off the chart, rt, with bounded positions. The shares n follow realised profit through the logit. At F the reflexive gain is μ·g·nc; past one the fundamental cannot hold.

  • The aether shiftobsolete
    Δ = 2Lλ · (v/c)2

    What a light wave in a medium at rest predicts: an interferometer moving through the aether at v sees its fringes slide by 2L(v/c)²/λ on turning a quarter — 0.4 of a fringe for Michelson and Morley’s 11 m. They found far less than a twentieth of it.

  • The alias foldthe illusion
    fperc = ffeat − fs ⌊ffeat/fs⌉N

    Anything above the limit folds back into (−fs/2, fs/2]. Divide by N and you get the spin you actually perceive — this is the line the instrument computes each frame.

  • The square–cube lawGalileo 1638
    A ∝ L², W ∝ L³ ⇒ σ = WA ∝ L

    Scale every length by L and a bone’s cross-section grows as L² while the weight it carries grows as L³. The stress — load per area — rises in proportion to size: a photocopied giant is a broken giant.

  • Galileo’s boned ∝ L^3/2
    d² ∝ L³ ⇒ d ∝ L3/2 ∝ M1/2

    To hold the stress constant the diameter must grow faster than the length — the femur of Galileo’s giant is stouter at every size. It works, and the bill is the skeleton: its share of the body grows as M1/3.

  • Rubner’s surface law2/3
    B ∝ A ∝ M2/3

    A warm body loses heat through its skin, so it should burn in proportion to its surface. The exponent Kleiber displaced — and that later re-analyses keep bringing back.

  • Kleiber’s law3/4
    B = B0 M3/4, B0 ≈ 3.4 W·kg−3/4

    Basal metabolic rate grows slower than mass and faster than surface: a whale is a hundred million times a mouse’s mass and burns a million times its power. B0 is Kleiber’s ≈ 70 kcal a day per kg3/4.

  • Heart-rate scaling−1/4
    f = QVstroke ∝M3/4M = M−1/4

    Cardiac output follows metabolism (M3/4); the stroke volume is a fixed share of blood volume (M). Their ratio is the heart rate: a mouse near 600 beats a minute, a whale near ten.

  • The power racea model
    PreqPavail ∝ W·VWe = W7/6 − e

    Power to fly grows as weight times speed, W7/6; flight muscle, a fixed share of the body, delivers ∝ We with e ≈ 1. The ratio creeps up as W1/6 and crosses one somewhere among the heaviest flapping birds.

  • Elastic similarityd ∝ h^3/2
    d ∝ h3/2, M ∝ d²h ∝ d8/3

    Grow so that every size sits the same fraction below its buckling height and the trunk must thicken as h3/2 — between the geometric twin (d ∝ h) that buckles and the stress rule (d ∝ h²) that never needs to.

  • Posture scalingBiewener 1989
    EMA ∝ M0.26 ⇒ σ ∝ M1 − 0.26 − 2·0.36 ≈ M0.02

    Larger mammals stand straighter, raising the effective mechanical advantage of their limbs. With near-isometric bones (d ∝ M0.36) that cancels almost all of the square–cube rise — until the legs are straight and EMA can grow no more.

  • Wing-loading scalingW^1/3
    W/S ∝ L³L² ∝ W1/3, W/S ≈ 47·W1/3

    Geometric similarity for flyers: weight is a volume, wing area a surface. Tennekes’s diagram runs one such line (N/m², W in N, approximate) from insects to airliners.

  • The WBE networkspace-filling
    β = n−a, γ = n−b ⇒ Ncap ∝ Vb1/(2a+b) → M3/4

    Branch n ways at every level; radii shrink by β, lengths by γ, down to a capillary that is the same in every mammal. Area-preserving (a = ½) and space-filling (b = ⅓) give 1/(2a+b) = 3/4: the number of capillaries — the power — grows as blood volume to the 3/4.

  • Allometrya power of size
    y = a·Mb ⇔ log y = log a + b·log M

    The Lemuel grammar: a quantity as a power of body size. On log–log axes every allometric law is a straight line whose slope is the exponent b — fitted, on every instrument of the edition, from the sizes you visit.

  • Cruise-speed scalingW^1/6
    V = (2 W/Sρ CL)1/2 ∝ W1/6

    The lift equation with wing loading ∝ W1/3: cruise speed grows as the sixth root of weight — a jumbo about a hundred times a fruit fly’s speed for 10¹¹ times its weight.

  • The disk integratorz = (1/a)∫y dx
    dz = (y/a) dx ⇒ z = (1/a) ∫ y dx

    A wheel resting on a turning disk at distance y from its centre is dragged round at y times the disk’s speed, so its shaft turns by the running integral of y against the disk’s rotation x. The radius is the integrand, the disk is the variable, the wheel is the answer.

  • Two integrators in a loopy″ = −ω²y
    dv = −ω y dt, dy = ω v dt ⇒ y″ = −ω2y

    Feed the second integrator’s output back, reversed, into the first one’s lead screw, and the shafts can only turn in the ratio the equation allows: y and v trace a circle. The machine does not compute the solution — the solution is the only motion it has.

  • Wardrop’s equilibriumequal cost on every used route
    xr > 0 ⇒ ℓr(x) = minq ℓq(x);  x100 + 45 = 2·N100 ⇒ 65 → 80

    At equilibrium every route anyone uses costs the same, and no unused route is cheaper — so no driver can gain by switching alone. On Braess’s diamond, 4 000 drivers split evenly at 20 + 45 = 65 minutes; open a free connector and the route through it (40 + 0 + 40 at full use) undercuts both old roads for each driver, until all 4 000 take it and every commute is 80. Both states are stable; the second is worse for everyone.

  • The price of anarchy≤ 4/3, linear latencies
    PoA = C(equilibrium)C(optimum) ≤ 43  (ℓ(x) = ax + b)

    The ratio of the total travel time at the selfish equilibrium to the least total any routing could achieve. With latencies linear in the flow it never exceeds 4/3, on any network. On the diamond it is 1 while the connector is used by everyone or by no one, rises with demand, and touches 4/3 exactly at 4 500 drivers (90 against 67.5 minutes); at 4 000 it is 80/64.69 = 1.237. A toll of 22.5 minutes on the connector alone — the extra delay a connector driver imposes on everyone else, beyond what an old-road driver does — moves the equilibrium onto the optimum (Pigou’s fix).

  • Choosing the delayfirst zero of ρ
    ρ(τ) = ⟨x(t)·x(t−τ)⟩σ²

    τ too small and the coordinates are near-identical — the cloud collapses onto the diagonal. τ too large and the fold tangles. The usable window is near the autocorrelation’s first zero (or the first minimum of mutual information): decorrelated, not yet folded.

  • Backpropagation∇L
    δl = (Wl+1ᵀ δl+1) ⊙ act’(zl)

    How to get ∇L cheaply: the chain rule, run backwards. The error at the output is propagated layer by layer, and each weight’s gradient is its incoming activation times the error flowing back through it.

    1986
  • One choice: the fullest binGonnet 1981
    max load ≈ Γ−1(n) − 3/2 ~ ln nln ln n

    Throw n balls into n bins uniformly at random and the fullest bin holds about this many — about 9 for a million bins. The asymptote is slow: at a million it says 5.3, and Γ⁻¹ is what the bins actually do.

  • The power of d choicesln ln n
    max load = ln ln nln d + O(1)

    Let each ball look at d bins and join the least loaded and the fullest falls from logarithmic to doubly logarithmic: 4 for a million bins with d = 2. Almost all the gain is in the second look; the third only divides by ln 3 instead of ln 2.

    1994
  • The supermarket taildoubly exponential
    si = λ(di − 1)/(d − 1), T = Σ si / λ

    Jobs arrive at rate λ per server, each joins the shortest of d sampled queues, service takes one unit on average. The fraction of queues with at least i jobs is si: geometric λi for one choice, doubly exponential for two. At λ = 0.9 the mean time in system falls from 10 to about 2.61.

  • Thompson samplingposterior sampling
    θa ∼ Beta(αa, βa); play argmaxa θa

    Keep a Bayesian belief (a Beta posterior) about each arm’s win rate, draw one sample from each, and play the arm whose sample is highest. Uncertainty makes wide posteriors likely to win a draw — exploration for free, from probability matching alone.

  • Regretcost of learning
    RT = T·μ* − 𝔼[ Σt μa(t) ]

    What you lose by not knowing the best arm from the start — the gap between always playing the true best (mean μ*) and what your policy actually earned. A good policy makes this grow like ln T, not T.

  • Upper confidence boundoptimism
    at = argmaxa [ Qa + c√(ln tna) ]

    Optimism in the face of uncertainty: judge each arm by the top of its confidence interval, not its mean. A rarely-pulled arm gets a big bonus, so it is tried; a well-known poor arm is dropped. No randomness, and provably logarithmic regret.

  • The run conditionwithdrawals beyond what the bank can pay
    fail ⇔ ℓ > ψ(θ) = m + (1 − m) θ

    ℓ is the share of deposits withdrawn early; the bank pays them from its liquid share m and then by selling loans worth θ on the unit. Under common knowledge every θ with 0 < ψ(θ) < 1 has two equilibria — nobody runs and it survives, everybody runs and it fails.

  • It is Bayesposterior
    p(x | z) ∝ p(z | x) · p(x)

    Why it’s optimal: for Gaussians, posterior ∝ prior × likelihood is again Gaussian, and its mean is precisely the Kalman update. The filter is Bayes’ rule, closed-form, once per tick.

  • Peak femur stressbeam model
    σ = 4Fπd²·[11 − ρ² + 8ε1 − ρ⁴], F = k M gEMA

    The femur as a cortical tube (inner/outer ratio ρ) loaded axially and bent by an eccentric load arm εd. The load is the muscle force that balances the ground reaction kMg through the limb’s mechanical advantage — so posture enters the stress directly.

  • The p-beauty contestclosest to p × the average
    x* = p · x̄,   xᵢ ∈ [0, 100] ⇒ Nash: x = 0 (p < 1), 100 (p > 1)

    Everyone names a number; the winner is closest to p times the average. For p < 1 no number above 100p can win, so none above 100p² can, and so on down: the only equilibrium is that everyone says zero.

  • Steps of reasoning50 · pᵏ
    xk = 50 · pk

    Start from a naive 50 and best-respond k times: 33, 22, 15 … for p = ⅔. A guess x implies a depth k = ln(x/50) / ln p. Nagel’s first rounds clustered near the first two steps.

  • Condensate fractionorder parameter
    N0/N = 1 − (T/Tc)3

    Below Tc the thermal density saturates, so every further-cooled atom is forced into the ground state. That exact curve is what the inset dot rides.

  • Critical temperatureharmonic trap
    kBTc = ℏω̄ (N/ζ(3))1/3

    In a harmonic trap the critical point follows the atom count and the trap stiffness — more atoms condense warmer. Slide N and the inset’s Tc line moves.

  • When it condensesdegeneracy
    n λdB3 = ζ(3/2) ≈ 2.612

    Condensation isn’t set by temperature alone but by phase-space density: when the de Broglie wave is as big as the spacing between atoms, nλ³ hits ζ(3/2).

    1925
  • One wavefunctionmacroscopic ψ
    ψ(r) = √n(r) eiφ

    The payoff: thousands of atoms share a single wavefunction — one amplitude, one phase. A quantum state grown to a size you can photograph.

  • CHSH inequality≤ 2 if local
    S = | E(a,b) − E(a,b′) + E(a′,b) + E(a′,b′) | ≤ 2

    Bell’s bound in testable form: four correlators at four angles. Every local-hidden-variable theory obeys S ≤ 2 — so one measured S above 2 rules them all out at once.

    1969
  • Tsirelson's bound2√2
    S ≤ 2√2 ≈ 2.828

    Quantum mechanics breaks the classical wall — but not infinitely. The most it reaches is 2√2, hit exactly at the optimal angles. Strange, yet still lawful.

    1980
  • The entangled statethe singlet
    |Ψ−⟩ = |↑↓⟩ − |↓↑⟩√2

    Two particles, one state — neither has a spin of its own, yet along any shared axis they always come out opposite. The whole is defined exactly where the parts are not.

  • Singlet correlationthe quantum curve
    E(a,b) = −cos(a − b)

    The average product of the two ±1 outcomes, fixed only by the angle between the detectors. It bows past every straight line a local model can draw — and that gap is the whole effect.

  • Bellman optimalitythe fixed point
    V*(s) = maxa [ r + γ Σs′ P(s′|s,a) V*(s′) ]

    The value of a state is the best action’s immediate reward plus the discounted value of wherever it lands you. Optimal play is exactly the solution of this self-referential equation — the present worth of the whole future, folded into one line.

  • Q-learning updatetemporal difference
    Q(s,a) ← Q(s,a) + α[ r + γ maxa′ Q(s′,a′) − Q(s,a) ]

    The one rule. Nudge the estimate toward the reward you just got plus the best value you now see ahead; the bracket is the TD error — the surprise. Learn the future from a better guess of the future, with no model of the world at all. This single line drives The Descent’s deep-RL descendants.

    1989
  • The action-value Qexpected return
    Q(s,a) = E[ Σt γt rt | s,a ]

    How good it is to take action a in state s and then play greedily forever after — the expected discounted reward. Learn this table and the policy is free: in every state, pick the action with the highest Q.

  • The bounce maparea-preserving
    (s, p) ↦ (s′, p′) , p = sin θ

    Birkhoff throws the flight between walls away: a billiard is fully told by two numbers per bounce — where it hit the wall and how steeply. The map preserves area, so what it draws in the (s, p) plane is an honest portrait of the dynamics: curves mean order, spray means chaos.

  • The mirror bouncethe whole rulebook
    θin = θout

    The one law on the table: the angle in is the angle out. Hero of Alexandria derived it around 60 AD from a least-path principle — light takes the shortest route via the wall — making this the oldest variational argument in physics. Everything else is the shape.

  • The hidden constantintegrable
    L₁ · L₂ = const

    In an ellipse the product of a chord’s angular momenta about the two foci survives every bounce. A trajectory can therefore never cross the confocal conic its constant pins — it wraps around it forever, and that untouchable curve is the glowing caustic the chords weave.

  • The shadow⌀ = √27 r_s
    bcrit = 3√3 GM/c2 ⇒ ⌀ = √27 rs

    The dark disc is bigger than the horizon itself: any ray aimed inside the critical impact parameter bcrit is swallowed. Lensing magnifies the rs horizon into a √27 ≈ 5.2 rs silhouette — the ring the Event Horizon Telescope photographed in 2019.

  • Nakamoto’s catch-up oddsNakamoto 2008 §11
    P = 1 − Σk=0z λk e−λk! (1 − (q/p)z−k), λ = z·q/p

    The attacker’s progress over the z confirmations taken as Poisson at the average pace, then the gambler’s ruin from wherever he is. 0.0009137 at q = 0.1, z = 5 — The Chain’s engine rebuilds the whole printed table. Run as an exact Poisson race the count is negative-binomial, and Rosenfeld’s form gives 0.00178 there.

  • The gambler’s ruin, as a catch-up(q/p)^z
    qz = 1 if p ≤ q · (q/p)z if p > q

    A random walk that steps toward the attacker with chance q and away with chance p ever closes a gap of z with probability (q/p)z — the gambler’s ruin. Below half the hash power the odds fall exponentially with every block; at half or above the attacker always catches up.

  • The selfish-mining thresholdEyal & Sirer 2013
    1 − γ3 − 2γ < α < ½

    A pool of share α, winning γ of the honest miners in a tie, earns more than its share by withholding blocks inside this range: above ⅓ when ties go against it, ¼ when they split, any size when they all go its way. The Chain simulates the state machine and lands on the revenue formula to the third decimal.

  • Three rulesReynolds boids
    v̂i ∝ wa·align + wc·cohere + ws·separate

    Three steering urges — match their heading, move toward them, but keep your distance — summed with weights. That is the whole brain; nothing in it knows the flock exists.

    1986
  • The Boltzmann factorheat to escape
    p(E) ∝ e−E / kBT

    In thermal equilibrium a state of energy E is occupied with weight e−E/kT. Temperature sets how readily the system pays energy to leave a comfortable valley — the physics annealing borrows and aims.

  • What we usually seeBorn rule
    P(x,y) = |ψ|² = R², ψ = R eiφ

    The cloud is the square of the wave. Squaring discards the phase — which is exactly the half this instrument paints back in as colour.

  • Leonardo’s ruleΣd² = d²
    dparentΔ = Σ dchildΔ, Δ = 2 (Leonardo) · Δ = 3 (Murray)

    With Δ = 2 the cross-section is conserved at every fork, so the branches of a tree squeezed together would make one trunk. Murray’s Δ = 3, the rule for vessels pumping fluid, thins the crown faster.

  • Brock–Hommes switchingintensity of choice β
    nh,t = eβUh,t−1Σk eβUk,t−1

    The share using predictor h follows its recent profit U — the rational predictor pays a cost C, the naive one is free. When the cobweb is unstable under naive expectations, raising β past β₁ = (2/C)·artanh(bd/bs) turns the steady price into a two-cycle, then a cascade of doublings into chaos: the crowd keeps leaving the expensive predictor while prices are calm, and coming back when they are not.

  • Brownian motionEinstein · Bachelier
    D = kBT / 6πηr

    Einstein 1905 tied this walk to Avogadro’s number — the jiggle of a pollen grain as proof that matter is atoms. Five years earlier Bachelier had modelled prices as the same walk: the seed of quantitative finance and every backtest.

  • Euler bucklingthe column
    Pcr = π² E I(KL)²

    A slender column carries its load straight until Pcr, then bows sideways all at once. Stiffness E·I grows as d⁴, so doubling the length quarters the load a column of the same girth can bear.

  • Greenhill’s heighth ∝ d^2/3
    hcr = (7.837 E d²16 ρ g)1/3 ∝ d2/3

    A column standing under its own weight w = ρgπd²/4 per metre buckles when w·h³/EI reaches 7.837. Tall and thin fails first: a trunk twice as thick may stand only 22/3 ≈ 1.6 times as tall.

  • Anchoring and adjustmentSterman 1989
    O = max(0, L̂ + α(S* − S) + αβ(SL* − SL)) · L̂ ← θ·L + (1 − θ)·L̂

    Anchor on the demand you expect, then adjust toward the stock you want (S*) and the supply line you should have (SL*). β is the weight a manager gives what is already on order: Sterman’s players gave it about a third. At β = 0 The Whip’s step game is still swinging at week 104; at β = 1 it settles.

  • The bullwhip boundChen et al. 2000
    Var(q)Var(D) ≥ 1 + (2Lp + 2L²p²)(1 − ρp)

    An order-up-to stage that forecasts by a p-period moving average and covers a lead time of L periods, under AR(1) demand of correlation ρ: equality when it keeps no safety stock. Long windows and short lead times damp the whip. In The Whip’s engine (L = order + shipping delay) the retailer measures 3.47 against 3.46 at L = 3, p = 4, ρ = 0.5.

    2000
  • The Byzantine boundn ≥ 3f + 1
    n ≥ 3f + 1 · with oral messages

    With oral messages, agreement among the loyal survives f traitors exactly when they are fewer than a third of the camp. The margin is a count: The Generals’ worst-case search breaks OM at n = 3, 4, 5, 6 for f ≤ 2 and finds no break at n ≥ 3f + 1. With unforgeable seals the third disappears: SM(m) handles any number of traitors up to m.

  • OM(m)’s courier countO(n^(m+1))
    M = (n−1) + (n−1)(n−2) + … + (n−1)(n−2)⋯(n−m−1) = O(nm+1)

    Round r of OM(m) relays every value along every path of r distinct generals, so the courier count is a falling product per round: 3 + 6 = 9 for OM(1) at n = 4, 6 + 30 + 120 = 156 for OM(2) at n = 7, 3 609 for OM(3) at n = 10 — all counted by the engine. Signed SM(m) relays each order once per lieutenant instead.

  • Four classesI · II · III · IV
    30 → III · 90 → II · 110 → IV

    The same eight bits sort into four behaviours — uniform, periodic, chaotic, complex. Rule 110 sits in the last and is Turing-complete: a universal computer hiding in three bits.

  • Local law3 inputs
    sit+1 = φ(si−1, si, si+1)

    Each cell sees only itself and its two neighbours, and a single local function decides its next state. No cell knows the whole world — all behaviour is local.

  • The rule spaceall of them
    φ : {0,1}3 → {0,1} ⇒ 28 = 256

    Three binary inputs make eight possible neighbourhoods, and a rule just picks an output bit for each — so the entire space of laws is countable and small.

  • The fall of caloricobsolete
    W = C · (Th − Tc) · Cin = Cout

    Carnot’s waterfall: caloric C, conserved, falls from the hot body to the cold one, and the motive power is the caloric times the height of its fall, as water turns a mill. False as physics — heat is not conserved — yet exact as bookkeeping, because the quantity his axioms define is Q/T.

  • Jurin’s lawh ∝ 1/r
    h = 2γ cos θρ g r

    A wetting liquid climbs a tube until the column’s weight matches the pull of its meniscus: halve the bore, double the rise (slope −1). Non-wetting mercury is pushed down instead.

  • Young’s equationthe contact angle
    γsv = γsl + γ cos θ

    The three surface tensions meeting at the contact line balance along the plate, fixing the angle the liquid meets it at — near 160° on a lotus leaf, near 10° on clean glass.

  • The Young–Laplace equationthe pressure jump
    Δp = γ(1R1 + 1R2)

    A curved liquid surface holds a pressure jump proportional to its surface tension and its total curvature. Add gravity’s hydrostatic pressure and integrate from the apex, and it is the whole shape of a drop, a meniscus or a soap film.

  • Capillary lengthBo = 1
    ℓc = (γ / ρg)1/2 ≈ 2.7 mm (water)

    The size at which surface tension and gravity are equal. Below it a drop is a spherical cap; above it, a puddle whose height stops at 2ℓc sin(θ/2).

  • Bond numbergravity ÷ tension
    Bo = ρ g R²γ = (R/ℓc)²

    Gravity on a volume against tension on a surface — the square–cube race for liquids. The crossover the drop’s fit may not cross is Bo = 1.

  • The Carnot boundthe limit
    η = 1 − TcTh

    The most any engine working between two temperatures can turn to work. Carnot showed in 1824 that it depends on the two temperatures alone; written on the absolute scale it is the second law’s first number, and it survived the caloric it was derived from.

  • The cat mapdet = 1
    (x, y) ↦ (2x + y, x + y) mod 1

    One integer matrix on the torus: shear, shear, wrap. The determinant is 1, so no information is ever lost — every pixel of the picture survives every step. It only stops looking like a picture.

  • The recurrence periodΠ(N) ≤ 3N
    MΠ(N) ≡ I (mod N)

    On an N×N pixel grid the map is a permutation, so the scrambled picture must come back exactly. Dyson & Falk bounded the wait at 3N steps — but the schedule is wild: 101 pixels return in 25 steps, 100 pixels take 150.

  • The stretchλ₊ = φ²
    λ± = 3 ± √52 , ln λ+ ≈ 0.9624

    The map’s eigenvalues are set by the golden ratio: along one irrational direction every gap is multiplied by φ² ≈ 2.618 per step, along the other it shrinks by the same factor. That one number is the Lyapunov exponent — the shredding rate — of the whole system.

  • The catenaryy = a cosh(x/a)
    y = a cosh(x/a),   a = H/w

    A chain of weight w per unit length, pulled horizontally by H, hangs in a hyperbolic cosine. Its slope grows with the weight of chain below each point — a parabola would need the weight spread along the ground instead.

    1691
  • Central Limit Theorem→ 𝒩(0,1)
    (Sn − nμ) / σ√n → 𝒩(0,1)

    Standardise the sum of any finite-variance step law and it converges to the same standard Gaussian. The bell is the universal attractor of sums — the reason it is everywhere.

    1810
  • Channel capacityBSC
    C = 1 − H(p)

    The bits per use you can send with vanishing error. Noise eats capacity; at p = ½ the channel carries nothing at all — input and output share no information.

  • Choreographiescentral configs
    three masses, one shared ∞-curve

    A measure-zero set of exact periodic orbits survives the chaos — Euler’s collinear line, Lagrange’s equilateral triangle (the origin of the L-points), and the figure-eight.

  • The four-timestamp offset±RTT/2
    θ = (T₂ − T₁) + (T₃ − T₄)2, |θ − θtrue| ≤ (T₄ − T₁) − (T₃ − T₂)2

    Send at T₁ by your clock, the server stamps T₂ and T₃, the reply lands at T₄. The midpoint estimate of the offset is wrong by at most half the round trip — whatever the split between the two directions — and Cristian’s version (T₂ = T₃) is the same bound.

  • The cobweb theoremwhich slope is steeper
    |bs / bd| < 1 ⇔ the spiral converges

    The spiral tightens when demand answers price more strongly than supply does, widens when supply is the steeper, and repeats a fixed square when they match. Nothing outside the market is needed for it to cycle for ever.

  • The cobweblast price plants this crop
    pt = a − c − bspetbd,   pet = pt−1 ⇒ pt − p* = −bsbd(pt−1 − p*)

    Demand a − bdp clears against supply c + bspe planted for the price farmers expect. With naive expectations — last season’s price — each deviation from the equilibrium p* comes back multiplied by −bs/bd: a glut, a shortage, a glut.

  • The cognitive hierarchyPoisson(τ) levels
    f(k) = e−ττkk!,   gk = p · Σh<k f(h) ghΣh<k f(h)

    Thinkers of level k believe everyone else thinks less, spread in proportion to the Poisson shares below k, and best-respond to that mix; level 0 picks at random. One parameter, τ, the mean number of steps. With τ = 1.5 and p = ⅔ the model crowd’s average is 33.5 and its winning number 22.3.

  • The double dixie cupm sets
    E[Tm] = N ∫0∞ [1 − (1 − Sm(t)e−t)N] dt,  Sm = Σj<m tj/j!

    The draws needed until every coupon has appeared m times — N ln N + (m − 1)·N ln ln N + O(N). The second complete set costs only the ln ln term, which is why two collectors who swap duplicates each finish at about two thirds of the solo cost.

  • The collector’s tailGumbel
    P(T ≤ N ln N + cN) → e−e−c

    Completion times are not bell-shaped: centred at N ln N with width N, they follow the double exponential of extreme values, with a long right tail. A thousand albums finish over a spread of several N, and one is always out near twice the mean.

  • The coupon collectorharmonic
    E[T] = N·HN ≈ N ln N + γN + ½  ·  gapk = NN − k + 1

    Drawing uniformly from N equally likely coupons, the k-th new one costs N/(N − k + 1) draws on average — one at the start, N at the end. The total is N times the harmonic number: the last coupon alone costs as much as the first N/2 together.

  • The elastic collisionnothing lost
    v′ = (m − M) v + 2M Vm + M

    Two conservation laws — momentum and energy — pin the outcome of a head-on elastic collision completely: no fitting, no friction, no choice. Huygens settled the rule in 1669; every clack in this instrument is this one line, applied exactly.

  • The countπ, in clacks
    N = ⌊π / θ⌋ , tan θ = √(m / M)

    Each clack turns the state around the circle by the same angle 2θ, so the total count is how many fit in a half-turn: π/θ. Set the mass ratio to 100^(d−1) and θ ≈ 10^(1−d) radians — and the count spells the first d digits of π.

  • The hidden circlephase space
    (√m v)² + (√M V)² = 2E

    Rescale each velocity by the square root of its mass and energy conservation draws a circle; momentum conservation draws a line. Every collision is then a reflection — the messy clatter of blocks becomes a polygon of equal chords inscribed in one circle.

  • The von Kries transformthree gains
    (L′, M′, S′) = (L · Lw′Lw, M · Mw′Mw, S · Sw′Sw)

    Discount a light by scaling each cone signal by the ratio of the target white to the assumed light’s white. Assume a different light and the same pixels become a different paint — the whole of the dress, in three ratios.

  • The daylight locusJudd et al. 1964
    yD = −3.000 xD² + 2.870 xD − 0.275

    The chromaticities of typical daylight, from warm late sun to blue north sky, with xD a cubic in 1/T (4000–25 000 K). The visual system’s guesses about the light live along this curve; D65 sits on it at 6504 K.

  • ℓ₁ recoverythe sparsest fit
    min ‖Ψᵀx‖1 subject to ‖Ax − y‖2 ≤ ε

    With fewer measurements than unknowns, infinitely many images fit y. Ask for the one whose coefficients in a basis Ψ have the smallest ℓ₁ norm: the corners of the ℓ₁ ball pick out sparse solutions, and if the scene really is sparse in Ψ, it is the scene. The prior is the basis — name it.

    2006
  • How many measurementsnot N
    M ≳ C · k · log(N / k)

    A k-sparse signal of length N is recovered from on the order of k log(N/k) random measurements — the scaling, not the constant. The Donoho–Tanner phase diagram draws where the bet actually pays.

  • Shrinkage step (ISTA / FISTA)gradient, then shrink
    c ← Sλ/L(c − 1L ΨᵀAᵀ(AΨc − y)), Sτ(v) = sign(v)·max(|v| − τ, 0)

    A gradient step on the misfit, then soft-thresholding: every coefficient shrinks toward zero by λ/L, and the small ones land on it exactly. FISTA adds a momentum term and converges as 1/k² rather than 1/k (the plain iteration is Daubechies, Defrise & De Mol, 2004).

  • The Turing machinea finite table
    δ : (state, symbol) → (symbol, move, state)

    A read/write head on an unbounded tape, steered by a finite table of rules — and that is the whole of what is computable. The universal machine reads another machine’s table off the tape and runs it: the first proof that one device can do anything any device can.

  • The halting problemundecidable
    halt(p, x)? — no algorithm decides it

    Will program p, run on input x, ever stop? Turing’s diagonal argument proves no single program can answer for all p — the first problem shown to have no algorithm. The busy beaver is its shadow: knowing S(n) would solve halting.

  • Wang tilesedge-matching
    a ▸ b ⟺ right-edge(a) = left-edge(b)

    Square tiles with coloured edges that may sit side by side only when the touching edges match. Wang asked the innocent question — can a given set tile the plane? — and conjectured that any set that can, can do so periodically. He was wrong, and the wrongness runs deep.

  • The busy beaveruncomputable
    S(n) = most steps a halting n-state machine runs

    Among all n-state machines that do halt, how long can the longest run? S(n) grows faster than any computable function — so no program can compute it. S(5) = 47,176,870 was settled only in 2024; S(6) already dwarfs every number in physics.

  • The domino problemundecidable
    “does this tile-set tile the plane?” — undecidable

    Berger proved no algorithm can decide it — by encoding a Turing machine into tiles, so “tiles forever” means “never halts”. The proof spits out the first aperiodic set: tiles that cover the plane but never repeat. Penrose shrank it to two; the hat, in 2023, to one.

  • Langton&rsquo;s anttwo rules
    white → turn right · black → turn left · flip, step

    On a white square turn right, on a black square turn left; either way, flip the square’s colour and step forward. That is the entire program — and from it comes ten thousand steps of apparent chaos, then, with nobody asking, an ordered “highway” that runs forever.

    1986
  • Computational irreducibilityno shortcut
    future ⟸ run every step

    For many simple programs there is no formula for the answer that is faster than just running them. The ant proves it cleanly: nothing in the rule foretells the highway — you can only watch. Wolfram’s claim is that this is the rule, not the exception.

    2002
  • Wave Function Collapsemin-entropy
    collapse argmin H · H = −Σ pi log pi

    Every cell starts as a superposition of all tiles. Collapse the most-constrained cell (lowest Shannon entropy) to one tile, propagate the adjacency rules to its neighbours, and repeat — a wave of certainty spreading by constraint, borrowing its name loosely from quantum measurement.

    2016
  • The turmitea 2-D Turing machine
    δ : (state, colour) → (colour, turn, state)

    Generalise the ant — more internal states, more colours — and you have a Turing machine whose tape is the plane. The Tape’s head crawling a line and the Ant crawling a grid are the same creature; some turmites are themselves universal computers.

  • The host’s update1/3 → 2/3
    P(prize behind the shut door | she opened o) = 1 · ⅓½ · ⅓ + 1 · ⅓ = 23

    Bayes’ rule on the host’s protocol, not on the empty door: if your door hides the prize she opens either of the other two (½); if the other closed door hides it she has exactly one legal door (1). That forced hand carries the closed door to 2/3. A host who opened at random and happened to show an empty door would leave ½ · ½ — a different game.

  • The improper priorno uniform on ℝ⁺
    E[other | A] = ½·2A + ½·A2 = 1.25 A  ⇐  P(A) = P(A/2) ∀A,  yet Σk P(2k) = ∞

    The swap argument needs “twice or half, equally likely” at every amount you could hold. That is a uniform weight on the doubling ladder, and its total is infinite: no probability distribution can give it. Under any real prior the odds of doubling move with the amount — and the 1.25× disappears.

  • Additive increase, multiplicative decreaseAIMD
    on ack: W ← W + 1W · on loss: W ← W2

    A sender keeps W packets in flight and grows W by one packet each round trip until the network drops one, then halves it. Chiu & Jain proved it the one linear rule that converges to a share both efficient and fair; the sawtooth is what it looks like.

  • The Mathis formularandom loss p
    B ≈ MSSRTT · C√p, C = √(3/2)

    A long-lived AIMD flow that loses a fraction p of its packets settles to this throughput: the window’s sawtooth between W/2 and W carries 1/p packets per cycle. Longer round trips and more loss both cost, whatever the wire’s speed.

  • Noether’s energy lawwhy it holds
    ∂L/∂t = 0  ⇔  dE/dt = 0

    Emmy Noether’s theorem: every continuous symmetry of the laws of motion yields a conserved quantity — and the symmetry of time (the laws are the same today as tomorrow) is exactly energy conservation. A wheel cannot gain energy from nothing because its physics does not know what time it is; run the same release a century later and the curves overlay. This is why the perpetual-motion wall stands, not merely that it does.

  • Conserved integralsonly 10 (Bruns)
    E, P, L = const

    Integrated by velocity-Verlet (leapfrog) — symplectic, so energy stays bounded over thousands of periods. The drift readout proves the orbits are the law’s, not the solver’s.

  • Barbier’s theoremP = πw
    P = ∮ ρ dθ = ∫02π h(θ) dθ = π w

    Every curve of constant width w — circle, Reuleaux polygon or anything stranger — has perimeter exactly π times its width. The perimeter is the integral of the support function, and only its constant term, w/2, survives a full turn.

  • The Blaschke–Lebesgue theoremleast area
    A ≥ ½ (π − √3) w2

    Of all curves with the same constant width, the Reuleaux triangle encloses the least area; the circle, π w²/4, encloses the most. Same width, same perimeter — about an 11% spread in area.

  • Constant width, by support functionodd harmonics
    h(θ) + h(θ + π) = w, ρ = h + h″ ≥ 0

    h(θ) is the distance from the centre to the tangent line facing direction θ; the caliper width is h(θ) + h(θ + π). It is the same in every direction exactly when h is w/2 plus odd harmonics only. The radius of curvature h + h″ must stay non-negative, or the curve folds over itself and is no roller at all.

  • The Watts drill’s squareswept area
    Aswept = (2√3 + π/6 − 3) w2 ≈ 0.98770 w2

    A Reuleaux triangle of width w can turn inside a square of side w touching all four sides at once. One full turn sweeps everything but four thin slivers at the corners — a hole that is square to within 1.23% of its area.

  • The constructible towerdeg = 2ᵏ
    [ℚ(construction) : ℚ] = 2k

    Every compass-and-straightedge step intersects lines and circles already drawn — at worst a quadratic — so each new point at most doubles the degree of the field it lives in. Every constructible number therefore sits at some finite floor 2ᵏ over ℚ; nothing built this way can ever reach degree infinity.

  • The cooling scheduleanneal, don't quench
    Tk = α·Tk−1,  0 < α < 1

    Temperature falls a little each sweep. Cool slowly (α → 1, anneal) and the search settles in the global minimum; cool fast (quench) and it freezes wherever it happens to be — a defect.

  • Cooperation fractionorder parameter
    fC = 1N Σi 𝟙[si = C]

    The share of the grid that cooperates — the population’s order parameter, exactly as magnetization is for The Threshold. On the spatial dilemma it settles near 0.318 across a band of b: neither all-C nor all-D, but a living mixture.

  • Holonomy — Gauss–BonnetΔψ = ∬K dA
    Δψ = ∬D K dA

    Carry a direction round a closed loop without ever turning it relative to the surface — parallel transport — and it comes home rotated by the total curvature inside the loop: zero on a plane, the enclosed area on a unit sphere, negative on a saddle.

  • Spherical excessA = R²(Σα − (n − 2)π)
    A = R2 (Σαi − (n − 2)π)

    A polygon of great-circle arcs on a sphere has angles that add up to more than a flat polygon’s; the excess, times R², is its area. The octant has three right angles: an excess of π/2, an eighth of the sphere.

  • The longest-cycle tail1 − ln 2
    P(longest = k) = 1k (k > n/2)  ⇒  P(survive) = 1 − Σk=51100 1k = 0.3118 → 1 − ln 2

    A prisoner who opens their own number’s box and then follows each slip walks exactly the cycle of the shuffle that contains them, and finds their slip within fifty opens if that cycle is no longer than fifty. So all hundred live or die together, on one property of one permutation: no cycle longer than half. A uniform shuffle of n has a cycle of length k > n/2 with probability exactly 1/k — at most one such cycle can exist — so the team fails with probability H₁₀₀ − H₅₀ and survives 31.18 % of the time, against (½)¹⁰⁰ ≈ 8 × 10⁻³¹ for a hundred independent guesses.

  • Matter is a wavede Broglie
    λ = hp

    Matter has a wavelength set by its momentum — the reason a particle can interfere at all. Heavier or faster ⇒ shorter λ ⇒ finer, harder-to-see fringes.

    1924
  • Matter wavede Broglie
    λdB = h√(2π m kT)

    Cold atoms have longer matter waves. Cool far enough and each atom’s wave grows until it overlaps its neighbours — the gas stops behaving like particles.

  • Wiener deconvolutionone division
    X̂ = H* · Y|H|² + Sn/Sx

    Divide by the blur where the signal is stronger than the noise and give up where it is not; the signal-to-noise per frequency is the prior, here estimated from the data’s own spectrum. Beyond the cutoff H = 0 and the filter returns exactly nothing.

  • Richardson–LucyPoisson, positive
    x ← x · Hᵀ[yHx + b]

    The maximum-likelihood iteration for photon counts: multiplicative, so it stays positive. Its misfit falls without end — past the noise, into fitting it — so it is stopped where χ²/N first reaches 1 (Morozov’s discrepancy principle). Positivity can extrapolate a little past the cutoff for point stars: that is the prior, not the sky.

  • The imaging equationblur + counts
    y = Poisson(PSF ∗ x + b), PSF = |F[P · ei2πW/λ]|²

    The image is the sky convolved with the point-spread function — the squared Fourier transform of the pupil P, with the wavefront error W as its phase — plus photon noise. For Hubble in 1990, W was a spherical aberration from a mirror ground to the wrong shape.

  • The optical cutoffthe null space
    OTF(f) = (P ⋆ P)(λf) = 0 for |f| > D/λ

    The transfer function is the pupil’s autocorrelation, and a pupil of diameter D has none beyond D/λ: finer detail reaches the detector multiplied by exactly zero. Deconvolution can divide by a small number; it cannot divide by zero.

  • The Chandrasekhar massα_G^−3/2 m_p
    MCh ≈ αG−3/2 mp = (ħcG)3/21mp² ≈ 1.9 M☉

    Relativistic electrons push back only as hard as gravity pulls, at every radius alike: above this mass a cold star has no size at which it can stand. With the full prefactor and two nucleons per electron, ≈ 1.44 solar masses.

  • The largest planet(α/α_G)^3/2
    a0N1/3 = a0ααGN−1/3 ⇒ Mpeak ≈ (ααG)3/2 mp

    Atomic matter keeps its density, so a ball of N proton masses grows as N1/3; squeezed by its own gravity, its electrons become degenerate and it shrinks as N−1/3. The two meet near the mass of Jupiter — no cold ball can be wider than about that.

  • The smallest starα^3/2
    Mmin ≈ α3/2 (mpme)3/4 αG−3/2 mp

    The hottest a ball of N protons can get before degeneracy stops the squeeze is kT ≈ αG²N4/3mec²; hydrogen burns near α²mpc². Lighter balls stay brown dwarfs. As an order of magnitude ≈ 0.3 M☉, against a measured ≈ 0.08.

  • Delay embeddingone signal
    v(t) = ( x(t), x(t−τ), …, x(t−(m−1)τ) )

    Stack a signal against delayed copies of itself to build an m-dimensional vector. As t runs, the point v(t) traces a curve — and that curve is a faithful copy of the hidden attractor you never measured directly.

    1980
  • The density fieldthe contrast
    δ(x) = ρ(x)/⟨ρ⟩ − 1

    Forget absolute density; what matters is the contrast — how much denser than average each patch is. The early universe was nearly smooth, δ ~ 10⁻⁵ everywhere. Gravity does the rest.

  • When η is too largeη ↑ → ∞
    overshoot → oscillate → diverge

    The same descent, pushed too hard, is the route to chaos. Raise η and the step overshoots the minimum, then overshoots back — a period-2 oscillation that doubles into divergence, exactly the bifurcation The Cascade draws.

  • The fixed pointself-reference
    G ⇔ ¬Prov(⌜G⌝)

    The same diagonal, aimed at proof. Arithmetize the syntax so that a formula can talk about provability; the fixed-point lemma then builds a sentence G that asserts “I have no proof.” If the system proved G it would prove a falsehood — so, if consistent, G is true and unprovable. The flipped diagonal, one register up.

  • The diagonalthe flip
    D(p) ≔ if H(p, p) halts, loop · else halt

    Cantor’s trick, aimed at machines. Suppose H decides halting for every program. Build D that asks H about itself and does the opposite. Then D(D) halts iff H says it doesn’t — H(D,D) = ¬H(D,D), a flat contradiction. No total H exists; every decider ships with its own nemesis.

  • The chariot’s differentialψ = (r/W)(θ_R − θ_L)
    Δψ = −(r/W)(ΔθR − ΔθL)

    When the cart turns, its outer wheel rolls further than its inner; the differential turns the figure back by the difference over the track W. On a flat floor that exactly cancels the cart’s turn, so the arm keeps its bearing.

  • The forward processnoising
    dx = −½ β(t) x dt + √(β(t)) dW

    Melt the data: drift each sample gently toward the origin while injecting Gaussian noise, step after step, until after enough time every trace of structure is gone and the cloud is indistinguishable standard noise. It is The Walk’s Brownian motion with a leash.

    1908
  • The reverse processAnderson 1982
    dx = [ −½βx − β ∇log pt(x) ] dt + √β dW̄

    Run time backward and the same diffusion becomes generative: the score term steers the noise uphill in probability, condensing formless static back into samples. Anderson proved a reversed diffusion is itself a diffusion — the theorem the whole field stands on.

    1982
  • The score∇ log p
    s(x,t) = ∇x log pt(x)

    The gradient of the log-density — at every point, the direction toward where data is denser. A real model learns it by denoising (gradient descent — The Descent); here we use the exact score of a known target, so the reverse flow is honest with nothing trained.

    2019
  • Probability-flow ODEdeterministic twin
    dxdt = −½βx − ½β ∇log pt(x)

    Drop the noise term and a plain ODE carries the same changing density: one fixed noise seed maps to one sample along a smooth path — the streamlines you watch the particles ride here. Same marginals, no randomness.

  • Discrete Laplacian3×3 stencil
    ∇²u ≈ 0.2·(N+S+E+W) + 0.05·(NE+NW+SE+SW) − u

    Diffusion on the grid is a weighted blur of the eight neighbours minus the centre — the only spatial coupling in the whole system.

  • Benford’s first-digit lawP(1) = 30.1 %
    P(d) = log10(1 + 1/d) · MAD = ⅑ Σd |p̂d − P(d)|

    Numbers made by multiplicative, scale-invariant growth spread evenly in the logarithm, so the leading digit d turns up with probability log₁₀(1 + 1/d) — a 1 thirty per cent of the time, a 9 under five. The mean absolute deviation of a ledger’s first-digit shares from the curve is Nigrini’s conformity measure; invented numbers fail it, because people make numbers by adding and spread them evenly.

    1938
  • The Clark–Evans ratiochance: R = 1
    R = d̄ / E[d] · E[d] = ½√(A/n) + (0.0514 + 0.041/√n) · P/n

    The mean distance from each point to its nearest neighbour, over what complete spatial randomness of the same density predicts — ½√(A/n) for n points in area A, plus Donnelly’s 1978 term for the points a boundary of perimeter P robs of neighbours. Chance sits at one; clumping pulls R below it, and anything that holds points apart — a rake, a planting, territory — pushes it above.

  • The switch biasfair: s = ½
    s = #(bi ≠ bi+1) / (n − 1) · (n − 1)s ~ Bin(n − 1, ½) for a fair process

    The share of neighbouring pairs that switch. Independent fair bits switch half the time, with an exact binomial spread; people asked to be random switch more — they under-repeat and avoid runs — and the bias survives being told about it. It is the simplest case of a Markov transition table, and the one a penny-matching machine needs.

  • The two-thirds power lawβ = −⅓
    v = γ · κ−1/3 ⇔ A = K · C2/3 · β = slope of log v on log κ

    A fluent hand slows in its bends by a fixed law: speed v goes as the path’s curvature κ to the −⅓ (angular velocity A as curvature C to the ⅔), the gain γ changing only from stroke to stroke. Fit the slope of log speed against log curvature and a written signature lands near −⅓; a tracer, creeping along the same line at one speed, lands near zero. The picture is the same — the binding of speed to shape is not.

    1983
  • Gradient kurtosisGaussian: κ̂ = 0
    κ̂ = m₄ / m₂² − 3 · over d = I(x + 1, y) − I(x, y)

    Take every step between horizontally neighbouring pixels and ask how heavy its tails are: the fourth moment over the square of the second, less the Gaussian’s three. A photograph steps gently almost everywhere and leaps at its edges, so the steps pile into a sharp peak with long tails (κ̂ ≫ 0). A field with the same spectrum but shuffled phases has nowhere to keep an edge; its steps fall into a bell.

  • The variance ratiowalk: VR = 1
    VR(q) = Var(rt + … + rt+q−1) / (q · Var(rt)) = 1 + 2 Σk<q (1 − k/q) ρk

    If steps are uncorrelated their variances add, and a q-step move has exactly q times the variance of one step — a ratio of one at every horizon. Persistence (positive autocorrelations ρk) compounds the ratio above one; mean reversion pulls it below. Lo and MacKinlay’s overlapping, bias-corrected estimator made it a test of the random walk.

  • The law of the longest runE[L] ≈ log₂ n
    P(Ln ≤ m) = Σj≤m qn(j) · qt+1(1) = ½ Σj qt(j), qt+1(j+1) = ½ qt(j)

    The longest streak of either face in n fair flips grows only like log₂ n — about eight in two hundred — and its law is exact by dynamic programming over the length j of the current run: each flip repeats the face (½, the run grows; past m the sequence is lost) or switches (½, back to one). A tape of two hundred that never reaches five happens to a fair coin less than once in a thousand.

  • The clustering coefficientfriends of friends
    Ci = 2Ti / (ki(ki − 1)) · C = (1/N) Σi Ci

    For each node, the share of its pairs of neighbours that are linked — Ti the triangles through it, ki its degree — averaged over the network. Growth through introductions leaves triangles behind; a degree-preserving shuffle keeps every degree and sinks C to the floor its hubs allow by accident, about (⟨k²⟩ − ⟨k⟩)² / (N⟨k⟩³).

  • Permutation entropyi.i.d.: H = 1
    H3 = −Σπ p(π) log p(π) / log 3! · π ∈ the six orderings of (xt, xt+1, xt+2)

    Sort every run of three consecutive samples into one of the six ways three numbers can be ordered, and take the entropy of the census. An i.i.d. sequence speaks all six equally and scores one, whatever its values; a deterministic one has a grammar — the logistic map at r = 4 never falls twice running, a word it cannot say.

  • The long-range echolag ≥ 6
    echo = maxℓ ≥ 6 (longest run of t with pt = pt+ℓ)

    Slide a line against itself and find the longest stretch where it lands on its own notes, at a distance of six or more. A composed melody returns to its phrases — the repeat is architecture; a Markov chain, however well it knows each step, forgets within a few and cannot walk back. Dependence in real music decays slowly with distance; in a chain it dies off exponentially.

  • The distance estimateHart 1989
    DE = ½ · r ln r|z′|, r = |z|

    A solid fractal has no formula for “distance to the surface”, but the escape-time potential gives a safe lower bound. The ray leaps by that bound each step — far away it strides, near the surface it creeps — so the whole bulb resolves in a few dozen steps per pixel.

  • Diffusion-limited aggregationstick on touch
    walker random-walks → freezes on contact

    Release one particle far away. It random-walks — Brownian, undirected — until it touches the cluster, then freezes there forever. Repeat a million times. A seed becomes a branching tree.

    1981
  • The screening instabilitytips win
    exposed tips ≫ shielded fjords

    A bump that pokes out intercepts more wandering walkers, so it grows faster and pokes out more — a runaway that sharpens tips and starves the bays between them. That instability is why the branches never fill in, and why they look the same at every magnification.

  • Terminal velocityweight = drag
    m g′ = ½ ρ Cd A vt² ⇒ vt ∝ √L

    A falling body speeds up until drag equals weight. Weight grows as L³, Newton’s drag as L²v², so the speed at which they balance grows only as √L — the mouse lands slowly, the horse does not. (g′ is gravity less buoyancy.)

  • Stokes dragRe ≪ 1
    F = 3πμLv ⇒ vt = Δρ g L²18 μ ∝ L²

    Below a Reynolds number of one, drag is viscous and linear in speed. Weight L³ against drag L·v gives v ∝ L²: halve a droplet and it falls four times slower — which is why fog hangs and bacteria never land.

  • The drag crossoverSchiller–Naumann
    Cd = 24Re(1 + 0.15 Re0.687), Re < 1000; ≈ 0.44 beyond

    One empirical curve joins Stokes (Cd = 24/Re) to Newton (Cd ≈ const). On the log–log plot of terminal speed it bends the slope from 2 down to ½ across a decade or two around Re = 1.

  • The truel recursionargmax = the air
    Pa first(a) = a1 − (1 − a)(1 − b)  ·  V(o₀) = h₀ + m₀h₁ + m₀m₁h₂1 − m₀m₁m₂

    Solve the endgames first: two left, the first shooter at accuracy a wins a / (1 − (1 − a)(1 − b)). With three standing, each turn either ends in one of those duels (hit) or passes the turn (a miss, or the air), so the three-standing states form a single cycle with an exact sum. At 0.3 · 0.8 · 1.0, firing at the Marksman leaves you 0.258, at the Gunner 0.237 — and into the air 729/2150 = 0.339, because the two strong shooters are each other’s largest threat and you enter every endgame holding the first shot.

  • The binomial bandexact 95 %, drawn before play
    K ~ Bin(n, p) · [k−, k+] : P(K < k−) ≤ 2.5 % ≥ P(K > k+)

    The exact central interval for K, the count of wins in n independent plays at win probability p — computed from the rules alone, before a single play. A record inside it is the rules working; a record outside it happens one time in twenty, and the band says so in advance. When the chance of each play is fixed before it is played but differs from play to play, the same band comes from the Poisson-binomial sum.

  • Field equationsEinstein 1915
    Gμν = 8πGc4 Tμν

    Matter tells spacetime how to curve; the curvature tells matter (and light) how to move. Gravity isn’t a force in the picture at all.

    1915
  • The El Farol rulego iff your best says < L
    goᵢ ⇔ Âi(t) < L,   ej ← λej + (1 − λ)|Âj − A|

    Agent i goes to the bar when the best predictor in its own bag — the one with the smallest running error e — forecasts attendance below the threshold L. Every predictor is rescored each week against what the week brought. A forecast everyone shares cannot survive the rule: believed below L it fills the bar, believed above it empties it.

  • Boltzmann entropyS = k ln Ω
    S = kB ln Ω

    Ω is the number of microstates that look like the current macrostate. The corner has vanishingly few; "spread out" has almost all of them — so that is where the gas goes. Engraved on Boltzmann’s tombstone.

  • InformationH(½) = 1 bit
    H(p) = −p log₂ p − (1 − p) log₂(1 − p)

    Surprise, in bits — the average doubt in one channel use. Maximal (one full bit) when p = ½, the channel a pure coin flip; zero when the outcome is certain. This is the thing a good code never wastes.

  • Landauer limitkT ln 2
    Emin = kBT ln 2

    Rolf Landauer’s 1961 result: erasing one bit of information in a thermal environment costs at least kBT ln 2 joules — irreversibly dissipated as heat. This is the bridge between information theory and thermodynamics: information is not just abstract, it is physical.

  • One equation, two worldsShannon ↔ Boltzmann
    H = −Σ pᵢ log pᵢ

    The same sum The Arrow carves as heat (S = k·ln Ω): there it counts microstates and points time forward, here it counts bits and bounds what a channel can carry. Entropy and information are one measure, read in two worlds.

  • Coarse-grainedGibbs / Shannon
    S = −kB Σ pi ln pi

    What the instrument plots: coarse-grain the box into cells, count the occupancy pᵢ, sum. It is maximal at a uniform spread — the Gibbs/Shannon form of the same quantity.

  • The equantobsolete
    C = e·â · Q = 2e·â · ∠AQP = λ̄ − λA (uniform)

    Ptolemy’s speed fix: the epicycle’s centre P rides a circle about C, e off the Earth, but turns uniformly about Q, twice as far — the motion is uniform about a point that is not the centre. On a sky of Mars it takes the worst error from six degrees to about one.

  • It is percolationT = β/(β+γ)
    outbreak ⇔ giant component of occupied edges

    The deep identity: an outbreak is the giant connected component of “occupied” edges, each present with probability T. Vaccination is site percolation; transmissibility is bond percolation. The epidemic threshold is a percolation threshold.

  • Explore vs. exploitε-greedy
    a = argmaxa Q(s,a)  w.p. 1−ε,  else random

    Most of the time take the best action you know; a fraction ε of the time gamble on a random one — because the only way to find a better path is to risk a worse one. Too little ε and the agent locks onto the first route it finds; too much and it never commits.

  • The two averagesensemble ≠ time
    ⟨Wn⟩ = ⟨r⟩n · Wtyp ≈ en ⟨ln r⟩

    The crowd’s average compounds the arithmetic mean of the factor; any one life compounds the geometric mean, which is always smaller. When ⟨r⟩ > 1 > e⟨ln r⟩ the same coin enriches the ensemble and ruins almost every member of it — both verdicts exact, both forever.

  • The multiplicative coinone line
    Wn+1 = rn · Wn

    Multiply your stake by a random factor, repeat — the whole law. No strategy, no memory, no choice. The line never says who takes the average of the outcomes — over copies of you, or over your own tosses — and that silence is worth a sign flip.

  • The all-pay laddersolved backward
    raise at h ⇔ (B − h)/s odd;  Vopen = −s;  x* = [(b − 1) mod (v − 1) + 1]·s

    Both top bidders pay. With fixed steps s and bankroll B, solve from the top rung down: the last affordable raise wins, so the verdict alternates rung by rung all the way to the first — and opening is worth −s. Allow jump bids and a single committed opening of x* (in units of s, with budget b and prize v) takes the note while the other rationally passes: O’Neill’s solution, commitment over reaction.

  • The nucleation barrierΔG* ∝ 1/ΔT²
    ΔG* = 16πσ33Δg2 · J ∝ e−ΔG*/kT

    A crystal embryo pays surface tension σ on its skin and earns Δg per volume; their sum is a hill in size whose crest — the critical nucleus — has height ΔG*. Undercooling raises Δg, so the crest falls as 1/ΔT² and the waiting time, e+ΔG*/kT, collapses exponentially. A dust grain lowers the hill by a contact factor f.

  • Euler’s identitythe hidden room
    eiπ + 1 = 0

    Five constants, one line. The Square’s hidden room runs it as a contradiction: if π sat on some floor of the tower it would be algebraic, so iπ would be algebraic too, so Lindemann–Weierstrass would force eiπ to be transcendental — but eiπ = −1, an integer, ground floor. Something has to give, and it is π.

  • Euler’s formulaalways 2
    V − E + F = 2

    For any connected drawing in the plane with no crossings, vertices minus edges plus faces is 2 — however the drawing writhes. The surface keeps this number, not the drawer.

  • The plane’s capacityE ≤ 2V − 4
    4F ≤ 2E, V − E + F = 2 ⇒ E ≤ 2V − 4

    Houses never wire to houses, so no face can be a triangle — every face spends at least four edge-sides and each edge has only two to give. Six vertices buy eight edges; the ninth has no face to live in.

  • The handle’s ledgerχ = 0
    V − E + F = 2 − 2g ⇒ E ≤ 2V on the torus

    Every handle a surface grows lowers its Euler number by two and raises the wire budget by four. One handle: the cap moves from 8 to 12, and K3,3’s nine wires land flat with room to spare.

  • Euler’s walk conditionodd ∈ {0, 2}
    trail crosses every edge once ⇔ #odd = 0 or 2

    A walk using every bridge exactly once may hold odd degree only at its two ends — every land mass it merely passes through spends its bridge-ends in one-in, one-out pairs. Zero or two odd shores, never four.

  • The handshake lemmaΣdeg = 2E
    Σv deg(v) = 2E ⇒ #odd is even

    Every edge lifts the degree at both of its ends by one, so the degrees sum to exactly twice the edge count — and an even sum can never be built from an odd number of odd terms. Odd-degree shores always come in pairs, so a count of four can never fall to two by routing.

  • Inside the barrierevanescent
    ψ ∼ e−κx, κ = √(2(V₀ − E))

    Inside a wall too tall (E < V₀) the wave doesn’t oscillate — it decays. Not to zero, just smaller: thin enough, and something survives to the far side.

  • FitzHugh–Nagumofast v, slow w
    v̇ = D∇²v + v − v³/3 − w  ·  ẇ = ε(v + β − γw)

    A fast excitable variable v with a cubic nullcline, coupled to a slow recovery w. A small kick decays; a kick past threshold fires a full spike, then a refractory pause — the engine of every wave here. (Simulated via Barkley’s numerically-robust cousin.)

    1961
  • Excitable kineticsall-or-none
    |kick| < θ : decay  ·  |kick| > θ : spike + refractory

    One stable rest state and a threshold. Below it perturbations die; above it the response is all-or-none and identical however hard you pushed, followed by a recovery window in which it cannot fire again. Threshold plus refractoriness is the whole secret.

  • Front speedc ∝ √(D·r)
    c ≈ 2√(D r)  ·  λ = c T

    Diffusion D and reaction rate r set how fast an excitation front marches — faster coupling, faster wave. Multiply by the recovery period T and you get the wavelength: the spacing between crests in a train, or the arms of a spiral.

  • The eikonal relationcurved fronts slow
    c(κ) = c0 − D κ

    A curved wavefront travels slower than a flat one by an amount set by its curvature κ. A broken end lags, the front pivots about the gap, and the lag winds it into a spiral that turns forever about its tip — why spirals are inevitable, not accidental.

  • Adaptive expectationscorrect a fraction of the error
    pet = pet−1 + w (pt−1 − pet−1)

    The expected price moves a fraction w of the way to the last price. The cobweb then converges whenever |1 − w(1 + bs/bd)| < 1 — enough caution tames a market that naive expectations (w = 1) would throw apart.

  • Jain’s fairness index1/n ≤ J ≤ 1
    J = (Σ xi)2n Σ xi2

    One number for how evenly n users share: 1 when all get the same, 1/n when one gets everything, and k/n when k share equally and the rest starve. It ignores how much there is — utilisation is a separate count.

  • Max–min fairnesswater-filling
    xi = min(di, Cleft / nopen), repeated

    Raise every flow’s rate together; freeze a flow when its demand is met and share what is left among the rest. No flow can gain without taking from one that has less. It needs every demand and the capacity in one place — an oracle, not a protocol.

  • Fatou–Julia dichotomyconnected ⟺ c ∈ M
    Jc connected ⟺ c ∈ M

    Fix c and iterate every starting point instead and you get its Julia set. It is one connected piece exactly when c lies in M — otherwise it shatters into a Cantor dust. M is the index of them all.

  • Feature ratewhat you sample
    ffeat = N ftrue

    An N-spoke wheel looks identical every 1/N turn, so the rate that matters isn’t its spin — it’s how often it returns to an identical pose: N × the spin.

  • Feigenbaum constantuniversal
    δ = limn→∞ rn − rn−1rn+1 − rn ≈ 4.6692

    The doubling intervals shrink by a fixed ratio that converges to δ ≈ 4.6692 — and Feigenbaum found the same number for any smooth one-hump map. The route into chaos is universal; it doesn’t care what the map is.

    1978
  • The Cascade, liftedδ ≈ 4.6692
    c ∈ [−2, ¼] ↔ r: 1 → 2 → 4 → 8 → …

    On the real axis z² + c is the logistic map in disguise: the period-1 cardioid, the period-2 disc, then 4, 8, … bulbs accumulating at Feigenbaum — the same cascade The Cascade (INST·02) walks, here as a map.

  • The flash-lagangle → time
    Δt = θPSEω

    The flash must be placed ahead of the moving ring by θPSE to be seen aligned with it; divided by the ring’s angular speed ω, that offset is a time — how far behind the present the flash is placed by the moving thing’s percept.

  • Kawasaki’s theoremalternate angles sum to π
    α1 − α2 + α3 − ⋯ − α2n = 0

    A single vertex folds flat only if the sector angles between its creases, taken alternately, add to 180° each — walk round the folded vertex and you must come back to where you started.

  • Maekawa’s theoremM − V = ±2
    M − V = ±2

    At every vertex of a flat fold the mountain creases outnumber the valleys, or the valleys the mountains, by exactly two — the turns of the paper round the vertex must close.

  • The lift equationlevel flight
    W = ½ ρ V² S CL

    In steady level flight the wing’s lift carries the weight: air density, the square of airspeed, wing area and a lift coefficient of order one. Solve it for V and size enters only through W/S.

  • The neighbourhoodlocal only
    Ni = { j : |xj − xi| < r }

    Each agent reads only the neighbours inside a radius r — never the whole flock. Everything below is computed over that handful, and that locality is the entire point.

  • Order parameterphase transition
    φ = (1/N)·|Σj eiθj| ∈ [0, 1]

    One number for "how much one mind": 0 when headings are random, 1 when unanimous. Raise the noise and φ collapses past a critical ηc — order is a phase, not a dial.

  • Fourier seriesthe claim
    f(θ) = a₀2 + Σn [ aₙ cos nθ + bₙ sin nθ ]

    Fourier claimed something outrageous in 1822: any periodic signal — even one with corners and jumps — is a sum of smooth sine and cosine waves at integer multiples of one fundamental frequency. The harmonics slider adds these terms one by one.

  • Complex form · epicyclese^{iθ} = cosθ + i·sinθ
    f(t) = Σn cₙ ei n ω t

    Euler’s formula turns each harmonic into a rotating vector of radius |cₙ| spinning at speed n. Stacked tip-to-tip they are the epicycle chain on the left; the pen’s height is the signal.

  • The transform (DFT)read the recipe out
    Xk = Σn xn e−2πi·kn/N

    Run it backwards: the discrete Fourier transform reads the amplitude of each harmonic out of the samples — the spectrum bars. The same operation lives inside MP3, JPEG and every spectrum analyser.

  • Fourier transformall frequency, no time
    X(f) = ∫ x(t) e−2πi f t dt

    The whole signal folded onto the frequency axis: perfect pitch resolution, and zero sense of when. The wavelet trades a sliver of this for all of the missing time.

  • Pluck = sum of modesFourier
    an = (2 / n²π² p(1−p)) · sin(nπp)

    Any shape is a sum of pure modes — a Fourier sine series. A triangular pluck at position p has these exact weights, so a centre pluck (p=½) cancels every even harmonic.

  • Correlation dimensionD₂, read from data
    C(r) ~ rD₂

    Count how many pairs of reconstructed points fall within distance r; it grows as a power of r whose exponent is the attractor’s fractal dimension. The deep payoff: an invariant of a system you only ever sampled through one wire — kin to The Set and The Dendrite.

    1983
  • Fractal dimensionD ≈ 1.71
    N(R) ~ RD

    The number of particles within radius R scales as a power of R, and the exponent D ≈ 1.71 is fractional — denser than a curve, sparser than the plane. The cluster is mostly empty space, at every scale.

  • Fringe spacingΔy
    Δy = λLd

    How far apart the bright bands sit. Wider wavelength or closer slits ⇒ wider fringes — the two spatial knobs on the panel move exactly this.

  • The funicular — a cable’s equilibriumH y″ = w(x)
    H y″(x) = w(x)

    A cable carries its load w(x) per horizontal length in pure tension, its horizontal pull H the same everywhere. A uniform deck gives a parabola; the cable’s own weight, spread along its length, gives the catenary. Inverted, the same line carries the same load in pure compression: the arch.

  • Selectionsurvival ∝ fitness
    P(choose i) = fi / Σk fk

    The only place the goal ever enters: the fitter a genome, the more likely it breeds the next generation. Everything else is blind — selection alone carries the information about what “better” means.

  • The schema theoremHolland 1975
    m(H, t+1) ≥ m(H, t) · f(H)f̄ · [ 1 − ε ]

    Holland’s “fundamental theorem”: short, above-average building blocks (schemata) get exponentially more trials each generation. It is the reason blind variation plus selection actually climbs instead of wandering.

  • Crossoverrecombination
    child = parentA[0:k] ⌢ parentB[k:]

    Two surviving parents cut their gene strings at a random point and swap tails. Recombination is what lets evolution mix two good half-solutions into one whole — building-block by building-block, faster than mutation could alone.

  • Mutationblind variation
    gene ← random, with probability μ

    Each gene has a small chance μ of flipping to a fresh random value — the one source of genuinely new material. Too little and the population stagnates; too much and selection can’t hold its gains.

  • The derived seriessolvable ⇔ reaches 1
    G ▷ G′ ▷ G″ ▷ …  →  1 ?

    Each G′ is the commutator subgroup of the one before it — the swaps a radical built on the previous layer can undo. A polynomial is solvable by radicals exactly when this ladder reaches the trivial group in finitely many rungs. S₂, S₃ and S₄ all ground out; S₅ does not.

  • A₅ is perfectthe wall in one line
    [A5, A5] = A5

    A₅ — the even permutations of five roots, order 60 — is simple and equal to its own commutator subgroup. Nest the ladder as deep as you like: it never shrinks, so no radical formula can ever name the roots of a general quintic. This is the rung that refuses to break.

  • Thick-barrier limitWKB · Gamow
    T ≈ exp(−2 ∫ κ(x) dx)

    For a fat or soft wall, the leak is set by the area of forbidden region the wave must cross — Gamow’s factor, the same integral that dates a rock by its alpha decay.

  • Geometric Brownian MotiondS
    dS = μ S dt + σ S dW

    The price earns a steady drift μ and is kicked by Brownian noise σ, both scaled by the price itself — so returns, not absolute moves, are what’s random. The multiplicative twin of The Walk’s additive diffusion.

    1965
  • Mean vs. medianE[S_t]
    E[St] = S0 eμt > median

    The skew made precise: the mean rides above the median by eσ²t/2. A symmetric shock to returns is an asymmetric shock to price — upside unbounded, downside floored at zero.

  • Exact solutionS_t
    St = S0 e(μ − σ²/2)t + σWt

    Itô’s lemma integrates the SDE exactly: log-price is a drifting Gaussian, so the price is log-normal. This is the increment each sample path multiplies by — no discretisation error.

  • Gibbs phenomenon≈ 8.95%
    overshoot → G − 12 ≈ 8.95% · G = (2/π) ∫₀π sin tt dt ≈ 1.179

    At a jump, no finite sum of smooth waves ever closes it — a ripple overshoots the edge and, as you add harmonics, only narrows, never flattening. Its height homes in on a fixed ≈ 8.95% of the jump: the square hits it almost at once, the sawtooth climbs toward it from below.

  • The stretched exponentialKohlrausch, β < 1
    φ(t) = e−(t/τ)^β,  0 < β < 1

    How a glass forgets: not one exponential but a lazy shoulder on a log-time axis — the sum of many exponentials with no single τ. The smaller β, the wider the spread of relaxation times. It breaks, on purpose, the one-τ decay a single well gives.

  • The Vogel–Fulcher lawa divergence above zero
    τ = τ0eA/(T − T₀)

    The empirical fit to a cooling liquid’s relaxation time: an Arrhenius law whose effective barrier blows up at a temperature T₀ above absolute zero. Fragile liquids follow it steeply, strong ones stay nearly Arrhenius (Angell’s map).

  • The trap modelaging below T_g
    ρ(E) ∝ e−E/T_g,  τ(E) = τ0eE/T ⇒ ⟨τ⟩ = ∞ for T < Tg

    Cages hop among traps whose depths are exponentially distributed; escape from a trap is Arrhenius. Below Tg the mean trapping time diverges, so the system keeps finding deeper traps: its relaxation time grows with its own age, and equilibrium recedes as it is approached.

  • The global-game thresholdone θ*, no multiplicity
    ψ(θ*) = 1 − c  ⇒  θ* = 1 − c − m1 − m

    Each depositor sees θ through private noise and withdraws when the chance the bank fails beats the cost ratio c. At the marginal depositor’s signal the share of others withdrawing is, to her, uniformly uncertain — so at the threshold the run is exactly 1 − c, whatever the noise. With ψ(θ) = θ this is Morris and Shin’s θ* = 1 − c.

  • Public against privatetoo precise a headline restores the run
    α²β ≤ 2π

    α is the precision of a public signal, β of each private one. While the inequality holds the threshold is unique; beyond it the public signal becomes common enough to coordinate, and a precise enough headline can start a run on a bank the private evidence calls sound. With the bank’s resilience slope ψ′ the bound is α ≤ ψ′√(2πβ).

  • The pseudorangerange plus a clock
    ρi = |si − x| + c·b + εi

    The delay between a satellite’s time stamp and the receiver’s clock, times c: the true range, plus the receiver clock’s error b — the same for every satellite — plus noise. Four of them fix x, y, z and b; a nanosecond of b is thirty centimetres of every range.

  • Geometric dilution of precisionthe sky’s gain
    GDOP = √tr (HTH)−1

    H holds, for each satellite, minus the unit vector toward it and a 1 for the clock. Range noise σ becomes position-and-time error σ·GDOP: satellites spread over the sky keep it near 2; bunched together, the rows of H nearly agree, HᵀH nearly loses its inverse, and the error ellipse swells.

  • Regressional Goodhartselection returns a fraction
    E[V | U = u] = u · σV²σV² + σX²,   U = V + X

    A proxy U is the value V plus an independent error X, both normal. Pick by U and you get back only this fraction of the pick in V — the rest of the selection was error. Galton’s regression toward the mean, read as a law of metrics.

    2018
  • Over-optimisation, best of nd(α − βd)
    R(d) = d(α − βd),   d = √KL,   KLbon = ln n − n − 1n

    Gao, Schulman & Hilton’s fitted form for the gold reward against the distance d a best-of-n policy moves from the initial one (d is the square root of its KL divergence, which for best of n is known in closed form). It rises, peaks at d = α/2β, and falls.

  • The best of n≈ √(2 ln n)
    E[maxi≤n Zᵢ] = ∫ x · nφ(x)Φ(x)n−1 dx ≈ √(2 ln n)

    The expected largest of n standard normals, exactly by the order-statistic density, and its slow asymptote. For a Gaussian proxy the true value of the best of n is σV²/σU times it — it keeps rising, ever more slowly.

  • The Google matrixcolumn-stochastic
    G = d·M + 1−dN 𝟙

    M is the link matrix (a column per node, its 1s split evenly down its out-links). The teleport term — jump to a random node with probability 1−d — makes G positive and connected, so a single, well-defined rank is guaranteed to exist.

  • The rumour’s residues ≈ 0.203
    s = e−(k+1)(1−s)

    A teller who finds the listener already knew stops with chance 1/k (Demers et al.’s feedback-and-coin rule). In the continuous-time limit the fraction never told solves this equation: ≈ 0.203 at k = 1, Daley & Kendall’s number, and falling exponentially in k.

  • Rounds of push gossiplog₂ n + ln n
    T(n) = log₂ n + ln n + O(1)

    Each informed node tells one random node a round. The informed set doubles until most know (log₂ n rounds), then the last stragglers are found like coupons — each round catches a fixed fraction of the rest (another ln n).

    1987
  • Pull’s last phaseu → u²
    ut+1 = ut2 (pull) · ut+1 = ut2 e−(1−ut) (push–pull)

    An uninformed node asks a random node and stays ignorant only if it asked another ignorant one: the ignorant share squares every round. Push spends its last rounds looking for the few who have not heard; pull lets them come looking.

  • The governed enginetwo equations
    J ω̇ = k cos φ − F, m l φ̈ = m l (nω)2 sin φ cos φ − m g sin φ − b φ̇

    The flywheel speeds up when the steam torque k cos φ beats the load F; the balls swing out when the spindle’s centrifugal pull beats gravity, and the sleeve’s friction b resists them. The throttle closes as they rise. Linearised about the working speed, the pair becomes one cubic.

  • Gradient descentthe rule
    θ ← θ − η ∇θL

    The whole of learning. The gradient ∇L points uphill, so step the opposite way, scaled by the learning rate η. Repeat. There is no cleverer secret underneath modern AI than this line.

  • Weisskopf’s mountainh ∝ 1/g
    hmax ≈ ξ RyA mp g ∝ 1g ∝ M−1/3

    A mountain’s base flows when the weight of one molecule’s column, A mpgh, matches the energy that melts it — a small fraction ξ of a Rydberg. Weight ∝ h³ on a base ∝ h²: the square–cube law of a bone, written on a planet. On a rocky body g ∝ M1/3; where hmax reaches the body’s own radius, the potato radius, gravity stops shaping it.

  • Gray–Scottthis engine
    ∂A/∂t = DA∇²A − AB² + F(1−A) ; ∂B/∂t = DB∇²B + AB² − (F+k)B

    B autocatalyses through the AB² term — it eats A to make more of itself; feed F tops up A, kill k removes B.

    1983
  • Green’s theorem, area formregion ↔ boundary
    A = ½ ∮ (x dy − y dx)

    The area inside a closed curve is a sum taken only along the curve. Walked counter-clockwise it is positive, clockwise negative; a figure-eight returns the difference of its lobes. On a polygon the sum is the shoelace formula.

  • The planimeter readingA = b·ℓ
    A = b · ℓ, dℓ = dW · n̂

    The wheel on the tracer arm rolls only with the part of its motion perpendicular to the arm. Summed round a closed trace, that sideways roll ℓ times the arm length b is the area the tracer enclosed — so long as the pole arm swung out and came back without going round.

  • The zero circlepole inside
    A = b · ℓ + π (a2 + b2 − 2bc)

    With the pole inside the outline both arms make a full turn, and the wheel misses the area of one special circle — the one the tracer could ride forever without the wheel turning. Here a is the pole arm, b the tracer arm and c the wheel’s place along the tracer arm from the hinge; a wheel set behind the hinge makes c negative.

  • Grover's iterationamplify
    G = (2|s⟩⟨s| − I) · Ow

    Two reflections per round: the oracle Ow flips the sign of the marked answer, then the diffusion operator reflects every amplitude about their mean — and the marked one grows while the rest shrink.

  • A rotation toward the answera rotation
    |ψk⟩ = sin((2k+1)θ)|w⟩ + cos((2k+1)θ)|s′⟩, sinθ = 1√N

    The whole algorithm is one rotation in the plane spanned by the answer |w⟩ and everything-else |s′⟩: each step turns the state by 2θ toward |w⟩. Run it too long and it sails past the answer.

  • The optimal stopping pointthe speedup
    k* ≈ π4√N

    About π4√N rounds land the amplitude almost entirely on the answer — a quadratic speedup over the N/2 a classical scan averages, and provably the best any quantum search can do.

  • H-theoremdH/dt ≤ 0
    H = ∫ f ln f d²v

    The velocity distribution relaxes to Maxwell–Boltzmann (the line in the histogram), and H falls monotonically: the H-theorem. The temperature never changes — only the shape does.

  • Hamming(7,4)dₘᵢₙ = 3
    4 data + 3 parity → 7 · corrects 1

    Codewords sit at least three flips apart, so any single error still lands nearest its original — sphere-packing in bit-space. A 3-bit syndrome names the flipped position outright. Hamming, 1950.

  • The Tower of Hanoi countexact
    M(n) = 2·M(n − 1) + 1 = 2n − 1

    To move n disks, move n − 1 out of the way, move the largest once, and move the n − 1 back on top: the count doubles and adds one, so it is 2ⁿ − 1 — and no solution is shorter. Each disk doubles the wait; sixty-four of them at one move a second outlast the stars.

  • The Frame–Stewart recursionfour pegs
    FS(n) = mink [2·FS(k) + 2n−k − 1]

    With a fourth peg, park the top k disks using all four, move the other n − k with three, and bring the k back: choose k to minimise. The wait stops doubling per disk — 64 disks fall from 1.8×10¹⁹ moves to 18,433. Optimal for four pegs (Bousch 2014); open beyond.

  • Laplacian growth∇²u = 0
    ∇²u = 0,  v ∝ ∂u/∂n

    The walker-arrival probability obeys Laplace’s equation — the same discrete Laplacian The Skin diffuses with — and the cluster grows fastest where that field is steepest. Drift biases the walkers and the tree leans into the wind.

  • Keys moved by a joinone arc
    ring: E[moved] = KS + 1 · mod N: ≈ K · SS + 1

    A new server on the ring takes only the keys between its predecessor and itself — on average one share in S + 1, the least any map could move. Under key mod N a key stays only if its remainder is the same for S and S + 1, about one key in S + 1; all the rest move.

    1997
  • Virtual nodes even the loadmax ÷ mean
    max/mean ≈ E[max of S Gamma(v, 1/v)] → ln S at v = 1

    The S·v positions are uniform points on the circle, so the arcs between them are nearly independent exponentials and a server’s share is a sum of v of them. With one position each, the busiest holds about ln S times the mean; with v, the spread shrinks like 1/√v.

  • Chord’s lookup length½ log₂ S
    finger i = succ(n + 2i) · E[hops] ≈ ½ log₂ S

    Each node keeps a finger to the node succeeding it at every power of two around the ring. Forwarding a lookup to the closest finger before the key at least halves the remaining distance, so a lookup takes at most log₂ S hops and about half that on average.

  • Borel’s family-size lawmean 1/(1 − n)
    P(s) = e−ns(ns)s−1s!,   E s = 11 − n

    The size of one immigrant’s family when every member has Poisson(n) children. Below n = 1 its tail is cut off exponentially; at n = 1 the cut-off vanishes and P(s) falls only as s−3/2.

  • The Hawkes intensityevery event raises the odds
    λ(t) = μ + Σtᵢ<t n·β e−β(t − tᵢ)

    The rate of events right now: an exogenous base μ, plus a decaying kick from every past event. Each kick integrates to n, the branching ratio — the mean number of events one event sets off — and fades on the time scale 1/β.

  • The stationary rateμ/(1 − n)
    Λ = μ1 − n,   n < 1

    Each immigrant brings 1 + n + n² + … events in all, so the long-run rate is the exogenous rate divided by 1 − n, and the endogenous share of activity is n. At n = 1 the series diverges: there is no finite rate, and no amount of news is needed to keep the market busy.

  • Diffusion / heat equation∂ρ/∂t = D ∂²ρ/∂x²
    ρ = (4πDt)−½ e−x²/4Dt

    The walker density obeys the heat equation; its point-source solution is the spreading Gaussian — the silver curve, of width σ = √(2Dt). One blind walker becomes a smooth, deterministic law.

  • The nodal setwhere sand settles
    { (x, y) : φ(x, y) = 0 }

    The curve where a mode never moves. Sand bounced off the violently vibrating antinodes drifts until it lands on these still lines and stays — so the silence, not the sound, is what you see. Chladni’s 1787 figures are these zero sets made visible.

  • Weyl's lawhear the area
    N(k) ≈ A4π k² − L4π k

    Count the modes below frequency k: they grow as the drum’s area times k², with a correction in its perimeter. So the spectrum betrays the area and the edge exactly — you can hear how big the drum is. Just not, Kac asked, its shape; in 1992 the answer came back: no.

  • The membrane eigenproblemfixed-edge modes
    ∇²φ + k²φ = 0,  φ|∂Ω = 0

    A drumhead’s standing waves are the eigenfunctions of the Laplacian on its domain Ω, pinned to zero at the rim. Each solution φ is a pure mode shape; its eigenvalue k² fixes the tone. The shape of Ω chooses the whole ladder of modes — that is the question.

  • Rectangular & disk modessin·sin · Bessel
    φmn = sinmπxa sinnπyb

    On a rectangle the modes are products of sines, frequencies following √((m/a)² + (n/b)²); on a disk they become Bessel functions Jl with circular and radial nodes. The sand settles wherever φ = 0.

  • Herd immunityp_c
    pc = 1 − 1 / R₀

    Immunise a fraction p and you thin the network of usable paths. Past p_c the infected component can no longer span the graph and the outbreak stops — even though a majority remain susceptible.

  • Use-inheritanceobsolete
    gchild = ½(ga + h·Δa + gb + h·Δb)

    Lamarck’s engine, made quantitative: the young inherit a share h of what their parents acquired, Δ. Under a standing use the trait climbs every generation with no selection at all.

  • The Weismann barriersoma ↛ germ
    gchild = ½(ga + gb) + ε · Δ ↛ g

    The germ line is set apart from the body: nothing acquired in a life (Δ) writes to what is passed on. Dock tails for five generations and every pup is born whole.

  • Hit probabilityP(S_T>K)
    Φ( [ln(S0/K) + (μ − σ²/2)T] / σ√T )

    A decision, in closed form: the chance of finishing above a target K is Φ(d₂). The fraction of sample endpoints above the line converges to it — the histogram checking the formula.

  • The Viterbi paththe best path
    δt(j) = maxi [ δt−1(i) aij ] bj(xt)

    Replace the forward sum with a max and keep a back-pointer: δ carries the score of the best path reaching each state, and the trace-back yields the single most-likely sequence of regimes — one global decode, not a string of per-step guesses.

  • Hidden Markov modelstates + emissions
    P(x, s) = πs₁ ∏t as→s′ bs(xt)

    A hidden chain of states you never see — transitions a — each state emitting an observation x through its own distribution b. Three things define it: the start π, the transition matrix, the emissions. Everything else is inference.

  • Forward–backwardsoft posterior
    γt(i) = αt(i) βt(i)P(x)

    Sweep forward accumulating α (the evidence up to t) and backward accumulating β (the evidence after it); their product, normalized, is the probability of each hidden state at every step — the soft regime bands painted under the price.

  • Baum–Welch (EM)learns itself
    θ ← argmaxθ E[ log P(x, s) | x, θ ]

    When the transitions and emissions are unknown, Expectation–Maximization re-estimates them from the very posteriors forward–backward computes, alternating until they stop moving. The chain learns its own parameters — the unsupervised cousin of The Lens’s given model.

  • Hodge projection∇·u = 0
    u = u* − ∇p  ·  ∇²p = ∇·u*

    Any vector field splits (Helmholtz) into a divergence-free flow plus the gradient of a scalar. We solve a Poisson equation ∇²p = ∇·u* for that scalar and subtract its gradient, leaving a mass-conserving field. This is the projection step every frame.

  • Recording a hologramthe twin is in it
    I = |R + O|² = |R|² + |O|² + R*O + RO*

    A plate records intensity only, but the cross term R*O carries the object wave’s phase, written into the spacing of the fringes. Its conjugate RO* — the twin — is recorded too. In-line the two share every spatial frequency; tilt the reference and they part.

  • Angular-spectrum propagationexact, by FFT
    U(z) = F−1[F[U(0)] · ei2πz√(1/λ² − f²)]

    A field is a sum of plane waves — its 2-D spectrum — and each plane wave advances by a known phase. Multiply, transform back, and the field is anywhere you like: run z negative and a plate refocuses to any depth of the scene that lit it.

  • Fresnel diffractionparaxial
    U(x, z) ∝ eikziλz ∫ U(ξ, 0) eik(x − ξ)²/2z dξ

    The near-axis limit of the same propagation: every point of the plate sends a spherical wavelet, approximated as a parabola. It is why one point’s hologram is a zone plate — rings that crowd outward — and why a small piece of plate still holds the point.

  • The linking number= 1
    lk = 14π ∮∮ (r₁ − r₂) · (dr₁ × dr₂)|r₁ − r₂|³ = 1

    Pick any two distinct fibres: Gauss’s integral, which counts how many times one loop threads the other, always returns exactly one. No two are unlinked — the whole weave is chain-mail with no free rings, kin to the trefoil of The Knot.

  • The Hopf mapS³ → S²
    (a,b,c,d) ↦ ( 2(ac+bd), 2(bc−ad), a²+b²−c²−d² )

    Every point of the 3-sphere (a,b,c,d) is sent to a point of the ordinary 2-sphere. The map forgets one circle’s worth of information — and the whole circle it forgets lands on the same spot below. That circle is a fibre.

  • A fibrea circle
    z1 = cosθ2 ei t , z2 = sinθ2 ei(t+φ)

    Fix a target (θ, φ) on the base sphere and let the shared phase t run once around: the pre-image is a full circle in the 3-sphere. Slide (θ, φ) and the circle sweeps out a torus — the whole weave is these circles, one per direction.

  • Stereographic projectionS³ → ℝ³
    (a,b,c,d) ↦ (a,b,c) / (1 − d)

    The 3-sphere lives in 4-D, so to see it we cast it into ordinary space from one pole — the same trick that flattens a globe to a map. It sends circles to circles, so every fibre stays a clean ring you can orbit.

  • Storing a memoryHebb
    wij = 1N Σμ ξiμ ξjμ

    Carve each pattern ξμ into the wiring: neurons that agree in a memory get a positive link, those that disagree a negative one. Memory is a property of the connections, never one address.

  • The energyLyapunov
    E = −12 Σi,j wij si sj

    The Ising energy of The Threshold, but with learned couplings wij. Every neuron flip can only lower it, so the state slides downhill and must stop — in a memory.

  • How much it holdsα_c ≈ 0.138
    pmax ≈ 0.138 N

    Store more than ≈0.138N patterns into N neurons and the basins overlap and collapse into spurious blends — the net forgets everything at once. The spin-glass price of a distributed memory.

    1985
  • Recalldescent
    si ← sgn(Σj wij sj)

    Show a corrupted cue, then let each neuron align with the majority vote of the rest. The state rolls down E into the nearest stored pattern — content-addressable memory: the cue is the address.

  • First incompletenessa true gap
    Con(T) ∧ arithmetic ⇒ ∃ G : (T ⊬ G) ∧ (T ⊬ ¬G)

    Any consistent system strong enough to do arithmetic leaves sentences it can neither prove nor refute — and one of them, G, is true. Add it as an axiom and the new system grows its own G′. The shoreline moves; the sea never empties.

  • Second incompletenessno self-audit
    T ⊬ Con(T)

    A consistent system cannot prove its own consistency — Con(T) is itself one of its dark stars. The auditor cannot audit herself. Gentzen (1936) proves it anyway, but only from a stronger vantage: induction up to ε₀.

  • The cascade conditionthe record outweighs one signal
    a = sign(λ + s·ℓ),   ℓ = ln p1 − p,   cascade ⇔ |λ| > ℓ

    λ is the public record’s log-odds, s = ±1 the private signal. Once the record leans by more than one signal’s worth, every arrival takes the same door whatever she saw; her action then tells the record nothing and λ stops moving. With binary signals that happens after the first two agreeing choices (a tie is a coin).

    1992
  • The odds of a right cascadeBHW, eq. 3
    P(right) = p(1 + p)2(1 − p + p²),   P(wrong) = (1 − p)(2 − p)2(1 − p + p²)

    The chance the crowd ends up locked on the better door, and on the worse one, for signals right with probability p. At p = 0.6 one crowd in three cascades the wrong way — with a hundred people still to come, each holding a signal that would have helped.

    1992
  • The van Cittert–Zernike theorema visibility is a Fourier coefficient
    V(u, v) = ∬ I(l, m) e−2πi(ul + vm) dl dm

    The correlation of the signals at two points — two dishes — is the Fourier transform of the source’s brightness, read at their separation in wavelengths. An interferometer never sees the sky; it hears one coefficient of its transform per baseline, per moment.

  • Högbom CLEANsubtract the beam
    Ires ← Ires − γ · max · B(x − xmax)

    The dirty image is the sky convolved with the dirty beam B. Find its brightest point, subtract a fraction γ of the beam there, record a delta-function component; repeat. Then convolve the components with a clean Gaussian beam and add back the residual. It is a prior — the sky is a sum of points — and it fills the unsampled uv plane with it.

  • Angular resolutionθ ≈ λ/B
    θ ≈ λB

    The finest detail a telescope separates is about the wavelength over its size — or, for an array, over its longest baseline. At 1.3 mm across the Earth that is some 25 micro-arcseconds: a grapefruit on the Moon. The shortest baselines set the largest scale it can see at all.

  • Inversion in a circleexact
    OP · OQ = ℓ2 − s2 = k2

    In the cell, P and Q stay on one ray from the fixed pivot O and the product of their distances never changes: Q is P inverted in the circle of radius k. Inversion sends every circle through O to a straight line — so drive P round a circle through O, and Q walks a line, exactly.

  • Critical pointexact
    kTc/J = 2ln(1 + √2) ≈ 2.269

    No fitting, no simulation — Onsager solved it on paper. Above it, noise wins; below it, order.

  • Lattice energyHamiltonian
    H = −J Σ⟨i,j⟩ σiσj − h Σi σi

    Each spin σ = ±1 only wants to match its four neighbours (coupling J > 0); the field h tilts the whole sheet one way.

  • Thermal samplingMetropolis
    P(σ) ∝ e−H/kT ; accept w.p. min(1, e−ΔE/kT), ΔE = 2σi(J Σj σj + h)

    Temperature is the only referee: cheap moves are always taken, costly ones only a fraction e−ΔE/kT of the time. Run it and the sheet is a draw from the thermal ensemble.

  • Order parameterβ = 1/8
    m = [1 − sinh−4(2J/kT)]1/8 ∼ (Tc − T)1/8

    The magnetisation leaves zero with a vertical tangent at Tc — that exact curve is what the inset dot rides.

  • Diverging scalesuniversality
    ξ ∼ |T − Tc|−ν, χ ∼ |T − Tc|−γ ; ν = 1, γ = 7/4

    At Tc the correlation length blows up and the lattice looks the same at every zoom. The kicker: these exponents are universal — the very same numbers describe water’s critical point and a ferromagnet.

  • The jury theoremindependent votes
    PN(majority right) = Σk>N/2 C(N, k) pk(1 − p)N−k → 1

    With every signal acted on and counted, a majority of N independent voters, each right with probability p > ½, is right with a probability that tends to one. The sealed street is a Condorcet jury; the published one is not.

  • Kalman gainK
    K = P⁻ Hᵀ S⁻¹

    The heart of it: the gain is the prediction’s share of the total uncertainty. Confident prediction (small P⁻) or noisy sensor (large R) → K small, trust the model. Loose model (large Q) → K → 1, trust the blip.

  • Predictx⁻, P⁻
    x⁻ = F x  ·  P⁻ = F P Fᵀ + Q

    The model pushes the belief forward — here a constant-velocity guess — and inflates the covariance by the process noise Q. With no correction this step alone is a forecast: the ellipse only grows.

  • Innovationy, S
    y = z − H x⁻  ·  S = H P⁻ Hᵀ + R

    How wrong the prediction was, and how much that surprise could be noise. S is the innovation’s covariance — the prediction’s uncertainty plus the sensor’s.

  • Correctx, P
    x = x⁻ + K y  ·  P = (I − K H) P⁻

    The new estimate is the prediction nudged toward the blip by K of the innovation, and the covariance shrinks — the ellipse closing. Information only ever sharpens the belief.

  • Blinded → forecastP↑
    predictn, no correct ⇒ P → ∞

    Run predict without correct, over and over, and P grows without bound: the covariance ellipse becomes a forecast cone. The Lens, blinded, is The Oracle.

  • Fractal dimensionKaplan–Yorke
    dKY = 2 + λ₁|λ₃| ≈ 2.06

    The attractor has a fractal (Hausdorff) dimension of ≈ 2.06 — more than a surface but less than a volume. Trajectories are confined to this set forever, yet are locally unstable in every direction. This is what "strange" means: attracting yet chaotic.

  • The exact solitonone pulse
    u(x, t) = c2 sech²(√(c/4)·(x − ct))

    The exact single-soliton solution to KdV. Amplitude = c/2, speed = c — taller solitons travel faster. The width narrows as amplitude rises, so a soliton’s energy is concentrated into ever-tighter pulses. When two solitons meet, they pass through each other, emerging with the same shape — only their positions shift (the phase shift).

    1895
  • The phase shiftelastic collision
    Δx = 1√c1 ln √c1 + √c2|√c1 − √c2|

    After two KdV solitons collide and separate, each soliton is displaced from where it would have been without the collision. This phase shift is the only trace of the interaction — their shapes are entirely preserved. Zabusky and Kruskal named these objects “solitons” in 1965, because they behave like particles.

    1965
  • Korteweg–de Vries equationPDE
    ∂u/∂t + 6u ∂u/∂x + ∂³u/∂x³ = 0

    The KdV equation governs shallow water waves, plasma waves, and ion-acoustic waves. The nonlinear term 6u∂u/∂x steepens tall waves (making them faster); the dispersive term ∂³u/∂x³ spreads them out. Their balance produces the soliton: a stable, shape-preserving pulse that travels at constant speed forever.

    1895
  • Maximise expected log-wealthgrowth-optimal
    maxf 𝔼[ ln Wn ] ⇔ max long-run growth

    Bernoulli’s idea, made a law by Breiman: value log-wealth, not wealth. Maximising expected log-wealth is exactly maximising the geometric growth rate — and it almost surely beats every other strategy over the long run, while a bigger expected dollar return can march the median to zero.

  • The Kelly fractionoptimal stake
    f* = p − qb = bp − qb

    Bet the fraction f* of your bankroll set by your edge over the odds — win probability p, loss q = 1 − p, payoff b to 1. No edge, no bet; a big edge at short odds, bet big. This one fraction maximises long-run growth.

  • The growth ratelog-growth
    g(f) = p ln(1 + bf) + q ln(1 − f)

    The exponential rate your wealth compounds at when you stake fraction f. It is a hill: zero at f = 0, peaking at f*, and back to zero at the no-growth point f = 2f* — past which a favourable game loses money to volatility drag.

  • Kelly for an assetcontinuous
    f* = μ − rσ²

    The same criterion for a lognormal asset (The Oracle’s GBM): stake the excess return μ − r over the variance σ². Thorp carried this from blackjack to the market — the growth-optimal portfolio is the Sharpe ratio, squared and scaled.

  • Kempe’s expansiona curve as cosines
    f(r cos α + r cos β, r sin α + r sin β) = C + Σ ck cos(mkα + nkβ + γk)

    Put a two-bar arm’s pen into the curve’s polynomial and products of cosines become sums of cosines. Each term is a rod at angle mα + nβ + γ, a point at radius c on it; add the points and the sum’s x-coordinate is f − C. Hold it on x = −C and the pen is on the curve.

  • Kempe’s reversorθ_E = 2θ_C − θ_A
    OABC ∼ OCDE  ⇒  ∠AOC = ∠COE,   θE = 2θC − θA

    Two similar contraparallelograms share a rod, the first’s corner pinned to the second’s long bar: the third rod is the mirror of the first in the second. Chained from a fixed rod it multiplies an angle; two of them through a bisecting rod add two.

  • Kepler’s ellipseone curve, one law
    r = a(1 − e²)1 + e cos ν · dA/dt = const

    The orbit is an ellipse with the Sun at a focus, and the line to the planet sweeps equal areas in equal times. Two curves fit the sky that a chain of circles only approaches — and the area law is a reason, the one Newton would derive.

  • The Jones polynomialskein relation
    t−1V(L+) − t V(L−) = (t½ − t−½) V(L0)

    A polynomial fingerprint built by switching one crossing at a time. It tells the trefoil from its mirror image — something no amount of wiggling can do by eye — and it fell out of the algebra of quantum observables, not geometry.

  • WritheGauss integral
    Wr = 14π ∮∮ (dr × dr′)·(r − r′)|r − r′|³

    A real number that measures how much the loop coils around itself, summed over every pair of points along the curve — the same double integral Gauss wrote for linking. Averaged over all viewing angles it is the signed crossing count you actually see.

  • Möbius energyself-repelling
    E(γ) = ∮∮ (1|r − r′|² − 1d(γ)²)

    Treat the curve as charged and let it push itself apart. Flowing downhill in this energy untangles every loop to its tidiest form — a plain circle for the unknot — yet it can never force a real knot open. That obstinacy is the invariant.

  • Kramers escape ratethermally activated
    rK ∝ exp(−ΔU/D)

    The rate at which thermal fluctuations knock a particle over a potential barrier ΔU, with D = σ²/2 the noise diffusion. In the threshold model of stochastic resonance, this sets how often the noisy signal crosses the threshold — and therefore the output firing rate that tracks the input signal.

  • Bifurcationr(K)
    r = √(1 − Kc/K)

    Above threshold the coherence rises continuously from zero — a second-order phase transition. Exact for a Lorentzian spread; the silver curve the live point tracks.

  • Critical couplingK_c
    Kc = 2 / (π g(0))

    The ignition point. Widen the frequency spread g and you raise the bar; sharpen it and the chorus locks for almost any nudge. The blue line on the bifurcation inset.

  • Locked vs drifting|ω| ≤ Kr
    locked ⟺ |ωi − Ω| ≤ K r

    Only oscillators slow enough to be captured join the cluster; |ωᵢ| > Kr keep drifting around the rim. The locked fraction is on the readout — and a bimodal spread can refuse to pick one rhythm at all.

  • Mean-field formexact
    θ̇i = ωi + K r sin(ψ − θi)

    Because the coupling is global, every oscillator feels a single pull toward the mean phase ψ, of strength Kr. The pull grows with the order it has already created — the feedback that makes the transition sharp.

  • Kuramoto modelθ̇ᵢ
    θ̇i = ωi + (K/N) Σj sin(θj − θi)

    Each oscillator runs at its own frequency ωᵢ and is pulled toward every other by the coupling K. With global coupling this is the whole model — Kuramoto’s solvable distillation of Winfree.

  • Order parameterr·e^{iψ}
    r eiψ = (1/N) Σj eiθⱼ

    The population’s mean phasor. r measures coherence: 0 when the phases are scattered, 1 when they march in lockstep. It is the arrow drawn on the circle and the number on the readout.

  • Kuratowski’s theoremtwo culprits
    planar ⇔ no K5, no K3,3

    Every graph that refuses to lie flat hides one of exactly two minimal culprits — the five-clique or the utilities board. This instrument is not an example of the wall; it is half of the wall’s entire statement.

  • The flashing ratchetrectified noise
    ẋ = −s(t) · V′(x) − F + √(2D) ξ(t)

    Overdamped particles in a sawtooth potential that switches on and off, s(t) ∈ {0, 1}. On: they slide into the teeth. Off: they spread as pure diffusion. On again: the asymmetry harvests the spread — net transport uphill against the load F, paid for entirely by the switching.

  • Snell&rsquo;s lawrefraction
    n1 sin θ1 = n2 sin θ2

    The bend at an interface — not a separate law but the price of least time: minimise T over the crossing point and the derivative vanishes exactly when n sin θ is conserved across the boundary. Past the critical angle there is no real θ2 and the light reflects entirely.

  • Fermat&rsquo;s principleleast time
    δ T = δ 1c ∫ n ds = 0

    Of all the paths light could take between two points, it takes the one whose travel time is stationary — almost always the least. The speed in a medium is c/n, so the ray bends to spend less length where it would crawl. Every law of ray optics falls out of this one line.

  • The brachistochronea cycloid
    x = r(φ − sin φ)   y = r(1 − cos φ)

    The fastest descent is a cycloid — the curve a chalk mark on a rolling wheel traces, run upside down. Bernoulli found it by a trick: a bead falling through ever-faster layers is a light ray refracting through ever-thinner glass, so sin θ / v = const — Snell’s law — gives the curve. Least time, twice.

    1696
  • The Euler–Lagrange equationδS = 0 ⇒ this
    ddt∂L∂q̇ − ∂L∂q = 0

    Demand that the action be stationary against every small wiggle of the path and this is what survives — the equation the extremal curve must obey. Feed it L = T − V and out comes Newton’s F = ma; feed it the descent-time integrand and out comes the cycloid. One machine, every path.

  • Stationary actionthe meta-law
    S = ∫t₁t₂ L dt   δS = 0

    Of all the ways a system could move between two fixed events, it takes the one whose action is stationary. This single sentence contains Newton, optics, relativity and — through Feynman — quantum mechanics. The rack’s deepest answer to what is the law under the laws?

    1834
  • The path integralFeynman 1948
    K = ∫ 𝒟q(t) · eiS/ℏ

    Why does the classical path win? A quantum particle takes every path, each carrying a phase e^(iS/ℏ). Away from the stationary path the phases churn and cancel; near it they march in step and add. “Least action” is the place where neighbouring histories stop interfering — classical mechanics as the loud silence of a quantum sum.

  • Optical path lengthwhat is minimised
    OPL = ∫ n ds = c · T

    Geometric length weighted by the index — the distance light “feels”. Fermat’s path is the one of stationary OPL; a perfect lens works by making every route to the focus take equal optical path, so they arrive in phase.

  • The descent-time functionalminimise this
    T = ∫ dsv = ∫ 1 + y′22gy1/2 dx

    The time for a bead to slide a curve y(x): arc length over speed, with v = √(2gy) from energy conservation. The brachistochrone is the y(x) that makes this integral least — a problem about a whole function, not a number, which is what the calculus of variations was invented to solve.

  • The Lagrangiankinetic − potential
    L = T − V

    Not energy (that is T + V) but its mirror — kinetic minus potential. Nature’s ledger: a path is cheap when it banks potential early and spends kinetic sparingly. Integrate L over the trip and you get the action, the thing the true motion makes least.

  • The updatecontinuous
    At+Δt = At + Δt · G(K ∗ At) , clipped to [0, 1]

    Every cell reads its neighbourhood through a kernel K, runs the total through a growth map G, and edges up or down by a small Δt. Take Δt → 1 and the kernel to Life’s 3×3 count and this is The Garden — Lenia just lets the grid, the clock and the state all go smooth.

  • The growth mapbirth ⇄ death
    G(u) = 2 e−(u−μ)² / 2σ² − 1

    A bell curve in [−1, 1]: a neighbourhood total near the sweet spot μ makes a cell grow, anything too sparse or too crowded makes it decay. Life’s hard “born on 3, survive on 2–3” became a smooth dial — and the width σ is how forgiving it is.

  • The potentialK ∗ A
    U(x) = Σy K(|y|) · A(x + y) , Σ K = 1

    The kernel is a soft ring — each cell weighs an annulus of neighbours, not just the eight touching it. That single ring, normalised to sum to one, is the whole difference between a pixel grid and a medium an organism can swim through.

  • Lens equationEinstein ring
    β = θ − θE2/θ

    Solve for the apparent position θ given the true source β. When the source sits dead behind the mass, every solution is a circle — a full ring of radius θE. (This is the line the shader inverts per pixel.)

  • When the bell failsLévy α-stable
    Var = ∞ ⇒ no 𝒩

    Remove the finite-variance condition (Cauchy / heavy tails) and the classical CLT is void: the sum stays a fat-tailed Lévy α-stable spike, no Gaussian ever forms — the wildness real markets live in (Mandelbrot vs Bachelier).

  • The ruleB3/S23
    born if n = 3 · survives if n ∈ {2, 3} · else dies

    One law for every cell: a dead cell with exactly three live neighbours is born, a live cell with two or three survives, anything else is death by loneliness or crowding. No parameters at all.

  • What a cell seestotalistic
    ni = Σj ∈ Moore(i) sj , s ∈ {0, 1}

    Each cell counts only the eight around it (the Moore neighbourhood) and cares only about the total — not which neighbours are alive. A local, totalistic rule, exactly like The Rule one dimension down.

  • Light deflectionmeasured 1919
    α = 4GMc2b

    A ray grazing a mass bends by this angle — exactly twice the Newtonian guess. Eddington’s eclipse confirmed it and made Einstein famous overnight.

    1919
  • The Thouless ratio0 = the screen · 1 = the paint
    TR = log M − log Slog R − log S

    M is the visitor’s match, S what a photometer would set (the luminance), R what full constancy would set (the paint). Zero is a luminance meter, one is perfect albedo recovery; the log form treats equal ratios as equal steps.

  • The anchoring rulethe brightest is white
    log ℓ = w · log0.9 LLmax, local + (1 − w) · log0.9 LLmax, global

    Within each framework of illumination the highest luminance is anchored to white (reflectance 0.9) and everything else scales from it; the reported lightness blends the local framework (the shadow) with the global one (the whole scene).

  • Light is paint times illuminationL = R·E
    L(x) = R(x) · E(x)

    What reaches the eye from a matte surface is its reflectance (the paint, R) times the light falling on it (E). One number, two causes: a dark paint in sun and a light paint in shade can send the same L, and the eye has to decide which product it is looking at.

  • Lindemann–Weierstrassπ is transcendental
    α algebraic, α ≠ 0  ⇒  eα transcendental

    For any nonzero algebraic number α, eα is transcendental — never the root of any polynomial with rational coefficients. Setting α = iπ forces π itself off every floor of the constructible tower: degree infinity, and no compass-and-straightedge square can ever equal the circle.

  • Gravitational growthδ ∝ D(t)
    δ(x, t) = D(t) δ0(x)

    A denser-than-average region pulls in more matter and grows denser still — runaway under its own gravity, slowed only by the cosmic expansion. While the contrast is small every patch grows by the same factor D(t): the rich get richer.

    1980
  • The four-bar looploop closure
    r1eiθ₁ + r2eiθ₂ = r0 + r3eiθ₃

    Walk round a four-bar linkage and you come home: crank and coupler, as complex vectors, add up to ground and rocker. Give the crank angle and the loop fixes the other two — one circle–circle intersection. Watt’s pen is the coupler’s midpoint; at his symmetric proportions its error grows only as the fifth power of the stroke.

  • Conserved, yet wildLiouville
    E = T + V = const ; dΓ/dt = 0

    No friction: energy is conserved and phase-space volume is preserved. The system is deterministic and reversible — and still unpredictable. Determinism ≠ predictability.

  • Happened-beforeLamport 1978
    a → b ⇔ same process and a first · or a sends what b receives · or a → c → b

    The only order a distributed system can know: a before b if a could have caused b — earlier on one machine, or linked by a chain of messages. Two events with no such chain either way are concurrent, a ∥ b, and no observer-independent order exists between them.

  • The Lamport clocka → b ⇒ C(a) < C(b)
    C ← C + 1 per event · on receive: C ← max(C, stamp) + 1

    One integer per process, ticked on every event and pushed past every stamp it receives. It never contradicts causality — but the converse fails: C(a) < C(b) does not mean a could have caused b.

  • The vector clocka → b ⇔ V(a) < V(b)
    V(a) < V(b) ⇔ ∀k V(a)k ≤ V(b)k and V(a) ≠ V(b)

    One counter per process instead of one in all. Each process ticks its own entry and takes the component-wise maximum of every stamp it receives; then a happened before b exactly when every entry of V(a) is at most V(b)’s — causality, captured both ways.

  • The logistic mapone line
    xn+1 = r · xn (1 − xn)

    Robert May’s 1976 toy ecology: a population is the previous one grown by r and damped by its own crowding (1 − x). One quadratic, iterated — and it breeds the full route to chaos.

  • Lorentz transformationboost
    t′ = γ(t − βx), x′ = γ(x − βt)

    The map between frames moving at relative speed β (units of c). It mixes space into time and time into space — which is why a boost shears the diagram instead of just sliding it.

    1905
  • The invariance of cpostulate
    c′ = c in every inertial frame

    Einstein’s second postulate: light’s speed is the same for every observer in uniform motion. No medium, no wind, no fringe to slide — the null the interferometer kept returning, made into a law.

  • Invariant intervals²
    s² = −(ct)² + x² = −(ct′)² + x′²

    The one number every observer measures the same. It is the metric of flat spacetime, and the reason the 45° cone is absolute while the axes are not.

    1908
  • Lorentz factorγ
    γ = 1 / √(1 − β²)

    The Lorentz factor. At β = 0 it is 1; as β → c it runs to infinity — moving clocks slow and moving rulers shrink without bound. The ember dot on the inset rides this curve.

  • Relativity of simultaneityΔt′
    Δt′ = γ(Δt − β Δx)

    Two events at the same time in one frame are not at the same time in another. "Now" is a frame-dependent slice — the root of every relativity paradox, drawn live in the Simultaneity scene.

  • Proper timetwin paradox
    τ = ∫ √(1 − β²) dt = t / γ

    The time a carried clock actually reads. The straight (inertial) path between two events racks up the most proper time, so the twin who turns around comes home younger — by exactly 1/γ.

  • Causal structurelight cone
    timelike s² < 0 · spacelike s² > 0

    Inside the cone (timelike, s² < 0) cause precedes effect for everyone — the order is absolute. Outside it (spacelike, s² > 0) no signal connects the events and their order is up for grabs.

  • The Lorentz forcebends every path
    F = q(E + v × B)

    The whole instrument in one line. E pushes a charge straight; the v × B term pushes it sideways, always perpendicular to motion — so it does no work and can only turn. Magnetism curves; electricity accelerates.

  • Cyclotron orbitgyration
    rL = mv⊥qB , ωc = qBm

    In a pure magnetic field the charge circles forever at the Larmor radius rL, turning at the cyclotron frequency ωc — which depends on the field and the charge-to-mass ratio, not on speed. Faster particles just trace bigger circles in the same time.

  • E × B driftcharge-blind
    vd = E × BB²

    Cross an electric field with a magnetic one and the gyrating charge slides sideways at a steady drift — and the drift velocity is the same for every charge, positive or negative, light or heavy. The plasma moves as one; that’s the “aha”.

  • Fixed points & bifurcationρ_c ≈ 24.74
    C± = (±√(β(ρ−1)), ±√(β(ρ−1)), ρ−1)

    Below ρ_c ≈ 24.74 the only attractors are the two fixed points C±. Above it they lose stability and the strange attractor is born — a fractal of dimension ≈ 2.06, never a surface, never a volume. This bifurcation is the instrument’s phase knob.

  • The Lorenz systemthree ODEs
    dxdt = σ(y − x) · dydt = x(ρ−z) − y · dzdt = xy − βz

    σ and ρ come from meteorology (Prandtl and Rayleigh numbers). Lorenz derived these equations in 1963 while modelling atmospheric convection rolls — and discovered chaos by accident when he re-ran a simulation from a rounded printout.

  • Lotka–Volterrapredator–prey
    ẋ = αx − βxy  ·  ẏ = δxy − γy

    Prey x breed freely (α) and are eaten on encounter (βxy); predators y grow only by eating (δxy) and starve otherwise (γ). Two coupled rates — the minimal ecology, and the first answer to “if conditions change, what happens to the populations?”

    1925
  • The balance pointa center, not a sink
    x* = γδ  ·  y* = αβ

    The one nonzero state where both rates cancel. But it is a center, not an attractor: nudge the system and it does not return — it circles. Counter-intuitively, x* depends only on the predator’s constants and y* only on the prey’s.

  • The conserved quantityclosed orbits
    V = δx − γ ln x + βy − α ln y = const

    Every trajectory rides a level set of V — Lotka–Volterra is conservative, like a frictionless pendulum or The Orbit’s two-body ellipse. The orbit you start on is the orbit you keep; the swing’s size is a memory of the initial state, not the parameters.

  • The production rulerewriting
    F → F[+F]F[−F]F

    The whole organism is one short rule that rewrites a symbol into a longer string. Here every F becomes a stem that sprouts two branches and carries on — the entire plant is folded into this one line, the way a seed folds the tree.

    1968
  • Turtle interpretationstring → form
    F: step · +/−: turn ±δ · [ ]: push/pop the branch

    A turtle reads the finished string like sheet music: F draws a step, +/− turn it by the angle δ, and the brackets [ ] save and restore its position — which is exactly how one stem remembers where to grow a side branch and then return.

    1990
  • Parallel rewritingall at once
    ωn+1 = P(ωn) — every symbol rewritten in parallel

    Lindenmayer’s twist on a grammar: each generation rewrites every symbol at once, not one at a time — modelling cells that all divide together. That single change turns a dry formal grammar into a growing, living thing.

  • Lyapunov exponentλ > 0 = chaos
    λ = lim 1N Σ ln | r (1 − 2xn) |

    Average stretching per step. λ < 0 → nearby states converge (a predictable cycle); λ > 0 → they separate exponentially (chaos). In the period-3 window near r ≈ 3.83 it dips back below zero — order hiding inside the disorder.

  • Sensitive dependenceλ > 0
    δ(t) ≈ δ0 eλt

    Any gap, however small, grows exponentially. λ is the slope of the climbing line in the inset — measured live by Benettin’s method.

  • Predictability wallLyapunov time
    τ = 1/λ ; δ0 ≤ ε e−λt

    Every extra Lyapunov time τ costs another factor of e in the precision you’d need. You can’t buy much future: the cost grows exponentially, the reach only logarithmically.

  • Lyapunov exponentBenettin 1980
    λ₁ = 1T Σ ln |δ(t)|δ₀

    Two nearby trajectories separate exponentially: δ(t) ≈ δ₀ eλ₁t. If λ₁ > 0 the system is chaotic — any uncertainty in the initial state doubles in 1/λ₁ ≈ 1.1 s. We measure λ₁ live by Benettin’s method: track a shadow orbit, renormalise, accumulate.

  • Self-similaritydim ∂M = 2
    ∂M fractal · mini-M at every scale

    The boundary is endlessly detailed — Hausdorff dimension 2 (Shishikura) — with whole mini-Mandelbrots buried at every depth, and the same shape governs the chaos onset of any smooth family. Infinite structure, five symbols of law.

  • The Mandelbrot setthe atlas
    M = { c : supn |zn| < ∞ }

    The black region: the set of c whose orbit of 0 never escapes. Outside, colour is the smooth escape rate — how many steps the point survives before it flies away.

    1980
  • The Mandelbulbz → zⁿ + c in 3D
    z → rn(sin nθ cos nφ, sin nθ sin nφ, cos nθ) + c

    No 3D multiplication behaves like the complex plane’s, so White & Nylander defined one by analogy: write the point in spherical coordinates (r, θ, φ), multiply the two angles by n and raise the radius to the n. The escape test |z| > 2 carries over unchanged.

    2009
  • The St. Petersburg meanfair(n) ≈ log₂ n
    E = Σk≥1 2−k·2k = ∞,  Snn log₂ n → 1

    A ticket that pays 2k with probability 2−k adds one coin to its expectation for every depth, forever. No single price is fair; Feller showed that n tickets together are fair at about n·log₂ n, so the running mean climbs like log₂ n with no ceiling — the more you buy, the more each is worth.

  • Optional stoppingno system moves the drift
    E[XT] = E[X0]  ⇒  E[P&LT] = −edge · E[stakedT]

    Stop a fair game by any rule that cannot see the next coin and its expected value does not move; give the house an edge and the expected loss is exactly the edge times the expected stake (Wald’s identity). A betting system chooses only when to stop and how much to stake — so it can change the shape of the outcome, never its mean.

  • Heyman’s safe theorema thrust line inside ⇒ safe
    ∃ (H, yc): Ri ≤ rj(H, yc) ≤ Ro for every joint j

    Stone is strong in compression, has no strength in tension and does not slide. Then if any line of thrust in equilibrium with the loads lies inside the masonry, the arch stands; it fails only when no line fits and hinges open where the line touches the edge.

  • Conservation of massclosed vessel
    Σ mbefore = Σ mafter (sealed)

    In a closed vessel nothing is made and nothing lost: weigh everything and the books close. Lavoisier’s sealed retort made it a measurement — the calx gains exactly what the air inside loses.

  • Maxwell's demonthought experiment
    ΔSsort + ΔSerase ≥ 0

    Maxwell’s 1867 thought experiment: a tiny intelligent being controls a trap-door and sorts fast from slow molecules, creating a temperature difference without doing work — apparently defeating the second law. The resolution, completed by Landauer and Bennett, is that the demon’s memory must be erased, and that erasure generates entropy that restores the ledger.

  • The fading iconletters available vs cue delay
    A(d) = A∞ + (A0 − A∞) · e−d/τ

    The letters available d ms after the flash — the cued row’s score times the number of rows — falling from A0 towards what whole report gives. Sperling plotted the curve; the exponential is the standard summary of it, fitted here by least squares.

  • Cowan’s Khow many you were holding
    K = N · (H − FA)

    For a change-detection test that probes one item of N: if K of them are held, a change to a held item is always caught and the rest are guessed, so the hit rate H exceeds the false-alarm rate FA by K/N. Below capacity K follows N; above it, K levels off — at about four.

  • Metropolis acceptancesometimes uphill
    Paccept = min(1, e−ΔE / T)

    Always take a downhill move; take an uphill move of cost ΔE only with probability e−ΔE/T. High T accepts almost anything; as T → 0 only improvements pass. This one rule lets a search climb out of a trap.

    1953
  • The catenoidr = c cosh(x/c)
    r = c cosh(x/c),   R = c cosh(h/2c)

    The film between two rings of radius R a distance h apart is a catenary spun round the axis. Its waist c solves R = c cosh(h/2c), which has a root only while h/2R < 0.6627 — then the film must snap.

  • Minimal surfaceH = 0
    H = ½(κ1 + κ2) = 0

    A film with the same air on both sides has no pressure across it, so by Young–Laplace its mean curvature is zero: curved one way exactly as much as the other, at every point — the surface of least area for its boundary.

  • Chebyshev equioscillationbest uniform
    f(xi) − p*(xi) = (−1)i E, i = 0 … n+1

    The degree-n polynomial closest to f in the worst case is the one whose error E = ‖f − p*‖∞ is reached n + 2 times with alternating sign. Born from trying to straighten Watt’s linkage: spread the error evenly across the stroke instead of making it tiny in the middle.

  • Minority-game volatilityone curve for every N
    σ²N = F(α),   α = 2mN,   A = Σ aᵢ

    N agents each pick a side, aᵢ = ±1, with the best of S strategies that read the last m minority sides; the minority wins. σ² is the variance of A. Coin tosses give σ²/N = 1. Measured for any N and m, σ²/N falls on one curve of α: far above 1 when α is small (the crowd herds), a minimum near α_c, back toward 1 when memory outgrows the crowd.

  • The lossL
    L = ½ ⟨(ŷ − y)²⟩

    How wrong the network is, averaged over every point. Squared so that big misses dominate and the surface is smooth — a landscape in weight space the optimiser can roll down.

  • Monte Carlo1/√N
    ŜT(1…N) → log-normal

    Monte Carlo is the other half: roll N independent futures and the empirical distribution converges to the closed form at rate 1/√N. Watch the ember cloud tighten onto the cone as you raise N.

    1949
  • The Reichardt detectordelay and multiply
    R(t) = s1(t − Δ) · s2(t) − s1(t) · s2(t − Δ)

    Two neighbouring inputs; each is delayed and multiplied with the other, and the two products subtracted. A pattern drifting from 1 to 2 makes the first product win, the opposite drift the second — an opponent motion signal from nothing but timing.

  • The muon zenith lawcos²θ
    I(θ) ≈ Iv cos2θ

    At sea level the muon intensity falls with zenith angle roughly as cos²θ — the slant path through the atmosphere is longer, and more of the parent pions and the muons themselves decay or stop on the way. A horizontal detector adds one more cos θ for its projected area.

  • Muon energy loss in matterdE/dX = −(a + bE)
    Emin(X) = ab (ebX − 1)

    Ionisation costs a muon about a = 2 MeV per g/cm²; above a few hundred GeV radiation (bremsstrahlung, pair production) adds bE with b ≈ 4 × 10⁻⁶ per g/cm². To cross a column density X it must start with Emin — about 46 GeV through 100 m of stone — and few do.

  • Counting significance(N − μ)/√μ
    z = N − μ√μ

    N counts observed where the model expects μ: Poisson noise is √μ, so a fixed fractional excess grows in significance as √t. A claim waits for 5σ — and for the excess to agree from more than one side, because the model of the stone is never perfect.

  • Nash equilibriumno solo gain
    ui(σ*i, σ*−i) ≥ ui(si, σ*−i)

    A profile where no player can do better by changing strategy alone. Every finite game has one. In the one-shot dilemma it is mutual defection — individually unbeatable, collectively the worst stable outcome.

  • Incompressible Navier–Stokesmomentum + mass
    ∂u∂t = −(u·∇)u − ∇pρ + ν∇²u + f  ·  ∇·u = 0

    The whole of classical fluid motion in two lines: a parcel accelerates from pressure gradients, viscous drag (ν∇²u) and any body force f. The constraint ∇·u = 0 says an incompressible fluid can swirl freely but never pile up.

  • The neuronforward
    a = act(W x + b)

    Each neuron takes a weighted sum of the layer below plus a bias, then bends it through a non-linearity. Stack two such layers and the network can carve curved, disconnected regions — not just a single straight cut.

  • Newton’s methodorder 2
    zn+1 = zn − f(zₙ)f′(zₙ)

    Stand on the curve, slide down the tangent to where it crosses zero, repeat. Near a root the error is squared every step — correct digits double — which is why this three-century-old update still sits inside every solver and every GPS fix.

  • Basins of attractionfractal border
    B(rk) = { z₀ : zn → rk }

    Colour every starting point by the root it falls to and the plane divides into basins. Two roots: a straight border, as Cayley proved. Three or more: the borders coincide — every boundary point touches all the basins at once, and the map of destiny turns fractal.

  • The damped step0 < a ≤ 1
    zn+1 = zn − a · f(zₙ)f′(zₙ)

    Take only a fraction a of each tangent step and the same roots keep their basins — but the borders between them swell, curl and bloom. The fractal isn’t decoration on the method; it is the method, rendered at the places where the next step is undecided.

  • Newtonian N-bodyinverse square
    ai = Σj≠i G mj(rj−ri)(|rj−ri|² + ε²)3/2

    Every mass pulls every other. ε is Plummer softening — it rounds off the 1/r² singularity so a near-collision integrates cleanly instead of exploding.

  • Noisy-channel theoremR < C ⇒ error → 0
    R < C(p*) reliable · R > C impossible

    The headline of 1948: codes exist that drive error to zero for any rate below capacity — and none above it. The wall p* is where R meets C. Our codes give out earlier; better ones reach all the way to it.

  • The dominance cycleno maximum
    A ≻ B ≻ C ≻ D ≻ A,   P(X beats next) = 23

    “Beats more often than not” is not an order: it can run in a circle. Efron’s four dice each beat the next two times in three, and the last beats the first — so whatever is chosen first, there is an answer that wins. The arrow refuses to be a straight line.

  • Conway’s leading numbersexact odds
    odds(B before A) = (AA − AB) : (BB − BA),   E[TA] = 2·AA

    XY reads as a binary number: its bit k is 1 when the last k symbols of X equal the first k of Y. Four such overlaps give the exact odds that B appears before A in a fair-coin stream — HHH loses to THH 7 : 1 — and the same self-overlap AA sets the mean wait (HHT 8 flips, HTH 10). Against b₁b₂b₃ the beater is b̄₂b₁b₂.

  • The light-bending equationnull geodesic
    u″ + u = 3(GM/c2) u2, u = 1/r

    A photon travels perfectly straight; spacetime is curved, so its path curves with it. This is the exact orbit equation the instrument integrates per ray — not a thin-lens shortcut like The Well. The extra u² term is general relativity; drop it and you get Newton’s half-strength, straight-ish light.

  • The number sensethe Weber fraction w
    P(correct) = Φ( |n1 − n2| / (w · √(n1² + n2²)) )

    Each numerosity n is represented as a noisy magnitude whose spread grows with n (scalar variability, Weber fraction w). The chance of picking the larger of two depends only on their ratio; w is fitted to the visitor’s answers by maximum likelihood.

    2008
  • Nyquist limitor it aliases
    fs > 2 ffeat

    Sample slower than twice the feature rate and the motion is lost — it can no longer be told apart from a slower one. This single inequality is the whole limit.

  • 1/f noiseflicker
    S(f) ~ 1 / fα

    Scale-free in time as well as size: the activity’s power spectrum runs as 1/f, the ubiquitous flicker noise found in resistors, rivers, heartbeats and starlight — the puzzle self-organized criticality was invented to explain.

  • Miura-ori kinematicsone degree of freedom
    S = b cos θ tan γ / √(1 + cos²θ tan²γ),   L = a √(1 − sin²θ sin²γ),   ν = −tan²γ cos²θ

    Parallelograms of sides a, b and angle γ, each tilted θ out of the plane: every vertex of the sheet follows from θ alone. Folding shrinks both the run S and the rise L at once, so the in-plane Poisson’s ratio is negative.

  • PageRankPage & Brin 1998
    r(i) = 1−dN + d · Σj→i r(j)L(j)

    A node’s rank is the rank poured into it along incoming links, each source splitting its own rank evenly across its out-links. Importance is recursive — you matter if things that matter point at you — plus a small flat share (1−d)/N so nothing starves.

    1998
  • The sign of a shuffle±1
    sgn(σ) = (−1)K, K = #inversions

    Count the pairs standing out of order and take minus-one to that power. Every arrangement is even or odd — a single swap flips the sign, so no chain of swaps can ever change a shuffle’s parity by stealth.

  • The locked productnever moves
    I = sgn(σ) · (−1)d(□) — conserved

    Both factors flip on every slide, so their product never does. Loyd’s deal carries I = −1; the solved board carries +1. The two values are two islands, and sliding never builds a bridge.

  • One slide, one swapflips both
    σ ← τt,□ ∘ σ ⇒ sgn ← −sgn, (−1)d ← −(−1)d

    A legal slide is exactly one transposition — the tile trades cells with the blank — and exactly one unit step of the blank. Two coins, flipped together on every move, with no third option on the board.

  • The two losing gamesboth losing
    A: p = ½ − ε  ·  B: p = 110 − ε if 3 | W, else 34 − ε

    Game A is a plain biased coin; game B picks its coin by whether your capital is divisible by 3 — a vicious coin on multiples of three, a generous one elsewhere. Played alone, each drifts your capital down. The law never mentions alternation.

  • The verdict of the chainexact
    g = π0(2pbad − 1) + (1 − π0)(2pgood − 1)

    The drift per round, exactly: π₀ is how often the capital sits on a multiple of 3 — the stationary state of a three-state Markov chain. Game B digs its own trap (π₀ high, the bad coin often); mixing in A shuffles π₀ down, and the same formula changes sign.

  • The small-angle liesin θ ≈ θ
    T0 = 2π √(l / g)

    Let sin θ ≈ θ and the amplitude vanishes from the period — every swing, wide or narrow, takes the same time. This is Galileo’s isochrony, the founding approximation of timekeeping, and it is only true in the limit of swings too small to see.

  • Equation of motion · θ₁″exact
    θ1″ = −g(2m1+m2)sinθ1 − m2g sin(θ1−2θ2) − 2sinΔ·m2(ω22l2+ω12l1cosΔ)l1(2m1+m2−m2cos2Δ)

    Straight from the Lagrangian L = T − V. No approximation, no randomness — every pendulum obeys this to the letter (Δ = θ₁ − θ₂).

  • Equation of motion · θ₂″Δ = θ₁ − θ₂
    θ2″ = 2sinΔ(ω12l1(m1+m2) + g(m1+m2)cosθ1 + ω22l2m2cosΔ)l2(2m1+m2−m2cos2Δ)

    The two arms are nonlinearly coupled — the source of the unpredictability. We integrate both with RK4 so energy doesn’t drift.

  • The pendulum lawone line
    θ″ = −gl sin θ

    A mass on a pivot under gravity — the most familiar law on Earth, and the whole of it. One line, no forcing, no friction, no dice. The law says nothing about how many joints the mass hangs from; that silence is where the two tempers come from.

  • The whole truthexact
    T = 4 √(l / g) · K(sin (θ0 / 2))

    The exact period, with K the complete elliptic integral — computable to machine precision by the arithmetic–geometric mean. The amplitude never left; it was hiding in K. At θ₀ = 90° the clock runs 18% slow, and as θ₀ → 180° the period grows without bound.

  • The forbidden linereleased at rest
    2 cos θ1 + cos θ2 = 1

    Energy alone decides who may flip: released at rest inside this line, a double pendulum (equal masses, equal arms) can never send its second arm over the top — however long you wait. Outside it, the flip is allowed, and the fractal decides how long you wait.

  • Cluster-size distribution at p_cτ = 187/91
    ns ~ s−τ

    At p_c the cluster-size distribution is a power law — there is no characteristic scale, clusters of all sizes coexist. This scale-free structure is why the critical point looks fractal and why the exponents are the same across wildly different physical realisations.

  • Correlation lengthν = 4/3
    ξ ~ |p − pc|−ν

    The correlation length — the typical size of the largest finite cluster — diverges as p approaches p_c from either side, meaning the system "sees" arbitrarily large distances just at the critical point.

  • Phase transitionp_c ≈ 0.5927
    P∞(p) ~ (p − pc)β, β = 5/36

    Below p_c only finite clusters exist. At p_c a giant cluster of infinite extent first appears (in the thermodynamic limit). Above p_c the giant component density P∞ grows as a power law — the order parameter of the transition.

  • Universalityd = 2 universality class
    β, ν independent of lattice

    The exponents β, ν, τ are universal — they do not depend on the lattice geometry (square, triangular, honeycomb) or the physical realisation. Only the spatial dimension d and the symmetry of the order parameter matter. That universality is why the same numbers describe water in rock, fire in a forest, and contagion in a network.

  • The strategic responsemove at a quadratic cost
    x′ = argmaxx [ θ·x − ‖x − x₀‖²2ε ] = x₀ + εθ

    A published rule invites a best response: change your features to improve your verdict, at a cost. With a reward linear in the rule and a quadratic cost the response is a shift of εθ — the location family in which the stable point and the optimum have closed forms.

    2016
  • Repeated risk minimisationretrain on what you caused
    θt+1 = argminθ EZ∼D(θₜ) ℓ(Z; θ)

    Deploy θₜ, collect the data D(θₜ) the deployment produces, refit on it, deploy again. With a γ-strongly convex, β-smooth loss and a distribution map that moves ε per unit of θ, each step contracts by εβ/γ.

    2020
  • Performative stabilityε < γ/β
    θPS = argminθ EZ∼D(θ_PS) ℓ(Z; θ),   ε < γβ

    A model that is optimal for the data it itself induces — calibrated against the future it causes. Retraining reaches it when ε < γ/β. It is not in general the performative optimum, which minimises the loss counting the data’s reaction.

    2020
  • Period-doubling cascader∞ ≈ 3.5699
    r₁=3 · r₂≈3.449 · r₃≈3.544 · r₄≈3.564 · … → r∞

    As r climbs, the stable cycle doubles: 1 → 2 → 4 → 8 → … The thresholds crowd together and accumulate at r∞ ≈ 3.5699, where the period becomes infinite — the door into chaos.

  • Perron–Frobeniusunique, positive
    G v = λ₁ v,  λ₁ = 1,  v > 0

    A positive matrix has one largest eigenvalue, and its eigenvector is all-positive and unique. That theorem is the guarantee under PageRank: exactly one rank vector, every entry a real fraction, and power iteration always finds it.

    1912
  • Diffraction of a helixthe X
    Il(R) ∝ Jl²(2π r R), Z = l / P

    A helix of radius r and pitch P diffracts onto layer lines spaced 1/P apart; on the l-th line the intensity follows the Bessel function Jl, whose first peak moves outward as l grows — so the brightest spots march away from the axis and draw an X.

  • Error-reductionproject, project
    g′ = F−1[|F|·ei arg G], g ← g′ on S where g′ ≥ 0, else 0

    Impose the measured modulus, then impose the object constraints — real, non-negative, inside the support S — and repeat. Each step cannot raise the error, which is why it so often stops at the wrong place: error-reduction stagnates.

  • Hybrid input-outputfeedback escapes
    g ← g′ where the constraints hold, else g − β·g′

    Fienup’s fix: where the new estimate violates the constraints, do not zero it — push back against it with feedback β (≈ 0.9). The iterate is no longer an estimate but a driving signal, and it climbs out of the traps error-reduction falls into; a few ER steps at the end polish.

  • The phase problemhalf the evidence
    I(k) = |F(k)|², F(k) = ∫ g(x) e−2πi k·x dx

    A detector records how brightly each spatial frequency lands, never its phase — and the phases carry where everything is. The object g(x), its 180° twin g(−x) and every translation share one |F|: those ambiguities are the null space of this evidence.

  • Phase vortexquantised circulation
    ∮ ∇φ · dl = 2πn, n ∈ ℤ

    Trapped at a node, the phase has nowhere to point — so it winds a whole turn around the point. A quantised vortex, its charge always an integer. The ones the rings mark.

  • The phlogiston balanceobsolete
    metal = calx + φ · air + φ → airphlogisticated

    Stahl’s ledger: burning, calcining and breathing all release phlogiston, φ, into the air, which can hold only so much — at saturation the fire dies. Read backward it is Lavoisier’s, calx = metal + oxygen, term for term: the arrow was wrong, the bookkeeping was not.

  • The photon spherer = 1.5 r_s
    rph = 3GM/c2 = 1.5 rs

    At exactly this radius gravity bends light into a circle. Rays that graze it wind around the hole many times before escaping, stacking into the thin bright photon ring at the rim of the shadow.

  • Trail sensingthree sensors
    sample FL, FC, FR → turn ±α toward the strongest

    Every agent reads the chemoattractant at three points ahead — left, centre, right — and rotates by the sensor angle α toward whichever is brightest, else holds course. No map, no goal: just follow the smell of where the swarm has already been.

  • Adaptive conductivitythe Tokyo rule
    dDijdt = f(|Qij|) − Dij

    Tero’s distillation: a tube’s conductivity D grows with the flux Q it carries and decays otherwise. Heavily-used routes thicken, idle ones wither — a feedback that prunes a tangle down to an efficient, fault-tolerant network with nobody designing it.

  • Deposit · diffuse · decaythe field
    T ← (1 − ρ)·(K ∗ T) + d·[agent here]

    The trail map T is the only shared memory: each agent drops a little d where it stands, the field blurs (K ∗) so neighbours can smell it, and a fraction ρ evaporates each step. Reinforcement plus evaporation is what lets unused paths fade and used ones thicken.

  • The pitch helixchroma round, height up
    (x, y, z) = (r sin 2πc/12, h · (o + c/12), r cos 2πc/12)

    Pitch as two dimensions at once: the chroma c (which note) goes round the circle, the height (octave o) goes up, one turn per octave. Tones an octave apart stand directly above one another; a Shepard tone lights the same chroma at every turn.

  • The Shepard tonea pitch class without an octave
    s(t) = Σk G(log₂ fk) · sin(2π fkt), fk = f0 · 2k + c/12

    Partials an octave apart, weighted by a fixed bell-shaped envelope G over log frequency. Raising the pitch class c by twelve semitones gives back exactly the same set of partials — so the tone has a chroma but no definite height.

  • Poincaré recurrenceτ ~ e^S
    P(return) ≈ 1 : 2N

    The gas will return arbitrarily close to its start — after a time that dwarfs the age of the universe. "Never," in practice, is a statement about how long you’d wait.

  • The Schläfli symbolrecipe
    {p, q, r}

    The whole recipe for a regular 4-polytope in three numbers: {p} faces, {q} around each vertex, {r} cells around each edge. The tesseract is {4,3,3} — cubes ({4,3}), three to an edge.

  • The 4-D double rotationtwo planes
    R = RXW(α) ∘ RYZ(β)

    A generic rotation in four dimensions turns in two completely independent planes at once. Drive the angles α and β at different rates and the motion never repeats — there is no axis to point at, and no 3-D analogue.

  • Perspective projectionthe shadow
    (x, y, z) ↦ (x, y, z)d − w

    Cast the 4-D point onto our space by dividing through its distance d − w from a 4-D eye. Cells nearer in the unseen axis w loom larger — which is why the inner cube can swell, pass through the outer one, and become it.

  • The 4-D Euler relationalways 0
    V − E + F − C = 0

    Euler’s V − E + F = 2 for solids, lifted a dimension: count vertices, edges, faces and cells with alternating sign and you always get zero, whichever of the six regular 4-polytopes you pick.

  • The swing equationsecond-order Kuramoto
    M θ̈i = Pi − D θ̇i − Σj Kij sin(θi − θj)

    Every generator and consumer a rotating mass in the frame turning at 50 Hz: power in, damping out, and each line pulling its two ends toward one phase. Kuramoto’s model with inertia. Two buses ±P on one line lock exactly when K ≥ |P| — The Grid’s engine locks at 1.02·P and slips at 0.98·P.

  • The DCB synchronisation conditionDörfler, Chertkov & Bullo 2013
    ‖BT L† ω‖∞ ≤ sin γ, γ < π/2 ⇒ a unique stable sync with |θi − θj| ≤ γ on every line

    L is the Laplacian weighted by the line couplings, ω the first-order natural frequencies, B the incidence matrix: each line’s |(L†ω)i − (L†ω)j| is the strain it will carry. Under 1 everywhere, the grid can lock. Over 150 random draws in The Grid it matched the swing equation 96.7% of the time, and never promised sync to a grid that slipped.

    2013
  • Droop controlprimary control
    Pi ← Pi0 − R · Δωi

    No dispatcher in the first seconds: a plant that sees its own frequency sag raises its output in proportion, and every plant does it at once, so a lost plant’s share is spread across the continent and the frequency settles slightly low instead of collapsing.

  • Power iterationr ← G r
    rk+1 = G rk → v₁

    Just apply G again and again: any starting guess is dragged onto G’s leading eigenvector. It converges at the spectral gap — about a factor d per sweep — which is why ~50 passes rank the whole web.

  • The avalanche lawτ ≈ 1.2
    P(s) ~ s−τ

    Avalanche sizes have no characteristic scale: a straight line on a log–log plot. Most are tiny, a few are system-spanning, and the same law fits earthquakes, forest fires, extinctions and market crashes — the rare giant is not an anomaly, it is the tail of the everyday.

  • The Euler productsum = product
    ζ(s) = Σn 1ns = Πp 11 − p−s

    A sum over all the integers equals a product over the primes — Euler’s bridge from analysis to arithmetic, and the reason the zeta function knows where the primes are.

  • The logarithmic integralbest estimate
    Li(x) = ∫2x dtln t

    Gauss’s sharper guess: integrate 1/ln t instead of taking it flat. Li(x) hugs the π(x) staircase far more tightly than x / ln x while obeying the same leading law.

  • The Riemann zerosthe music
    π(x) = li(x) − Σρ li(xρ) − … , ζ(ρ) = 0

    Riemann’s explicit formula writes the prime count exactly as a smooth term corrected by a sum over the zeros ρ of ζ — the primes are the spectrum of the zeta function. The Hypothesis: every ρ has Re(ρ) = ½.

  • The Prime Number Theoremthe law
    π(x) ~ xln x , so π(x) ln xx → 1

    The count of primes up to x grows like x / ln x — they thin out, but only logarithmically. Gauss guessed it at fifteen, Legendre published it, and it went unproven for a century.

    1896
  • The Prisoner's DilemmaT > R > P > S
    T > R > P > S,  2R > T+S

    Reward R for mutual cooperation, punishment P for mutual defection, temptation T for defecting on a cooperator who gets the sucker payoff S. Whatever the other does, defecting scores more — yet two defectors both do worse than two cooperators. The trap of rational self-interest.

  • Probability conequantile
    S0 e(μ − σ²/2)t ± z·σ√t

    Because log S is Gaussian, every quantile is a formula: the silver cone. Its half-width grows like σ√t — uncertainty widens with the root of time, the same law as diffusion.

  • The hidden halfprobability current
    j = |ψ|² ∇φ = Im(ψ* ∇ψ)

    Probability actually flows, and its current is set by the gradient of the phase φ. Flat phase, no flow; a winding phase is a current going somewhere.

  • Weber’s lawa constant ratio
    ΔI / I = k

    The smallest change ΔI that can be told apart grows in proportion to the intensity I: a fixed fraction k, the Weber fraction. For number, the same law means two quantities are discriminated by their ratio.

  • Fechner’s lawsensation ∝ log stimulus
    S = k · ln(I / I0)

    Count just-noticeable differences up from the threshold I0: if each is a constant fraction of I (Weber), the count — the sensation S — grows with the logarithm of the stimulus. The logarithmic number line is the same idea laid out in space.

  • The transformed staircase2-down / 1-up → 70.7 %
    P(+)n = ½ ⇒ P(+) = 0.5 · 0.707 · 0.794 (n = 1 · 2 · 3)

    Step down after n positive answers in a row, up after one negative: the track can only stop moving where a run of n positives is as likely as not. Two interleaved tracks — one rising, one its mirror — bracket the PSE from both sides and spend the trials where the curve is steep.

  • The virtual chinrestdistance from the blind spot
    d = s / tan α, α ≈ 13.5° · px/° = 2·d·tan 0.5° · (px/mm)

    The dot vanishes at the blind spot’s entry, a distance s from the fixation cross on the screen (in mm, from the card match). Taking that entry at α ≈ 13.5° gives the viewing distance d, and d gives the pixels in a degree. Measure d with a ruler instead and α becomes a measurement of your own.

    2020
  • The psychometric functionevery Iris trial
    Ψ(x) = γ + (1 − γ − λ) · F(x; α, β)

    The probability of one answer as the stimulus x changes: F a sigmoid (logistic or cumulative normal) centred on α, the point of subjective equality, with spread β; γ the guess rate, λ the lapse rate. Every Iris instrument fits it to the visitor’s own trials by maximum likelihood, with the PSE’s error from the profile likelihood.

  • Which story the trials favoura likelihood ratio
    Λ = max L(postdiction | trials)max L(extrapolation | trials)

    Each model predicts the three conditions’ PSEs from one lag: extrapolation (L, 0, L), postdiction (L, L, 0) for continuous, flash-initiated, flash-terminated. Each is fitted to every trial; the ratio of their best likelihoods says which the data favour — and how strongly.

  • A prior on the light’s directionvon Mises
    p(θ) ∝ eκ cos(θ − μ) ⇒ P(bump | φ) = σ(2κ cos(φ − μ) + logit c)

    A circular normal over where the light comes from, centred on μ with concentration κ. The image allows only two stories — a bump lit from φ or a dent lit from φ + 180° — so the posterior odds are the prior’s ratio at those two angles, times the odds c of convexity.

  • The aftereffect’s decayτ
    vnull(t) = A · e−t/τ

    The real motion needed to cancel the aftereffect, t seconds after the adapting motion stops: strength A, falling by a factor e every τ — the time an adapted channel takes to win back its gain.

  • Sensitivity d′signal detection
    d′ = z(H) − z(FA)

    The distance between the “changed” and “same” evidence distributions in units of their spread, from the hit and false-alarm rates through the inverse normal z. Unlike percent correct it does not move when the observer simply says “changed” more often. A log-linear correction keeps it finite at 0 and 1.

  • The visual anglesize as the eye sees it
    θ = 2 · atan(s / 2d)

    An object of size s at distance d fills the angle θ at the eye — the unit every Iris stimulus is specified in. At 60 cm a CSS pixel of 0.26 mm is about 1/40 of a degree; the same picture on a phone held at 30 cm is twice as large in the eye.

  • Brouwer’s fixed-point theoremcontinuous ⇒ a fixed point
    g : Bⁿ → Bⁿ continuous ⇒ ∃ x* : g(x*) = x*

    A continuous map of a closed ball into itself fixes some point. On the interval [0, 1] it is the intermediate value theorem: a continuous curve from the bottom edge to the top must cross the diagonal. A jump can step over it.

  • The public prediction that comes trueR(f*) = f*
    R : [0, 1] → [0, 1] continuous ⇒ ∃ f* : R(f*) = f*,   stable iff |R′(f*)| < 1

    R is what happens when the forecast f is published. A forecast that accounts for its own effect is a fixed point of R; Brouwer says one exists when R is continuous. Iterating publication — publishing what happened last time — reaches it only where R is a contraction.

    1954
  • The bandwagon responseb > 0 bandwagon · b < 0 underdog
    R(f) = 1N Σᵢ σ(zᵢ + b·(f − ½)),   R′(f) = bN Σᵢ σ′(zᵢ + b·(f − ½))

    The Prophecy’s model electorate: voter i leans zᵢ, and a published poll f moves every voter toward the leader it shows by the gain b (against it when b < 0). Simon’s argument needs only that R be continuous; this logistic form is the instrument’s choice, smooth and with R′ in closed form.

  • Womersley numberα = 1
    α = r·(ωρ/μ)1/2

    Vessel radius against the viscous skin depth of one heartbeat. Above one, blood travels as a reflected-free pulse wave and area-preserving branches match impedance; below one it is viscous Poiseuille flow, and radii should shrink as n−1/3 instead.

  • The Hadamard gatesuperpose
    H : |0⟩ ↦ |0⟩ + |1⟩√2, |1⟩ ↦ |0⟩ − |1⟩√2

    The gate that makes superposition: it turns a definite bit into an even 50/50 mix — but with a sign. Apply it twice and the signs interfere to undo it, which is the whole trick of quantum algorithms.

  • The CNOT gateentangle
    CNOT : |a, b⟩ ↦ |a, a ⊕ b⟩

    Flip the second bit only if the first is 1. Fed a superposed control it can’t be written as “this qubit times that qubit” any more — the two become one inseparable state. This is where entanglement enters.

  • The Bell pairone circuit
    |Φ⁺⟩ = |00⟩ + |11⟩√2

    A Hadamard then a CNOT, and two qubits are perfectly correlated: each alone is a coin-flip, yet measure one and the other’s answer is fixed — the computational cousin of The Pact’s singlet.

  • The quadratic mapz₀ = 0
    zn+1 = zn2 + c

    Feed the rule its own output, starting from zero. If |z| ever passes 2 the orbit escapes to infinity; otherwise it is trapped forever. That single yes/no, asked at every point, draws everything here.

  • The Bloch sphereon the sphere
    |ψ⟩ = cosθ2|0⟩ + eiφ sinθ2|1⟩

    Drop the invisible overall phase and every pure qubit state is a point on a sphere: θ tips it between the poles |0⟩ and |1⟩, φ spins it round the equator. The whole one-qubit world is its surface.

  • The qubita unit vector
    |ψ⟩ = α|0⟩ + β|1⟩, |α|² + |β|² = 1

    Not a 0 or a 1 but a blend of both — two complex amplitudes whose squares must sum to one. A classical bit is the two ends; the qubit is the whole line between them.

  • The M/M/1 waitdiverges at ρ = 1
    W = 1/μ1 − ρ  ·  ρ = λ/μ

    One server, Poisson arrivals at rate λ, exponential service at rate μ: the mean time in the system is the service time divided by the idle share 1 − ρ. Ten times the service at 90% load, a hundred at 99% — the wait diverges while the work stays finite.

  • Little’s lawany queue
    L = λ·W

    The mean number in a system equals its arrival rate times the mean time each one spends inside — for any queue in steady state, whatever the laws of arrival and service. Count heads, or time visitors: two independent measurements of one quantity, an audit built into the subject.

  • Kingman’s formulavariance × utilization × time
    E[Wq] ≈ ρ1 − ρ · Ca2 + Cs22 · 1μ

    The heavy-traffic wait in line is three levers multiplied: load, variability (the squared coefficients of variation of arrivals and service) and speed. Clockwork arrivals halve the variance term and so the wait — but nothing touches the 1/(1 − ρ) wall.

  • The Bring–Jerrard quinticreduced form
    x5 + a x + b = 0

    Every general quintic reduces, by a Tschirnhaus transformation, to this two-parameter form — the object the braid floor actually drags: two numbers, five roots, and a coefficient plane whose loops are the whole instrument.

  • On a networkNewman 2002
    R₀ = T · (⟨k²⟩ − ⟨k⟩) / ⟨k⟩

    On a contact network the right number uses the excess degree — who your neighbours’ neighbours are. A heavy-tailed (scale-free) degree distribution makes ⟨k²⟩ blow up, so hubs push R₀ far above the well-mixed estimate.

  • The Stefan–Boltzmann lawthe whole oven
    u = aT4 · j = σT4

    The total energy in a cavity, and the power from each square metre of a black surface, grow as the fourth power of the temperature — finite, and measured, before anyone could derive the spectrum it sums.

  • Wien’s displacement lawthe peak
    λmax · T = b ≈ 2.898 mm·K

    The hotter the body, the bluer its peak: a kiln at 1600 K peaks in the infrared near 1.8 µm, the Sun near 500 nm.

  • Planck’s lawthe continuum
    Bλ(T) = 2hc²λ⁵ · 1ehc/λkT − 1

    The light of a body in thermal equilibrium at temperature T. Its peak slides blueward as T rises — the colour of a star — and to derive it Planck had to let energy come in quanta.

  • The Rayleigh–Jeans lawexact · low ν
    u(ν, T) = 8πν²c³ · kT

    Equipartition, applied without error: a cavity holds 8πν²/c³ standing waves per unit volume and frequency, and each carries kT. Exact where hν ≪ kT; summed over every mode it diverges — the ultraviolet catastrophe.

  • Random walk⟨x²⟩ = 2Dt
    Sn = Σ Xi ; Var = nσ²

    A sum of i.i.d. steps goes nowhere on average, yet its spread grows as the square root of the number of steps — the signature √t law of every diffusion.

  • Lateral inhibitionHartline–Ratliff
    rp = ep − Σj kpj (rj − r⁰pj)

    Each unit’s response r is its own excitation e minus its neighbours’ responses above a threshold r⁰, weighted by k — measured in the horseshoe crab’s eye. A network of such units exaggerates every edge it sees.

  • Difference of Gaussianscentre − surround
    B = (Gσc − w·Gσs) * L1 − w

    A narrow excitatory centre minus a broad inhibitory surround of weight w, normalised so a uniform field comes through unchanged. On a ramp it is silent in the middle and overshoots at the knees — a bright band where the ramp meets the bright plateau, a dark one where it meets the dark.

  • Kac’s lemmathe mean return
    ⟨TA⟩ = Δ / μ(A)

    Look at a stationary system every Δ; the mean time between looks that find it in a set A is exactly Δ over A’s share of the states. Poincaré says the return is certain; Kac says how long it takes — and for a gas, μ(A) shrinks geometrically with every particle.

  • Gravitational redshifttime slows
    1 + z = 1√(1 − rs/r)

    Clocks tick slower deeper in the well; light climbing out loses energy and reddens, diverging as r → rs. Together with the disk’s orbital Doppler shift it makes the approaching side blaze blue-white and the receding side fall dim and red.

  • Harmonics = particlesRegge
    α′M² = n − 1

    The string-theory punchline: the same string at oscillator level n is a different particle, its mass climbing the Regge trajectory. n=1 is massless; the spectrum above is the particle zoo.

    1968
  • The inspection paradoxlength-biased
    E[W] = E[X2]2·E[X] = (1 + C2)·E[X]2

    Arrive at a random moment and you land in a gap with probability in proportion to its length — so you sample long gaps more often than they occur, and wait, on average, E[X²]/2E[X]. Clockwork gaps (C = 0) give half the mean gap; Poisson gaps (C = 1) the whole of it; any burstier schedule more than the gap it promised.

  • Reproduction numberR₀
    R₀ = β / γ  (·⟨k⟩ on a network)

    The single most important number in epidemiology: the average number of new infections one case produces in a fully susceptible population. Above 1 the outbreak grows; below 1 it dies.

  • The reservoir updateweights frozen
    xt+1 = tanh( W xt + Win ut )

    A big pool of neurons, wired to each other by a fixed random matrix W and fed the input through random Win. Its state x is a living, fading echo of everything it has heard — a high-dimensional memory that is never trained, only listened to.

  • The linear readoutthe only trained part
    yt = Wout xt

    The answer is just a weighted sum of the reservoir’s state. All the learning lives in this one output layer Wout — the recurrent tangle stays random, so training a recurrent network collapses to fitting a single line.

  • Ridge regressionclosed form
    Wout = Y X⊤( X X⊤ + βI )−1

    Collect the reservoir states X against the targets Y and solve one regularised least-squares system — no backpropagation, no epochs, one matrix inverse. The ridge term β keeps the readout from chasing noise in the echo.

  • The echo state propertyfading memory
    ρ(W) < 1

    For the network to be a reservoir at all, it must forget where it started: scale W so its spectral radius — its largest eigenvalue — sits just below one. Below, old inputs decay and the state is a clean echo of recent history; above, the activity never settles and prediction falls apart.

  • Retinexratios across edges
    log ℓB − log ℓA = Σi T[log Li+1 − log Li]

    Walk from A to B summing the log luminance ratios between neighbours, but pass each through a threshold T that zeroes gentle steps: sharp steps are paint edges, gradual ones are illumination. A soft shadow edge is discounted; a hard one is read as paint.

  • Reversibility paradoxLoschmidt 1876
    v → −v ; ΔS ≥ 0

    Loschmidt’s objection: reversible micro-laws can’t logically give an irreversible world. The resolution is the Past Hypothesis — the universe began in a low-entropy state — plus sensitive dependence, which is why real reversal fails.

    1876
  • Reynolds numberRe = UL / ν
    Re = U Lν

    One dimensionless number decides everything: the ratio of inertia to viscosity. Low Re, the fluid moves in smooth sheets; past a critical value (~2000 in a pipe) it shatters into eddies on every scale. The readout above is this number, in the instrument’s own units.

  • Euler's equationstorque-free
    I1ω̇1 = (I2 − I3)ω2ω3 , & cyclic

    In the body frame, with no torque, the spin obeys three coupled equations. The sign of (I2 − I3) and its cyclic kin decides whether a small wobble decays or grows.

  • Two locked invariantssphere ∩ ellipsoid
    |L|² = Σ (Iiωi)² , 2T = Σ Iiωi²

    Angular momentum and energy are both conserved; in ω-space one is a sphere, the other an ellipsoid, and the motion must ride their intersection — the polhode of Poinsot’s rolling construction.

  • The tennis-racket theoremunstable
    I1 < I2 < I3 ⇒ axis 2 unstable

    Spin about the largest or smallest principal axis is stable; spin about the middle one is not — the body flips 180° over and over, the Dzhanibekov effect filmed on a wing-nut in orbit.

    1991
  • Holling type IIsaturation
    f(x) = a x1 + a h x

    A real predator can only catch and handle so much: at high prey density its kill rate saturates (handling time h). Swapping the cartoon βx for this curve is the step from the textbook cycle to a model that can settle — or crash.

  • Rosenzweig–MacArthurlogistic + saturation
    ẋ = rx(1 − x/K) − y f(x)  ·  ẏ = e·y f(x) − γy

    Give the prey a carrying capacity K (à la The Cascade’s logistic growth) and the predator a saturating appetite f(x), and the neutral cycles resolve — into a stable equilibrium, or a self-sustaining limit cycle born in a Hopf bifurcation.

    1963
  • Paradox of enrichmentmore food → crash
    K ↑  ⇒  Hopf  ⇒  extinction

    Raise the prey’s carrying capacity — enrich the system — and the stable point goes unstable: the populations swing wider and wider until a trough touches zero and someone goes extinct. Generosity destabilises. The same fragility The Cascade finds when a rate is pushed too hard.

    1971
  • Hot-potato routing, learned backwardsthe postman
    next(d) = arg mink live [hk(d) + α(t − τk(d))] · hk(o) ← a

    Pass every block at once to the live link whose table entry toward its destination is lowest. The table is learned backwards: a block from station o that arrives over link k after a hops says k reaches o in a. Entries age by α a round until refreshed, so routes that stop working are forgotten.

  • The toppling rulez_c = 4
    z(x) ≥ 4 ⟹ z(x) − 4, each neighbour + 1

    Drop one grain on a cell. The instant a cell holds four, it spills one grain to each of its four neighbours — which may themselves topple. One grain can set off a chain reaction: an avalanche.

  • The random graphno hubs
    P(k) = e−⟨k⟩ ⟨k⟩kk!

    Wire the same nodes at random and degrees follow a Poisson bell tight around the mean — almost every node near-average, hubs essentially impossible. The null model that makes the scale-free tail look as strange as it is.

    1959
  • Preferential attachmentrich get richer
    Π(ki) = kiΣj kj

    A new node links to an existing one with probability proportional to how many links it already has. Growth plus this one bias — the already-popular get more popular — is enough to breed hubs from nothing.

  • Scale-free degreesγ = 3
    P(k) ≈ 2m²k³ ∼ k−3

    The resulting degree distribution is a power law with no typical scale — most nodes have a few links, a rare few are giant hubs. The airline map, the web, the citation graph and the protein network all wear this tail.

    1999
  • Why hubs matter⟨k²⟩ → ∞
    scale-free ⇒ threshold → 0

    On scale-free networks the threshold can vanish entirely: with a diverging ⟨k²⟩ any transmissibility spreads. The cure is not more vaccine but smarter targeting — immunise the hubs and the giant component shatters.

  • The pudding’s nudgeobsolete
    θatom ≲ zZe²4πε₀ · E · R ≈ 10⁻⁴ rad · θfoil ≈ θatom√N

    A charge smeared through a sphere of radius R can only nudge a fast, heavy alpha — about a hundredth of a degree per atom — and a foil’s N random nudges add to a Gaussian cone a few degrees wide. Past 90° it allows, in effect, nothing.

  • Rutherford scattering1/sin⁴(θ/2)
    dσdΩ = (zZe² / 16πε₀E)2sin⁴(θ/2)

    A point charge deflects an alpha along a hyperbola, and the rate of large deflections falls as 1/sin⁴(θ/2) — rare but never zero. The counted fraction past 90° fixes how close the alpha came, and so how small the charge must be.

  • Energy levelsn only
    En = −13.6 eVn²

    Bohr’s ladder, re-derived exactly by Schrödinger. The energy depends on n alone — so 2s and 2p share a rung, a whole shell of shapes at one energy. That accidental degeneracy is the scaffold the periodic table hangs on.

  • Equation of motionSchrödinger
    i ∂ψ∂t = [ −½ ∇² + V(x,y) ] ψ

    The wavefunction’s whole future, set by its curvature and the potential V. It is complex — amplitude and phase — and only the pair together is the state.

  • Equation of motionSchrödinger
    i ∂ψ∂t = [ −½ ∂²/∂x² + V(x) ] ψ

    The wave function’s whole future, set by its curvature and the potential V. We solve it exactly in time by splitting it into a potential kick and a kinetic drift.

    1926
  • Hydrogen wavefunctionseparable
    ψnℓm = Rnℓ(r) · Yℓm(θ, φ)

    The bound state factors cleanly into a radial part Rnℓ — the shells, fixed by n and ℓ — and an angular part, the spherical harmonic Yℓm that draws the lobes. Square it and you have the cloud.

    1926
  • Angular momentumquantised
    |L| = √(ℓ(ℓ+1)) ℏ , Lz = m ℏ

    The length of L and one projection are sharp together; the other two never are. ℓ names the shape (s, p, d, f…), m the tilt — and in the complex basis, m is the visible circulation winding the phase.

  • The metricSchwarzschild
    ds2 = −(1 − rs/r)c2dt2 + … , rs = 2GM/c2

    Clocks deeper in the well tick slower (the −(1−rs/r) term). The event horizon is at rs = 2GM/c², where the escape speed reaches light.

  • The happiness ruleτ
    happy ⇔ nsamenoccupied ≥ τ

    An agent stays put if at least a fraction τ of its occupied neighbours share its colour, else it moves to a random empty cell. That is the whole rule — no agent wants segregation, each only dislikes being a small local minority.

  • Threshold modeltipping
    act ⇔ (fraction already acting) ≥ θi

    The same logic beyond space: each person joins in once enough others have. Whether a crowd tips depends on the distribution of thresholds, not the average — one stubborn holdout can stall a cascade a lower-threshold crowd would ignite.

    1978
  • Segregation indexorder param
    s = ⟨ nsame / noccupied ⟩

    The mean like-neighbour fraction across all agents — the order parameter. A random mix sits near 0.5; run the rule with τ as low as 0.3 and s climbs past 0.8, the grid coarsening into single-colour domains exactly like The Threshold’s magnet.

  • Herglotz–Wiechert inversionv(r) from T(Δ)
    ln Rr1 = 1π ∫0Δ1 cosh−1(p(Δ)/p1) dΔ

    The exact inverse of the travel-time curve: for the ray landing at Δ₁ with parameter p₁, the integral gives the radius r₁ where it turned, and v(r₁) = r₁ / p₁. It needs p to fall steadily with distance — a slow layer, which no ray turns in, breaks it, and that layer is in the null space.

  • The ray parameterSnell on a sphere
    p = r sin iv(r) = const

    In a spherically layered Earth a ray keeps r sin i / v all the way down and back — Snell’s law with the radius folded in. It turns where r / v(r) = p; and where the velocity drops, as under the mantle, no ray can turn at all: those depths leave a shadow at the surface.

  • The travel-time slopedT/dΔ = p
    dTdΔ = p

    The slope of the travel-time curve at a distance Δ is the ray parameter of the ray that lands there — the record section measures p directly. A shadow is a stretch of Δ where no ray lands; a kink is a layer where the velocity jumps.

  • Self-organized criticalityno knob
    slow drive + fast relaxation ⟹ ⟨z⟩ → z_c

    Add grains slowly, let avalanches run fast, and the pile climbs to the critical slope and stays there — balanced at the edge with no parameter tuned. The critical point is an attractor, not a setting. This is what The Threshold, The Percolation and The Cascade had to be dialled to by hand.

    1987
  • Semi-Lagrangian advectionunconditionally stable
    u(x, t+Δt) = u(x − u Δt, t)

    Stam’s insight: don’t step the velocity forward (which blows up), trace it backward — ask where each parcel came from and copy that value. This semi-Lagrangian advection is unconditionally stable, which is why it runs in real time in your browser.

  • Moore’s law×2 / 2 yr
    N(t) = N0 · 2(t − t₀)/T, T ≈ 2 yr

    Not a law of nature but a schedule: components per chip doubling every year (1965), then every two (1975). Each doubling is a shrink of about √2 in feature size.

  • Dennard scaling1/k
    L, V, C → 1k, f → k f, PA → PA

    Shrink every length and the voltage by 1/k: each transistor switches k times faster on k² less power, and k² more of them fit — so the power per area holds still. It held until the voltage could fall no further.

  • Switching powerC V² f
    P = α C V2 f, PA ∝ C V² fL²

    Every switch charges a capacitance C to V and dumps it: energy CV², f times a second. With V pinned and f stalled, shrinking L by k raises the power density by k² — the power wall.

  • The subthreshold swing60 mV/dec
    S = k Tq ln 10 ≈ 59.5 mV/decade (300 K)

    Below threshold a transistor’s current is a Boltzmann tail: it can rise no faster than tenfold per S of gate voltage. Turning it off by five decades costs ≈ 0.3 V of swing — a floor the voltage cannot pass.

  • Lambert shadingthe cosine law
    L = ρ · (Ea + Ed · max(0, n̂·l̂))

    A matte surface’s brightness: its paint ρ times ambient light plus the direct light scaled by the cosine between the surface normal and the light. Flip the surface’s tilt and the light’s side together and n̂·l̂ is unchanged — a bump and a dent send the same picture.

  • SIR dynamicswell-mixed
    Ṡ = −βSI   İ = βSI − γI   Ṙ = γI

    The classical well-mixed model: susceptibles become infected in proportion to contacts between S and I, and infecteds recover at rate γ. The whole epidemic is a flow S → I → R.

  • Site percolationp ∈ [0, 1]
    P(site open) = p

    Each site opens independently with probability p. Two open sites belong to the same cluster iff a path of open nearest-neighbours (up/down/left/right) connects them — the von Neumann neighbourhood.

  • The Condorcet cyclethe storm
    α ≻ β ≻ γ ≻ α

    Perfectly consistent individuals can build a society that is not: three voters ranking α>β>γ, β>γ>α, γ>α>β give a majority for α over β, for β over γ, and for γ over α — a preference that eats its own tail. There is no Condorcet winner; the circle was there before any rule.

  • Arrow’s theoremthe wall
    unanimity ∧ IIA ∧ ¬dictator ⇒ ⊥  (n ≥ 3)

    For three or more options, no rule that turns individual rankings into a social ranking can honour all three of unanimity (if all prefer x to y, so does society), independence of irrelevant alternatives (the x-vs-y verdict depends only on x-vs-y opinions), and non-dictatorship. Every rule surrenders one — or crowns a dictator. Not hard: impossible.

  • May’s theoremthe door
    n = 2 ⇒ majority is unique

    Shrink the choice to two candidates and the wall opens exactly: simple majority rule is the one and only method that is anonymous (voters interchangeable), neutral (candidates interchangeable) and monotone (more support never hurts). With two options fairness is not just possible but forced — the door Arrow’s wall leaves standing.

  • Imitate the bestcopy your winner
    si ← argmaxj ∈ N(i)∪{i} Πj

    Each cell plays its eight neighbours, sums its payoff, then adopts the strategy of whoever scored highest nearby (itself included). No memory, no foresight, no referee — and yet cooperation organizes in space.

  • How it integratesStrang split-step
    ψ ← e−iVΔt/2 ψ ; ψ̂ ← e−i|k|²Δt/2 ψ̂

    We solve it exactly per step by splitting: a potential kick in real space, a kinetic drift in Fourier space (where −½∇² is just ×|k|²/2). Unitary — |ψ|² is conserved.

  • The sRGB transfer curvecode → light
    L = ((V + 0.055) / 1.055)2.4 (V > 0.04045; else V / 12.92)

    A pixel value V is not a luminance. Every Iris stimulus is specified in linear light and encoded through this curve (IEC 61966-2-1), and ramps are dithered in linear light — so a ramp on the screen is a ramp of light, and two equal codes are two equal luminances.

  • Vyshnegradsky’s criterionXY > 1
    λ3 + Xλ2 + Yλ + 1 = 0 stable ⇔ XY > 1

    Scale the governed engine’s cubic until its outer coefficients are 1, and stability is the region above a hyperbola in the (X, Y) plane. In this model it reads b J ω₀ > 2 m l F: friction and inertia on one side, ball mass and load on the other.

  • Routh–Hurwitz for a cubica₁a₂ > a₀a₃
    a0s3 + a1s2 + a2s + a3 stable ⇔ aᵢ > 0, a1a2 > a0a3

    Every root of the cubic has a negative real part — every small disturbance dies away — exactly when the coefficients are positive and the middle pair’s product beats the outer pair’s. No root need ever be found.

  • Standing modesn = 1, 2, 3…
    yn = sin(nπx/L) cos(ωnt), ωn = nω1

    Only a whole number of half-waves fits, so the solutions are a harmonic series: mode n has n loops and frequency exactly n times the fundamental — quantisation straight out of geometry.

  • Stationary distributionπ = πP
    π = πP,  Σi πi = 1

    The rank is the equilibrium of a Markov chain — the fraction of eternity a blind random surfer spends at each node. Run one walker forever, or solve for the fixed point: the same vector. Kin to The Walk, set on a graph.

  • Plateau’s junction120°
    ê1 + ê2 + ê3 = 0

    Three films meeting along a line pull with equal tension, so their unit directions sum to zero: they meet at 120°. A shortest network between points has only such junctions — at most n − 2 of them.

  • The Saha equationionisation
    Ni+1 neNi = 2Zi+1Zi (2πmekT/h²)3/2 e−χi/kT

    How much of an element is ionised, from temperature and electron density alone. Hydrogen gives up its electron between about 8 000 and 12 000 K in a stellar photosphere — which is why its lines fade in the hottest stars.

  • Boltzmann excitationwho is up a level
    NbNa = gbga e−(Eb − Ea)/kT

    The fraction of atoms in an excited level. A Balmer line needs hydrogen in n = 2, 10.2 eV up: rare in a cool star, and — by Saha — ionised away in a hot one. The product of the two peaks near 10 000 K.

  • Stochastic resonanceoptimal noise
    SNRout peaks at σopt > 0

    A subthreshold signal — too weak to cross a detection threshold on its own — can become detectable by adding the right amount of noise. The output signal-to-noise ratio rises with noise from zero, peaks at an optimal amplitude, then falls again. Noise is not just a nuisance; here it is the carrier.

    1981
  • The Stokes equationsRe → 0
    ∇p = μ∇²u, ∇·u = 0

    Navier–Stokes with inertia struck out. No time appears: the flow is fixed by the boundary motion at this instant, linear in it, and reversed exactly when the boundary retraces its path.

  • Resistive-force theoryGray–Hancock
    f = −ζ∥(v·t)t − ζ⊥(v − (v·t)t), ζ⊥ ≈ 2ζ∥

    Each piece of a slender body feels local drag, twice as stiff sideways as lengthways. Demand zero total force and torque and the body’s velocity follows from its shape change alone — the engine behind the scallop, Purcell’s swimmer and the bacterial helix here.

  • The scallop theoremPurcell 1977
    Δx = ∮ G(θ) dθ = 0 for a stroke θ(t) with one degree of freedom

    Under Stokes flow the body’s velocity is linear in the shape’s rate of change, so displacement depends only on the path through shape space, never on speed. A one-hinge stroke goes out and back along the same path — and the displacement cancels exactly.

  • The look-then-leap rulestop optimally
    reject first r = ⌊N/e⌋, then take first > all seen

    Candidates arrive one at a time in random order; you must accept or reject each on the spot, and you want the single best. The rule: watch the first r ≈ N/e to calibrate, then grab the next one that beats everyone so far.

  • The win probabilitywin rate
    P(win | r) = rN Σi=r+1N 1i − 1

    The chance the look-then-leap rule lands the very best, as a function of the cutoff r. Too small and you leap before you have calibrated; too large and the best has likely already walked past. The maximum sits at r/N ≈ 1/e.

  • The 1/e lawthe 37% law
    rN → 1e, P(win) → 1e ≈ 0.368

    The striking part: as the field grows, both the optimal cutoff fraction and the probability of success converge to the same number, 1/e ≈ 37% — whether there are ten candidates or ten million. Skill cannot push a random order past this wall.

  • Bruss’ odds theoremone-line stop
    stop at the last index where Σ ri ≥ 1, ri = pi1 − pi

    A single rule for a whole family of stopping problems: sum the odds of the remaining “records” and stop the moment that running sum reaches one. The secretary’s 1/e is just this theorem applied to a random permutation.

  • Freeze & reversek = 1, 2, …
    fs = N ftrue/k ⇒ frozen

    When the strobe lands exactly on a spoke-pass the fold hits zero — frozen. Sample a touch faster and it tips negative: the wheel runs backward while it truly turns forward.

  • Hill’s tail indexthe largest k only
    α̂ = kΣi≤k ln X(i) − k ln X(k+1)

    For a tail P(|r| > x) ∝ x−α, the mean log-excess of the k largest moves over the next one estimates 1/α. For returns of real markets Cont (2001) puts it between two and five; a Gaussian sample gives a number that keeps climbing as the sample grows.

  • Fat tailskurtosis > 3
    κ = E[(r − r̄)⁴](E[(r − r̄)²])²,   κGauss = 3

    The fourth moment against the square of the second: 3 for a Gaussian, far more when large moves come too often. With it, the autocorrelation of r near zero and of |r| decaying slowly — calm and storm come in spells — are the stylised facts of real returns.

  • Subthreshold signalinvisible without noise
    x(t) = A sin(2π f0t) + σ ξ(t), A < θ

    A periodic signal of amplitude A sits below detection threshold θ — by itself it never fires the detector. Adding Gaussian white noise σ ξ(t) creates random barrier crossings; at the right σ, crossings cluster at the signal’s peaks and the output becomes correlated with the input. Biological sensory neurons, inner-ear hair cells, and crayfish mechanoreceptors all appear to operate near this optimal noise level.

  • Redundancy levelR = links ÷ stations
    R = LN (R = 1: the minimum-span net)

    Baran’s measure of how meshed a network is: the link-to-node ratio of a large array, with the fewest links that connect everyone (a tree) as level one. A square grid is R = 2, eight neighbours R = 4. Redundancy is what a strike has to cut through.

  • Survivabilityconnected pairs
    S = Σ ck(ck − 1)ns(ns − 1)

    After a strike, the share of surviving station pairs that can still reach each other: ck are the sizes of the surviving components, ns the survivors. Baran’s own criterion is its cousin — the share of survivors still in the largest group. Both are counted by breadth-first search, which no station can run.

  • Szilard engineW = kT ln 2
    Wmax = kBT ln 2 per measurement

    Leó Szilárd’s one-molecule heat engine: a demon who measures which half of a box the molecule is in can extract kBT ln 2 of work per cycle. But the demon must erase that single bit afterward — and that erasure, Landauer showed three decades later, costs exactly what was gained.

    1929
  • Takens' theoremm ≥ 2d+1
    Φτ : M ↪ ℝm  (m ≥ 2d+1)

    For almost any delay τ and any m more than twice the attractor’s dimension d, the delay map is a diffeomorphism — the reconstruction has the same topology and the same invariants (Lyapunov λ, fractal dimension) as the real thing. The geometry survives the shadow.

  • Galileo’s tideobsolete
    u = Vorbit + ω × r · h ∝ −d|u|dt = ω²(V · r⊥)/|u|

    A point of the Earth moves with the orbit plus the spin: faster by night, slower by day. The basin is sped up and slowed once a day, and the water piles where it slows — at dawn. One high a solar day, 24h00m, at every port.

  • The tidal force · two bulgesthe Moon’s gradient
    h(ψ) = 32mMRd3 R (cos²ψ − ⅓)

    The Moon pulls the near side more than the Earth’s centre and the far side less: the difference falls as 1/d³ and stretches the sea into two bulges, beneath the Moon and opposite it. The Earth turns each port through both on the Moon’s day — two highs, 12h25m apart. About half a metre at equilibrium; the Sun adds 46%.

  • The harmonic tidea harbor’s curve
    h(t) = Σk Ak cos(ωkt + Vk − gk)

    A port’s tide as a sum of constituents, each turning at a speed set by the Moon and the Sun — M2 (12h25m), S2 (12h00m), N2 — with amplitudes and lags fitted to the port’s own gauge. Kelvin’s machine summed them with pulleys.

  • The colour count1B + 1W
    domino = 1B + 1W ⇒ tiling needs B = W

    On a chequered board every domino, however laid, covers exactly one black and one white square — so any region that tiles must hold equal counts. Opposite corners share a colour; cutting both leaves 32 to 30, and no arrangement closes a gap of two.

  • Gomory’s circuitthe door
    cut 1B + 1W ⇒ ∃ tiling

    Thread all 64 squares onto one closed rook’s tour, colours alternating like beads. Remove one black and one white and the necklace falls into two arcs of even length — each tiles trivially along its own strand. Balance the colours and the wall opens, wherever the knife falls.

  • Tit-for-TatAxelrod's champion
    at = a′t−1,  a0 = C

    Open by cooperating, then echo your opponent’s last move. Nice, retaliatory, forgiving and clear — it won Axelrod’s tournament against every craftier rule. The spatial grid trades this memory for neighbourhood structure: clustering does the work repetition once did.

  • The Titius–Bode laddernumerology
    a = 0.4 + 0.3 · 2n AU · n = −∞, 0, 1, 2, …

    A rule, not a law: no force in it, only a pattern in seven numbers. It found Ceres in the gap at n = 3 and missed Neptune by 29% at n = 7 — a fit to a short series, read as if it were dynamics.

  • The Radon transformthe sinogram
    p(θ, t) = ∫ f(t cos θ − s sin θ, t sin θ + s cos θ) ds

    A beam at angle θ and offset t never sees the body — it returns one number, the integral of attenuation along its line. Stack those numbers over every θ and t and you have the sinogram; a single point of the body traces one sinusoid through it.

  • Kaczmarz row-action (ART)one ray at a time
    x ← x + λ pi − ⟨ai, x⟩‖ai‖² ai

    Treat every ray as one linear equation and correct the estimate onto each in turn. Renamed the algebraic reconstruction technique by Gordon, Bender & Herman in 1970; with positivity and a support as priors, it fills more of a missing wedge than FBP can — but never the frequencies the evidence did not sample.

  • Projection-slice theoremthe ring
    P̂θ(ω) = F̂(ω cos θ, ω sin θ)

    The 1-D Fourier transform of one projection is a radial line through the 2-D Fourier transform of the body. Every angle samples one spoke of Fourier space; the spokes you never measured — a missing wedge, the gaps past the Crowther radius — are the null space. The same theorem runs a CT scanner and an interferometer.

  • Filtered back-projectionramp |ω|
    f(x, y) = ∫0π (pθ ∗ h)(x cos θ + y sin θ) dθ, ĥ(ω) = |ω|

    Smear each projection back along its rays and sum: plain back-projection blurs every point into a 1/r halo, because polar Fourier samples crowd the centre. Filtering each row by |ω| first undoes that exactly (Bracewell & Riddle 1967; Ramachandran & Lakshminarayanan 1971). Shepp–Logan and Hann windows trade sharpness for noise.

  • The fundamental diagramflow–density
    q(ρ) = ρ · v(ρ)

    Flow (cars per minute past a point) is density times speed. As density ρ climbs, flow rises to a peak at the critical density ρc, then falls back to zero in gridlock — the inverted-U every road obeys, and the proof that more cars can mean less throughput.

    1935
  • The LWR equationa jam is a shock
    ∂ρ∂t + ∂q∂x = 0

    Cars are conserved, so density flows like a fluid — and where fast, sparse traffic meets a dense pack the solution forms a shock: the sharp upstream edge of a jam, a discontinuity that drifts backward through the cars at a speed set by the diagram, not by the drivers.

    1955
  • The Nagel–Schreckenberg rulesfour steps
    v ← min(v + 1, g, vmax); v ← v − 1 (prob p); x ← x + v

    Each car, each tick: speed up by one, but never faster than the gap g ahead or the limit vmax; then with probability p tap the brake for no reason; then move. That one random tap is the whole story — set p = 0 and traffic never jams on its own.

    1992
  • Transmissionexact (E < V₀)
    T = [ 1 + V₀² sinh²(κa)4E(V₀ − E) ]⁻¹

    The exact leak for one sharp energy. The sinh makes it fall exponentially with width a and with √(V₀−E) — why a hair more wall kills it.

  • Conservationunitary
    T + R = 1

    Probability is conserved: what doesn’t get through comes back. The instrument measures T by integrating |ψ|² past the wall — its own answer for a packet of many energies.

  • Why it patternsTuring instability
    Dv / Du > 1 ⟹ uniform state unstable → λ* selected

    Short-range activation, long-range inhibition: when the inhibitor diffuses faster (Dv ≫ Du), one wavelength λ* wins and tiles the plane.

  • Turing’s mechanism1952
    ∂u/∂t = Du∇²u + f(u,v) ; ∂v/∂t = Dv∇²v + g(u,v)

    Two morphogens, each diffusing and reacting. Turing’s shock: diffusion, usually a smoother, can destabilise a uniform state and manufacture pattern.

  • The correlation functionhow clustered
    ξ(r) = ⟨δ(x) δ(x + r)⟩

    The honest measure of clustering: how much more likely than chance you are to find a second clump a distance r from a first. It climbs from zero as the web sharpens — read live as the cloud collapses.

  • Two-slit intensityfar field
    I = I0 · sinc²(πa·u/λ) · cos²(πd·u/λ)

    A single-slit diffraction envelope times a two-slit interference fringe. The cos² term is the interference; kill it and only the smooth envelope is left.

  • Universal approximationCybenko 1989
    any continuous f, to any ε

    Why one hidden layer is, in principle, enough: a wide enough network of these neurons can approximate any continuous function. Capacity is the catch — Spiral needs more units to bend the boundary far enough.

  • Urban scalingY = Y₀ N^β
    Y = Y0 Nβ, β ≈ 1.15 (output) · β ≈ 0.85 (roads)

    A city of twice the people makes about 21.15 ≈ 2.2 times the wages, patents and crime, on about 20.85 ≈ 1.8 times the road — more per head and less per head, from one doubling.

  • The mixing exponents1 ± 1/6
    A ∝ Nα, roads ∝ N(1+α)/2, contacts ∝ N(3−α)/2; α = 2/3 → 56, 76

    Settled area grows as Nα; roads fill it as √(area·nodes); everyone reached within a day’s travel along those roads gives the contacts. While a day crosses the whole town, everyone meets everyone and contacts grow as N².

  • Vicsek updateminimal model
    θit+1 = ⟨θj⟩Ni + η·ξ

    Replace steering with "copy the average heading of your neighbours, plus a kick of noise η". This minimal model is the one with an exact order parameter and a phase transition.

  • The Voronoi cellnearest seed
    V(pi) = { x : ‖x − pi‖ ≤ ‖x − pj‖ ∀ j }

    Every point in the plane belongs to its nearest seed. That one rule tiles space into cells whose walls are the perpendicular bisectors between seeds — the partition nature reaches for whenever things compete for the space around them.

  • The Delaunay dualthe mesh
    Del(P) = dual of Vor(P) · empty-circumcircle

    Join the seeds of every pair of touching cells and you get the Delaunay triangulation — the Voronoi diagram’s twin. It is the “best” mesh on a set of points: no point sits inside any triangle’s circumcircle, so it avoids thin slivers.

  • Lloyd relaxationmove to centroid
    pi ← centroid( V(pi) ), repeat

    Let each seed hop to the centre of mass of its own cell, recompute the cells, and repeat. A random scatter relaxes into an even honeycomb — the identical assign-then-average loop as k-means clustering.

  • Quantization energyminimised
    E = Σi ∫V(pi) ‖x − pi‖² dx

    The total squared distance from every point to its seed. Lloyd’s step is gradient descent on this energy, and its minimum — a centroidal Voronoi tessellation — is the most uniform scatter there is: blue noise, no clumps and no holes.

  • Conservation of circulationthe true kernel
    Γ = ∮ v · dl · dΓdt = 0

    In an ideal fluid the circulation round a vortex never changes and its lines move with the fluid (Helmholtz 1858; Kelvin’s theorem 1869) — so a knotted vortex ring can never untie. True, and the seed of Kelvin’s atoms.

  • A vortex ring’s speedthin core
    U = Γ4πR (ln 8Ra − 14)

    How fast a smoke ring of radius R and core a carries itself forward — Kelvin’s formula, the motion the vortex atom was built on. It sets a speed from free numbers; nothing in it is a wavelength.

  • The exposure tradepopulation is time
    ⟨events⟩ = ε·N·Tτ  ·  τ > ε·N·T2.303 (none seen, 90%)

    A crowd of N memoryless things, each with mean life τ, watched for a time T ≪ τ, loses about N·T/τ of them — so N protons for one year are worth one proton for N years. You cannot wait 10³⁴ years; you can watch 10³⁴ protons. And when nothing happens, the nothing is a measurement: zero events bound τ above ε·N·T/2.303. Every waitclock in the edition prints this exchange rate; here physics signs it.

  • Exponential survival — memorylessnessthe spine law
    P(T > t) = e−t/τ  ·  P(T > s + t | T > s) = P(T > t)

    A process with a constant hazard survives a span t with probability e−t/τ — and having already waited s changes nothing about the wait ahead. The bus you just missed never existed; a queue’s arrivals and services forget; a metastable state has a half-life and no age. The one law under every Zeno wait.

  • Wave equationclamped ends
    ∂²y/∂t² = c² ∂²y/∂x²

    A clamped string obeys the wave equation with both ends pinned. Those boundary conditions are the whole story — they force the continuous string into a discrete set of allowed shapes.

    1747
  • Wave–particle dualitythe law
    V2 + D2 ≤ 1

    The exact trade. V is fringe contrast (wave), D is how well the path is known (particle). You can buy one only by selling the other — never both at once.

    1995
  • Time–frequency limitthe wall
    Δt · Δf ≥ 14π

    No signal is sharp in both time and frequency: their spreads multiply to a fixed floor. The Gaussian alone meets it — the classical twin of Heisenberg’s Δx·Δp ≥ ℏ/2.

  • Continuous wavelet transformtime × scale
    W(a,b) = 1√a ∫ x(t) ψ*(t − ba) dt

    Correlate the signal against one little wave ψ, scaled by a (≈ 1/frequency) and slid to time b. |W|² is the scalogram — energy laid across time and frequency at once.

    1984
  • The Morlet waveleta ripple in a bell
    ψ(t) = π−¼ eiω₀t e−t²/2

    A complex sinusoid wrapped in a Gaussian. The width ω₀ is its wiggle-count: more of them sharpens frequency and blurs time, fewer does the reverse — the one knob on the trade.

  • Admissibilityzero mean
    Cψ = ∫ |Ψ(ω)|²|ω| dω < ∞ ⇒ ∫ ψ dt = 0

    For the transform to invert, the wavelet must average to zero — a true little wave, not a lump. That’s why ω₀ can’t fall too low: too few wiggles and the mean stops vanishing.

  • Why it spreadsdispersion
    σ(t) = σ0 √(1 + (t/2σ0²)²)

    A free packet doesn’t hold together — the faster components outrun the slow, so σ grows with time. The spread you watch is the price of ever having localised it.

  • Pareto’s tailpower law
    P(m > x) ∝ x−α,  α ≈ 1–2

    The richest tail is not exponential but a power law: a straight line on log–log axes, no characteristic scale. The famous 80/20 — roughly a fifth of the agents hold four-fifths of the wealth — and the same fat tail The Sandpile and The Web draw in other keys.

  • Money as energyexponential
    P(m) ∝ e−m/T,  T = ⟨m⟩

    Conserve a fixed pot of money and let agents swap it at random, and wealth settles into the same exponential law that energy takes among gas molecules — the money temperature T is just the average balance. Fairness at the microscale, a heavy skew at the macroscale.

    2000
  • The transactionconserved
    mi′ = mi + Δ,   mj′ = mj − Δ

    Every trade only moves money between two agents — the total is exactly conserved, like energy. A fair Δ relaxes to the exponential; a biased Δ (the loser always pays a slice of the poorer purse) condenses nearly all of it onto one agent — Bouchaud–Mézard oligarchy.

  • The Gini coefficient0…1
    G = Σij |mi − mj|2N²⟨m⟩

    One number for the whole distribution: 0 is perfect equality, 1 is one agent owning everything. A fair-exchange economy sits at exactly G = ½; the biased rule drives it toward 1 as the money condenses.

  • The digit enginex ↦ {1/x}
    a = ⌊1/x⌋ , x ↦ 1x − a

    Flip the fraction, keep the fractional part, write down the whole number that fell out — repeat, and any ω dictates its continued fraction digit by digit. Euclid’s algorithm wearing dynamical clothes; the golden number is the fixed point of the first branch, x = 1/(1+x), so its tape reads 1, 1, 1, … forever.

  • The floor√5 is sharp
    |ω − p/q| < 1√5 · q²

    Every irrational is caught this close by infinitely many fractions — and no constant larger than √5 works, because one number sits exactly on the floor: for φ the normalized error q²·|φ − p/q| converges to 1/√5 ≈ 0.4472 and never below. The worst-served number in arithmetic, named by its distance from every fraction.

  • The armor conditionDiophantine
    |ω − p/q| ≥ γqτ ⇒ the torus survives

    KAM’s survival clause is an arithmetic condition: a torus outlives a small kick exactly when its winding number keeps this far from every fraction. Being badly approximable — the crime Hurwitz convicted φ of — is here the armor. Same property, opposite loyalty.

  • The kicked windingarea-preserving
    p ← p + k2π sin 2πx , x ← x + p

    The winding law with one kick bolted on. At k = 0 every horizontal line p = ω is an invariant circle — the pure winding, stacked. Turn k up and the kick feeds on rational windings first: each p/q grows an island chain, the chains grow a chaotic sea, and the only circles left standing are the ones the kick cannot resonate with.

  • The residue dialR → ¼
    R = 2 − Tr M4 , R(Fn/Fn+1) → ¼ at kc

    Greene’s dial: linearize the map around the periodic orbits whose rotation numbers are the golden winding’s own convergents — the Fibonacci fractions besieging it — and read the residue R. Falling toward 0: the circle holds. Rising past ¼ up the whole ladder: it is gone. The ladder pinches flat at kc ≈ 0.971635.

  • The windingthe whole law
    x ← x + ω (mod 1)

    Turn by the same twist ω, forever — the pair’s entire law. Arithmetic reads its near-returns as fractions; dynamics stacks one such circle for every ω and asks which survive a shake. The law itself never says what the number is for.

  • The ladderpₙ/qₙ
    pn = anpn−1 + pn−2 , qn = anqn−1 + qn−2

    Each digit buys one rung: the convergents pₙ/qₙ, the only fractions that beat every smaller denominator. Big digits are windfalls — π’s a₄ = 292 turns 355/113 into seven digits of π. All-ones digits are the slowest ladder arithmetic allows: the golden ratio’s rungs are the Fibonacci numbers.

  • Number = behaviourWolfram code
    sit+1 = ⌊R / 2k⌋ mod 2

    Index the neighbourhood by its bits, k = 4si−1 + 2si + si+1, and the rule number R is literally the truth table — bit k of R is the next state.

  • Zel'dovich's approximationthe engine
    x(q, t) = q − D(t) ∇ψ(q)

    Instead of summing forces between millions of particles, give each one a fixed shove from the initial field and slide it that way, scaled by the growth D(t). Crude once paths cross — but it draws the whole cosmic web of sheets, filaments and knots almost for free. This instrument runs exactly this.

    1970
  • The quantum Zeno effectthe watched pot
    Pn(t) = [P(t/n)]n ≈ [1 − (t/nτZ)2]n → 1

    Every quantum decay starts flat: the survival is 1 − (t/τ_Z)² before it is ever e−t/τ. Look n times in a time t and each look catches the state still at the top of the parabola — the product of n near-ones climbs to one as n grows. Asked often enough, it never decays.

  • The decay rate under looksZeno or anti-Zeno
    Γeff(τ) = −ln P(τ)/τ = 2π∫G(ω)·Fτ(ω − ω0) dω

    Looks every τ broaden the state’s line to a width ~1/τ (F_τ), and the decay rate is the reservoir’s spectrum G(ω) seen through that window. A very short τ smears it thin — Zeno, Γ_eff → 0; a τ that slides the window onto a peak of G lifts Γ_eff above the unwatched rate — anti-Zeno. Which one you get is the environment’s choice.

Locator — Scan the archive
esc
Indexing archive…

Indexing the archive…

↑↓ select↵ openesc close