Dynamics

A point that pulls, and a point that pushes

Every crossing of a curve with the diagonal is a value the rule leaves alone. Whether anything ever arrives there is decided by one number — the slope at the crossing — and the picture makes the reason obvious.

Worth reading first: The staircase that shows the whole orbit · Something always stays put.

A fixed point is a value the rule returns unchanged: f(x)=xf(x) = x. On the graph it is a crossing of the curve with the diagonal, and finding one is a matter of solving an equation.

Whether an orbit ever gets to one is a completely different question, and the answer is not in the equation. It is in the slope.

A fixed point that attracts, and one that does notThe same map at two parameters, with the staircase walking towards the crossing in one and away in the other.r = 2.6slope in (−1, 1) — attractingr = 3.3slope outside (−1, 1) — repellingat r = 2.6 the slope at the crossing is -0.60 and the staircase walks inat r = 3.3 it is -1.30 and the staircase walks out — the crossing has not moved, its steepness has
Fig. 1 The same map at two parameters. Both have a fixed point where the curve meets the diagonal; in the left panel the staircase walks into it and in the right panel it walks away. The generator computes the slope at each crossing and checks that the orbit’s behaviour matches it.

The crossing has not moved between the panels — it is at 11/r1 - 1/r in both, and it exists for every rr above one. What has changed is how steeply the curve passes through it, and that is what decides everything.

Why the slope is the whole answer

Zoom in on a crossing until the curve looks straight. Near a fixed point x\*x^\* the map is approximately

f(x)x\*+m(xx\*),f(x) \approx x^\* + m\,(x - x^\*),

where m=f(x\*)m = f'(x^\*) is the slope there. Write dd for the distance from the fixed point. One step multiplies dd by mm; two steps multiply it by m2m^2; nn steps by mnm^n.

So the whole question is what powers of mm do. If m<1|m| < 1 they shrink to nothing and the orbit converges. If m>1|m| > 1 they grow and the orbit leaves. The threshold is exactly one and there is nothing arbitrary about it: one is the slope of the diagonal, and the diagonal is what the map is being compared against.

That is the entire theory of stability for a one-dimensional map, and it takes four lines because the linearisation is honest — near enough to the fixed point, a smooth map really is its own tangent line.

Reading it off the picture

The algebra is short and the picture is shorter.

the logistic map at 2.6, iterated from 0.2A map drawn as a curve with the diagonal across it, and the staircase that iterating it produces.xf(x)x ↦ 2.6x(1 − x), started at 0.2the orbit never repeats — no cycle of eight or fewer closes it
Fig. 2 At r=2.6r = 2.6 the slope at the crossing is 2r=0.62 - r = -0.6. Negative, so the staircase alternates sides; smaller than one in size, so each alternation is closer in. The result is the inward spiral.

The sign of the slope is visible too, and it says something the size does not. A positive slope gives a staircase that approaches from one side, monotonically. A negative slope gives one that overshoots each time and comes back — the spiral above. Both converge when the size is under one; they simply converge differently, and an orbit that alternates about its limit is a different experience from one that creeps up on it.

the logistic map at 3.3, iterated from 0.2A map drawn as a curve with the diagonal across it, and the staircase that iterating it produces.xf(x)x ↦ 3.3x(1 − x), started at 0.2the orbit settles into a cycle of 2 points
Fig. 3 At r=3.3r = 3.3 the slope is 1.3-1.3. The alternation is still there and the overshoot now grows: each swing is thirty percent wider than the last until the staircase gives up on the crossing entirely and finds the rectangle instead.

Watching those two side by side is the argument. Nothing about the crossing changed. The map’s steepness there crossed a threshold, and a point that was pulling started pushing.

The logistic map’s own arithmetic

For the map f(x)=rx(1x)f(x) = rx(1-x) the whole thing is computable in closed form, which is unusual and worth using.

The fixed points are 00 and x\*=11/rx^\* = 1 - 1/r. Differentiating gives f(x)=r(12x)f'(x) = r(1 - 2x), so the slope at the origin is rr and the slope at the other crossing is

f(x\*)=r(12(11r))=2r.f'(x^\*) = r\left(1 - 2\left(1 - \tfrac1r\right)\right) = 2 - r.

Both facts fall out immediately. The origin is attracting for r<1r < 1 — a population that dies out — and repelling above it. The other fixed point is attracting when 2r<1|2 - r| < 1, which is 1<r<31 < r < 3, and repelling above three. Every figure in this essay checks both against the drawn orbit.

That is why r=3r = 3 appears where it does. It is not a value anybody chose; it is where 2r2 - r reaches 1-1, and it is the last parameter at which this map has anything to settle on.

What happens exactly at the threshold

The interesting case is the one the inequality leaves out: m=1|m| = 1 exactly.

Linearisation says nothing there. One step multiplies the distance by something of size one, so the tangent line predicts an orbit that neither approaches nor leaves — and the tangent line is not the map. What actually happens is decided by the next term in the expansion, the curvature, and it can go either way.

the logistic map at 3, iterated from 0.2A map drawn as a curve with the diagonal across it, and the staircase that iterating it produces.xf(x)x ↦ 3x(1 − x), started at 0.2the orbit never repeats — no cycle of eight or fewer closes it
Fig. 4 Exactly at r=3r = 3, where the slope is 1-1. The orbit does converge, and it converges so slowly that after two hundred steps it is still visibly short of the crossing — the approach is like 1/n1/\sqrt{n} rather than exponential.

This matters practically as much as theoretically. Near a stability threshold, convergence becomes arbitrarily slow, so a system observed for a fixed length of time looks unsettled long before it becomes unstable. “It has not settled yet” and “it will never settle” are different claims, and near the threshold no finite observation separates them.

The name for a parameter where a fixed point changes stability is a bifurcation, and this one — where the slope passes through 1-1 and a two-cycle appears — is a period-doubling bifurcation. It is the first of infinitely many, which is the next rung of this ladder.

The cycle has a slope too

A two-cycle is a fixed point of the map applied twice, so the same test applies to it — and applying it says where the next threshold will be.

If aa and bb are the two points of the cycle, the slope of fff \circ f at aa is, by the chain rule, f(b)f(a)f'(b)\,f'(a): the product of the slopes at the two points. The cycle attracts exactly when that product is smaller than one in size.

the logistic map at 3.4, iterated from 0.2A map drawn as a curve with the diagonal across it, and the staircase that iterating it produces.xf(x)x ↦ 3.4x(1 − x), started at 0.2the orbit settles into a cycle of 4 points
Fig. 5 At r=3.4r = 3.4 the two-cycle is attracting: the staircase spirals onto the rectangle rather than merely reaching it. The fixed point at the centre is still there and still repelling, which is why the orbit passes it and keeps going.

The product starts near zero when the cycle first appears, grows as rr rises, and eventually reaches 1-1 — at which point the two-cycle stops attracting and a four-cycle appears. Then the four-cycle’s slope, a product of four numbers, does the same thing.

So the mechanism repeats itself, and it repeats on an object that is twice as complicated each time. That is the reason the doublings keep coming, and it is also the reason they come faster and faster: a product of many slopes passes through 1-1 over a shorter range of rr than a product of few.

The same criterion, everywhere

The slope test is not about the logistic map or about this field. It is the criterion for any iterative process, and this collection has already relied on it twice without naming it.

Newton’s method is the map xxf(x)/f(x)x \mapsto x - f(x)/f'(x). Its fixed points are the roots of ff, and its derivative at a root is zero — which is a slope not merely under one but as favourable as it can be. That is precisely why Newton’s method converges so fast: the distance is not multiplied by a constant, it is squared, so the number of correct digits doubles each step.

Euclid’s algorithm is the same story in whole numbers. Each step strictly shrinks the smaller side, which is a contraction, and a contraction on a set with no infinite descending chain must terminate. The geometric version — the leftover strip always smaller than the one before — is the slope test with the analysis removed.

The general statement is the contraction mapping theorem: a map that shrinks every distance by a fixed factor has exactly one fixed point, and every orbit converges to it. It is one of the few theorems that gives existence, uniqueness and a method of finding it in the same sentence.

A fixed point worth computing

The test is not only diagnostic. Turning an equation into a map whose fixed point is its solution, and then iterating, is a method — and the collection already has its best example.

The golden ratio satisfies φ=1+1/φ\varphi = 1 + 1/\varphi. Read the right-hand side as a map, g(x)=1+1/xg(x) = 1 + 1/x, and φ\varphi is its fixed point. The slope there is 1/φ2=0.382-1/\varphi^2 = -0.382, comfortably inside the interval, so iterating from almost any positive start converges — and the successive values are 1,2,3/2,5/3,8/51, 2, 3/2, 5/3, 8/5, which are the Fibonacci fractions.

That is the same sequence the continued fraction produces, arrived at as an orbit rather than as a peeling, and the slope of 0.382-0.382 is a quantitative statement about it: each step cuts the error to about two fifths, so the convergents gain roughly one bit of accuracy per step. The slowness of the golden ratio’s approximation, which the number field measured with Hurwitz’s theorem, is this slope seen from the other side.

The general recipe is worth naming because it is how a great deal of numerical work is actually done: to solve F(x)=0F(x) = 0, invent a map whose fixed point is the root, check its slope, and iterate. Newton’s method is the choice of map that makes the slope zero.

What it does not give

The contrast with Brouwer’s theorem is the useful part of this essay, and it is a contrast between two kinds of guarantee.

Brouwer says a continuous map of a disc to itself has a fixed point. It needs no slope, no smoothness, no contraction — only continuity and the right shape of domain. And it hands over nothing: the proof is topological, it produces no point, and it says nothing about whether iterating the map would find one.

The slope test says the opposite thing in both directions. It needs smoothness and a fixed point already in hand, and in exchange it says whether iteration will get there and how fast. One theorem is about existence and gives no method; the other is about method and assumes existence.

Most of the useful work in this field lives in the gap. A system can have a fixed point that nothing ever reaches — the right panel of the first figure is exactly that — and calling such a point an equilibrium is true and misleading, because nothing observed will ever be found at it.

Several fixed points at once

A map with more than one crossing has more than one fixed point, and each gets its own verdict from its own slope.

the tent map at 1.6, iterated from 0.2A map drawn as a curve with the diagonal across it, and the staircase that iterating it produces.xf(x)x ↦ 1.6 min(x, 1 − x), started at 0.2the orbit never repeats — no cycle of eight or fewer closes it
Fig. 6 The tent map at 1.61.6 has fixed points at 00 and at r/(1+r)r/(1+r). The slope is ±1.6\pm1.6 at both, so both repel, and the orbit settles on neither — it has nowhere to go and keeps moving.

When every fixed point repels, the orbit is forced into whatever else the map offers: a cycle, or nothing periodic at all. That is not a rare configuration. It is what “chaotic” looks like from the point of view of this essay’s criterion — every fixed point and every cycle repelling, so the orbit is pushed off each of them in turn and never lands.

Which is worth stating plainly, because it inverts the usual intuition. Chaos is not the absence of structure. It is a great deal of structure, all of it unstable.

What a physicist calls it, and what changes

The same criterion is called several things, and the names are worth knowing because they are how the idea reaches other subjects.

An engineer calls a fixed point with every slope under one asymptotically stable, and the analysis of whether a control system settles is exactly this test applied to the loop’s Jacobian. An economist calls it a stable equilibrium and asks whether a market perturbed slightly returns. A biologist modelling a population asks whether a species knocked off its carrying capacity comes back or oscillates — which is the logistic map’s own original setting, and the reason rr in these figures is a growth rate.

What none of those renamings changes is the mathematics: a number is compared against one. What they change is what the number means, and therefore what the threshold costs. Passing r=3r = 3 in a fisheries model is a population that starts alternating between good years and bad ones for reasons internal to the arithmetic rather than external to the fishery, which is a genuinely counter-intuitive thing to have to explain to somebody looking for a cause.

That is the standing use of this field: it supplies causes that are not events. A system can change its behaviour completely with nothing happening to it, because a parameter drifted past a threshold that nothing marks.

The basin, and its edges

Attracting is a local property, and the natural next question is how far it reaches. The set of starting points whose orbits converge to a given fixed point is its basin of attraction, and for a one-dimensional map with a single hump the basins are usually intervals with tidy ends.

Usually. The basin’s boundary consists of points that go to neither — typically the repelling fixed points and their preimages — and those form a set that can be very far from tidy.

In one dimension the mess is limited, because there is only one direction to be pushed in. In the plane it is not limited at all: the basin boundary of a perfectly reasonable map can be a curve of infinite length, and every point on it can have every basin arbitrarily close by. That is what Newton’s method on a complex cubic does, and it is one of the places where a picture is the only reasonable way to state the result.

More than one dimension

In the plane the same criterion holds with the slope replaced by a matrix, and the replacement brings in a piece of machinery this collection has already built.

Near a fixed point, a smooth map of the plane is approximately its Jacobian — the matrix of partial derivatives — and iterating the map is approximately iterating that matrix. What a matrix does under repetition is decided by the directions it leaves alone and the factors it stretches them by: its eigenvectors and eigenvalues.

So the test becomes: the fixed point attracts when every eigenvalue is smaller than one in size. If one exceeds one and another does not, the point attracts along one direction and repels along another — a saddle, which has no one-dimensional analogue and is the configuration most of the interesting behaviour in the plane is built from.

A linear map redrawing the planeThe integer grid before and after a linear transformation; the shaded unit square becomes a parallelogram whose area is the determinant.beforeafter · area × 0.820.60.20.11.4
Fig. 7 A linear map with one eigenvalue under one and one over it. Repeating it pulls the plane in along one direction and pushes it out along another; a fixed point of a nonlinear map with this Jacobian behaves the same way near enough to itself.

That is also where the tidy one-dimensional picture ends. A saddle’s stable and unstable directions are curves rather than lines once the map is nonlinear, and when those curves cross each other transversally the result is a tangle that Poincaré found in 1890 and declined to draw, writing that he would not even attempt it. The Lorenz system is what one looks like when a computer draws it instead.

Where the ladder goes next

This is the third rung of the fixed-point ladder. Brouwer gave existence with no method; the cobweb gave a way to see an orbit; this rung says which fixed points an orbit will actually find.

The next question is what happens after the slope passes 1-1. The fixed point stops attracting and a two-cycle appears — and the two-cycle has a slope of its own, which eventually passes 1-1 too. Following that all the way is the road paved with doublings, and it ends somewhere finite, which is the surprising part.

What links here

Computed from the collection, not written here: the essays that point at this one.

Reads more easily once this is understood

Essays that name this one as worth reading first.

Shares its objects with

Essays that name at least two of the same things, and that neither author linked.

Named objects

A dashed tag is an object no other essay names yet.

CobwebContractionConvergenceDerivativeFixed pointIterationLinearisationLogistic mapOrbitStability