Chaos explorer: logistic map cobweb and bifurcation diagram, Lorenz attractor in 3D, Hénon map, Lyapunov exponents

Explore chaos in your browser: the logistic map (or any typed map) with cobweb, orbit diagram, Lyapunov exponent and Feigenbaum period-doubling points; the Lorenz and Rössler attractors in a 3D view you can turn, with their Lyapunov spectrum; and the Hénon map.

Keep it one tap away for everyday use.

Bookmark Chaos explorer?

Press Ctrl + D now. Your browser opens its own bookmark box, where you choose the name and folder.

  1. Choose r
  2. Cobweb
  3. Diagram
  4. Lyapunov
Nothing is uploadedWorks offline once open

One-dimensional map

The logistic map r·x·(1 − x) is the classic. Try r*sin(pi*x) (the sine map) or r*x*exp(-x) (Ricker).

Orbit diagram and Lyapunov exponent

For 600 values of r: 500 steps are dropped, then the next 150 values of x are drawn. One dot means a fixed point, two a 2-cycle, a smear chaos. The dashed line is your r.

The Lyapunov exponent λ is the average of ln|f′(x)| along the orbit (20,000 steps): below 0 nearby starts come together, above 0 they separate exponentially, which is chaos.

Period-doubling points

Report

The report has every part you have opened.

Download

What it does

Chaos means a simple, fully determined rule whose future still cannot be predicted in practice, because tiny differences in the start grow exponentially. This tool lets you see it in three classic systems and measure it with the Lyapunov exponent, the average rate at which nearby starts separate.

  • Logistic map, x → r·x·(1 − x), or any map you type in x and r. You get the cobweb diagram, the orbit step by step, the orbit (bifurcation) diagram over a range of r, the Lyapunov exponent against r, and the period-doubling points with the ratios that tend to Feigenbaum’s constant δ = 4.6692.
  • Lorenz and Rössler systems, the famous strange attractors, drawn in a 3D view you can turn by dragging. A button works out all three Lyapunov exponents.
  • Hénon map, a two-dimensional map with a fractal attractor, and both of its exponents.

How to use it

  1. Logistic map: type r or move the slider. The big line says whether the orbit settles, cycles (with its period) or is chaotic. The orbit diagram and the λ curve show the whole range of r at once, with your r marked by a dashed line. Change “r from” and “r to” to zoom in, for example 3.54 to 3.58 to see the doublings crowd together.
  2. Lorenz and Rössler: pick the system and its three parameters. The path is solved with fourth-order Runge-Kutta. Drag the picture or use the Turn and Tilt sliders. Press Work out the exponents for the Lyapunov spectrum.
  3. Hénon map: change a and b; up to 200,000 points are drawn.

Each part has a CSV download, and the report (image, PDF or Word) holds every part you have opened.

Worked example

At r = 3.7 the logistic map is chaotic: the page shows λ = 0.3528 (averaged over 20,000 steps; longer runs give 0.3535). At r = 3.5 it settles on a cycle of period 4 with λ = −0.8725, and at r = 4 the exponent is ln 2 = 0.6931, the exact answer. Over r = 2.5 to 4 the tool finds seven period doublings: r = 3, 3.449489743 (= 1 + √6), 3.544090360, 3.564407266, 3.568759420, 3.569691610 and 3.569891259. The ratios of the gaps are 4.7514, 4.6563, 4.6682, 4.6687 and 4.6691, closing in on δ = 4.6692, and the cascade ends near r∞ ≈ 3.569946.

For the Lorenz system with σ = 10, ρ = 28, β = 8/3, averaging over 1000 time units gives λ₁ = 0.9069, λ₂ = −0.0007 and λ₃ = −14.573, close to the known 0.9056, 0, −14.5723. Their sum is −13.667 = −(σ + 1 + β), as it must be, and the Kaplan-Yorke dimension is about 2.06. The Hénon map with a = 1.4, b = 0.3 gives λ₁ = 0.41801 and λ₂ = −1.6220 from 100,000 steps (a million steps give 0.4195 and −1.6235), always summing to ln 0.3 = −1.2040; the Kaplan-Yorke dimension is 1.258.

How the maths was checked

The map exponents, period-doubling points and Hénon exponents were compared with independent NumPy and SciPy code (exact derivatives, fsolve and QR), agreeing to 10⁻⁷ or better; the Lorenz and Rössler spectra with SciPy’s solve_ivp and NumPy QR, and with the known published values.

Good to know

  • A finite average is an estimate: the Lorenz λ₁ wobbles by about ±0.01 for 1000 time units. λ₂ of a flow should be 0; how close it gets shows the accuracy.
  • The period is found by checking when the orbit repeats to 10⁻⁷, after 5000 steps, up to period 64. Right at a bifurcation the orbit settles very slowly and the period can be misread.
  • The 3D view is a plain projection, so depth is shown by turning it, not by perspective.

Frequently asked questions

Is anything uploaded? No. Everything is worked out in your browser and nothing is stored.

What does a negative λ mean? Nearby starts come together, so the motion is predictable: a fixed point or a cycle. Positive λ is chaos.

Can I study another map? Yes: type any formula in x and r, such as rsin(pix). The period-doubling search works for any smooth map with one hump.

Where the numbers come from

  • Methods: Lyapunov exponents of maps as the mean of ln|f′|; Benettin et al. (1980) and Shimada and Nagashima (1979) for flow spectra; period-doubling points solved from f^p(x) = x and (f^p)′(x) = −1. The systems and their classic parameter values are standard textbook mathematics.

Standard scientific values, such as constants, atomic weights and melting points, are facts found in every textbook; we name where ours were checked. If you are an author or publisher and want a value, a credit or a table changed or removed, tell us and we will act quickly.

Last updated 2026-11-09.

Related tools

Light curve period finder #48 Paste times and brightness measurements, even with gaps and uneven spacing, and find the period with a Lomb-Scargle periodogram, false alarm probabilities, alias checks and a folded light curve. Publication plot maker #27 Paste columns or open a CSV or Excel file and get a clean scientific figure with error bars, sized for a journal column. Download vector SVG or PDF, or PNG at 300 or 600 dpi. Your data stays in your browser. Fractal explorer #547 Zoom deep into the Mandelbrot set and Julia sets, drawn in a background worker with smooth colours, or grow the Barnsley fern, Sierpinski triangle, carpet and Heighway dragon by the chaos game, with their fractal dimensions. Download a 2400 × 1800 PNG. ODE solver with phase portraits #548 Type a system of differential equations and solve it in your browser with RK4 or adaptive RK45 (as SciPy). Time plots, phase portrait with direction field and nullclines, equilibria with eigenvalues, a bifurcation diagram, and a heat and wave equation playground. Cellular automata explorer #589 Explore Wolfram's 256 elementary cellular automata as space-time pictures, and run Conway's Game of Life or any Life-like B/S rule on a grid you can draw on, with famous patterns, step and run, wrap or dead edges, cycle detection and RLE import and export. Fourier series and FFT visualiser #549 See how sines and cosines build a square, sawtooth or any typed wave, watch the Gibbs overshoot, draw shapes with epicycles, and take the FFT of a formula or your own data with Hann, Hamming, Blackman or flat-top windows. All in your browser.