Rack // Equations

Every Law on the Rack

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

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

  • The alias foldthe illusion
    fperc = ffeatfsffeat/fsN

    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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    1925
  • One wavefunctionmacroscopic ψ
    ψ(r) = √n(r) e

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

  • The mirror bouncethe whole rulebook
    θin = θout

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

  • The hidden constantintegrable
    L₁ · L₂ = const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

  • Central Limit Theorem→ 𝒩(0,1)
    (Sn) / σn → 𝒩(0,1)

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

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

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

  • Choreographiescentral configs
    three masses, one shared ∞-curve

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

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

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

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

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

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

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

  • The 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 universalmachine reads another machine’s table off the tape and runs it: the first proof that one device can do anything any device can.

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

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

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

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

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

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

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

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

  • Langton’s anttwo rules
    white → turn right · black → turn left · flip, step

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

    1986
  • Computational irreducibilityno shortcut
    future ⟸ run every step

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

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

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

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

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

  • 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.

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

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

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

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

  • 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 longermatter waves. Cool far enough and each atom’s wave grows until it overlaps its neighbours — the gas stops behaving like particles.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    1981
  • The screening instabilitytips win
    exposed tips ≫ shielded fjords

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

  • Field equationsEinstein 1915
    Gμν = 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
  • Boltzmann entropyS = k ln Ω
    S = kB ln Ω

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

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

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

  • Landauer limitkT ln 2
    Emin = kBT ln 2

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

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

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

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

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

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

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

  • Explore vs. exploitε-greedy
    a = argmaxaQ(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 geometricmean, which is always smaller. When ⟨r⟩ > 1 > e⟨ln r⟩ the same coin enriches the ensemble and ruins almost every member of it — both verdicts exact, both forever.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

  • Feature ratewhat you sample
    ffeat = N ftrue

    An N-spoke wheel looks identical every 1/Nturn, 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 samenumber for any smooth one-hump map. The route into chaos is universal; it doesn’t care what the map is.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

  • Fringe spacingΔy
    Δy = λLd

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

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

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

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

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

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

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

  • Mutationblind variation
    gene ← random, with probability μ

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

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

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

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

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

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

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

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

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

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

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

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

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

  • 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.

  • Gray–Scottthis engine
    A/∂t = DA∇²AAB² + 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
  • Grover's iterationamplify
    G = (2|s⟩⟨s| − I) · Ow

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

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

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

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

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

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

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

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

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

  • 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.

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

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

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

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

  • Weyl's lawhear the area
    N(k) ≈ A k² − L k

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

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

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

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

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

  • Herd immunityp_c
    pc = 1 − 1 / R₀

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

  • Diverging scalesuniversality
    ξ ∼ |TTc|−ν, χ ∼ |TTc|−γ ; ν = 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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

  • WritheGauss integral
    Wr = 1 ∮∮ (dr × dr′)·(rr′)|rr′|³

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

  • The Lagrangiankinetic − potential
    L = T − V

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

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

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

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

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

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

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

  • Lens equationEinstein ring
    β = θθE2/θ

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

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

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

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

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

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

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

  • Light deflectionmeasured 1919
    α = 4GMc2b

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

    1919
  • 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
  • Conserved, yet wildLiouville
    E = T + V = const ; dΓ/dt = 0

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

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

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

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

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

    1905
  • Invariant interval
    = −(ct)² + = −(ct′)² + x′²

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

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

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

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

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

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

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

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

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

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

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

  • Cyclotron orbitgyration
    rL = mvqB , ωc = qBm

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    1980
  • The Mandelbulbz → zⁿ + c in 3D
    zrn(sin cos , sin sin , cos ) + 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
  • 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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

  • Nyquist limitor it aliases
    fs > 2 ffeat

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

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

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

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

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

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

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

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

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

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

    A legal slide is exactly one transposition — the tile trades cells with the blank — and exactly one unit step of the blank. Two coins, flipped together on every move, with no third option on the board.

  • The two losing gamesboth losing
    A: p = ½ − ε  ·  B: p = 110 − ε if 3 | W, else 34 − ε

    Game A is a plain biased coin; game B picks its coin by whether your capital is divisible by 3 — a vicious coin on multiples of three, a generous one elsewhere. Played alone, each drifts your capital down. The law never mentions alternation.

  • The verdict of the chainexact
    g = π0(2pbad − 1) + (1 − π0)(2pgood − 1)

    The drift per round, exactly: π₀ is how often the capital sits on a multiple of 3 — the stationary state of a three-state Markov chain. Game B digs its own trap (π₀ high, the bad coin often); mixing in A shuffles π₀ down, and the same formula changes sign.

  • The small-angle liesin θ ≈ θ
    T0 = 2π √(l / g)

    Let sin θ ≈ θ and the amplitude vanishes from the period — every swing, wide or narrow, takes the same time. This is Galileo’s isochrony, the founding approximation of timekeeping, and it is only true in the limit of swings too small to see.

  • Equation of motion · θ₁″exact
    θ1″ = g(2m1+m2)sinθ1m2g sin(θ1−2θ2) − 2sinΔ·m222l212l1cosΔ)l1(2m1+m2m2cos2Δ)

    Straight from the Lagrangian L = TV. 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+m2m2cos2Δ)

    The two arms are nonlinearly coupled — the source of the unpredictability. We integrate both with RK4 so energy doesn’t drift.

  • The pendulum lawone line
    θ″ = −gl sin θ

    A mass on a pivot under gravity — the most familiar law on Earth, and the whole of it. One line, no forcing, no friction, no dice. The law says nothing about how many joints the mass hangs from; that silence is where the two tempers come from.

  • The whole truthexact
    T = 4 √(l / g) · K(sin (θ0 / 2))

    The exact period, with K the complete elliptic integral — computable to machine precision by the arithmetic–geometric mean. The amplitude never left; it was hiding in K. At θ₀ = 90° the clock runs 18% slow, and as θ₀ → 180° the period grows without bound.

  • The forbidden linereleased at rest
    2 cos θ1 + cos θ2 = 1

    Energy alone decides who may flip: released at rest inside this line, a double pendulum (equal masses, equal arms) can never send its second arm over the top — however long you wait. Outside it, the flip is allowed, and the fractal decides how long you wait.

  • Cluster-size distribution at p_cτ = 187/91
    ns ~ s−τ

    At p_c the cluster-size distribution is a power law — there is no characteristic scale, clusters of all sizes coexist. This scale-free structure is why the critical point looks fractal and why the exponents are the same across wildly different physical realisations.

  • Correlation lengthν = 4/3
    ξ ~ |p − pc|−ν

    The correlation length — the typical size of the largest finite cluster — diverges as p approaches p_c from either side, meaning the system "sees" arbitrarily large distances just at the critical point.

  • Phase transitionp_c ≈ 0.5927
    P∞(p) ~ (p − pc)β, β = 5/36

    Below p_c only finite clusters exist. At p_c a giant cluster of infinite extent first appears (in the thermodynamic limit). Above p_c the giant component density P∞ grows as a power law — the order parameter of the transition.

  • Universalityd = 2 universality class
    β, ν independent of lattice

    The exponents β, ν, τ are universal — they do not depend on the lattice geometry (square, triangular, honeycomb) or the physical realisation. Only the spatial dimension d and the symmetry of the order parameter matter. That universality is why the same numbers describe water in rock, fire in a forest, and contagion in a network.

  • Period-doubling cascader∞ ≈ 3.5699
    r₁=3 · r₂≈3.449 · r₃≈3.544 · r₄≈3.564 · … → r

    As r climbs, the stable cycle doubles: 1 → 2 → 4 → 8 → … The thresholds crowd together and accumulate at r∞ ≈ 3.5699, where the period becomes infinite — the door into chaos.

  • Perron–Frobeniusunique, positive
    G v = λ₁ v,  λ₁ = 1,  v > 0

    A positive matrix has one largest eigenvalue, and its eigenvector is all-positive and unique. That theorem is the guarantee under PageRank: exactly one rank vector, every entry a real fraction, and power iteration always finds it.

    1912
  • 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 photon spherer = 1.5 r_s
    rph = 3GM/c2 = 1.5 rs

    At exactly this radius gravity bends light into a circle. Rays that graze it wind around the hole many times before escaping, stacking into the thin bright photon ring at the rim of the shadow.

  • Trail sensingthree sensors
    sample FL, FC, FR → turn ±α toward the strongest

    Every agent reads the chemoattractant at three points ahead — left, centre, right — and rotates by the sensor angle α toward whichever is brightest, else holds course. No map, no goal: just follow the smell of where the swarm has already been.

  • Adaptive conductivitythe Tokyo rule
    dDijdt = f(|Qij|) − Dij

    Tero’s distillation: a tube’s conductivity D grows with the flux Q it carries and decays otherwise. Heavily-used routes thicken, idle ones wither — a feedback that prunes a tangle down to an efficient, fault-tolerant network with nobody designing it.

  • Deposit · diffuse · decaythe field
    T ← (1 − ρ)·(K ∗ T) + d·[agent here]

    The trail map T is the only shared memory: each agent drops a little d where it stands, the field blurs (K ∗) so neighbours can smell it, and a fraction ρ evaporates each step. Reinforcement plus evaporation is what lets unused paths fade and used ones thicken.

  • Poincaré recurrenceτ ~ e^S
    P(return) ≈ 1 : 2N

    The gas willreturn arbitrarily close to its start — after a time that dwarfs the age of the universe. "Never," in practice, is a statement about how long you’d wait.

  • The Schläfli symbolrecipe
    {p, q, r}

    The whole recipe for a regular 4-polytope in three numbers: {p} faces, {q} around each vertex, {r} cells around each edge. The tesseract is {4,3,3} — cubes ({4,3}), three to an edge.

  • The 4-D double rotationtwo planes
    R = RXW(α) ∘ RYZ(β)

    A generic rotation in four dimensions turns in two completely independent planes at once. Drive the angles α and β at different rates and the motion never repeats — there is no axis to point at, and no 3-D analogue.

  • Perspective projectionthe shadow
    (x, y, z) ↦ (x, y, z)d − w

    Cast the 4-D point onto our space by dividing through its distance d − w from a 4-D eye. Cells nearer in the unseen axis w loom larger — which is why the inner cube can swell, pass through the outer one, and become it.

  • The 4-D Euler relationalways 0
    V − E + F − C = 0

    Euler’s V − E + F = 2 for solids, lifted a dimension: count vertices, edges, faces and cells with alternating sign and you always get zero, whichever of the six regular 4-polytopes you pick.

  • Power iterationr ← G r
    rk+1 = G rk → v₁

    Just apply G again and again: any starting guess is dragged onto G’s leading eigenvector. It converges at the spectral gap — about a factor d per sweep — which is why ~50 passes rank the whole web.

  • The avalanche lawτ ≈ 1.2
    P(s) ~ s−τ

    Avalanche sizes have no characteristic scale: a straight line on a log–log plot. Most are tiny, a few are system-spanning, and the same law fits earthquakes, forest fires, extinctions and market crashes — the rare giant is not an anomaly, it is the tail of the everyday.

  • The Euler productsum = product
    ζ(s) = Σn 1ns = Πp 11 − p−s

    A sum over all the integers equals a product over the primes— Euler’s bridge from analysis to arithmetic, and the reason the zeta function knows where the primes are.

  • The logarithmic integralbest estimate
    Li(x) = ∫2x dtln t

    Gauss’s sharper guess: integrate 1/ln t instead of taking it flat. Li(x) hugs the π(x) staircase far more tightly than x / ln x while obeying the same leading law.

  • The Riemann zerosthe music
    π(x) = li(x) − Σρ li(xρ) − … , ζ(ρ) = 0

    Riemann’s explicit formula writes the prime count exactly as a smooth term corrected by a sum over the zeros ρ of ζ — the primes are the spectrum of the zeta function. The Hypothesis: every ρ has Re(ρ) = ½.

  • The Prime Number Theoremthe law
    π(x) ~ xln x , so π(x) ln xx → 1

    The count of primes up to x grows like x / ln x — they thin out, but only logarithmically. Gauss guessed it at fifteen, Legendre published it, and it went unproven for a century.

    1896
  • The Prisoner's DilemmaT > R > P > S
    T > R > P > S,  2R > T+S

    Reward R for mutual cooperation, punishment P for mutual defection, temptation T for defecting on a cooperator who gets the sucker payoff S. Whatever the other does, defecting scores more — yet two defectors both do worse than two cooperators. The trap of rational self-interest.

  • Probability conequantile
    S0 e(μ − σ²/2)t ± z·σ√t

    Because log S is Gaussian, every quantile is a formula: the silver cone. Its half-width grows like σ√t — uncertainty widens with the root of time, the same law as diffusion.

  • The hidden halfprobability current
    j = |ψ|² ∇φ = Im(ψ* ∇ψ)

    Probability actually flows, and its current is set by the gradient of the phase φ. Flat phase, no flow; a winding phase is a current going somewhere.

  • The Hadamard gatesuperpose
    H : |0⟩ ↦ |0⟩ + |1⟩√2, |1⟩ ↦ |0⟩ − |1⟩√2

    The gate that makes superposition: it turns a definite bit into an even 50/50 mix — but with a sign. Apply it twice and the signs interfere to undo it, which is the whole trick of quantum algorithms.

  • The CNOT gateentangle
    CNOT : |a, b⟩ ↦ |a, a ⊕ b⟩

    Flip the second bit only if the first is 1. Fed a superposed control it can’t be written as “this qubit times that qubit” any more — the two become one inseparable state. This is where entanglement enters.

  • The Bell pairone circuit
    |Φ⁺⟩ = |00⟩ + |11⟩√2

    A Hadamard then a CNOT, and two qubits are perfectly correlated: each alone is a coin-flip, yet measure one and the other’s answer is fixed — the computational cousin of The Pact’s singlet.

  • The quadratic mapz₀ = 0
    zn+1 = zn2 + c

    Feed the rule its own output, starting from zero. If |z| ever passes 2 the orbit escapes to infinity; otherwise it is trapped forever. That single yes/no, asked at every point, draws everything here.

  • The Bloch sphereon the sphere
    |ψ⟩ = cosθ2|0⟩ + e 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.

  • On a networkNewman 2002
    R₀ = T · (⟨k²⟩ − ⟨k⟩) / ⟨k⟩

    On a contact network the right number uses the excess degree — who your neighbours’ neighbours are. A heavy-tailed (scale-free) degree distribution makes ⟨k²⟩ blow up, so hubs push R₀ far above the well-mixed estimate.

  • Random walk⟨x²⟩ = 2Dt
    Sn = Σ Xi ; Var = nσ²

    A sum of i.i.d. steps goes nowhere on average, yet its spread grows as the square root of the number of steps — the signature √t law of every diffusion.

  • Gravitational redshifttime slows
    1 + z = 1√(1 − rs/r)

    Clocks tick slower deeper in the well; light climbing out loses energy and reddens, diverging as rrs. Together with the disk’s orbital Doppler shift it makes the approaching side blaze blue-white and the receding side fall dim and red.

  • Harmonics = particlesRegge
    α′M² = n − 1

    The string-theory punchline: the same string at oscillator level n is a different particle, its mass climbing the Regge trajectory. n=1 is massless; the spectrum above is the particle zoo.

    1968
  • Reproduction numberR₀
    R₀ = β / γ  (·⟨k⟩ on a network)

    The single most important number in epidemiology: the average number of new infections one case produces in a fully susceptible population. Above 1 the outbreak grows; below 1 it dies.

  • The reservoir updateweights frozen
    xt+1 = tanh( W xt + Win ut )

    A big pool of neurons, wired to each other by a fixed random matrix W and fed the input through random Win. Its state x is a living, fading echo of everything it has heard — a high-dimensional memory that is never trained, only listened to.

  • The linear readoutthe only trained part
    yt = Wout xt

    The answer is just a weighted sum of the reservoir’s state. All the learning lives in this one output layer Wout — the recurrent tangle stays random, so training a recurrent network collapses to fitting a single line.

  • Ridge regressionclosed form
    Wout = Y X( X X + βI )−1

    Collect the reservoir states X against the targets Y and solve one regularised least-squares system — no backpropagation, no epochs, one matrix inverse. The ridge term β keeps the readout from chasing noise in the echo.

  • The echo state propertyfading memory
    ρ(W) < 1

    For the network to be a reservoir at all, it must forget where it started: scale W so its spectral radius — its largest eigenvalue — sits just below one. Below, old inputs decay and the state is a clean echo of recent history; above, the activity never settles and prediction falls apart.

  • Reversibility paradoxLoschmidt 1876
    v → −v ; ΔS ≥ 0

    Loschmidt’s objection: reversible micro-laws can’t logically give an irreversible world. The resolution is the Past Hypothesis — the universe began in a low-entropy state — plus sensitive dependence, which is why real reversal fails.

    1876
  • Reynolds numberRe = UL / ν
    Re = U Lν

    One dimensionless number decides everything: the ratio of inertia to viscosity. Low Re, the fluid moves in smooth sheets; past a critical value (~2000 in a pipe) it shatters into eddies on every scale. The readout above is this number, in the instrument’s own units.

    1883
  • Euler's equationstorque-free
    I1ω̇1 = (I2 − I32ω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 polhodeof Poinsot’s rolling construction.

  • The tennis-racket theoremunstable
    I1 < I2 < I3 ⇒ axis 2 unstable

    Spin about the largest or smallest principal axis is stable; spin about the middle one is not — the body flips 180° over and over, the Dzhanibekov effect filmed on a wing-nut in orbit.

    1991
  • Holling type IIsaturation
    f(x) = a x1 + a h x

    A real predator can only catch and handle so much: at high prey density its kill rate saturates (handling time h). Swapping the cartoon βx for this curve is the step from the textbook cycle to a model that can settle — or crash.

  • Rosenzweig–MacArthurlogistic + saturation
    ẋ = rx(1 − x/K) − y f(x)  ·  ẏ = e·y f(x) − γy

    Give the prey a carrying capacity K (à la The Cascade’s logistic growth) and the predator a saturating appetite f(x), and the neutral cycles resolve — into a stable equilibrium, or a self-sustaining limit cycle born in a Hopf bifurcation.

    1963
  • Paradox of enrichmentmore food → crash
    K ↑  ⇒  Hopf  ⇒  extinction

    Raise the prey’s carrying capacity — enrich the system — and the stable point goes unstable: the populations swing wider and wider until a trough touches zero and someone goes extinct. Generosity destabilises. The same fragility The Cascade finds when a rate is pushed too hard.

    1971
  • The toppling rulez_c = 4
    z(x) ≥ 4 ⟹ z(x) − 4, each neighbour + 1

    Drop one grain on a cell. The instant a cell holds four, it spills one grain to each of its four neighbours — which may themselves topple. One grain can set off a chain reaction: an avalanche.

  • The random graphno hubs
    P(k) = e−⟨k⟩ ⟨k⟩kk!

    Wire the same nodes at random and degrees follow a Poisson bell tight around the mean — almost every node near-average, hubs essentially impossible. The null model that makes the scale-free tail look as strange as it is.

    1959
  • Preferential attachmentrich get richer
    Π(ki) = kiΣj kj

    A new node links to an existing one with probability proportional to how many links it already has. Growth plus this one bias — the already-popular get more popular — is enough to breed hubs from nothing.

    1999
  • Scale-free degreesγ = 3
    P(k) ≈ 2m² ∼ 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.

  • Energy levelsn only
    En = −13.6 eVn²

    Bohr’s ladder, re-derived exactly by Schrödinger. The energy depends on n alone — so 2s and 2p share a rung, a whole shell of shapes at one energy. That accidental degeneracy is the scaffold the periodic table hangs on.

    1913
  • Equation of motionSchrödinger
    i ψt = [ −½ ∇² + V(x,y) ] ψ

    The wavefunction’s whole future, set by its curvature and the potential V. It is complex — amplitude and phase — and only the pair together is the state.

  • Equation of motionSchrödinger
    i ψt = [ −½ ∂²/∂x² + V(x) ] ψ

    The wave function’s whole future, set by its curvature and the potential V. We solve it exactly in time by splitting it into a potential kick and a kinetic drift.

    1926
  • Hydrogen wavefunctionseparable
    ψnℓm = Rnℓ(r) · Ym(θ, φ)

    The bound state factors cleanly into a radial part Rnℓ — the shells, fixed by n and ℓ — and an angular part, the spherical harmonic Ym that draws the lobes. Square it and you have the cloud.

    1926
  • Angular momentumquantised
    |L| = √(ℓ(ℓ+1)) ℏ , Lz = m ℏ

    The length of L and one projection are sharp together; the other two never are. ℓ names the shape (s, p, d, f…), m the tilt — and in the complex basis, m is the visible circulation winding the phase.

  • The metricSchwarzschild
    ds2 = −(1 − rs/r)c2dt2 + … , rs = 2GM/c2

    Clocks deeper in the well tick slower (the −(1−rs/r) term). The event horizon is at rs = 2GM/c², where the escape speed reaches light.

  • The happiness ruleτ
    happy ⇔ nsamenoccupied ≥ τ

    An agent stays put if at least a fraction τ of its occupied neighbours share its colour, else it moves to a random empty cell. That is the whole rule — no agent wants segregation, each only dislikes being a small local minority.

  • Threshold modeltipping
    act ⇔ (fraction already acting) ≥ θi

    The same logic beyond space: each person joins in once enough others have. Whether a crowd tips depends on the distribution of thresholds, not the average — one stubborn holdout can stall a cascade a lower-threshold crowd would ignite.

    1978
  • Segregation indexorder param
    s = ⟨ nsame / noccupied

    The mean like-neighbour fraction across all agents — the order parameter. A random mix sits near 0.5; run the rule with τ as low as 0.3 and s climbs past 0.8, the grid coarsening into single-colour domains exactly like The Threshold’s magnet.

  • Self-organized criticalityno knob
    slow drive + fast relaxation ⟹ ⟨z⟩ → z_c

    Add grains slowly, let avalanches run fast, and the pile climbs to the critical slope and stays there — balanced at the edge with no parameter tuned. The critical point is an attractor, not a setting. This is what The Threshold, The Percolation and The Cascade had to be dialled to by hand.

    1987
  • Semi-Lagrangian advectionunconditionally stable
    u(x, t+Δt) = u(x − u Δt, t)

    Stam’s insight: don’t step the velocity forward (which blows up), trace it backward — ask where each parcel came from and copy that value. This semi-Lagrangian advection is unconditionally stable, which is why it runs in real time in your browser.

    1999
  • 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.

    1927
  • 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.

    1957
  • Imitate the bestcopy your winner
    si ← argmaxj ∈ N(i)∪{i} Πj

    Each cell plays its eight neighbours, sums its payoff, then adopts the strategy of whoever scored highest nearby (itself included). No memory, no foresight, no referee — and yet cooperation organizes in space.

  • How it integratesStrang split-step
    ψeiVΔt/2 ψ ; ψ̂ei|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.

  • 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.

  • Stochastic resonanceoptimal noise
    SNRout peaks at σopt > 0

    A subthreshold signal — too weak to cross a detection threshold on its own — can become detectable by adding the right amount of noise. The output signal-to-noise ratio rises with noise from zero, peaks at an optimal amplitude, then falls again. Noise is not just a nuisance; here it is the carrier.

    1981
  • The look-then-leap rulestop optimally
    reject first r = ⌊N/e⌋, then take first > all seen

    Candidates arrive one at a time in random order; you must accept or reject each on the spot, and you want the single best. The rule: watch the first r ≈ N/e to calibrate, then grab the next one that beats everyone so far.

  • The win probabilitywin rate
    P(win | r) = rN Σi=r+1N 1i − 1

    The chance the look-then-leap rule lands the very best, as a function of the cutoff r. Too small and you leap before you have calibrated; too large and the best has likely already walked past. The maximum sits at r/N ≈ 1/e.

  • The 1/e lawthe 37% law
    rN1e, P(win) → 1e ≈ 0.368

    The striking part: as the field grows, both the optimal cutoff fraction and the probability of success converge to the same number, 1/e ≈ 37% — whether there are ten candidates or ten million. Skill cannot push a random order past this wall.

  • Bruss’ odds theoremone-line stop
    stop at the last index where Σ ri ≥ 1, ri = pi1 − pi

    A single rule for a whole family of stopping problems: sum the odds of the remaining “records” and stop the moment that running sum reaches one. The secretary’s 1/e is just this theorem applied to a random permutation.

  • Freeze & reversek = 1, 2, …
    fs = N ftrue/k ⇒ frozen

    When the strobe lands exactly on a spoke-pass the fold hits zero — frozen. Sample a touch faster and it tips negative: the wheel runs backward while it truly turns forward.

  • Subthreshold signalinvisible without noise
    x(t) = A sin(2π f0t) + σ ξ(t), A < θ

    A periodic signal of amplitude A sits below detection threshold θ — by itself it never fires the detector. Adding Gaussian white noise σ ξ(t) creates random barrier crossings; at the right σ, crossings cluster at the signal’s peaks and the output becomes correlated with the input. Biological sensory neurons, inner-ear hair cells, and crayfish mechanoreceptors all appear to operate near this optimal noise level.

  • Szilard engineW = kT ln 2
    Wmax = kBT ln 2 per measurement

    Leó Szilárd’s one-molecule heat engine: a demon who measures which half of a box the molecule is in can extract kBT ln 2 of work per cycle. But the demon must erase that single bit afterward — and that erasure, Landauer showed three decades later, costs exactly what was gained.

    1929
  • Takens' theoremm ≥ 2d+1
    Φτ : M ↪ ℝm (m ≥ 2d+1)

    For almost any delay τ and any m more than twice the attractor’s dimension d, the delay map is a diffeomorphism — the reconstruction has the same topology and the same invariants (Lyapunov λ, fractal dimension) as the real thing. The geometry survives the shadow.

    1981
  • The colour count1B + 1W
    domino = 1B + 1W ⇒ tiling needs B = W

    On a chequered board every domino, however laid, covers exactly one black and one white square — so any region that tiles must hold equal counts. Opposite corners share a colour; cutting both leaves 32 to 30, and no arrangement closes a gap of two.

  • Gomory’s circuitthe door
    cut 1B + 1W ⇒ ∃ tiling

    Thread all 64 squares onto one closed rook’s tour, colours alternating like beads. Remove one black and one white and the necklace falls into two arcs of even length — each tiles trivially along its own strand. Balance the colours and the wall opens, wherever the knife falls.

  • Tit-for-TatAxelrod's champion
    at = a′t−1,  a0 = C

    Open by cooperating, then echo your opponent’s last move. Nice, retaliatory, forgiving and clear — it won Axelrod’s tournament against every craftier rule. The spatial grid trades this memory for neighbourhood structure: clustering does the work repetition once did.

  • The fundamental diagramflow–density
    q(ρ) = ρ · v(ρ)

    Flow (cars per minute past a point) is density times speed. As density ρ climbs, flow rises to a peak at the critical density ρc, then falls back to zero in gridlock — the inverted-U every road obeys, and the proof that more cars can mean less throughput.

    1935
  • The LWR equationa jam is a shock
    ∂ρ∂t + ∂q∂x = 0

    Cars are conserved, so density flows like a fluid — and where fast, sparse traffic meets a dense pack the solution forms a shock: the sharp upstream edge of a jam, a discontinuity that drifts backward through the cars at a speed set by the diagram, not by the drivers.

    1955
  • The Nagel–Schreckenberg rulesfour steps
    v ← min(v + 1, g, vmax); v ← v − 1 (prob p); x ← x + v

    Each car, each tick: speed up by one, but never faster than the gap g ahead or the limit vmax; then with probability p tap the brake for no reason; then move. That one random tap is the whole story — set p = 0 and traffic never jams on its own.

    1992
  • Transmissionexact (E < V₀)
    T = [ 1 + V₀² sinh²(κa)4E(V₀E) ]⁻¹

    The exact leak for one sharp energy. The sinh makes it fall exponentially with width a and with √(V₀−E) — why a hair more wall kills it.

  • Conservationunitary
    T + R = 1

    Probability is conserved: what doesn’t get through comes back. The instrument measures T by integrating |ψ|² past the wall — its own answer for a packet of many energies.

  • Why it patternsTuring instability
    Dv / Du> 1 ⟹ uniform state unstable → λ* selected

    Short-range activation, long-range inhibition: when the inhibitor diffuses faster (DvDu), one wavelength λ* wins and tiles the plane.

  • Turing’s mechanism1952
    u/∂t = Du∇²u + f(u,v) ; ∂v/∂t = Dv∇²v + g(u,v)

    Two morphogens, each diffusing and reacting. Turing’s shock: diffusion, usually a smoother, can destabilise a uniform state and manufacture pattern.

  • The correlation functionhow clustered
    ξ(r) = ⟨δ(x) δ(x + r)⟩

    The honest measure of clustering: how much more likely than chance you are to find a second clump a distance r from a first. It climbs from zero as the web sharpens — read live as the cloud collapses.

  • Two-slit intensityfar field
    I = I0 · sinc²(πa·u/λ) · cos²(πd·u/λ)

    A single-slit diffraction envelope times a two-slit interference fringe. The cos² term is the interference; kill it and only the smooth envelope is left.

  • Universal approximationCybenko 1989
    any continuous f, to any ε

    Why one hidden layer is, in principle, enough: a wide enough network of these neurons can approximate any continuous function. Capacity is the catch — Spiral needs more units to bend the boundary far enough.

  • Vicsek updateminimal model
    θit+1 = ⟨θjNi + η·ξ

    Replace steering with "copy the average heading of your neighbours, plus a kick of noise η". This minimal model is the one with an exact order parameter and a phase transition.

  • The Voronoi cellnearest seed
    V(pi) = { x : ‖x − pi‖ ≤ ‖x − pj‖ ∀ j }

    Every point in the plane belongs to its nearest seed. That one rule tiles space into cells whose walls are the perpendicular bisectors between seeds — the partition nature reaches for whenever things compete for the space around them.

  • The Delaunay dualthe mesh
    Del(P) = dual of Vor(P) · empty-circumcircle

    Join the seeds of every pair of touching cells and you get the Delaunay triangulation — the Voronoi diagram’s twin. It is the “best” mesh on a set of points: no point sits inside any triangle’s circumcircle, so it avoids thin slivers.

  • Lloyd relaxationmove to centroid
    pi ← centroid( V(pi) ), repeat

    Let each seed hop to the centre of mass of its own cell, recompute the cells, and repeat. A random scatter relaxes into an even honeycomb — the identical assign-then-average loop as k-means clustering.

  • Quantization energyminimised
    E = ΣiV(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.

  • Wave equationclamped ends
    ∂²y/∂t² = c² ∂²y/∂x²

    A clamped string obeys the wave equation with both ends pinned. Those boundary conditions are the whole story — they force the continuous string into a discrete set of allowed shapes.

    1747
  • Wave–particle dualitythe law
    V2 + D2 ≤ 1

    The exact trade. V is fringe contrast (wave), D is how well the path is known (particle). You can buy one only by selling the other — never both at once.

    1995
  • Time–frequency limitthe wall
    Δt · Δf ≥ 1

    No signal is sharp in both time and frequency: their spreads multiply to a fixed floor. The Gaussian alone meets it — the classical twin of Heisenberg’s Δx·Δp ≥ ℏ/2.

  • Continuous wavelet transformtime × scale
    W(a,b) = 1√a ∫ x(t) ψ*(t − ba) dt

    Correlate the signal against one little wave ψ, scaled by a (≈ 1/frequency) and slid to time b. |W|² is the scalogram — energy laid across time and frequency at once.

    1984
  • The Morlet waveleta ripple in a bell
    ψ(t) = π−¼ eiω₀t e−t²/2

    A complex sinusoid wrapped in a Gaussian. The width ω₀ is its wiggle-count: more of them sharpens frequency and blurs time, fewer does the reverse — the one knob on the trade.

  • Admissibilityzero mean
    Cψ = ∫ |Ψ(ω)|²|ω|dω < ∞ ⇒ ∫ ψ dt = 0

    For the transform to invert, the wavelet must average to zero — a true little wave, not a lump. That’s why ω₀ can’t fall too low: too few wiggles and the mean stops vanishing.

  • Why it spreadsdispersion
    σ(t) = σ0 √(1 + (t/2σ0²)²)

    A free packet doesn’t hold together — the faster components outrun the slow, so σ grows with time. The spread you watch is the price of ever having localised it.

  • Pareto’s tailpower law
    P(m > x) ∝ x−α,  α ≈ 1–2

    The richest tail is not exponential but a power law: a straight line on log–log axes, no characteristic scale. The famous 80/20 — roughly a fifth of the agents hold four-fifths of the wealth — and the same fat tail The Sandpile and The Web draw in other keys.

  • Money as energyexponential
    P(m) ∝ e−m/T,  T = ⟨m⟩

    Conserve a fixed pot of money and let agents swap it at random, and wealth settles into the same exponential law that energy takes among gas molecules — the money temperature T is just the average balance. Fairness at the microscale, a heavy skew at the macroscale.

    2000
  • The transactionconserved
    mi′ = mi + Δ,   mj′ = mj − Δ

    Every trade only moves money between two agents — the total is exactly conserved, like energy. A fair Δ relaxes to the exponential; a biased Δ (the loser always pays a slice of the poorer purse) condenses nearly all of it onto one agent — Bouchaud–Mézard oligarchy.

  • The Gini coefficient0…1
    G = Σij |mi − mj|2N²⟨m⟩

    One number for the whole distribution: 0 is perfect equality, 1 is one agent owning everything. A fair-exchange economy sits at exactly G = ½; the biased rule drives it toward 1 as the money condenses.

  • The digit enginex ↦ {1/x}
    a = ⌊1/x⌋ , x ↦ 1x − a

    Flip the fraction, keep the fractional part, write down the whole number that fell out — repeat, and any ω dictates its continued fraction digit by digit. Euclid’s algorithm wearing dynamical clothes; the golden number is the fixed point of the first branch, x = 1/(1+x), so its tape reads 1, 1, 1, … forever.

  • The floor√5 is sharp
    |ω − p/q| < 1√5 · q²

    Every irrational is caught this close by infinitely many fractions — and no constant larger than √5 works, because one number sits exactly on the floor: for φ the normalized error q²·|φ − p/q| converges to 1/√5 ≈ 0.4472 and never below. The worst-served number in arithmetic, named by its distance from every fraction.

  • The armor conditionDiophantine
    |ω − p/q| ≥ γqτ ⇒ the torus survives

    KAM’s survival clause is an arithmetic condition: a torus outlives a small kick exactly when its winding number keeps this far from every fraction. Being badly approximable — the crime Hurwitz convicted φ of — is here the armor. Same property, opposite loyalty.

  • The kicked windingarea-preserving
    p ← p + k sin 2πx , x ← x + p

    The winding law with one kick bolted on. At k = 0 every horizontal line p = ω is an invariant circle — the pure winding, stacked. Turn k up and the kick feeds on rational windings first: each p/q grows an island chain, the chains grow a chaotic sea, and the only circles left standing are the ones the kick cannot resonate with.

  • The residue dialR → ¼
    R = 2 − Tr M4 , R(Fn/Fn+1) → ¼ at kc

    Greene’s dial: linearize the map around the periodic orbits whose rotation numbers are the golden winding’s own convergents — the Fibonacci fractions besieging it — and read the residue R. Falling toward 0: the circle holds. Rising past ¼ up the whole ladder: it is gone. The ladder pinches flat at kc ≈ 0.971635.

  • The windingthe whole law
    x ← x + ω (mod 1)

    Turn by the same twist ω, forever — the pair’s entire law. Arithmetic reads its near-returns as fractions; dynamics stacks one such circle for every ω and asks which survive a shake. The law itself never says what the number is for.

  • The ladderpₙ/qₙ
    pn = anpn−1 + pn−2 , qn = anqn−1 + qn−2

    Each digit buys one rung: the convergents pₙ/qₙ, the only fractions that beat every smaller denominator. Big digits are windfalls — π’s a₄ = 292 turns 355/113 into seven digits of π. All-ones digits are the slowest ladder arithmetic allows: the golden ratio’s rungs are the Fibonacci numbers.

  • Number = behaviourWolfram code
    sit+1 = ⌊R / 2k⌋ mod 2

    Index the neighbourhood by its bits, k = 4si−1 + 2si + si+1, and the rule number R is literally the truth table — bit k of R is the next state.

  • Zel'dovich's approximationthe engine
    x(q, t) = qD(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
Locator — Scan the archive
esc
Indexing archive…

Indexing the archive…

select openesc close