Approximating Areas Under Curves and Riemann Sum Foundations

Foundations of Calculus and Approximating Area

  • Calculus is built upon three foundational geometric and physical questions:
    • The behavior of points infinitely close to each other or infinitely far away, which establishes the concept of limits.
    • The definition of a tangent line to a curve or the instantaneous rate of change, which establishes the concept of derivatives.
    • The definition of the area of curved shapes and the accumulation of changing quantities over an interval, which establishes the concept of integration (explored in Chapter 5).
  • Development strategy for integration:
    • Rather than computing areas under curved curves directly from the outset, approximations are constructed using simple geometric polygons (specifically rectangles).
    • The approximation is subsequently transformed into an exact equality by evaluating the limit as the number of approximating rectangles approaches infinity.
  • Section 5.1 focuses entirely on these rectangular approximations.
  • Rectangles are selected because their geometric area is straightforward to compute:
    • Area=base×height\text{Area} = \text{base} \times \text{height}
  • Area under the curve:
    • Defined as the area of the region bounded between the graph of a function and the xx-axis over a specified interval.
  • Mechanics of basic rectangular approximation:
    • Rectangular approximations are not exact because rectangles inevitably overshoot (extending above the curve and including unwanted area) or undershoot (failing to cover the entirety of the area beneath the curve).
    • Over a partition of the domain, the total approximate area is the sum of the individual rectangular areas:
    • Total Area=(base1×height1)+(base2×height2)+⋯+(basen×heightn)\text{Total Area} = (\text{base}_1 \times \text{height}_1) + (\text{base}_2 \times \text{height}_2) + \dots + (\text{base}_n \times \text{height}_n)
    • The base corresponds to the width of the interval along the xx-axis (Δx\Delta x), and the height corresponds to the yy-coordinate of the function.

Example 1: Arbitrary Rectangular Approximation

  • Approximation of the area under a curve using five distinct rectangles with varying base widths and heights:
    • Rectangle 1:
    • Base width: 22
    • Height: 66
    • Area: 2×6=122 \times 6 = 12
    • Rectangle 2:
    • Spans from x=2x = 2 to x=4x = 4, giving a base width of 4−2=24 - 2 = 2
    • Height: 55
    • Area: 2×5=102 \times 5 = 10
    • Rectangle 3:
    • Spans from x=4x = 4 to x=8x = 8, giving a base width of 8−4=48 - 4 = 4
    • Height: 33
    • Area: 4×3=124 \times 3 = 12
    • Rectangle 4:
    • Spans a base width of 44
    • Height: 55
    • Area: 4×5=204 \times 5 = 20
    • Rectangle 5:
    • Spans a base width of 44
    • Height: 77
    • Area: 4×7=284 \times 7 = 28
  • Total area computation:
    • Total Area=12+10+12+20+28=82 square units\text{Total Area} = 12 + 10 + 12 + 20 + 28 = 82\,\text{square units}
  • While base widths do not strictly need to be identical, setting all bases to an equal width significantly streamlines systematic algebraic computations.

Right Endpoint Approximation (RnR_n)

  • Subdividing an interval [a,b][a, b] into nn equal subintervals:
    • The uniform base width (denoted Δx\Delta x) is computed by taking the total length of the interval and dividing by the desired number of rectangles nn:
    • Δx=b−an\Delta x = \frac{b - a}{n}
    • The interval [a,b][a, b] is partitioned into nn subintervals of equal width: [xi−1,xi][x_{i-1}, x_i] for i=1,2,…,ni = 1, 2, \dots, n.
    • Endpoints are computed iteratively by shifting to the right from aa by multiples of Δx\Delta x:
    • x0=ax_0 = a
    • x1=a+Δxx_1 = a + \Delta x
    • x2=a+2 Δxx_2 = a + 2\,\Delta x
    • xi=a+i Δxx_i = a + i\,\Delta x
    • xn=a+n Δx=bx_n = a + n\,\Delta x = b
  • Height selection for the right endpoint approximation:
    • Rectangles are constructed such that their top-right corners intersect the curve.
    • For each subinterval [xi−1,xi][x_{i-1}, x_i], the height of the ii-th rectangle is determined by evaluating the function at the right endpoint xix_i:
    • Heighti=f(xi)\text{Height}_i = f(x_i)
  • Formula for the nn-th Right Endpoint Approximation (RnR_n):
    • Factoring out the uniform base width Δx\Delta x yields:
    • Rn=Δx⋅[f(x1)+f(x2)+⋯+f(xn)]R_n = \Delta x \cdot [f(x_1) + f(x_2) + \dots + f(x_n)]
  • Behavior as nn increases:
    • Making rectangles thinner by increasing nn diminishes the magnitude of both undershooting and overshooting gaps.
    • As n→∞n \to \infty, the approximation converges to the true area under the curve.

Example 2: Approximating the Area Under f(x)=x2f(x) = x^2 on [1,3][1, 3]

  • Problem objective: Estimate the area under f(x)=x2f(x) = x^2 on the interval [1,3][1, 3] using right endpoint approximations with n=4n = 4 and n=6n = 6.
  • Part (a): Four rectangles (n=4n = 4):
    • Calculate the uniform base width Δx\Delta x:
    • Δx=3−14=24=12\Delta x = \frac{3 - 1}{4} = \frac{2}{4} = \frac{1}{2}
    • Identify the subinterval endpoints by starting at a=1a = 1 and repeatedly adding Δx=12\Delta x = \frac{1}{2}:
    • x0=1x_0 = 1
    • x1=1+12=32x_1 = 1 + \frac{1}{2} = \frac{3}{2}
    • x2=32+12=42=2x_2 = \frac{3}{2} + \frac{1}{2} = \frac{4}{2} = 2
    • x3=42+12=52x_3 = \frac{4}{2} + \frac{1}{2} = \frac{5}{2}
    • x4=52+12=62=3x_4 = \frac{5}{2} + \frac{1}{2} = \frac{6}{2} = 3
    • Calculate rectangular heights using f(x)=x2f(x) = x^2:
    • f(x1)=f(32)=(32)2=94f(x_1) = f\left(\frac{3}{2}\right) = \left(\frac{3}{2}\right)^2 = \frac{9}{4}
    • f(x2)=f(42)=(42)2=164f(x_2) = f\left(\frac{4}{2}\right) = \left(\frac{4}{2}\right)^2 = \frac{16}{4}
    • f(x3)=f(52)=(52)2=254f(x_3) = f\left(\frac{5}{2}\right) = \left(\frac{5}{2}\right)^2 = \frac{25}{4}
    • f(x4)=f(62)=(62)2=364f(x_4) = f\left(\frac{6}{2}\right) = \left(\frac{6}{2}\right)^2 = \frac{36}{4}
    • Compute R4R_4:
    • R4=Δx⋅[f(x1)+f(x2)+f(x3)+f(x4)]R_4 = \Delta x \cdot [f(x_1) + f(x_2) + f(x_3) + f(x_4)]
    • R4=12[94+164+254+364]R_4 = \frac{1}{2} \left[ \frac{9}{4} + \frac{16}{4} + \frac{25}{4} + \frac{36}{4} \right]
    • Sum of numerators: 9+16+25+36=869 + 16 + 25 + 36 = 86
    • R4=12(864)=868=10.75R_4 = \frac{1}{2} \left( \frac{86}{4} \right) = \frac{86}{8} = 10.75
  • Part (b): Six rectangles (n=6n = 6):
    • Calculate the uniform base width Δx\Delta x:
    • Δx=3−16=26=13\Delta x = \frac{3 - 1}{6} = \frac{2}{6} = \frac{1}{3}
    • Identify the subinterval endpoints by starting at a=1a = 1 and repeatedly adding Δx=13\Delta x = \frac{1}{3}:
    • x0=1x_0 = 1
    • x1=1+13=43x_1 = 1 + \frac{1}{3} = \frac{4}{3}
    • x2=43+13=53x_2 = \frac{4}{3} + \frac{1}{3} = \frac{5}{3}
    • x3=53+13=63=2x_3 = \frac{5}{3} + \frac{1}{3} = \frac{6}{3} = 2
    • x4=63+13=73x_4 = \frac{6}{3} + \frac{1}{3} = \frac{7}{3}
    • x5=73+13=83x_5 = \frac{7}{3} + \frac{1}{3} = \frac{8}{3}
    • x6=83+13=93=3x_6 = \frac{8}{3} + \frac{1}{3} = \frac{9}{3} = 3
    • Calculate rectangular heights using f(x)=x2f(x) = x^2:
    • f(x1)=(43)2=169f(x_1) = \left(\frac{4}{3}\right)^2 = \frac{16}{9}
    • f(x2)=(53)2=259f(x_2) = \left(\frac{5}{3}\right)^2 = \frac{25}{9}
    • f(x3)=(63)2=369f(x_3) = \left(\frac{6}{3}\right)^2 = \frac{36}{9}
    • f(x4)=(73)2=499f(x_4) = \left(\frac{7}{3}\right)^2 = \frac{49}{9}
    • f(x5)=(83)2=649f(x_5) = \left(\frac{8}{3}\right)^2 = \frac{64}{9}
    • f(x6)=(93)2=819f(x_6) = \left(\frac{9}{3}\right)^2 = \frac{81}{9}
    • Compute R6R_6:
    • R6=13[169+259+369+499+649+819]R_6 = \frac{1}{3} \left[ \frac{16}{9} + \frac{25}{9} + \frac{36}{9} + \frac{49}{9} + \frac{64}{9} + \frac{81}{9} \right]
    • Sum of numerators: 16+25+36+49+64+81=27116 + 25 + 36 + 49 + 64 + 81 = 271
    • R6=13(2719)=27127≈10.04R_6 = \frac{1}{3} \left( \frac{271}{9} \right) = \frac{271}{27} \approx 10.04
  • Analysis of results:
    • The approximation decreased from 10.7510.75 to 10.0410.04.
    • Because f(x)=x2f(x) = x^2 is an increasing function on [1,3][1, 3], right endpoint rectangles consistently overshoot the actual area under the graph.
    • Increasing nn from 44 to 66 reduces the overcounting margin, driving the estimate closer to the actual value.

Summation (Sigma) Notation

  • Definition and structure:
    • Summation notation provides a concise mathematical shorthand for expressing long sums of indexed numbers am,am+1,…,ana_m, a_{m+1}, \dots, a_n:
    • ∑i=mnai=am+am+1+am+2+⋯+an\sum_{i=m}^n a_i = a_m + a_{m+1} + a_{m+2} + \dots + a_n
    • Read aloud as: "the sum as ii goes from mm to nn of aia_i".
  • Components of the notation:
    • Capital Greek letter Sigma (Σ\Sigma): Represents "sum" (originating from the English letter S).
    • Summation index (ii): A "dummy variable" that keeps track of the current term. Any symbol can serve as the index.
    • Lower limit (mm): The starting integer value substituted into the expression.
    • Upper limit (nn): The stopping integer value substituted into the expression.
  • Algorithmic analogy:
    • Sigma notation is equivalent to a programming for loop initialized at i=mi = m, running up to i=ni = n, incrementing by 11 each iteration, and accumulating the term aia_i into a total sum.
  • Illustrative evaluation example:
    • ∑i=35(i2+i)\sum_{i=3}^5 (i^2 + i)
    • Term for i=3i = 3: 32+3=9+3=123^2 + 3 = 9 + 3 = 12
    • Term for i=4i = 4: 42+4=16+4=204^2 + 4 = 16 + 4 = 20
    • Term for i=5i = 5: 52+5=25+5=305^2 + 5 = 25 + 5 = 30
    • Sum evaluation:
    • ∑i=35(i2+i)=(32+3)+(42+4)+(52+5)=12+20+30=62\sum_{i=3}^5 (i^2 + i) = (3^2 + 3) + (4^2 + 4) + (5^2 + 5) = 12 + 20 + 30 = 62
  • Compact representation of Right Endpoint Approximation (RnR_n):
    • Rn=Δx∑i=1nf(xi)R_n = \Delta x \sum_{i=1}^n f(x_i)

Left Endpoint and Midpoint Approximations

  • Left Endpoint Approximation (LnL_n):
    • Rectangles are constructed such that their top-left corners touch the graph of f(x)f(x).
    • On each subinterval [xi−1,xi][x_{i-1}, x_i], the height is evaluated at the left boundary xi−1x_{i-1}:
    • Heighti=f(xi−1)\text{Height}_i = f(x_{i-1})
    • Expressed using an index starting at 11:
    • Ln=Δx∑i=1nf(xi−1)L_n = \Delta x \sum_{i=1}^n f(x_{i-1})
    • Expressed alternatively by re-indexing to start at 00 and stop at n−1n - 1:
    • Ln=Δx∑i=0n−1f(xi)L_n = \Delta x \sum_{i=0}^{n-1} f(x_i)
    • The index begins at x0=ax_0 = a and stops at xn−1x_{n-1}, omitting the right boundary xn=bx_n = b.
  • Midpoint Approximation (MnM_n):
    • Rectangles are constructed such that their heights match the value of the function at the exact midpoint of each subinterval.
    • The midpoint of subinterval [xi−1,xi][x_{i-1}, x_i] is defined as the arithmetic average of the two endpoints:
    • xi∗=xi−1+xi2x_i^* = \frac{x_{i-1} + x_i}{2}
    • The height of the ii-th rectangle is:
    • Heighti=f(xi−1+xi2)\text{Height}_i = f\left(\frac{x_{i-1} + x_i}{2}\right)
    • Formula for MnM_n:
    • Mn=Δx∑i=1nf(xi−1+xi2)M_n = \Delta x \sum_{i=1}^n f\left(\frac{x_{i-1} + x_i}{2}\right)
  • Accuracy comparison across approximations:
    • For a fixed number of rectangles nn, the midpoint approximation MnM_n typically yields a substantially more accurate approximation than either LnL_n or RnR_n.
    • Because the midpoint sample lies in the interior of each subinterval, the rectangle overshoots the curve on one side and undershoots it on the other side. These two opposite errors partially cancel each other out across each subinterval.

Physical Application: Estimating Distance from Velocity Data

  • Physical foundations:
    • At a constant speed vv, distance traveled across elapsed time tt is given by:
    • Distance=speed×time=v⋅t\text{Distance} = \text{speed} \times \text{time} = v \cdot t
    • When velocity is non-constant, this product is no longer directly applicable over the entire duration.
    • Rectangular approximations operate on the premise that over sufficiently small time increments (subintervals), velocity is approximately constant.
    • Total distance is approximated by accumulating the estimated distances across all subintervals.
  • Example 3 setup and data:
    • An object's velocity v(t)v(t) was recorded at half-second intervals over the span t∈[0,3]t \in [0, 3]:
    • At t=0.0 st = 0.0\,\text{s}: v(0.0)=0 m/sv(0.0) = 0\,\text{m/s}
    • At t=0.5 st = 0.5\,\text{s}: v(0.5)=12 m/sv(0.5) = 12\,\text{m/s}
    • At t=1.0 st = 1.0\,\text{s}: v(1.0)=18 m/sv(1.0) = 18\,\text{m/s}
    • At t=1.5 st = 1.5\,\text{s}: v(1.5)=25 m/sv(1.5) = 25\,\text{m/s}
    • At t=2.0 st = 2.0\,\text{s}: v(2.0)=20 m/sv(2.0) = 20\,\text{m/s}
    • At t=2.5 st = 2.5\,\text{s}: v(2.5)=14 m/sv(2.5) = 14\,\text{m/s}
    • At t=3.0 st = 3.0\,\text{s}: v(3.0)=20 m/sv(3.0) = 20\,\text{m/s}
  • Part (a): Right endpoint approximation with six rectangles (R6R_6):
    • Calculate time increment Δt\Delta t for n=6n = 6:
    • Δt=3−06=0.5 s\Delta t = \frac{3 - 0}{6} = 0.5\,\text{s}
    • Right endpoint sample times: t1=0.5t_1 = 0.5, t2=1.0t_2 = 1.0, t3=1.5t_3 = 1.5, t4=2.0t_4 = 2.0, t5=2.5t_5 = 2.5, t6=3.0t_6 = 3.0
    • Computation:
    • R6=Δt∑i=16v(ti)R_6 = \Delta t \sum_{i=1}^6 v(t_i)
    • R6=0.5⋅[v(0.5)+v(1.0)+v(1.5)+v(2.0)+v(2.5)+v(3.0)]R_6 = 0.5 \cdot [v(0.5) + v(1.0) + v(1.5) + v(2.0) + v(2.5) + v(3.0)]
    • R6=0.5⋅[12+18+25+20+14+20]R_6 = 0.5 \cdot [12 + 18 + 25 + 20 + 14 + 20]
    • R6=0.5⋅[109]=54.5 metersR_6 = 0.5 \cdot [109] = 54.5\,\text{meters}
  • Part (b): Left endpoint approximation with six rectangles (L6L_6):
    • The subinterval width remains Δt=0.5 s\Delta t = 0.5\,\text{s}.
    • Left endpoint sample times: t0=0.0t_0 = 0.0, t1=0.5t_1 = 0.5, t2=1.0t_2 = 1.0, t3=1.5t_3 = 1.5, t4=2.0t_4 = 2.0, t5=2.5t_5 = 2.5
    • Computation:
    • L6=Δt∑i=05v(ti)L_6 = \Delta t \sum_{i=0}^5 v(t_i)
    • L6=0.5⋅[v(0.0)+v(0.5)+v(1.0)+v(1.5)+v(2.0)+v(2.5)]L_6 = 0.5 \cdot [v(0.0) + v(0.5) + v(1.0) + v(1.5) + v(2.0) + v(2.5)]
    • L6=0.5⋅[0+12+18+25+20+14]L_6 = 0.5 \cdot [0 + 12 + 18 + 25 + 20 + 14]
    • L6=0.5⋅[89]=44.5 metersL_6 = 0.5 \cdot [89] = 44.5\,\text{meters}
    • With L6=44.5 mL_6 = 44.5\,\text{m} and R6=54.5 mR_6 = 54.5\,\text{m}, one method undershoots and the other overshoots; the actual total distance is estimated to lie near the center of these bounds (around 50 m50\,\text{m}.
  • Part (c): Midpoint approximation with three rectangles (M3M_3):
    • Because n=3n = 3, Δt\Delta t must be recalculated:
    • Δt=3−03=1.0 s\Delta t = \frac{3 - 0}{3} = 1.0\,\text{s}
    • Subintervals: [0,1][0, 1], [1,2][1, 2], and [2,3][2, 3]
    • Subinterval midpoints:
    • For [0,1][0, 1]: 0+12=0.5 s\frac{0 + 1}{2} = 0.5\,\text{s}
    • For [1,2][1, 2]: 1+22=1.5 s\frac{1 + 2}{2} = 1.5\,\text{s}
    • For [2,3][2, 3]: 2+32=2.5 s\frac{2 + 3}{2} = 2.5\,\text{s}
    • Computation:
    • M3=1.0⋅[v(0.5)+v(1.5)+v(2.5)]M_3 = 1.0 \cdot [v(0.5) + v(1.5) + v(2.5)]
    • M3=1.0⋅[12+25+14]=51 metersM_3 = 1.0 \cdot [12 + 25 + 14] = 51\,\text{meters}
    • Result analysis:
    • The value 51 meters51\,\text{meters} lands directly between L6L_6 (44.5 m44.5\,\text{m}) and R6R_6 (54.5 m54.5\,\text{m}), demonstrating that the midpoint approximation achieves superior balance and accuracy even with half the number of rectangles (n=3n = 3 vs. n=6n = 6).

Convergence and the Definition of Area Under a Curve

  • Equivalence of approximations in the limit:
    • Turning discrete geometric approximations into exact values requires taking the limit as the number of subintervals approaches infinity (n→∞n \to \infty).
    • For well-behaved functions, the choice of evaluation point within each subinterval (right endpoint, left endpoint, or midpoint) becomes irrelevant in the infinite limit.
  • Formal Theorem for Continuous Functions:
    • Let ff be a continuous function defined on the closed interval [a,b][a, b].
    • There exists a unique real number AA such that:
    • lim⁡n→∞Rn=lim⁡n→∞Ln=lim⁡n→∞Mn=A\lim_{n \to \infty} R_n = \lim_{n \to \infty} L_n = \lim_{n \to \infty} M_n = A
    • If f(x)≥0f(x) \ge 0 for all x∈[a,b]x \in [a, b], then A≥0A \ge 0.
    • This common limit AA is formally defined as the exact area under the curve y=f(x)y = f(x) on [a,b][a, b].
  • Theoretical role of continuity in convergence:
    • The width of each subinterval is Δx=b−an\Delta x = \frac{b - a}{n}.
    • As n→∞n \to \infty, Δx→0\Delta x \to 0, meaning the endpoints xi−1x_{i-1}, xix_i, and the midpoint xi−1+xi2\frac{x_{i-1} + x_i}{2} converge toward the exact same value.
    • By the fundamental definition of continuity, if input coordinates get arbitrarily close to one another, their corresponding function outputs f(xi−1)f(x_{i-1}), f(xi)f(x_i), and f(xi−1+xi2)f\left(\frac{x_{i-1} + x_i}{2}\right) must also get arbitrarily close to one another.
    • Consequently, the evaluations across all three methods converge to the identical sum in the limit.

Questions & Discussion

  • Summing an infinite collection of rectangles:
    • Question: How do we add up an infinite amount of rectangles?
    • Answer: By evaluating limits. The discrete sum with nn rectangles is evaluated under the limit process as n→∞n \to \infty.
  • Upper bound in sigma notation:
    • Question: The nn on the top is basically the number we stop at?
    • Answer: Yes. The variable mm at the base of the sigma denotes the starting index, and the variable nn at the top denotes the stopping index. It is important to emphasize that these limits are not the numbers being summed; they are indices that specify which evaluated terms enter the summation.
  • Dummy variables and indices:
    • Question: Are the ii and the nn like normal variables?
    • Answer: The variable ii is a dummy variable/index that keeps track of the progression through the sum, and any letter can be substituted in its place. The variable nn represents the total number of rectangles into which the interval is partitioned. Because right endpoint approximations evaluate after shifting over by Δx\Delta x, they initialize at i=1i = 1.
  • Starting index for left endpoint approximations:
    • Question: And left endpoint thing was a thing when we start with zero?
    • Answer: Yes. The left endpoint approximation begins evaluating at the far left boundary x0=ax_0 = a, which corresponds to starting the summation index at i=0i = 0 and terminating at n−1n - 1.
  • Behavior when functions reach zero or negative values:
    • Question: What if our line goes to zero?
    • Answer: That condition is addressed in subsequent sections. For introductory approximations, the function is assumed to satisfy f(x)≥0f(x) \ge 0 to keep the geometric concept of area straightforward.
  • Derivation of evaluation points in midpoint approximation:
    • Question: How do we know to derive the equation for tt? Like, where am I getting these numbers from?… The line under that M3M_3 equals summation…
    • Answer: Those values are obtained because the midpoint approximation requires evaluating the function at the center of each subinterval, which is the average of successive endpoints. For subintervals [0,1][0, 1], [1,2][1, 2], and [2,3][2, 3], the midpoints are 0+12=0.5\frac{0 + 1}{2} = 0.5, 1+22=1.5\frac{1 + 2}{2} = 1.5, and 2+32=2.5\frac{2 + 3}{2} = 2.5. These times are then referenced directly in the table to retrieve their respective velocities.