Worked solution: Tucker's rigorous computer-assisted proof via interval arithmetic (1999)
The eigenvalues of the Lorenz system linearised at the origin are three specific real numbers: one positive (, an outward direction) and two negative (, inward directions), with . Because these three numbers do not satisfy any accidental integer-linear relationship (a "resonance"), a classical theorem guarantees a smooth change of coordinates that makes the nonlinear system near the origin look exactly like three uncoupled, trivially solvable equations , , .
In these new coordinates, exactly how long a trajectory takes to enter and leave the small cube , and exactly where it exits, can be written down with explicit formulas rather than approximated — turning the single most dangerous part of the whole system into the only part solved with pencil and paper instead of a computer.
The Lorenz vector field linearised at the origin has three real eigenvalues satisfying for the classical parameters (with , , ), making the origin a saddle with a two-dimensional stable manifold and a one-dimensional unstable manifold. Because these eigenvalues satisfy a non-resonance condition (no small-integer linear combination of them vanishes, other than the trivial one), the Sternberg linearisation theorem guarantees a smooth (in fact, sufficiently differentiable) change of coordinates in a neighbourhood of the origin conjugating the nonlinear flow to its linear part .
Tucker makes this classical existence theorem effective: working within a small cube around the origin in these linearising coordinates, entry and exit times and positions through the faces of are computed by direct, closed-form integration of the linear system, with all approximation errors (from the higher-order terms neglected by the linearisation, and from the coordinate change itself) rigorously bounded using interval arithmetic on the relevant Taylor expansions.
This produces explicit, rigorously verified formulas describing exactly how the flow carries a small patch of initial conditions entering near the stable manifold of the origin and re-injects it, split into two pieces by the direction of the unstable manifold, back out into the region where the numerical Poincaré map of the next steps takes over.