Integrability of Continuous Functions

Suppose you want to compute the area under a wildly wiggling curve. An example is f(x)=sin⁡(1/x)f(x) = \sin(1/x) patched to be continuous at 00. The interval considered is [a,b][a,b]. Can you always trust that "the area" is a well-defined number? It should be computable as a limit of rectangle sums. This should hold no matter how erratic the curve looks. This question sits at the very foundation of integral calculus. Before one can compute an integral with the Fundamental Theorem of Calculus, one must know it exists. This existence must hold in a rigorous sense. The Cauchy Integrability Criterion gives a clean test. A bounded function is Riemann integrable on [a,b][a,b] exactly when the following holds. For every ε>0\varepsilon>0 there is a partition making a certain difference small. That difference is between upper and lower Riemann sums, and it must be smaller than ε\varepsilon. This criterion is powerful but abstract. It does not, by itself, tell us which functions actually satisfy it. The theorem in this article closes that gap. It does so for the single most important class of functions in analysis: continuous ones. The key insight involves a closed, bounded interval. On such an interval, continuity upgrades automatically to Uniform Continuity on Compact Sets. Uniform continuity is precisely the tool needed to control the oscillation of ff. It controls this oscillation simultaneously across an entire partition. This is stronger than controlling it just near each individual point. This is the mechanism that converts "ff is continuous" into something stronger. It becomes the fact that "the upper and lower sums can be squeezed together." This article gives a complete, self-contained proof. It shows that every continuous function on [a,b][a,b] is Riemann integrable. It builds directly on the Cauchy criterion and uniform continuity. It also sets the stage for deeper theorems that follow. These include the Mean Value Theorem for Integrals. They also include the First Fundamental Theorem of Calculus. Both of these presuppose something important. They assume the integral of a continuous function is a well-defined real number to begin with.

Key takeaways
  • Every f∈C([a,b])f \in C([a,b]) is bounded (by the Extreme Value Theorem) and hence a candidate for Riemann integration.
  • The Cauchy Integrability Criterion reduces integrability to making U(f,P)−L(f,P)U(f,P) - L(f,P) arbitrarily small.
  • Continuity on the compact interval [a,b][a,b] is automatically uniform continuity.
  • Uniform continuity lets one choose a single δ\delta controlling oscillation of ff over every subinterval of a fine-enough partition.
  • Taking a partition with mesh less than that δ\delta forces U(f,P)−L(f,P)<εU(f,P) - L(f,P) < \varepsilon, proving integrability.
  • Key idea: uniform continuity turns local wiggle-control into global, partition-wide wiggle-control, which is exactly what the Cauchy criterion demands.

Where this lives today

The integrability of continuous functions is the load-bearing beam beneath the entire calculus curriculum: the First Fundamental Theorem of Calculus and the Mean Value Theorem for Integrals both presuppose that ∫abf\int_a^b f is a well-defined real number before differentiating or estimating it. In numerical analysis, quadrature methods (Simpson's rule, trapezoidal rule) are justified by exactly this kind of uniform-continuity error control — the same δ\delta-vs-mesh argument bounds numerical quadrature error. In probability and statistics, continuity of density functions on bounded supports guarantees that expectations E[g(X)]=∫g(x)fX(x) dx\mathbb{E}[g(X)] = \int g(x) f_X(x)\,dx are well-posed integrals rather than formal expressions. In physics and engineering, continuity of physical fields (temperature, potential, stress) over bounded regions is routinely invoked to justify that line and volume integrals used in energy or work computations actually exist. In theoretical computer science, uniform continuity estimates underlie rigorous error bounds for algorithms that approximate integrals or continuous functions by discretized samples, connecting this 19th-century theorem directly to modern computational guarantees.

:::quiz[Which property of [a,b][a,b] is essential to upgrade continuity to uniform continuity?]

  • [a,b][a,b] is an interval of real numbers
  • [a,b][a,b] is compact (closed and bounded) ✓
  • [a,b][a,b] contains rational numbers
  • [a,b][a,b] has positive length

:::quiz[What breaks in the theorem if ff is continuous on (0,1](0,1] but domain is not compact, as with f(x)=sin⁡(1/x)f(x) = \sin(1/x)?]

  • Nothing breaks; the theorem still applies
  • ff fails to be bounded near 00
  • Continuity no longer guarantees uniform continuity, so the oscillation on subintervals near 00 cannot be uniformly controlled ✓
  • The Cauchy criterion becomes false in general

:::faq[Does uniform continuity alone guarantee integrability, or is compactness of [a,b][a,b] still doing work elsewhere in the proof?] Compactness of [a,b][a,b] is used twice: first, to upgrade continuity to uniform continuity (Step 2), and second, to guarantee that ff attains its supremum and infimum on each compact subinterval [xi−1,xi][x_{i-1},x_i] via the Extreme Value Theorem (Step 5), which is what lets us write Mi−miM_i - m_i as an actual difference of function values f(ui)−f(vi)f(u_i)-f(v_i) rather than merely a difference of suprema/infima that might not be attained. Both uses trace back to the same underlying compactness of closed bounded intervals in R\mathbb{R}.

Ask the tutor: How does this integrability proof change if ff is allowed finitely many jump discontinuities on [a,b][a,b]?

DefinitionRiemann Integrability of a Bounded Function

Let f:[a,b]→Rf:[a,b]\to\mathbb{R} be bounded. For a partition P={a=x0<x1<⋯<xn=b}P = \{a = x_0 < x_1 < \cdots < x_n = b\} of [a,b][a,b], define

Mi=sup⁡x∈[xi−1,xi]f(x),mi=inf⁡x∈[xi−1,xi]f(x),M_i = \sup_{x \in [x_{i-1},x_i]} f(x), \qquad m_i = \inf_{x \in [x_{i-1},x_i]} f(x),

and the upper and lower Riemann sums

U(f,P)=∑i=1nMi(xi−xi−1),L(f,P)=∑i=1nmi(xi−xi−1).U(f,P) = \sum_{i=1}^n M_i (x_i - x_{i-1}), \qquad L(f,P) = \sum_{i=1}^n m_i (x_i - x_{i-1}).

Define the upper integral ∫ab‾f=inf⁡PU(f,P)\overline{\int_a^b} f = \inf_P U(f,P) and the lower integral ∫ab‾f=sup⁡PL(f,P)\underline{\int_a^b} f = \sup_P L(f,P), infimum and supremum taken over all partitions PP of [a,b][a,b]. The function ff is Riemann integrable on [a,b][a,b] if

∫ab‾f=∫ab‾f,\overline{\int_a^b} f = \underline{\int_a^b} f,

and this common value is denoted ∫abf(x) dx\int_a^b f(x)\,dx.

Intuition

Picture the graph of a continuous function over [a,b][a,b], and imagine slicing [a,b][a,b] into thin vertical strips. Over each strip [xi−1,xi][x_{i-1}, x_i], the function's values range between a minimum mim_i and maximum MiM_i. The upper sum U(f,P)U(f,P) stacks up rectangles that overestimate the area (using the tallest value on each strip); the lower sum L(f,P)L(f,P) underestimates (using the shortest value). The gap U(f,P)−L(f,P)U(f,P) - L(f,P) measures how much "slop" the partition leaves. Now here is the crucial fact that continuity supplies: because ff is continuous on a closed, bounded interval, it is uniformly continuous. This means that if two points x,yx,y are close together (within some δ\delta, chosen once and for all, not depending on where x,yx,y sit in [a,b][a,b]), then f(x)f(x) and f(y)f(y) are guaranteed to be close together (within ε\varepsilon). If every strip of the partition has width less than this δ\delta, then ff cannot vary by more than ε\varepsilon within any single strip — so Mi−mi≤εM_i - m_i \le \varepsilon for every ii simultaneously. This is exactly the uniform control needed to make U(f,P)−L(f,P)U(f,P)-L(f,P) small, strip by strip, all at once. Without uniform continuity — if δ\delta were allowed to depend on the location — some strip somewhere might harbor huge oscillation no matter how fine the partition, and the gap between upper and lower sums could refuse to shrink. Uniform continuity is what forecloses that possibility.

Figure
Upper and lower sums squeezing together as the partition refines
f(x) x a b M_i − m_i ≤ ε on every strip of width < δ U(f,P) (overlaid) L(f,P)

:::theorem[Continuous Functions on [a,b][a,b] are Riemann Integrable] Let f:[a,b]→Rf:[a,b] \to \mathbb{R} be continuous on the closed, bounded interval [a,b][a,b], where a<ba < b. Then ff is Riemann integrable on [a,b][a,b].

Proof

Technique: direct construction of a partition via uniform continuity, verified against the Cauchy Integrability Criterion.

Step 1: ff is bounded, so upper and lower sums are well-defined. Since [a,b][a,b] is closed and bounded, it is compact. A continuous real-valued function on a compact set is bounded (the Extreme Value Theorem, itself a standard consequence of compactness). Hence there exists B>0B > 0 with ∣f(x)∣≤B|f(x)| \le B for all x∈[a,b]x \in [a,b]. Boundedness guarantees that for every partition PP, each supremum Mi=sup⁡[xi−1,xi]fM_i = \sup_{[x_{i-1},x_i]} f and infimum mi=inf⁡[xi−1,xi]fm_i = \inf_{[x_{i-1},x_i]} f exists as a finite real number, so U(f,P)U(f,P) and L(f,P)L(f,P) are well-defined finite sums, and the upper and lower integrals exist as finite real numbers (they are, respectively, the infimum of a set of reals bounded below and the supremum of a set of reals bounded above, since L(f,P)≤U(f,P′)L(f,P) \le U(f,P') for all partitions P,P′P, P' — a standard fact underlying the definition of upper/lower integrals).

Step 2: Continuity on [a,b][a,b] upgrades to uniform continuity. By the theorem on Uniform Continuity on Compact Sets, since [a,b][a,b] is compact and ff is continuous at every point of [a,b][a,b], ff is uniformly continuous on [a,b][a,b]. Explicitly: for every ε>0\varepsilon > 0 there exists δ>0\delta > 0 such that for all x,y∈[a,b]x, y \in [a,b],

∣x−y∣<δ  ⟹  ∣f(x)−f(y)∣<ε.|x - y| < \delta \implies |f(x) - f(y)| < \varepsilon.

Crucially, this single δ\delta works for every pair x,yx,y in [a,b][a,b] — it does not depend on the location of xx or yy within the interval.

Step 3: Invoke the Cauchy Integrability Criterion. By the Cauchy Integrability Criterion, to show ff is Riemann integrable it suffices to show:

∀ε>0, ∃ a partition P of [a,b] such that U(f,P)−L(f,P)<ε.\forall \varepsilon > 0,\ \exists \text{ a partition } P \text{ of } [a,b] \text{ such that } U(f,P) - L(f,P) < \varepsilon.

Fix ε>0\varepsilon > 0. We now construct such a partition.

Step 4: Choose δ\delta and build a partition with mesh less than δ\delta. Apply Step 2 with tolerance ε′=εb−a\varepsilon' = \dfrac{\varepsilon}{b-a} (a valid choice since ε′>0\varepsilon' > 0, because b>ab > a). Obtain δ>0\delta > 0 such that

x,y∈[a,b], ∣x−y∣<δ  ⟹  ∣f(x)−f(y)∣<εb−a.x, y \in [a,b],\ |x-y| < \delta \implies |f(x) - f(y)| < \frac{\varepsilon}{b-a}.

Let n∈Nn \in \mathbb{N} be large enough that b−an<δ\frac{b-a}{n} < \delta (such nn exists by the Archimedean property). Let P={x0,x1,…,xn}P = \{x_0, x_1, \ldots, x_n\} be the partition of [a,b][a,b] into nn equal subintervals:

xi=a+i⋅b−an,i=0,1,…,n.x_i = a + i \cdot \frac{b-a}{n}, \qquad i = 0, 1, \ldots, n.

Each subinterval [xi−1,xi][x_{i-1}, x_i] has length b−an<δ\frac{b-a}{n} < \delta.

Step 5: Bound Mi−miM_i - m_i on each subinterval. Fix i∈{1,…,n}i \in \{1, \ldots, n\}. The subinterval [xi−1,xi][x_{i-1}, x_i] is closed and bounded, hence compact, and ff restricted to it is continuous, so ff attains a maximum and a minimum on [xi−1,xi][x_{i-1},x_i] (Extreme Value Theorem): there exist ui,vi∈[xi−1,xi]u_i, v_i \in [x_{i-1}, x_i] with

f(ui)=Mi=sup⁡[xi−1,xi]f,f(vi)=mi=inf⁡[xi−1,xi]f.f(u_i) = M_i = \sup_{[x_{i-1},x_i]} f, \qquad f(v_i) = m_i = \inf_{[x_{i-1},x_i]} f.

Since ui,vi∈[xi−1,xi]u_i, v_i \in [x_{i-1}, x_i], we have ∣ui−vi∣≤xi−xi−1=b−an<δ|u_i - v_i| \le x_i - x_{i-1} = \frac{b-a}{n} < \delta. By the uniform continuity bound from Step 4,

Mi−mi=f(ui)−f(vi)=∣f(ui)−f(vi)∣<εb−a.M_i - m_i = f(u_i) - f(v_i) = |f(u_i) - f(v_i)| < \frac{\varepsilon}{b-a}.

(The middle equality holds because Mi≥miM_i \ge m_i, so f(ui)−f(vi)≥0f(u_i) - f(v_i) \ge 0, and thus it equals its own absolute value.)

Step 6: Sum the estimate over all subintervals.

U(f,P)−L(f,P)=∑i=1n(Mi−mi)(xi−xi−1)<∑i=1nεb−a⋅b−an=∑i=1nεn=ε.U(f,P) - L(f,P) = \sum_{i=1}^n (M_i - m_i)(x_i - x_{i-1}) < \sum_{i=1}^n \frac{\varepsilon}{b-a} \cdot \frac{b-a}{n} = \sum_{i=1}^n \frac{\varepsilon}{n} = \varepsilon.

Step 7: Conclude. We have exhibited, for arbitrary ε>0\varepsilon > 0, a partition PP of [a,b][a,b] with U(f,P)−L(f,P)<εU(f,P) - L(f,P) < \varepsilon. By the Cauchy Integrability Criterion, ff is Riemann integrable on [a,b][a,b]. ■\blacksquare

∎
CorollaryThe Integral is Sandwiched by Every Refining Partition Sequence

If (Pn)(P_n) is any sequence of partitions of [a,b][a,b] with mesh ∥Pn∥→0\|P_n\| \to 0, and ff is continuous on [a,b][a,b], then

U(f,Pn)→∫abf(x) dxandL(f,Pn)→∫abf(x) dx.U(f, P_n) \to \int_a^b f(x)\,dx \quad \text{and} \quad L(f, P_n) \to \int_a^b f(x)\,dx.

This follows immediately by rerunning Step 4–6 of the proof with δ\delta chosen relative to any ε>0\varepsilon > 0, and noting that once ∥Pn∥<δ\|P_n\| < \delta, the same oscillation bound Mi−mi<ε/(b−a)M_i - m_i < \varepsilon/(b-a) applies to every subinterval of PnP_n (not just equal-width ones), giving U(f,Pn)−L(f,Pn)<εU(f,P_n) - L(f,P_n) < \varepsilon for all sufficiently large nn.

Worked examples

Example 1. Let f(x)=x3−4xf(x) = x^3 - 4x on [−1,3][-1, 3]. This is a polynomial, hence continuous on R\mathbb{R}, hence continuous on the compact interval [−1,3][-1,3]; by the theorem, ff is Riemann integrable on [−1,3][-1,3]. Let us carry out the ε\varepsilon-construction explicitly for ε=0.6\varepsilon = 0.6, to see the machine work on an asymmetric, non-monotonic function. First bound the derivative to get an explicit modulus of continuity: for x∈[−1,3]x \in [-1,3],

∣f′(x)∣=∣3x2−4∣≤max⁡(∣3(−1)2−4∣, ∣3(3)2−4∣)=max⁡(1,23)=23,|f'(x)| = |3x^2 - 4| \le \max(|3(-1)^2 - 4|,\ |3(3)^2-4|) = \max(1, 23) = 23,

since 3x2−43x^2 - 4 is increasing in ∣x∣|x| on this interval and ∣−1∣<∣3∣|{-1}| < |3|. So by the Mean Value Theorem, for any x,y∈[−1,3]x,y \in [-1,3],

∣f(x)−f(y)∣≤23∣x−y∣.|f(x) - f(y)| \le 23|x-y|.

(This is a Lipschitz bound, a stronger and more explicit version of uniform continuity available here because ff is differentiable with bounded derivative; it is not needed in general but illustrates the mechanism concretely.) To force ∣f(x)−f(y)∣<ε′=0.63−(−1)=0.15|f(x)-f(y)| < \varepsilon' = \frac{0.6}{3-(-1)} = 0.15, it suffices to take δ=0.1523≈0.006522\delta = \frac{0.15}{23} \approx 0.006522. The interval has length b−a=4b - a = 4. Choose nn so that 4n<δ\frac{4}{n} < \delta, i.e. n>40.006522≈613.3n > \frac{4}{0.006522} \approx 613.3, so take n=614n = 614. Partitioning [−1,3][-1,3] into 614614 equal pieces of width 4614≈0.006515<δ\frac{4}{614} \approx 0.006515 < \delta guarantees, by Step 5 of the proof,

Mi−mi<0.15for every i,M_i - m_i < 0.15 \quad \text{for every } i,

and hence

U(f,P)−L(f,P)<∑i=16140.15⋅4614=0.15⋅4=0.6=ε.U(f,P) - L(f,P) < \sum_{i=1}^{614} 0.15 \cdot \frac{4}{614} = 0.15 \cdot 4 = 0.6 = \varepsilon.

This confirms concretely that a fine enough uniform partition squeezes the upper and lower sums together, exactly as the general theorem promises — no cleverness about where the partition points sit is required, only that the mesh is small enough.

Example 2. Let g(x)=xg(x) = \sqrt{x} on [0,5][0, 5]. This function is continuous on [0,5][0,5] (in particular, right-continuous at 00, since lim⁡x→0+x=0=0\lim_{x\to 0^+}\sqrt x = 0 = \sqrt 0), so the theorem applies and gg is Riemann integrable on [0,5][0,5], even though gg is not differentiable at x=0x=0 and has unbounded derivative near 00 (so the Lipschitz shortcut of Example 1 fails here). This shows why the theorem's proof must genuinely use uniform continuity rather than a derivative bound. To see uniform continuity directly: given ε>0\varepsilon > 0, we want δ>0\delta > 0 with ∣x−y∣<ε|\sqrt x - \sqrt y| < \varepsilon whenever x,y∈[0,5]x, y \in [0,5], ∣x−y∣<δ|x-y| < \delta. Using the identity x−y=x−yx+y\sqrt{x} - \sqrt{y} = \frac{x - y}{\sqrt x + \sqrt y} for x,y≥0x, y \ge 0 not both zero, and separately checking the case where one of them is 00:

  • If x,y∈[0,5]x, y \in [0,5] and (say) y=0y = 0: ∣x−0∣=x<ε  ⟺  x<ε2|\sqrt x - \sqrt 0| = \sqrt x < \varepsilon \iff x < \varepsilon^2.
  • If x,y>0x, y > 0: one shows directly ∣x−y∣≤∣x−y∣|\sqrt x - \sqrt y| \le \sqrt{|x-y|} (squaring both sides of x≤y+x−y\sqrt x \le \sqrt y + \sqrt{x-y} when x≥yx \ge y, which holds since (y+x−y)2=y+(x−y)+2y(x−y)≥x(\sqrt y + \sqrt{x-y})^2 = y + (x-y) + 2\sqrt{y(x-y)} \ge x).

So in all cases ∣x−y∣≤∣x−y∣|\sqrt x - \sqrt y| \le \sqrt{|x - y|} for x,y≥0x, y \ge 0. Given ε=0.4\varepsilon = 0.4, choose δ=ε2=0.16\delta = \varepsilon^2 = 0.16; then ∣x−y∣<0.16  ⟹  ∣x−y∣≤0.16=0.4|x - y| < 0.16 \implies |\sqrt x - \sqrt y| \le \sqrt{0.16} = 0.4. Take ε′=0.45=0.08\varepsilon' = \frac{0.4}{5} = 0.08, needing δ=0.082=0.0064\delta = 0.08^2 = 0.0064. Choosing nn with 5n<0.0064\frac{5}{n} < 0.0064, i.e. n>781.25n > 781.25, so n=782n = 782 subintervals of width 5782≈0.006394<δ\frac{5}{782} \approx 0.006394 < \delta, forces

U(g,P)−L(g,P)<∑i=17820.08⋅5782=0.08⋅5=0.4.U(g,P) - L(g,P) < \sum_{i=1}^{782} 0.08 \cdot \frac{5}{782} = 0.08 \cdot 5 = 0.4.

Again the Cauchy criterion is met, confirming integrability at the level of an explicit partition — and demonstrating that the square-root modulus of continuity (rather than a linear Lipschitz bound) still suffices, because the theorem only ever needs some uniform continuity estimate, of whatever shape.

CounterexampleDropping compactness of the domain breaks the theorem

The theorem requires [a,b][a,b] to be a closed, bounded interval. Consider f:(0,1]→Rf:(0,1] \to \mathbb{R}, f(x)=1xf(x) = \frac{1}{x}, which is continuous on (0,1](0,1] but (0,1](0,1] is not compact (not closed). Then:

  • ff is unbounded on (0,1](0,1]: for any B>0B > 0, taking x=min⁡(1,12B)x = \min(1, \frac{1}{2B}) gives f(x)≥Bf(x) \ge B. Since Riemann integrability requires boundedness by definition, ff cannot be Riemann integrable on (0,1](0,1] in the ordinary sense — for any partition P={0<x1<⋯<xn=1}P = \{0 < x_1 < \cdots < x_n = 1\} attempting to include the endpoint 00, we'd have sup⁡[0,x1]f=+∞\sup_{[0,x_1]} f = +\infty, so U(f,P)=+∞U(f,P) = +\infty for every partition, and the upper integral is not a finite real number.
  • Uniform continuity also genuinely fails here (not merely boundedness): even if one instead considers a bounded but still non-uniformly-continuous example, such as h(x)=sin⁡(1/x)h(x) = \sin(1/x) on (0,1](0,1] (bounded, continuous, but not uniformly continuous, since near 00 arbitrarily close points x,yx,y can have h(x),h(y)h(x), h(y) far apart), the proof's Step 4–5 construction breaks: no single δ\delta controls oscillation on subintervals near 00, however the partition is refined elsewhere. Concretely, for xn=12nπ+π/2x_n = \frac{1}{2n\pi + \pi/2} and yn=12nπ+3π/2y_n = \frac{1}{2n\pi + 3\pi/2}, we get xn,yn→0x_n, y_n \to 0 (so ∣xn−yn∣→0|x_n - y_n| \to 0) while h(xn)=1h(x_n) = 1 and h(yn)=−1h(y_n) = -1 for all nn, so ∣h(xn)−h(yn)∣=2|h(x_n) - h(y_n)| = 2 never shrinks — violating the uniform continuity conclusion the proof relies on.

This shows both hypotheses "closed" and "bounded" (packaged together as compactness) are load-bearing: remove either, and continuity no longer forces uniform continuity, and the proof's central mechanism collapses.

Scope, edge cases & common traps

SituationDoes the theorem apply?Why
ff continuous on [a,b][a,b], a<ba<bYesExactly the hypothesis of the theorem
ff continuous on (a,b)(a,b) openNot directly(a,b)(a,b) not compact; uniform continuity can fail (e.g. f(x)=1/xf(x)=1/x on (0,1)(0,1))
ff continuous on [a,b][a,b] except at finitely many jump pointsYes, but needs a separate theoremHandled by a refinement (piecewise-continuous / bounded-with-finite-discontinuities integrability), not this theorem alone
ff continuous but domain is [a,∞)[a,\infty)NoUnbounded domain is not compact; improper integrals require separate limiting definitions
ff Lipschitz on [a,b][a,b]YesLipschitz   ⟹  \implies uniformly continuous   ⟹  \implies theorem applies (stronger hypothesis, same conclusion)
ff merely bounded, not continuous, on [a,b][a,b]Not covered by this theoremNeed the Cauchy criterion checked directly, or another sufficient condition (e.g. monotonicity)

The mathematician behind it

Origins · 1 minBernhard Riemann and how this was discovered
BR
Bernhard Riemann
1826–1866 · Göttingen, mid-19th century analysis

Bernhard Riemann formalized the modern notion of the definite integral in his 1854 habilitation work, seeking rigorous foundations for representing functions by trigonometric (Fourier) series. His construction — via upper and lower sums over partitions — replaced looser, geometrically-motivated notions of "area under a curve" with an arithmetic criterion, later sharpened by Darboux into the upper/lower sum formulation used throughout this article. Riemann's broader legacy spans complex analysis, differential geometry (Riemannian manifolds), and the Riemann Hypothesis in number theory, but the integral bearing his name remains the entry point through which most students first encounter rigorous analysis.

History
- 1854 — Riemann introduces a rigorous definition of the definite integral in his habilitation thesis on trigonometric series.
- 1875 — Gaston Darboux reformulates Riemann's criterion using upper and lower sums, giving the version used in modern textbooks.
- 1854–1870s — Uniform continuity is isolated and named as a distinct, stronger notion than pointwise continuity, notably in the work of Eduard Heine.
- 1872 — Heine proves that continuous functions on compact intervals are uniformly continuous, supplying the exact tool this theorem depends on.
- 1900s–present — The Riemann integral becomes the standard first rigorous integral in undergraduate analysis, later supplemented (not replaced) by the more general Lebesgue integral.
Why isn't ordinary (pointwise) continuity enough to prove integrability?

Pointwise continuity only guarantees that for each fixed point x0x_0 and each ε>0\varepsilon>0, there is some δ\delta (possibly depending on x0x_0) controlling nearby values. To build one partition that controls oscillation everywhere, we need a single δ\delta that works uniformly across the whole interval, so that a mesh smaller than δ\delta controls every subinterval at once. That upgrade — from a δ\delta depending on the point to a δ\delta depending only on ε\varepsilon — is precisely uniform continuity, and it holds automatically on compact intervals like [a,b][a,b], but can fail on open or unbounded domains, as the counterexample sin⁡(1/x)\sin(1/x) on (0,1](0,1] shows.

Does this theorem also apply to functions continuous except at a few points?

Not directly — this theorem's hypothesis is continuity at every point of [a,b][a,b]. However, a closely related and provable extension states that any bounded function with only finitely many (or even countably many, in measure-zero settings) discontinuities on [a,b][a,b] is still Riemann integrable. The proof strategy is similar in spirit: isolate small intervals around the bad points to control their contribution, and use uniform continuity of ff away from those points to control the rest. That extension is a separate theorem, not automatically covered by the pure continuous case proven here.

Is boundedness of ff a separate hypothesis, or does continuity give it for free?

On a compact interval [a,b][a,b], boundedness is automatic: continuity plus compactness of the domain implies boundedness by the Extreme Value Theorem, so it need not be assumed separately. This is a special feature of compact domains — on a non-compact domain like (0,1](0,1], a continuous function such as f(x)=1/xf(x) = 1/x can be unbounded, which is exactly why the Riemann integral (which requires bounded integrands by definition) cannot apply there without further care, such as an improper integral.

How does this theorem relate to the Fundamental Theorem of Calculus?

The Fundamental Theorem of Calculus asserts that if ff is continuous on [a,b][a,b], then F(x)=∫axf(t) dtF(x) = \int_a^x f(t)\,dt is differentiable with F′(x)=f(x)F'(x) = f(x). But this statement is meaningless unless ∫axf(t) dt\int_a^x f(t)\,dt is already known to exist as a well-defined real number for every x∈[a,b]x \in [a,b] — which is precisely what this integrability theorem supplies. In this sense, the integrability theorem is a logical prerequisite that must be established before the Fundamental Theorem of Calculus can even be stated meaningfully, let alone proved.

Can the proof's partition be non-uniform (unequal subinterval widths)?

Yes — the proof only needs the mesh (the maximum subinterval width) to be smaller than δ\delta; it does not require equal widths. Using an equal-width partition is a convenient, concrete choice that makes the arithmetic in Step 6 clean, but any partition, however irregular, with all subintervals shorter than δ\delta achieves the same bound U(f,P)−L(f,P)<εU(f,P) - L(f,P) < \varepsilon. This flexibility is used implicitly in the corollary about refining sequences of partitions with mesh tending to zero.

What is the connecting idea to the Mean Value Theorem for Integrals?

Once integrability is established, the Mean Value Theorem for Integrals states that for ff continuous on [a,b][a,b], there exists c∈[a,b]c \in [a,b] with ∫abf(x) dx=f(c)(b−a)\int_a^b f(x)\,dx = f(c)(b-a). Its proof leans on the Extreme Value Theorem and Intermediate Value Theorem applied to ff, but it fundamentally needs ∫abf(x) dx\int_a^b f(x)\,dx to already be a well-defined real number — again supplied by the integrability theorem proven here. The two results form a natural pair: existence of the integral, then a precise description of its average value.

Before this

Where this leads