A logbook for curious waters
← All Reckonings

The Daily Reckoning

How Does Mathematics Make ‘Close’ Exact?

The epsilon–delta definition turns an intuitive approach toward a limit into a guarantee that survives every requested tolerance.

Difficulty
Intermediate
Area
Calculus and analysis
Reading time
9 minutes
Prerequisites
Algebra · Absolute value as distance · Basic function notation
A luminous blue curve crosses nested acid-green horizontal tolerance bands and magenta vertical corridors on a dark maritime plotting grid.

Enter numbers increasingly close to 22 into the function f(x)=x2f(x)=x^2:

Input xxOutput f(x)f(x)
1.91.93.613.61
1.991.993.96013.9601
1.9991.9993.9960013.996001

The outputs appear to approach 44. A table can suggest that behavior, but no finite list can test every input near 22. How could a proof cover infinitely many nearby inputs at once?

The answer is a challenge-and-response game. A challenger names any acceptable output error, however tiny. The proof responds with an input radius that guarantees the requested accuracy. The Greek letters epsilon, written ε\varepsilon, and delta, written δ\delta, traditionally name those two positive numbers.

A guarantee with two distances

First consider a specific challenge. Suppose we want x2x^2 to land within 0.50.5 of 44. Choosing xx within 0.10.1 of 22 works: the condition 1.9<x<2.11.9<x<2.1 forces

3.61<x2<4.41,3.61<x^2<4.41,

so ∣x2−4∣<0.5|x^2-4|<0.5. Here the output tolerance is ε=0.5\varepsilon=0.5, and the input radius is δ=0.1\delta=0.1.

A limit proof must answer every positive tolerance, including tolerances far smaller than 0.50.5. For f(x)=x2f(x)=x^2, a dependable response is

δ=min⁡(1,ε5).\delta=\min\left(1,\frac{\varepsilon}{5}\right).

The notation min⁡(u,v)\min(u,v) means the smaller of the two numbers uu and vv. The fixed choice 11 keeps xx in a region where the curve cannot grow too steeply; the fraction ε/5\varepsilon/5 adjusts the input radius to the challenger’s requested accuracy.

The interactive graph below shows both distances. The green horizontal band contains outputs within ε\varepsilon of 44. The magenta vertical corridor contains inputs within δ\delta of 22. Move the test point inside and outside the corridor, or tighten the tolerance. A striped boundary and geometric shape accompany each color cue.

Interactive diagram

Turn an output tolerance into an input tolerance

ε output bandδ input windowf(x) = x²
Epsilon and delta windows for the limit of x squared at x equals twoA blue parabola crosses a horizontal epsilon band around y equals four and a vertical delta window around x equals two. A movable test point reports whether the epsilon-delta guarantee applies.
Proof's responseδ = min(1, ε/5) = 0.2
Input distance|x − 2| = 0.16
Output distance|x² − 4| = 0.666

Guarantee active: the point lies inside both tolerance windows.

What to notice: reducing ε tightens the striped horizontal band. The rule δ = min(1, ε/5) responds by tightening the vertical window. Every point on the blue curve inside that vertical window must land inside the horizontal band.

The precise definition

Let ff be a function. Let aa be the input that nearby values approach, and let LL be the output that f(x)f(x) appears to approach. We write

lim⁡x→af(x)=L\lim_{x\to a}f(x)=L

when the following statement is true:

For every ε>0\varepsilon>0, there exists a δ>0\delta>0 such that, for every input xx, if 0<∣x−a∣<δ0<|x-a|<\delta, then ∣f(x)−L∣<ε|f(x)-L|<\varepsilon.

OpenStax gives this same precise definition and develops its geometric meaning. Every symbol in the statement measures something specific:

  • ∣f(x)−L∣|f(x)-L| is the distance from the output f(x)f(x) to the proposed limit LL.
  • ε\varepsilon is the largest output error the challenger will accept.
  • ∣x−a∣|x-a| is the distance from the input xx to the approached input aa.
  • δ\delta is the input radius supplied by the proof.
  • The condition 0<∣x−a∣0<|x-a| excludes x=ax=a itself.

The order of the words matters. The challenger chooses ε\varepsilon first. The proof may choose δ\delta after seeing ε\varepsilon, so δ\delta may depend on ε\varepsilon, the function ff, the input aa, and the proposed limit LL. The proof must then handle every xx in the punctured interval 0<∣x−a∣<δ0<|x-a|<\delta; it cannot choose a different δ\delta after seeing an inconvenient xx.

Excluding the center lets a limit describe nearby behavior even when the function has a hole or an unusual value there. Define

g(x)={x2,x≠2,100,x=2.g(x)= \begin{cases} x^2, & x\ne 2,\\ 100, & x=2. \end{cases}

Although g(2)=100g(2)=100, every nearby value g(x)g(x) with x≠2x\ne2 equals x2x^2. Therefore

lim⁡x→2g(x)=4.\lim_{x\to2}g(x)=4.

A limit records where nearby outputs go. It does not, by itself, dictate the output at the approached input.

Worked example: prove that x2x^2 approaches 44

We will prove

lim⁡x→2x2=4.\lim_{x\to2}x^2=4.

Begin with an arbitrary number ε>0\varepsilon>0. We need to find a positive δ\delta that turns 0<∣x−2∣<δ0<|x-2|<\delta into ∣x2−4∣<ε|x^2-4|<\varepsilon.

Factoring the output error reveals two factors:

∣x2−4∣=∣(x−2)(x+2)∣=∣x−2∣ ∣x+2∣.|x^2-4|=|(x-2)(x+2)|=|x-2|\,|x+2|.

The first factor is exactly the input distance that we can control. To control the second factor, require δ≤1\delta\le1. Then ∣x−2∣<δ≤1|x-2|<\delta\le1 implies

−1<x−2<1,-1<x-2<1,

and adding 22 to all three parts gives

1<x<3.1<x<3.

Adding 22 again yields 3<x+2<53<x+2<5, so ∣x+2∣<5|x+2|<5. Now choose

δ=min⁡(1,ε5).\delta=\min\left(1,\frac{\varepsilon}{5}\right).

This choice is positive because ε>0\varepsilon>0. If 0<∣x−2∣<δ0<|x-2|<\delta, then the earlier bound ∣x+2∣<5|x+2|<5 applies, and

∣x2−4∣=∣x−2∣ ∣x+2∣<5∣x−2∣<5δ≤5(ε5)=ε.\begin{aligned} |x^2-4| &=|x-2|\,|x+2|\\ &<5|x-2|\\ &<5\delta\\ &\le5\left(\frac{\varepsilon}{5}\right)\\ &=\varepsilon. \end{aligned}

The final inequality uses δ≤ε/5\delta\le\varepsilon/5, which follows from the definition of min⁡\min. Because the argument works for every positive ε\varepsilon, the limit is proved.

The response δ=min⁡(1,ε/5)\delta=\min(1,\varepsilon/5) is sufficient rather than uniquely best. A smaller positive radius would also work. Epsilon–delta proofs ask for a guarantee, not the widest possible interval.

Try the challenge

Problem 1

Use the epsilon–delta definition to prove

lim⁡x→3(2x+1)=7.\lim_{x\to3}(2x+1)=7.
Reveal the solution

Let ε>0\varepsilon>0 be arbitrary. The output error is

∣(2x+1)−7∣=∣2x−6∣=2∣x−3∣.|(2x+1)-7|=|2x-6|=2|x-3|.

Choose δ=ε/2\delta=\varepsilon/2. This number is positive. If 0<∣x−3∣<δ0<|x-3|<\delta, then

∣(2x+1)−7∣=2∣x−3∣<2δ=2(ε2)=ε.|(2x+1)-7|=2|x-3|<2\delta=2\left(\frac{\varepsilon}{2}\right)=\varepsilon.

Therefore lim⁡x→3(2x+1)=7\lim_{x\to3}(2x+1)=7.

Problem 2

Use the epsilon–delta definition to prove

lim⁡x→3x2=9.\lim_{x\to3}x^2=9.
Reveal the solution

Let ε>0\varepsilon>0 be arbitrary. Factor the output error:

∣x2−9∣=∣x−3∣ ∣x+3∣.|x^2-9|=|x-3|\,|x+3|.

Choose

δ=min⁡(1,ε7).\delta=\min\left(1,\frac{\varepsilon}{7}\right).

If 0<∣x−3∣<δ0<|x-3|<\delta, then ∣x−3∣<1|x-3|<1, so 2<x<42<x<4. Consequently, 5<x+3<75<x+3<7 and ∣x+3∣<7|x+3|<7. Therefore

∣x2−9∣=∣x−3∣ ∣x+3∣<7∣x−3∣<7δ≤7(ε7)=ε.\begin{aligned} |x^2-9| &=|x-3|\,|x+3|\\ &<7|x-3|\\ &<7\delta\\ &\le7\left(\frac{\varepsilon}{7}\right)\\ &=\varepsilon. \end{aligned}

The argument proves lim⁡x→3x2=9\lim_{x\to3}x^2=9.

How “as close as desired” became mathematics

The modern definition did not appear in one person’s notebook on one date. It emerged as mathematicians tried to replace persuasive geometric pictures and infinitesimal quantities with statements that could survive every exceptional case.

In 1817, the Prague mathematician and philosopher Bernard Bolzano published a “purely analytic” proof of a theorem about continuous functions. His work pursued rigor without relying on geometric intuition, but it circulated poorly and was probably unknown to the mathematicians who soon made similar advances. The MacTutor biography of Bolzano places his mathematics inside Habsburg rule in Bohemia: after Bolzano advocated pacifism, social equality, and intellectual freedom, the Austrian government removed him from his university post in 1819, censored his mail, and restricted his publication.

Augustin-Louis Cauchy’s 1821 textbook Cours d’analyse defined a limit in verbal terms, through variable values approaching a fixed value without end. A public-domain scan of the 1821 book preserves the primary text. Cauchy used both ε\varepsilon and δ\delta in an 1823 proof, although he did not yet express their dependency with the modern quantifier structure. Karl Weierstrass systematized that structure in his 1861 lectures. Historian G. I. Sinkevich traces the transition from verbal limits to full “epsilontics”, while Judith V. Grabiner’s study reconstructs the longer route to Cauchy’s rigorous calculus.

Politics touched Cauchy’s career too. He gained an Académie des Sciences seat in 1816 after the restored Bourbon government removed Lazare Carnot and Gaspard Monge. Following the July Revolution of 1830, Cauchy refused to swear allegiance to the new regime and lost academic posts that required the oath. The MacTutor biography of Cauchy documents both turns. Bolzano’s suppression and Cauchy’s shifting access to institutions did not determine whether their proofs were correct. Those events did affect careers, audiences, and which work traveled easily enough to be read.

Where exact tolerances matter

Epsilon–delta language anchors the modern theory of continuity, derivatives, and integrals. The same pattern also appears whenever a calculation must meet a declared error budget: state the permitted output error, then determine what accuracy the inputs or intermediate steps require.

Numerical simulation software uses this pattern operationally. An adaptive time-stepping program estimates the error produced by one simulated step and compares that estimate with a specified tolerance. The program accepts, rejects, or resizes the step in response. The National Institute of Standards and Technology guide to numerical implementation describes that error-control loop.

Iterative algorithms create a sequence of improving approximations. For example, Newton’s method estimates 2\sqrt2 by starting with a positive number x0x_0 and repeatedly calculating

xn+1=12(xn+2xn),x_{n+1}=\frac12\left(x_n+\frac{2}{x_n}\right),

where nn counts completed steps. Starting from x0=1.5x_0=1.5, three repetitions produce an estimate whose error is about 1.6×10−121.6\times10^{-12}. The four number lines below use a different scale in each row, so each row can show both the current estimate and a much tighter permitted error.

Numerical convergence

Each calculation earns a tighter error budget

Four Newton-method estimates approaching the square root of twoFour separately magnified number lines show successive estimates moving closer to the square root of two while the allowed error shrinks from one tenth to one billionth.Startx = 1.5√2tolerance: 0.1error: 0.0858After 1 stepx = 1.416666666667√2tolerance: 0.01error: 0.00245After 2 stepsx = 1.414215686275√2tolerance: 0.00001error: 0.00000212After 3 stepsx = 1.414213562375√2tolerance: 1 × 10⁻⁹error: 1.59e-12
What to notice: each row magnifies a different interval centered on √2. The pink dot marks the current estimate; the green center mark is the target. Three Newton steps reduce the error from about 0.086 to roughly 1.6 × 10−12, comfortably inside the final tolerance.

Production solvers turn the visual idea into a stopping rule. For example, the National Institute of Standards and Technology documentation for a conjugate-gradient solver describes repeated approximations that stop when a relative error falls below a chosen tolerance.

Speculative application: a future formally verified artificial-intelligence controller could attach mathematical error certificates to decisions made near a safety boundary. Such a certificate might guarantee that sensor uncertainty below a computed δ\delta changes a risk score by less than a regulator’s chosen ε\varepsilon. Building a trustworthy physical model would remain the hard part; epsilon–delta reasoning would supply the form of the guarantee.

One idea to remember

To prove that f(x)f(x) approaches LL as xx approaches aa, accept every positive output tolerance ε\varepsilon and answer it with a positive input radius δ\delta that works for every nearby xx.

Yesterday’s Reckoning, “How Can Motion Have a Speed at One Instant?”, used a limit to define instantaneous velocity. Today’s guarantee supplies the missing logical machinery. Natural directions from here are continuity, which asks whether a function’s value agrees with its limit, and sequential limits, which test approach using every convergent sequence of inputs.

Sources