the contraction-mapping principle
The contraction-mapping principle is a clean abstract theorem with a surprising reach: it is the machine that produces and pins down the solution behind Picard-Lindelof. A map T is called a contraction if it always pulls points closer together — applying T to any two points shrinks the distance between them by at least a fixed factor k < 1: distance(T(x), T(y)) <= k times distance(x, y). The principle says that in a complete space, every contraction has exactly one fixed point, a point p with T(p) = p, and you reach it by iterating T from anywhere.
The picture is irresistible. Pick any starting point and apply T over and over: x, T(x), T(T(x)), and so on. Because each step cuts distances by the factor k, consecutive points get geometrically closer, the sequence cannot help but converge (this is where completeness — no holes in the space — is needed), and the limit must be a fixed point because T is continuous. Uniqueness is just as quick: if there were two fixed points, T would have to shrink the distance between them, but a fixed point keeps its own distance, so the only escape is that the distance was zero all along. One fixed point, reached constructively, error falling off geometrically.
The leap that makes this central to differential equations is that the 'points' need not be numbers — they can be entire functions, with distance measured as the largest gap between two functions over an interval. The Picard integral operator, y goes to y0 plus the integral of f(t, y(t)) dt, is exactly such a map, and a short enough interval makes it a contraction thanks to the Lipschitz condition. So existence, uniqueness, and the iterative construction all fall out of this one principle. The same theorem also underwrites the inverse and implicit function theorems and many numerical iterations — it is a genuine cornerstone, not an ODE-only trick.
The map T(x) = (x + 2/x)/2 on numbers near sqrt(2) is a contraction, and iterating it from x0 = 1 gives 1.5, 1.4167, 1.41422, ... rushing to the fixed point sqrt(2) where T(x) = x. The same principle, applied to functions instead of numbers, is what makes Picard iteration converge to the unique solution.
A contraction drags every start to one fixed point; with functions as the 'points,' that fixed point is the ODE solution.
Two hypotheses are doing the work and both are essential: the space must be complete (no missing limit points), and the factor k must be strictly less than 1. A map with k = 1 (an isometry) can fail to have any fixed point at all.