Riemann sums

Let f(x) be defined for axb (Figure 1). Assume at pleasure n1 points between a and b a=x0<x1<x2<<xn=b, thus dividing the interval I=[a,b] into n subintervals: [a,x1],[x1,x2],,[xn1,b]. In each of these subintervals, take a value of x=xk where xk1xkxk, and form the sum (8.1.1)Sn=f(x1)(x1x0)+f(x2)(x2x1)++f(xn)(xnxn1) This sum is called a Riemann sum for f on the interval [a,b], after German mathematician George Friedrich Bernhard Riemann (1826–1866). Now let’s see what the geometrical interpretation of Sn is. First suppose f(x)0 for all x[a,b]. In this case

f(x1)(x1x0)= the area of the rectangle AP1R1M1,

f(x2)(x2x1)= the area of the rectangle M1P2R2M2,

f(xn1)(xn1xn2)= the area of the rectangle Mn2Pn1Rn1Mn1,

f(xn)(xnxn1)= the area of the rectangle Mn1PnRnB.

Figure 1.
Figure 1.

Thus Sn represents the total area of all rectangles in Figure 1, and approximates the area of the region that lies between the curve y=f(x) and the x-axis between x=a and x=b.

If f(x) takes both positive and negative values on [a,b], then Sn is an approximation for the net signed area; that is, the area below the curve y=f(x) and above the x-axis minus the area above the curve y=f(x) and below the x-axis from x=a to x=b, as suggested in Figure 2.

 

Figure 2
Figure 2. The Riemann sum is an approximation to the net area if f takes on both positive and negative values.

 

The Riemann sum is often shortened as Sn=k=1nf(xk)(xkxk1), where Σ (sigma), Greek form of the letter S, stands for the word “summation.” The whole expression indicates that the sum is to be taken of all terms obtained from f(xk)(xkxk1) by giving to k in succession the values of 1,2,3,,n.

By introducing the symbol Δkx=xkxk1, we can shorten it further and write: (8.1.2)Sn=k=1nf(xk)Δkx. Please note that here Δk is not a factor, but denotes the difference (here the length of the k-th subinterval).

  • Some books, write Δxk instead of Δkx.

 


Read about some specific types of Riemann sums


 

 

Limit of the Riemann Sums as n

As shown in Figure 8, if the number of subintervals n increases such that the subintervals lengths (or widths) get smaller and the rectangles become thinner, then the Riemann sum Sn will be a better approximation to the net area between the curve y=f(x) and the x-axis.

n=11 n=20 n=40

Figure 8. As the number of subintervals n increases, Sn becomes a better approximation to the net area that lie between the curve y=f(x) and the x-axis.

 

Now let the number of subintervals n increase indefinitely (n), while the width of the widest subinterval approaches zero. If the sum (8.1.1) (or (8.1.2)) approaches a number which is independent of the choice of the xk’s and of the xk’s, its limit is called the definite integral of f(x) between a and b and is denoted by abf(x)dx.

Definition 1. The limit of the sum (8.1.2) provided that the number of subintervals n tends to infinity and at the same time the length of each subinterval approaches zero is called the definite integral of f(x) from a to b and is denoted by abf(x)dx. If this limit exists, the function is called integrable on the interval [a,b]. abf(x)dx=limmaxΔkx0k=1nf(xk)Δkx.

We note that

  • The use of the word “integral” and of the symbol suggests a connection with the indefinite integrals in the previous Chapter. This connection will be shown in Section on the Fundamental Theorem of Calculus.

 

  • The symbol is a modified form of S (that stands for summation).

 

  • We call a the lower limit and b the upper limit of integration (or of the internal), and the function f(x) the integrand.

 

  • If the integral abf(x)dx exists, f(x) is called integrable between a and b.

 

  • If f is continuous on the interval (a,b) then abf(x)dx exists. The proof of this theorem is beyond the scope of an elementary course.

 


The variable of integration is a dummy variable


 

In the previous discussion, we assumed a<b. However, if we omit this condition and assume b<a, we can still retain our definition of integral; the only change is that we have to number the subdivision points from the right side, which make the differences Δkx negative. We are thus led to the relation baf(x)dx=abf(x)dx. In conformity, we define aaf(x)dx=0.

Definition 2.

(a and b interchanged)abf(x)dx=baf(x)dx. and (same upper and lower limits)aaf(x)dx=0.

abf(x)dx= area above the x-axis area below the x-axis   (when a<b)

 

Example 1

Find 13(2x)dx.

Solution 1

Let divide the interval [1,3] into n subintervals of equal width Δx=(3(1))/n=4/n. In each subinterval [xk,xk+1], let use xk=xk as the sample point x0=1x1=1+Δx=1+4nx2=1+2Δx=1+2(4n)xk=1+kΔx=1+k(4n)xn=1+nΔx=1+n(4n)=1+4=3. Therefore, f(xk)=f(xk)=2xk=2[1+k(4n)]=34kn and k=1nf(xk)Δkx=k=1n(34kn)4n=k=1n(12n16kn2)=12nk=1n116n2k=1nk=12n(n)16n2n(n+1)2=128(1+1n) Notice that in the above summation, n is a fixed number (say n=10) and not a function of the summation index k. That is why we could move n outside the sigma sign. For these calculations, we have also used the formulas k=1n1=n and k=1nk=n(n+1)/2 that we reviewed in the Section on Sigma Notation.

Because the width of each subinterval is the same maxΔkx=4/n0 is equivalent to n, and we can conclude 13(2x)dx=limnk=1nf(xk)Δx=limn[128(1+1n)]=128(1+0)=4.

You can easily show that if we use xk=xk+1 (instead of xk=xk) as the sample point in each subinterval [xk,xk+1], we will arrive at the same result.

Because the graph of y=2x is a straight line, we can easily check our answer. We should subtract the area of the part that lies below the x-axis from the area of the above that is above the x-axis (see Figure 9). It follows from the familiar triangle area formula A=base×height/2 that 13(2x)dx=3×321×12=4.

Figure 9
Figure 9. The integral of 2x from 1 to 3 is 13(2x)dx=A1A2=9/21/2=4.
Example 2

Find 13x2dx.

Solution 2

Let’s divide the interval [1,3] into n subintervals of equal width Δkx=Δx=(31)n=2n: x0=1x1=1+Δx=1+2nx2=1+2Δx=1+2(2n)xk=1+kΔx=1+k(2n)xn=1+nΔx=1+n(2n)=3. Assume f(xk)=f(xk)=(1+2kn)2=1+4kn+4k2n2. Then k=1nf(xk)Δx=k=1n[2n+8kn2+8k2n3]=2nk=1n1+8n2k=1nk+8n3k=1nk2=2nn+8n2n(n+1)2+8n3n(n+1)(2n+1)6=2+4n+1n+43(n+1n)(2n+1n)=2+4(1+1n)+43(1+1n)(2+1n) The width of every subinterval is the same and maxΔkx0 is equivalent to n. As n, 1n0 and hence 13x2dx=limnk=1nf(xk)Δx=limn[2+4(1+1n)+86(1+1n)(2+1n)]=2+4(1+0)+43(1+0)(2+0)=263.

There are no formulas from elementary geometry that can help us verify our answer (see Figure 10), but we will soon learn a theorem, called the Fundamental Theorem of Calculus, that enables us to evaluate such definite integrals easily.

Figure 10
Figure 10.