Get AI summaries of any video or article — Sign up free
Newton’s fractal (which Newton knew nothing about) thumbnail

Newton’s fractal (which Newton knew nothing about)

3Blue1Brown·
6 min read

Based on 3Blue1Brown's video on YouTube. If you like this content, support the original creators by watching, liking and subscribing to their content.

TL;DR

Newton’s method updates a guess using x_{n+1}=x_n−p(x_n)/p′(x_n), repeatedly applying a local linear approximation to p(x).

Briefing

Newton’s method turns a simple root-finding rule into an endlessly intricate fractal when it’s run over the complex plane. Starting from a seed value, the iteration repeatedly replaces the polynomial by its local linear (tangent) approximation and jumps to where that line hits zero. Most seeds quickly converge to one of the polynomial’s true roots, but the “map” of which root each seed reaches changes violently under tiny perturbations—so the boundaries between outcomes become fractal and never settle into smooth curves.

The setup begins with a practical problem: solving polynomial equations. In engineering and computer graphics, polynomials show up constantly—for instance, Bezier curves used to render fonts. To decide how pixels relate to a curve, systems often reduce distance-to-curve questions to calculus: the minimum of a polynomial-derived function is found by setting its derivative (another polynomial) equal to zero. That naturally raises the computational question of how to find polynomial roots beyond the quadratic and cubic cases where closed-form formulas exist.

For higher degrees, exact formulas become unreliable or nonexistent in the usual sense. The transcript notes the classical result that general quintic polynomials cannot be solved by a universal formula built from standard operations and radicals, pushing practice toward numerical methods. Newton’s method is one such workhorse: pick an initial guess x0, compute x1 = x0 − p(x0)/p′(x0), then repeat. Convergence is fast when the starting point lies near a root, but it can bounce around or fail when it doesn’t.

The fractal emerges when the same iteration is applied to complex starting points. Even if a polynomial has only one real root, allowing complex inputs guarantees additional roots in the complex factorization (via the fundamental theorem of algebra). When seeds are colored by which root they approach, large regions become uniform—each region corresponds to seeds that reliably converge to a particular root. The chaos concentrates on the boundaries between these regions: two seeds that are extremely close can end up converging to different roots.

A key structural reason for the boundary’s roughness is described through a shared-boundary property: for any one color region (seeds converging to a specific root), the boundary of that region matches the boundary of every other color region. This implies that no small neighborhood can contain exactly two colors; any tiny circle around a boundary point must intersect all colors or only one. As a result, the boundary cannot include smooth segments, because smoothness would typically separate only two outcomes locally. The transcript connects this to fractal dimension, citing a boundary dimension around 1.44 for the shown example.

The complexity is therefore not an accident of one polynomial. It also explains why quadratic cases look comparatively tame: with only two roots, a “two-color” boundary is the only possibility, so the constraint doesn’t force the same kind of visible fractal boundary behavior. The deeper mathematical framework behind the shared-boundary phenomenon is attributed to holomorphic dynamics.

Finally, the transcript closes by emphasizing the irony of naming: Newton had no access to these complex-plane visualizations, yet his method generates them. The broader point is that simple ideas can hide future relevance—so questions about Newton’s method (like whether iterations can cycle) connect to other major objects such as the Mandelbrot set, setting up a follow-up discussion.

Cornell Notes

Newton’s method iteratively finds polynomial roots by jumping from a seed x to where the tangent-line (linear approximation) to p(x) crosses zero: x_{n+1}=x_n−p(x_n)/p′(x_n). In the complex plane, coloring each seed by which root it converges to produces large uniform regions, but the boundaries between regions are fractal and infinitely detailed. The transcript highlights a structural theorem-like property: the boundary of the basin of attraction for any one root is the same as the boundary for every other root, meaning tiny circles near the boundary intersect either all colors or only one. This forbids smooth boundary segments and forces roughness at every zoom level, consistent with a measured fractal dimension (~1.44 in the example). The behavior is framed as a consequence of holomorphic dynamics rather than a quirk of one polynomial.

Why does Newton’s method sometimes converge quickly and sometimes behave chaotically?

Newton’s method works best when the initial guess is near a true root, because the tangent-line approximation to p(x) is then a good local model of where p(x)=0. When the guess is far from any root, the linear approximation can point toward the wrong region, producing bouncing or wandering before (or without) eventually landing near a root. The transcript illustrates this with a shifted polynomial on the real line: starting from the same x0, changing the polynomial upward makes the iterates bounce around a local minimum instead of homing in on the true root.

How does the complex-plane version of Newton’s method create a “map” with sharp boundaries?

For complex seeds, the iteration still uses the same update rule based on p(x) and p′(x), even though tangent lines and x-axis crossings no longer have a direct geometric meaning. Each seed is colored by which of the polynomial’s complex roots it approaches. Most seeds fall into basins of attraction (uniform colors), but seeds near basin boundaries can converge to different roots. The transcript emphasizes that tiny changes to the seed can flip the final root, concentrating unpredictability on the boundary.

What is the “shared boundary” property, and why does it force fractal roughness?

Pick one root’s basin (one color). The boundary of that basin turns out to be exactly the same set as the boundary of every other basin. Using the formal definition of boundary (every arbitrarily small circle around a boundary point contains points both inside and outside the set), the transcript argues that near such points, tiny neighborhoods cannot contain only two colors. Therefore, smooth boundary segments—which would typically separate just two outcomes locally—cannot persist at any scale. The boundary must remain jagged under zoom, producing infinite detail.

Why is the quadratic case comparatively “boring” in these pictures?

Quadratic polynomials have only two roots, so the coloring has only two colors. The shared-boundary constraint still allows fractal boundaries between two colors, but the transcript notes that the Newton fractal doesn’t become more visually complicated than necessary under the constraint. With only two outcomes, there’s no requirement for a boundary that simultaneously separates three or more basins, so the pattern looks much simpler than for degree 3 and higher.

What role does the unsolvability of the quintic play here?

The transcript mentions the unsolvability of the quintic to motivate why numerical methods matter for degree 5 and above. However, it stresses that the fractal complexity of Newton’s method is essentially unrelated to the quintic’s algebraic unsolvability. The fractal arises from iterative dynamics and basin boundaries, not from the impossibility of expressing roots by radicals.

How do Voronoi diagrams relate to Newton fractals?

If each point is colored by the closest root without iterating (zero Newton steps), the result is a Voronoi diagram with straight-line boundaries. After one Newton step, the coloring changes because each seed moves according to the Newton update before being assigned to the nearest root; with two steps, the pattern becomes more intricate. As the number of steps grows, these approximations converge toward the full Newton fractal boundary structure.

Review Questions

  1. In Newton’s method, what exact quantity determines the next iterate x_{n+1}, and how does p′(x) influence step size?
  2. What does it mean for the boundary of one color region to equal the boundary of all other color regions, and how does that prevent smooth boundary segments?
  3. Why does increasing the polynomial degree from 2 to 3 change the qualitative behavior of the Newton fractal’s boundaries?

Key Points

  1. 1

    Newton’s method updates a guess using x_{n+1}=x_n−p(x_n)/p′(x_n), repeatedly applying a local linear approximation to p(x).

  2. 2

    In the complex plane, coloring seeds by which root they converge to produces large uniform basins separated by fractal boundaries.

  3. 3

    Tiny perturbations to initial complex seeds can switch the final root, concentrating sensitivity on basin boundaries.

  4. 4

    A shared-boundary property implies the boundary of one root’s basin matches the boundary of every other basin, preventing neighborhoods from containing only two colors near boundary points.

  5. 5

    Because smooth boundary segments would typically separate only two outcomes locally, the shared-boundary constraint forces rough, infinitely detailed boundaries at every zoom level.

  6. 6

    Quadratic polynomials (two roots) look simpler because there are only two colors, so the shared-boundary constraint doesn’t force the same multi-color boundary complexity.

  7. 7

    The transcript attributes the underlying boundary behavior to holomorphic dynamics rather than to algebraic issues like quintic unsolvability.

Highlights

Newton’s method can look tame on a graph, but in the complex plane it generates a root-attraction “map” with infinitely detailed boundaries.
The fractal roughness isn’t just visual flair: the boundary of each basin is the same set for every root, ruling out smooth segments.
Voronoi diagrams appear as the zero-iteration starting point; each additional Newton step increases boundary complexity toward the full fractal.
The quintic’s unsolvability motivates numerical root-finding, but the Newton fractal’s complexity comes from iterative dynamics, not from radicals vs. no radicals.
Naming irony: Newton had no way to see these complex-plane fractals, yet his method produces them.

Topics

  • Newton’s Method
  • Complex Roots
  • Fractal Boundaries
  • Holomorphic Dynamics
  • Bezier Curves

Mentioned

  • Joseph Raphson