5.1: Approximating Areas
- Page ID
- 2511
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)- Use sigma (summation) notation to calculate sums and powers of integers.
- Use the sum of rectangular areas to approximate the area under a curve.
- Use Riemann sums to approximate area.
Archimedes was fascinated with calculating the areas of various shapes—in other words, the amount of space enclosed by the shape. He used a process that has come to be known as the method of exhaustion, which used smaller and smaller shapes, the areas of which could be calculated exactly, to fill an irregular region and thereby obtain closer and closer approximations to the total area. In this process, an area bounded by curves is filled with rectangles, triangles, and shapes with exact area formulas. These areas are then summed to approximate the area of the curved region.
In this section, we develop techniques to approximate the area between a curve, defined by a function \(f(x),\) and the x-axis on a closed interval \([a,b].\) Like Archimedes, we first approximate the area under the curve using shapes of known area (namely, rectangles). By using smaller and smaller rectangles, we get closer and closer approximations to the area. Taking a limit allows us to calculate the exact area under the curve.
Let’s start by introducing some notation to make the calculations easier. We then consider the case when \(f(x)\) is continuous and nonnegative. Later in the chapter, we relax some of these restrictions and develop techniques that apply in more general cases.
Sigma (Summation) Notation
As mentioned, we will use shapes of known area to approximate the area of an irregular region bounded by curves. This process often requires adding up long strings of numbers. To make it easier to write down these lengthy sums, we look at some new notation here, called sigma notation (also known as summation notation). The Greek capital letter \(Σ\), sigma, is used to express long sums of values in a compact form. For example, if we want to add all the integers from 1 to 20 without sigma notation, we have to write
\[1+2+3+4+5+6+7+8+9+10+11+12+13+14+15+16+17+18+19+20. \nonumber \]
We could probably skip writing a couple of terms and write
\[1+2+3+4+⋯+19+20, \nonumber \]
which is better, but still cumbersome. With sigma notation, we write this sum as
\[\sum_{i=1}^{20}i \nonumber \]
which is much more compact. Typically, sigma notation is presented in the form
\[\sum_{i=1}^{n}a_i \nonumber \]
where \(a_i\) describes the terms to be added, and the \(i\) is called the \(index\). Each term is evaluated, then we sum all the values, beginning with the value when \(i=1\) and ending with the value when \(i=n.\) For example, an expression like \(\displaystyle \sum_{i=2}^{7}s_i\) is interpreted as \(s_2+s_3+s_4+s_5+s_6+s_7\). Note that the index is used only to keep track of the terms to be added; it does not factor into the calculation of the sum itself. The index is therefore called a dummy variable. We can use any letter we like for the index. Typically, mathematicians use \(i, \,j, \,k, \,m\), and \(n\) for indices.
Let’s try a couple of examples of using sigma notation.
- Write in sigma notation and evaluate the sum of terms \(3^i\) for \(i=1,2,3,4,5.\)
- Write the sum in sigma notation:
\[1+\dfrac{1}{4}+\dfrac{1}{9}+\dfrac{1}{16}+\dfrac{1}{25}. \nonumber \]
Solution
- Write \[\sum_{i=1}^{5}3^i=3+3^2+3^3+3^4+3^5=363. \nonumber \]
- The denominator of each term is a perfect square. Using sigma notation, this sum can be written as \(\displaystyle \sum_{i=1}^5\dfrac{1}{i^2}\).
Write in sigma notation and evaluate the sum of terms \(2^i\) for \(i=3,4,5,6.\)
- Hint
-
Use the solving steps in Example \(\PageIndex{1}\) as a guide.
- Answer
-
\(\displaystyle \sum_{i=3}^{6}2^i=2^3+2^4+2^5+2^6=120\)
The properties associated with the summation process are given in the following rule.
Let \(a_1,a_2,…,a_n\) and \(b_1,b_2,…,b_n\) represent two sequences of terms and let \(c\) be a constant. The following properties hold for all positive integers \(n\) and for integers \(m\), with \(1≤m≤n.\)
- \(\displaystyle \sum_{i=1}^n c=nc\)
- \(\displaystyle \sum_{i=1}^n ca_i=c\sum_{i=1}^na_i\)
- \(\displaystyle \sum_{i=1}^n(a_i+b_i)=\sum_{i=1}^na_i+\sum_{i=1}^nb_i\)
- \(\displaystyle \sum_{i=1}^n(a_i−b_i)=\sum_{i=1}^na_i−\sum_{i=1}^nb_i\)
- \(\displaystyle \sum_{i=1}^na_i=\sum_{i=1}^ma_i+\sum_{i=m+1}^na_i\)
We prove properties (ii.) and (iii.) here, and leave proof of the other properties to the Exercises.
(ii.) We have
\[\sum_{i=1}^nca_i=ca_1+ca_2+ca_3+⋯+ca_n=c(a_1+a_2+a_3+⋯+a_n)=c\sum_{i=1}^na_i. \nonumber \]
(iii.) We have
\[ \begin{align} \sum_{i=1}^{n}(a_i+b_i) &=(a_1+b_1)+(a_2+b_2)+(a_3+b_3)+⋯+(a_n+b_n) \\[4pt] &=(a_1+a_2+a_3+⋯+a_n)+(b_1+b_2+b_3+⋯+b_n) \\[4pt] &=\sum_{i=1}^na_i+\sum_{i=1}^nb_i. \end {align} \nonumber \]
□
A few more formulas for frequently found functions simplify the summation process further. These are shown in the next rule, for sums and powers of integers, and we use them in the next set of examples.
1. The sum of \(n\) integers is given by
\[\sum_{i=1}^n i=1+2+⋯+n=\dfrac{n(n+1)}{2}. \label{sum1} \]
2. The sum of consecutive integers squared is given by
\[\sum_{i=1}^n i^2=1^2+2^2+⋯+n^2=\dfrac{n(n+1)(2n+1)}{6}. \label{sum2} \]
3. The sum of consecutive integers cubed is given by
\[\sum_{i=1}^n i^3=1^3+2^3+⋯+n^3=\dfrac{n^2(n+1)^2}{4}. \label{sum3} \]
Write using sigma notation and evaluate:
- The sum of the terms \((i−3)^2\) for \(i=1,2,…,200.\)
- The sum of the terms \((i^3−i^2)\) for \(i=1,2,3,4,5,6\)
Solution
a. Multiplying out \((i−3)^2\), we can break the expression into three terms.
\[\begin{align*} \sum_{i=1}^{200}(i−3)^2 &=\sum_{i=1}^{200}(i^2−6i+9) \\[4pt]
&=\sum_{i=1}^{200}i^2−\sum_{i=1}^{200}6i+\sum_{i=1}^{200}9 \\[4pt]
&=\sum_{i=1}^{200}i^2−6\sum_{i=1}^{200}i+\sum_{i=1}^{200}9 \\[4pt]
&=\dfrac{200(200+1)(400+1)}{6}−6 \left[\dfrac{200(200+1)}{2}\right]+9(200) \\[4pt]
&=2,686,700−120,600+1800 \\[4pt]
&=2,567,900 \end{align*}\]
b. Use sigma notation property iv. and the rules for the sum of squared terms and the sum of cubed terms.
\[\begin{align*} \sum_{i=1}^{6}(i^3−i^2) &=\sum_{i=1}^6 i^3−\sum_{i=1}^6 i^2 \\[4pt]
&=\dfrac{6^2(6+1)^2}{4}−\dfrac{6(6+1)(2(6)+1)}{6} \\[4pt]
&=\dfrac{1764}{4}−\dfrac{546}{6} \\[4pt]
&=350 \end{align*} \nonumber \]
Find the sum of the values of \(4+3i\) for \(i=1,2,…,100.\)
- Hint
-
Use the properties of sigma notation to solve the problem.
- Answer
-
\(15,550\)
Find the sum of the values of \(f(x)=x^3\) over the integers \(1,2,3,…,10.\)
Solution
Using Equation \ref{sum3}, we have
\[\sum_{i=0}^{10}i^3=\dfrac{(10)^2(10+1)^2}{4}=\dfrac{100(121)}{4}=3025 \nonumber \]
Evaluate the sum indicated by the notation \(\displaystyle \sum_{k=1}^{20}(2k+1)\).
- Hint
-
Use the rule on sum and powers of integers (Equations \ref{sum1}-\ref{sum3}).
- Answer
-
\(440\)
Approximating Area
Now that we have the necessary notation, we return to the problem at hand: approximating the area under a curve. Let \(f(x)\) be a continuous, nonnegative function defined on the closed interval \([a,b]\). We want to approximate the area \(A\) bounded by \(f(x)\) above, the \(x\)-axis below, the line \(x=a\) on the left, and the line \(x=b\) on the right (Figure \(\PageIndex{1}\)).
How do we approximate the area under this curve? The approach is a geometric one. By dividing a region into many small shapes that have known area formulas, we can sum these areas and obtain a reasonable estimate of the true area. We begin by dividing the interval \([a,b]\) into \(n\) subintervals of equal width, \(\dfrac{b−a}{n}\). We do this by selecting equally spaced points \(x_0,x_1,x_2,…,x_n\) with \(x_0=a,x_n=b,\) and
\[x_i−x_{i−1}=\dfrac{b−a}{n} \nonumber \]
for \(i=1,2,3,…,n.\)
We denote the width of each subinterval with the notation \(Δx,\) so \(Δx=\frac{b−a}{n}\) and
\[x_i=x_0+iΔx \nonumber \]
for \(i=1,2,3,…,n.\) This notion of dividing an interval \([a,b]\) into subintervals by selecting points from within the interval is used quite often in approximating the area under a curve, so let’s define some relevant terminology.
A set of points \(P={x_i}\) for \(i=0,1,2,…,n\) with \(a=x_0 < x_1 < x_2 < ... < x_n = b\), which divides the interval \([a, b]\) into subintervals of the form \([x_0, x_1], [x_1, x_2], ..., [x_{n-1}, x_n]\) is called a
We can use this regular partition as the basis of a method for estimating the area under the curve. We next examine two methods: the left-endpoint approximation and the right-endpoint approximation.
On each subinterval \([x_{i−1},x_i]\) (for \(i=1,2,3,…,n\)), construct a rectangle with width \(Δx\) and height equal to \(f(x_{i−1})\), which is the function value at the left endpoint of the subinterval. Then the area of this rectangle is \(f(x_{i−1})Δx\). Adding the areas of all these rectangles, we get an approximate value for \(A\) (Figure \(\PageIndex{2}\)). We use the notation \(L_n\) to denote that this is a left-endpoint approximation of \(A\) using \(n\) subintervals.
\[A≈L_n=f(x_0)Δx+f(x_1)Δx+⋯+f(x_{n−1})Δx=\sum_{i=1}^nf(x_{i−1})Δx \nonumber \]
The second method for approximating area under a curve is the right-endpoint approximation. It is almost the same as the left-endpoint approximation, but now the heights of the rectangles are determined by the function values at the right of each subinterval.
Construct a rectangle on each subinterval \([x_{i−1},x_i]\), only this time the height of the rectangle is determined by the function value \(f(x_i)\) at the right endpoint of the subinterval. Then, the area of each rectangle is \(f(x_i)\,Δx\) and the approximation for \(A\) is given by
\[A≈R_n=f(x_1)Δx+f(x_2)Δx+⋯+f(x_n)Δx=\sum_{i=1}^nf(x_i)Δx. \nonumber \]
The notation \(R_n\) indicates this is a right-endpoint approximation for \(A\) (Figure \(\PageIndex{3}\)).
The graphs in Figure \(\PageIndex{4}\) represent the curve \(f(x)=\dfrac{x^2}{2}\). In Figure \(\PageIndex{4b}\) we divide the region represented by the interval \([0,3]\) into six subintervals, each of width \(0.5\). Thus, \(Δx=0.5\). We then form six rectangles by drawing vertical lines perpendicular to \(x_{i−1}\), the left endpoint of each subinterval. We determine the height of each rectangle by calculating \(f(x_{i−1})\) for \(i=1,2,3,4,5,6.\) The intervals are \([0,0.5],[0.5,1],[1,1.5],[1.5,2],[2,2.5],[2.5,3]\). We find the area of each rectangle by multiplying the height by the width. Then, the sum of the rectangular areas approximates the area between \(f(x)\) and the \(x\)-axis. When the left endpoints are used to calculate height, we have a left-endpoint approximation. Thus,
\[\begin{align*} A≈L_6 &=\sum_{i=1}^6f(x_{i−1})Δx =f(x_0)Δx+f(x_1)Δx+f(x_2)Δx+f(x_3)Δx+f(x_4)Δx+f(x_5)Δx \\[4pt]
&=f(0)0.5+f(0.5)0.5+f(1)0.5+f(1.5)0.5+f(2)0.5+f(2.5)0.5 \\[4pt]
&=(0)0.5+(0.125)0.5+(0.5)0.5+(1.125)0.5+(2)0.5+(3.125)0.5 \\[4pt]
&=0+0.0625+0.25+0.5625+1+1.5625 \\[4pt]
&=3.4375 \,\text{units}^2\end{align*} \nonumber \]
In Figure \(\PageIndex{4b}\), we draw vertical lines perpendicular to \(x_i\) such that \(x_i\) is the right endpoint of each subinterval, and calculate \(f(x_i)\) for \(i=1,2,3,4,5,6\). We multiply each \(f(x_i)\) by \(Δx\) to find the rectangular areas, and then add them. This is a right-endpoint approximation of the area under \(f(x)\). Thus,
\[ \begin{align*} A≈R_6 &=\sum_{i=1}^6f(x_i)Δx=f(x_1)Δx+f(x_2)Δx+f(x_3)Δx+f(x_4)Δx+f(x_5)Δx+f(x_6)Δx\\[4pt]
&=f(0.5)0.5+f(1)0.5+f(1.5)0.5+f(2)0.5+f(2.5)0.5+f(3)0.5 \\[4pt]
&=(0.125)0.5+(0.5)0.5+(1.125)0.5+(2)0.5+(3.125)0.5+(4.5)0.5 \\[4pt]
&=0.0625+0.25+0.5625+1+1.5625+2.25 \\[4pt]
&=5.6875 \,\text{units}^2.\end{align*} \nonumber \]
Use both left-endpoint and right-endpoint approximations to approximate the area under the curve of \(f(x)=x^2\) on the interval \([0,2]\); use \(n=4\).
Solution
First, divide the interval \([0,2]\) into \(n\) equal subintervals. Using \(n=4,\, Δx=\dfrac{(2−0)}{4}=0.5\). This is the width of each rectangle. The intervals \([0,0.5],[0.5,1],[1,1.5],[1.5,2]\) are shown in Figure \(\PageIndex{5}\). Using a left-endpoint approximation, the heights are \(f(0)=0,\,f(0.5)=0.25,\,f(1)=1,\) and \(f(1.5)=2.25.\) Then,
\[ \begin{align*} L_4 &=f(x_0)Δx+f(x_1)Δx+f(x_2)Δx+f(x_3)Δx \\[4pt] &=0(0.5)+0.25(0.5)+1(0.5)+2.25(0.5) \\[4pt] &=1.75 \,\text{units}^2 \end{align*} \nonumber \]
The right-endpoint approximation is shown in Figure \(\PageIndex{6}\). The intervals are the same, \(Δx=0.5,\) but now use the right endpoint to calculate the height of the rectangles. We have
\[ \begin{align*} R_4 &=f(x_1)Δx+f(x_2)Δx+f(x_3)Δx+f(x_4)Δx \\[4pt] &=0.25(0.5)+1(0.5)+2.25(0.5)+4(0.5) \\[4pt] &=3.75 \,\text{units}^2 \end{align*} \nonumber \]
The left-endpoint approximation is \(1.75\,\text{units}^2\); the right-endpoint approximation is \(3.75 \,\text{units}^2\).
Sketch left-endpoint and right-endpoint approximations for \(f(x)=\dfrac{1}{x}\) on \([1,2]\); use \(n=4\). Approximate the area using both methods.
- Hint
-
Follow the solving strategy in Example \(\PageIndex{4}\) step-by-step.
- Answer
-
The left-endpoint approximation is \(0.7595 \,\text{units}^2\). The right-endpoint approximation is \(0.6345 \,\text{units}^2\). See the below Media.
Looking at Figure \(\PageIndex{4}\) and the graphs in Example \(\PageIndex{4}\), we can see that when we use a small number of intervals, neither the left-endpoint approximation nor the right-endpoint approximation is a particularly accurate estimate of the area under the curve. However, it seems logical that if we increase the number of points in our partition, our estimate of \(A\) will improve. We will have more rectangles, but each rectangle will be thinner, so we will be able to fit the rectangles to the curve more precisely.
We can demonstrate the improved approximation obtained through smaller intervals with an example. Let’s explore the idea of increasing \(n\), first in a left-endpoint approximation with four rectangles, then eight rectangles, and finally \(32\) rectangles. Then, let’s do the same thing in a right-endpoint approximation, using the same sets of intervals, of the same curved region. Figure \(\PageIndex{7}\) shows the area of the region under the curve \(f(x)=(x−1)^3+4\) on the interval \([0,2]\) using a left-endpoint approximation where \(n=4.\) The width of each rectangle is
\[Δx=\dfrac{2−0}{4}=\dfrac{1}{2}.\nonumber \]
The area is approximated by the summed areas of the rectangles, or
\[L_4=f(0)(0.5)+f(0.5)(0.5)+f(1)(0.5)+f(1.5)0.5=7.5 \,\text{units}^2\nonumber \]
Figure \(\PageIndex{8}\) shows the same curve divided into eight subintervals. Comparing the graph with four rectangles in Figure \(\PageIndex{7}\) with this graph with eight rectangles, we can see there appears to be less white space under the curve when \(n=8.\) This white space is area under the curve we are unable to include using our approximation. The area of the rectangles is
\[L_8=f(0)(0.25)+f(0.25)(0.25)+f(0.5)(0.25)+f(0.75)(0.25)+f(1)(0.25)+f(1.25)(0.25)+f(1.5)(0.25)+f(1.75)(0.25)=7.75 \,\text{units}^2\nonumber \]
The graph in Figure \(\PageIndex{9}\) shows the same function with \(32\) rectangles inscribed under the curve. There appears to be little white space left. The area occupied by the rectangles is
\[L_{32}=f(0)(0.0625)+f(0.0625)(0.0625)+f(0.125)(0.0625)+⋯+f(1.9375)(0.0625)=7.9375 \,\text{units}^2.\nonumber \]
We can carry out a similar process for the right-endpoint approximation method. A right-endpoint approximation of the same curve, using four rectangles (Figure \(\PageIndex{10}\)), yields an area
\[R_4=f(0.5)(0.5)+f(1)(0.5)+f(1.5)(0.5)+f(2)(0.5)=8.5 \,\text{units}^2.\nonumber \]
Dividing the region over the interval \([0,2]\) into eight rectangles results in \(Δx=\dfrac{2−0}{8}=0.25.\) The graph is shown in Figure \(\PageIndex{11}\). The area is
\[R_8=f(0.25)(0.25)+f(0.5)(0.25)+f(0.75)(0.25)+f(1)(0.25)+f(1.25)(0.25)+f(1.5)(0.25)+f(1.75)(0.25)+f(2)(0.25)=8.25 \,\text{units}^2\nonumber \]
Last, the right-endpoint approximation with \(n=32\) is close to the actual area (Figure \(\PageIndex{12}\)). The area is approximately
\[R_{32}=f(0.0625)(0.0625)+f(0.125)(0.0625)+f(0.1875)(0.0625)+⋯+f(2)(0.0625)=8.0625 \,\text{units}^2\nonumber \]
Based on these figures and calculations, it appears we are on the right track; the rectangles appear to approximate the area under the curve better as \(n\) gets larger. Furthermore, as \(n\) increases, both the left-endpoint and right-endpoint approximations appear to approach an area of \(8\) square units. Table \(\PageIndex{15}\) shows a numerical comparison of the left- and right-endpoint methods. The idea that the approximations of the area under the curve get better and better as \(n\) gets larger and larger is very important, and we now explore this idea in more detail.
Value of \(n\) | Approximate Area \(L_n\) | Approximate Area \(R_n\) |
---|---|---|
\(n=4\) | \(7.5\) | \(8.5\) |
\(n=8\) | \(7.75\) | \(8.25\) |
\(n=32\) | \(7.94\) | \(8.06\) |
Forming Riemann Sums
So far we have been using rectangles to approximate the area under a curve. The heights of these rectangles have been determined by evaluating the function at either the right or left endpoints of the subinterval \([x_{i−1},x_i]\). In reality, there is no reason to restrict evaluation of the function to one of these two points only. We could evaluate the function at any point \(x^∗_i\) in the subinterval \([x_{i−1},x_i]\), and use \(f(x^∗_i)\) as the height of our rectangle. This gives us an estimate for the area of the form
\[A≈\sum_{i=1}^nf(x^∗_i)\,Δx. \nonumber \]
A sum of this form is called a Riemann sum, named for the 19th-century mathematician Bernhard Riemann, who developed the idea.
Let \(f(x)\) be defined on a closed interval \([a,b]\) and let \(P\) be any partition of \([a,b]\). Let \(Δx_i\) be the width of each subinterval \([x_{i−1},x_i]\) and for each \(i\), let \(x^∗_i\) be any point in \([x_{i−1},\,x_i]\). A Riemann sum is defined for \(f(x)\) as
\[\sum_{i=1}^nf(x^∗_i)\,Δx_i. \nonumber \]
At this point, we'll choose a regular partition \(P\), as we have in our examples above. This forces all \(Δx_i\) to be equal to \(Δx = \dfrac{b-a}{n}\) for any natural number of intervals \(n\).
Recall that with the left- and right-endpoint approximations, the estimates seem to get better and better as \(n\) get larger and larger. The same thing happens with Riemann sums. Riemann sums give better approximations for larger values of \(n\). We are now ready to define the area under a curve in terms of Riemann sums.
Let \(f(x)\) be a continuous, nonnegative function on an interval \([a,b]\), and let \(\displaystyle \sum_{i=1}^nf(x^∗_i)\,Δx\) be a Riemann sum for \(f(x)\) with a regular partition \(P\). Then, the area under the curve \(y=f(x)\) on \([a,b]\) is given by
\[A=\lim_{n→∞}\sum_{i=1}^nf(x^∗_i)\,Δx. \nonumber \]
See a graphical demonstration of the construction of a Riemann sum.
Some subtleties here are worth discussing. First, note that taking the limit of a sum is a little different from taking the limit of a function \(f(x)\) as \(x\) goes to infinity. Limits of sums are discussed in detail in the chapter on Sequences and Series; however, for now we can assume that the computational techniques we used to compute limits of functions can also be used to calculate limits of sums.
Second, we must consider what to do if the expression converges to different limits for different choices of \({x^∗_i}.\) Fortunately, this does not happen. Although the proof is beyond the scope of this text, it can be shown that if \(f(x)\) is continuous on the closed interval \([a,b]\), then \(\displaystyle \lim_{n→∞}\sum_{i=1}^nf(x^∗_i)Δx\) exists and is unique (in other words, it does not depend on the choice of \({x^∗_i}\)).
We look at some examples shortly. But, before we do, let’s take a moment and talk about some specific choices for \({x^∗_i}\). Although any choice for \({x^∗_i}\) gives us an estimate of the area under the curve, we don’t necessarily know whether that estimate is too high (overestimate) or too low (underestimate). If it is important to know whether our estimate is high or low, we can select our value for \({x^∗_i}\) to guarantee one result or the other.
If we want an overestimate, for example, we can choose \({x^∗_i}\) such that for \(i=1,2,3,…,n,\) \(f(x^∗_i)≥f(x)\) for all \(x∈[x_i−1,x_i]\). In other words, we choose \({x^∗_i}\) so that for \(i=1,2,3,…,n,\) \(f(x^∗_i)\) is the maximum function value on the interval \([x_{i−1},x_i]\). If we select \({x^∗_i}\) in this way, then the Riemann sum \(\displaystyle \sum_{i=1}^nf(x^∗_i)Δx\) is called an upper sum. Similarly, if we want an underestimate, we can choose \({x^∗_i}\) so that for \(i=1,2,3,…,n,\) \(f(x^∗_i)\) is the minimum function value on the interval \([x_{i−1},x_i]\). In this case, the associated Riemann sum is called a lower sum. Note that if \(f(x)\) is either increasing or decreasing throughout the interval \([a,b]\), then the maximum and minimum values of the function occur at the endpoints of the subintervals, so the upper and lower sums are just the same as the left- and right-endpoint approximations.
Find a lower sum for \(f(x)=10−x^2\) on \([1,2]\); let \(n=4\) subintervals.
Solution
With \(n=4\) over the interval \([1,2], \,Δx=\dfrac{1}{4}\). We can list the intervals as \([1,1.25],\,[1.25,1.5],\,[1.5,1.75],\) and \([1.75,2]\). Because the function is decreasing over the interval \([1,2],\) Figure shows that a lower sum is obtained by using the right endpoints.
The Riemann sum is
\[\begin{align*} \sum_{k=1}^4(10−x^2)(0.25) &=0.25[10−(1.25)^2+10−(1.5)^2+10−(1.75)^2+10−(2)^2] \\[4pt]
&=0.25[8.4375+7.75+6.9375+6] \\[4pt]
&=7.28 \,\text{units}^2.\end{align*}\]
The area of \(7.28\) \(\text{units}^2\) is a lower sum and an underestimate.
- Find an upper sum for \(f(x)=10−x^2\) on \([1,2]\); let \(n=4.\)
- Sketch the approximation.
- Hint
-
\(f(x)\) is decreasing on \([1,2]\), so the maximum function values occur at the left endpoints of the subintervals.
- Answer
-
a. Upper sum=\(8.0313 \,\text{units}^2.\)
b.
Find a lower sum for \(f(x)=\sin x\) over the interval \([a,b]=\left[0,\frac{π}{2} \right]\); let \(n=6.\)
Solution
Let’s first look at the graph in Figure \(\PageIndex{14}\) to get a better idea of the area of interest.
The intervals are \(\left[0,\frac{π}{12}\right],\,\left[\frac{π}{12},\frac{π}{6}\right],\,\left[\frac{π}{6},\frac{π}{4}\right],\,\left[\frac{π}{4},\frac{π}{3}\right],\,\left[\frac{π}{3},\frac{5π}{12}\right]\), and \(\left[\frac{5π}{12},\frac{π}{2}\right]\). Note that \(f(x)=\sin x\) is increasing on the interval \(\left[0,\frac{π}{2}\right]\), so a left-endpoint approximation gives us the lower sum. A left-endpoint approximation is the Riemann sum \(\sum_{i=0}^5\sin x_i\left(\tfrac{π}{12}\right)\).We have
\[A≈\sin(0)\left(\tfrac{π}{12}\right)+\sin\left(\tfrac{π}{12}\right)\left(\tfrac{π}{12}\right)+\sin\left(\tfrac{π}{6}\right)\left(\tfrac{π}{12}\right)+\sin\left(\tfrac{π}{4}\right)\left(\tfrac{π}{12}\right)+\sin\left(\tfrac{π}{3}\right)\left(\tfrac{π}{12}\right)+\sin\left(\tfrac{5π}{12}\right)\left(\tfrac{π}{12}\right)\approx 0.863 \,\text{units}^2. \nonumber \]
Using the function \(f(x)=\sin x\) over the interval \(\left[0,\frac{π}{2}\right],\) find an upper sum; let \(n=6.\)
- Hint
-
Follow the steps from Example \(\PageIndex{6}\).
- Answer
-
\(A≈1.125 \,\text{units}^2\)
Key Concepts
- The use of sigma (summation) notation of the form \(\displaystyle \sum_{i=1}^na_i\) is useful for expressing long sums of values in compact form.
- For a continuous function defined over an interval \([a,b],\) the process of dividing the interval into \(n\) equal parts, extending a rectangle to the graph of the function, calculating the areas of the series of rectangles, and then summing the areas yields an approximation of the area of that region.
- When using a regular partition, the width of each rectangle is \(Δx=\dfrac{b−a}{n}\).
- Riemann sums are expressions of the form \(\displaystyle \sum_{i=1}^nf(x^∗_i)Δx,\) and can be used to estimate the area under the curve \(y=f(x).\) Left- and right-endpoint approximations are special kinds of Riemann sums where the values of \({x^∗_i}\) are chosen to be the left or right endpoints of the subintervals, respectively.
- Riemann sums allow for much flexibility in choosing the set of points \({x^∗_i}\) at which the function is evaluated, often with an eye to obtaining a lower sum or an upper sum.
Key Equations
- Properties of Sigma Notation
\[\begin{align*} \sum_{i=1}^nc&=nc \\[4pt]
\sum_{i=1}^nca_i &=c\sum_{i=1}^na_i \\[4pt]
\sum_{i=1}^n(a_i+b_i) &=\sum_{i=1}^na_i+\sum_{i=1}^nb_i \\[4pt]
\sum_{i=1}^n(a_i−b_i) &=\sum_{i=1}^na_i−\sum_{i=1}^nb_i \\[4pt]
\sum_{i=1}^na_i&=\sum_{i=1}^ma_i+\sum_{i=m+1}^na_i \end{align*}\]
- Sums and Powers of Integers
\[\sum_{i=1}^ni=1+2+⋯+n=\dfrac{n(n+1)}{2} \nonumber \]
\[\sum_{i=1}^ni^2=1^2+2^2+⋯+n^2=\dfrac{n(n+1)(2n+1)}{6} \nonumber \]
\[\sum_{i=0}^ni^3=1^3+2^3+⋯+n^3=\dfrac{n^2(n+1)^2}{4} \nonumber \]
- Left-Endpoint Approximation
\(A≈L_n=f(x_0)Δx+f(x_1)Δx+⋯+f(x_{n−1})Δx=\displaystyle \sum_{i=1}^nf(x_{i−1})Δx\)
- Right-Endpoint Approximation
\(A≈R_n=f(x_1)Δx+f(x_2)Δx+⋯+f(x_n)Δx=\displaystyle \sum_{i=1}^nf(x_i)Δx\)
Glossary
- left-endpoint approximation
- an approximation of the area under a curve computed by using the left endpoint of each subinterval to calculate the height of the vertical sides of each rectangle
- lower sum
- a sum obtained by using the minimum value of \(f(x)\) on each subinterval
- partition
- a set of points that divides an interval into subintervals
- regular partition
- a partition in which the subintervals all have the same width
- riemann sum
- an estimate of the area under the curve of the form \(A≈\displaystyle \sum_{i=1}^nf(x^∗_i)Δx\)
- right-endpoint approximation
- the right-endpoint approximation is an approximation of the area of the rectangles under a curve using the right endpoint of each subinterval to construct the vertical sides of each rectangle
- sigma notation
- (also, summation notation) the Greek letter sigma (\(Σ\)) indicates addition of the values; the values of the index above and below the sigma indicate where to begin the summation and where to end it
- upper sum
- a sum obtained by using the maximum value of \(f(x)\) on each subinterval