Approximations of the Riemann Integral in R

Subjects: Real Analysis
Links: Riemann Integral in R

Equal Partitions

Let Pn be uniform partitions of n elements. Evaluating the Reimann sums at left or right endpoints is an important construction. They receive a particular notation, where hn=ban:

Ln(f)=hnk=0n1f(a+khn)Rn(f)=hnk=1nf(a+khn)

If we don’t have a reason to prefer either one of them, it’s reasonable to take the mean we get:

Tn(f)=Ln(f)+Rn(f)2=hn(f(a)2+k=1n1f(a+khn)+f(b)2)

In particular, if f is monotonic then the integral will be trapped between the two, and we will get that the average will differ of the integral by a small amount:

|abfTn(f)|ban|f(b)f(a)2|

Trapezoidal Rule

Continue to work with equal partitions, Pn and a width of hn=ban . Let’s consider the function gn is a piecewise function, of a linear approximation of f at the interval [a+(k1)hn,a+khn]=Ik. In other words:

Ikgn=12hn[f(a+(k1)hn)+f(a+khn)]

Then the integral along all [a,b] will be equal to:

abgn=k=1nIkgn=12hnk=1nf(a+(k1)hn)+f(a+khn)=hn(f(a)2+k=1nf(a+khn)+f(b)2)=Tn(f)

Tn(f) is the n-th Trapezoidal Approximation.

Theorem: If fC2[a,b], then there’s a c[a,b]:

Tn(f)abf=(ba)hn212f(c)

Corollary: If fC2[a,b] and let |f|B2, then:

|abfTn(f)|(ba)312n2B2

Midpoint Rule

Mn(f)=hnk=1nf(a+(k12)hn)

Theorem: if fC2[a,b], then there’s a c[a,b] that:

abfMn(f)=(ba)hn224f(c)

Corollary: If fC2[a,b] and let |f|B2, then:

|abfMn(f)|(ba)324n2B2

How to transform knowledge of Tn(f) and Mn(f) into a type of Tn(f):

T2n(f)=12Mn(f)+12Tn(f)

Simpson’s Rule

Given a n2N, the Simpson Sum is defined as follows:

Sn(f)=hn3(k=1n/2f(a+(2k2)hn)+4f(a+(2k1)hn)+f(a+2khn))

Theorem: If fC4[a,b], there’s a c[a,b] that:

Sn(f)abf=(ba)5180n4f(4)(c)

Corollary: If fC4[a,b], and |f(4)|B4, then:

|abfSn(f)|(ba)5180n4B4

How to transform knowledge of Tn(f) and Mn(f) into a type of Sn(f):

S2n(f)=23Mn(f)+13Tn(f)