Every Law on the Rack
616 lawsConstellation Gallery- The abelian propertyDhar 1990order-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 chartistsrt+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 illusionfperc = 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 1638A ∝ 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/2d² ∝ 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/3B ∝ 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/4B = 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/4f = 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 modelPreqPavail ∝ 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/2d ∝ 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 1989EMA ∝ 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/3W/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 sizey = 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/6V = (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 dxdz = (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″ = −ω²ydv = −ω 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 routexr > 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 latenciesPoA = 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.
- One choice: the fullest binGonnet 1981max 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 nmax 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.
- The supermarket taildoubly exponentialsi = λ(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 learningRT = 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 boundoptimismat = 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 payfail ⇔ ℓ > ψ(θ) = 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 Bayesposteriorp(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 averagex* = 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 parameterN0/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 trapkBTc = ℏω̄ (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 condensesdegeneracyn λ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).
- 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 localS = | 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√2S ≤ 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 curveE(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 pointV*(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.
1957 - Q-learning updatetemporal differenceQ(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 returnQ(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 constantintegrableL₁ · 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_sbcrit = 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 §11P = 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)^zqz = 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 20131 − γ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 boidsv̂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.
- The Boltzmann factorheat to escapep(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 ruleP(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 · BachelierD = 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 columnPcr = π² 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/3hcr = (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 1989O = 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. 2000Var(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.
- The Byzantine boundn ≥ 3f + 1n ≥ 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 · IV30 → 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 inputssit+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 caloricobsoleteW = 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/rh = 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 ÷ tensionBo = ρ 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) ≤ 3NMΠ(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.
- 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 capacityBSCC = 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.
1948 - Choreographiescentral configsthree 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 croppt = 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(τ) levelsf(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 setsE[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 tailGumbelP(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 collectorharmonicE[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 lostv′ = (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 clacksN = ⌊π / θ⌋ , 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. 1964yD = −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 fitmin ‖Ψᵀ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.
- How many measurementsnot NM ≳ 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 shrinkc ← 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 problemundecidablehalt(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-matchinga ▸ 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.
1961 - The busy beaveruncomputableS(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.
1962 - 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.
1966 - Langton’s anttwo ruleswhite → 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 shortcutfuture ⟸ 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-entropycollapse 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/3P(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 decreaseAIMDon 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 pB ≈ 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 = πwP = ∮ ρ 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 areaA ≥ ½ (π − √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 harmonicsh(θ) + 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 areaAswept = (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 quenchTk = α·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 parameterfC = 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 2P(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 divisionX̂ = 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, positivex ← 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 + countsy = 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 spaceOTF(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_pMCh ≈ α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/2a0N1/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/2Mmin ≈ α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 signalv(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.
- 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-referenceG ⇔ ¬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 flipD(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 processnoisingdx = −½ β(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.
- The reverse processAnderson 1982dx = [ −½β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.
- The score∇ log ps(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.
- Probability-flow ODEdeterministic twindxdt = −½β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.
- The Clark–Evans ratiochance: R = 1R = 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.
- 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 = 1VR(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₂ nP(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 friendsCi = 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 = 1H3 = −Σπ 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 ≥ 6echo = 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 1989DE = ½ · 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.
1989 - Diffusion-limited aggregationstick on touchwalker 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.
- The screening instabilitytips winexposed 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 = dragm 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 ≪ 1F = 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–NaumannCd = 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 airPa 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 playK ~ 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.
INST·140◆Vega Edition The DiceINST·141◆Vega Edition The CallINST·142◆Vega Edition The HostINST·143◆Vega Edition The EnvelopeINST·144◆Vega Edition The MartingaleINST·145◆Vega Edition The TicketINST·146◆Vega Edition The AuctionINST·147◆Vega Edition The LineupINST·148◆Vega Edition The DetourINST·149◆Vega Edition The TruelINST·170◆Vera Edition The WireINST·171◆Vera Edition The BooksINST·172◆Vera Edition The KeyboardINST·173◆Vera Edition The CanvasINST·174◆Vera Edition The DriftINST·175◆Vera Edition The PulseINST·176◆Vera Edition The CrowdINST·177◆Vera Edition The ScatterINST·178◆Vera Edition The MelodyINST·179◆Vera Edition The Hand - Field equationsEinstein 1915Gμν = 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 < Lgoᵢ ⇔ Â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 bitH(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.
1948 - Landauer limitkT ln 2Emin = 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 ↔ BoltzmannH = −Σ 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 / ShannonS = −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 equantobsoleteC = 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ε-greedya = 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 lineWn+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 backwardraise 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 roomeiπ + 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 2V − 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 − 44F ≤ 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χ = 0V − 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 wv̇ = 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.)
- 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 slowc(κ) = 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 errorpet = 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 ≤ 1J = (Σ 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-fillingxi = 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 ∈ MJc 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 sampleffeat = 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.
- The Cascade, liftedδ ≈ 4.6692c ∈ [−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 = ±2M − 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 flightW = ½ ρ 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 onlyNi = { 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 claimf(θ) = 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 outXk = Σ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 timeX(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 modesFourieran = (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 dataC(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.
- Fractal dimensionD ≈ 1.71N(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 ∝ fitnessP(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 1975m(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.
- Crossoverrecombinationchild = 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 variationgene ← 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 1G ▷ 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 · GamowT ≈ 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 MotiondSdS = μ 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.
- 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_tSt = 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 gainGDOP = √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 fractionE[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.
- 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-stochasticG = 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.203s = 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 nT(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).
- 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 equationsJ ω̇ = 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/ghmax ≈ ξ 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 ↔ boundaryA = ½ ∮ (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 insideA = 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 iterationamplifyG = (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 speedupk* ≈ π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 ≤ 0H = ∫ 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ₘᵢₙ = 34 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.
1950 - The Tower of Hanoi countexactM(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 pegsFS(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 arcring: 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.
- Virtual nodes even the loadmax ÷ meanmax/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₂ Sfinger 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.
1787 - Weyl's lawhear the areaN(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.
1911 - 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_cpc = 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-inheritanceobsoletegchild = ½(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 ↛ germgchild = ½(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 + emissionsP(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 = 0u = 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 itI = |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 FFTU(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 diffractionparaxialU(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= 1lk = 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 circlez1 = 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 memoryHebbwij = 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 energyLyapunovE = −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.138pmax ≈ 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.
- Recalldescentsi ← 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 gapCon(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-auditT ⊬ 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 signala = 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).
- The odds of a right cascadeBHW, eq. 3P(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.
- The van Cittert–Zernike theorema visibility is a Fourier coefficientV(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 beamIres ← 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 circleexactOP · 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 pointexactkTc/J = 2ln(1 + √2) ≈ 2.269
No fitting, no simulation — Onsager solved it on paper. Above it, noise wins; below it, order.
- Lattice energyHamiltonianH = −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 samplingMetropolisP(σ) ∝ 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/8m = [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 votesPN(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 gainKK = 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.
1960 - 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, Sy = 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, Px = 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–YorkedKY = 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 pulseu(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).
- 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.
- 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.
- Maximise expected log-wealthgrowth-optimalmaxf 𝔼[ 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 stakef* = 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-growthg(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 assetcontinuousf* = μ − 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 cosinesf(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 − θ_AOABC ∼ 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 lawr = 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 relationt−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 integralWr = 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-repellingE(γ) = ∮∮ (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 activatedrK ∝ 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_cKc = 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|ω| ≤ Krlocked ⟺ |ω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 culpritsplanar ⇔ 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’s lawrefractionn1 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’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.
1662 - The brachistochronea cycloidx = 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 ⇒ thisddt∂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-lawS = ∫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?
- The path integralFeynman 1948K = ∫ 𝒟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 minimisedOPL = ∫ 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 thisT = ∫ 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 − potentialL = 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 updatecontinuousAt+Δ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 ⇄ deathG(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 ∗ AU(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 α-stableVar = ∞ ⇒ 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/S23born 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 seestotalisticni = Σ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 paintTR = 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 whitelog ℓ = 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·EL(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 closurer1eiθ₁ + 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 wildLiouvilleE = 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 1978a → 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 linexn+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 transformationboostt′ = γ(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 cpostulatec′ = 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 conetimelike 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 pathF = 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 orbitgyrationrL = 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-blindvd = 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.74C± = (±√(β(ρ−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 ODEsdxdt = σ(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 sinkx* = γδ · 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 orbitsV = δ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 rulerewritingF → 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 → formF: 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 atlasM = { 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 3Dz → 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₂ nE = Σ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 driftE[XT] = E[X0] ⇒ E[P<] = −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 delayA(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 holdingK = 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 uphillPaccept = 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.
- 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 = 0H = ½(κ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 uniformf(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 lossLL = ½ ⟨(ŷ − 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.
- The Reichardt detectordelay and multiplyR(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 gainui(σ*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 neuronforwarda = 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 2zn+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 borderB(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 ≤ 1zn+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 squareai = Σ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 → 0R < 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.
1948 - The dominance cycleno maximumA ≻ 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 oddsodds(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 geodesicu″ + 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 wP(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.
- Nyquist limitor it aliasesfs > 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 noiseflickerS(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 freedomS = 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 1998r(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±1sgn(σ) = (−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 movesI = 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 losingA: 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 chainexactg = π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 truthexactT = 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 rest2 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/91ns ~ 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.5927P∞(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 costx′ = 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.
- 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 εβ/γ.
- 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.
- Period-doubling cascader∞ ≈ 3.5699r₁=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, positiveG 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 XIl(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, projectg′ = 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 escapesg ← 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 evidenceI(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 balanceobsoletemetal = 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_srph = 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 sensorssample 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 ruledDijdt = 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 fieldT ← (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 octaves(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^SP(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 planesR = 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 0V − 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 KuramotoM θ̈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.
- Droop controlprimary controlPi ← 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 rrk+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.2P(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 estimateLi(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.
- The Prisoner's DilemmaT > R > P > ST > 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 conequantileS0 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 currentj = |ψ|² ∇φ = 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 stimulusS = 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 spotd = 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.
- 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 Misesp(θ) ∝ 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 detectiond′ = 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 pointg : 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.
- The bandwagon responseb > 0 bandwagon · b < 0 underdogR(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 gatesuperposeH : |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 gateentangleCNOT : |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₀ = 0zn+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 ρ = 1W = 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 queueL = λ·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 × timeE[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 formx5 + 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 2002R₀ = 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 ovenu = 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 continuumBλ(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²⟩ = 2DtSn = Σ 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–Ratliffrp = 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 − surroundB = (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 slows1 + 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.
- The inspection paradoxlength-biasedE[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 frozenxt+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.
2001 - The linear readoutthe only trained partyt = 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 formWout = 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 edgeslog ℓ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 1876v → −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.
- 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-freeI1ω̇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 theoremunstableI1 < 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.
- Holling type IIsaturationf(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.
1959 - 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 → crashK ↑ ⇒ 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 postmannext(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 = 4z(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 hubsP(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γ = 3P(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 onlyEn = −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ödingeri ∂ψ∂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ödingeri ∂ψ∂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.
- 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.
- 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 metricSchwarzschildds2 = −(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 modeltippingact ⇔ (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.
- Segregation indexorder params = ⟨ 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 spherep = 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Δ = pdTdΔ = 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 knobslow 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.
- Semi-Lagrangian advectionunconditionally stableu(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 yrN(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/kL, 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² fP = α 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/decS = 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 lawL = ρ · (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 wallunanimity ∧ 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 doorn = 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 winnersi ← 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 → lightL = ((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.
1906 - 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 equationionisationNi+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 levelNbNa = 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 noiseSNRout 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.
- 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–Hancockf = −ζ∥(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 optimallyreject 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 rateP(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% lawrN → 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 stopstop 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 noisex(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 ÷ stationsR = 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 pairsS = Σ 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 2Wmax = 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.
- 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 tideobsoleteu = 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 gradienth(ψ) = 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 curveh(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 + 1Wdomino = 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 doorcut 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 championat = 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 laddernumerologya = 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 sinogramp(θ, 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 timex ← 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 ringP̂θ(ω) = 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–densityq(ρ) = ρ · 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 stepsv ← 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.
- ConservationunitaryT + 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 instabilityDv / 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.
1952 - 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 fieldI = 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 1989any 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/6A ∝ 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 seedV(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 meshDel(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 centroidpi ← 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 energyminimisedE = Σ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 coreU = Γ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 lawP(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.
- Wave–particle dualitythe lawV2 + 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.
1946 - Continuous wavelet transformtime × scaleW(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.
1984 - Admissibilityzero meanCψ = ∫ |Ψ(ω)|²|ω| 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 lawP(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 energyexponentialP(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.
- The transactionconservedmi′ = 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…1G = Σ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-preservingp ← 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 lawx ← 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 codesit+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.
1983 - Zel'dovich's approximationthe enginex(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 potPn(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.