Mc34, Mc35 - Trapezium rule and area bounds

Mc34, Mc35 - Trapezium rule and area bounds

Could a zero estimate hide a positive area? Under Mc34, you will construct and interpret trapezium estimates and justify their error direction on an interval of consistent concavity. Under Mc35, you will use rectangle sums for bounds, then read an unfamiliar expression for the same trapezium method.

Build the composite trapezium rule

The curve 16x2(1−x)216x^2(1-x)^2 has height zero at both ends of [0,1][0,1], so a single trapezium gives zero. Yet the curve is positive inside and its exact integral is 8/158/15. Sampling can miss an entire hill!

You already know

A trapezium of width hh and parallel heights p,qp,q contributes h(p+q)/2h(p+q)/2. Definite integrals are signed; geometric area needs magnitudes if the curve crosses the axis.

Use nn equal strips of width h=(b−a)/nh=(b-a)/n. There are n+1n+1 ordinates yi=f(a+ih)y_i=f(a+ih), including both endpoints.

Trapezium rule

For equally spaced ordinates,

Tn=h2[y0+2(y1+⋯+yn−1)+yn]≈∫abf(x) dx.T_n=\frac h2[y_0+2(y_1+\cdots+y_{n-1})+y_n] \approx\int_a^b f(x)\,\mathrm dx.

In the formula sheet.

Every interior ordinate belongs to two adjacent strips, explaining its weight 22.

Worked example 1: applying the weights

Question. Estimate ∫01ex dx\int_0^1\mathrm e^x\,\mathrm dx with two equal strips.

Plan. Calculate three heights, then form the weighted sum.

h=12,(y0,y1,y2)=(1,e1/2,e)T2=14[1+2e1/2+e]=1.753931092…\begin{aligned} h&=\frac12,\qquad (y_0,y_1,y_2)=(1,\mathrm e^{1/2},\mathrm e)\\ T_2&=\frac14[1+2\mathrm e^{1/2}+\mathrm e]\\ &=1.753931092\ldots \end{aligned}

Why three heights? Two strips share the middle ordinate but each needs two boundary heights.

Answer. The estimate is 1.753931.75393 to five decimal places.

Check. The exact integral is e−1≈1.71828\mathrm e-1\approx1.71828, so this is an overestimate.

Change only the strip count to four.

Worked example 2: more strips

Question. Estimate the same integral using four equal strips.

Plan. Insert the quarter-step heights and retain their full values.

h=14T4=18[1+2(e1/4+e1/2+e3/4)+e]=1.727221905…\begin{aligned} h&=\frac14\\ T_4&=\frac18[1+2(\mathrm e^{1/4}+\mathrm e^{1/2}+\mathrm e^{3/4})+\mathrm e]\\ &=1.727221905\ldots \end{aligned}

Why retain the middle height? Refining adds ordinates; it does not discard those already sampled.

Answer. The estimate is 1.727221.72722 to five decimal places.

Check. It is closer to e−1\mathrm e-1 than the two-strip estimate.

For the polynomial x2+1x^2+1 on [0,2][0,2], five heights 1,1.25,2,3.25,51,1.25,2,3.25,5 give T4=4.75T_4=4.75; the individual trapezia 0.5625,0.8125,1.3125,2.06250.5625,0.8125,1.3125,2.0625 give the same total. Another table, 2,3,6,112,3,6,11 at 0,1,2,30,1,2,3, gives 15.515.5.

Exam tip

If instructed to use the trapezium rule, show a substituted expression with every ordinate. “Determine” requires justified working, and four strips require five heights.

Refine the exponential example once more.

Use technology and refine the estimate

For e−x2\mathrm e^{-x^2}, elementary antiderivatives do not supply a formula, so a calculator table is especially useful. A spreadsheet can store inputs, ordinates, weights and weighted values; use full stored precision and round only the final total.

Worked example 3: a table without an elementary antiderivative

Question. Estimate ∫01e−x2 dx\int_0^1\mathrm e^{-x^2}\,\mathrm dx with two equal strips.

Plan. Evaluate the function at 0,0.5,10,0.5,1, then apply endpoint/interior weights.

(y0,y1,y2)=(1,e−1/4,e−1)T2=14[1+2e−1/4+e−1]=0.731370252…\begin{aligned} (y_0,y_1,y_2)&=(1,\mathrm e^{-1/4},\mathrm e^{-1})\\ T_2&=\frac14[1+2\mathrm e^{-1/4}+\mathrm e^{-1}]\\ &=0.731370252\ldots \end{aligned}

Why evaluate the function, rather than use the inputs as heights? The inputs give strip positions; the outputs supply areas.

Answer. T2≈0.731370T_2\approx0.731370.

Check. Adding the two individual trapezia gives the same total.

Change only the strip count to four.

Worked example 4: refining the table

Question. Repeat with four equal strips.

Plan. Add values at 0.250.25 and 0.750.75 without rounding the stored ordinates.

T4=18[1+2(e−1/16+e−1/4+e−9/16)+e−1]=0.742984098…\begin{aligned} T_4&=\frac18[1+2(\mathrm e^{-1/16}+\mathrm e^{-1/4}+\mathrm e^{-9/16})+\mathrm e^{-1}]\\ &=0.742984098\ldots \end{aligned}

Why does the changed answer not alone give an error bound? It compares two approximations; neither has yet been proved to enclose the integral.

Answer. T4≈0.742984T_4\approx0.742984.

Check. Eight strips give 0.745865615…0.745865615\ldots, showing further movement rather than exact agreement.

For unequal spacings, use separate widths. Readings (0,4),(1,6),(3,2)(0,4),(1,6),(3,2) give 12(4+6)(1)+12(6+2)(2)=13\frac12(4+6)(1)+\frac12(6+2)(2)=13. A common-width formula is invalid there. Do not invent new measured heights when refining a sparse data table: take more readings or state an interpolation assumption.

Watch out

More printed digits do not mean more reliable data. A narrow peak between sampled inputs can remain hidden, as the opening zero-height example shows.

Keep the same Gaussian function and refine once more.

Decide the direction of trapezium error

The curve's concavity determines which side of each chord it lies on. Increasing or decreasing behaviour alone does not decide trapezium bias.

Trapezium chords above and below curves

The diagram shows x2+1x^2+1, whose chords lie above it, and 4−x24-x^2, whose chords lie below it. Four strips give 4.75>14/34.75>14/3 and 5.25<16/35.25<16/3, respectively.

Concavity and the error direction

On an interval of consistent concavity:

f′′(x)≥0⇒Tn≥Iconcave upwardsf′′(x)≤0⇒Tn≤Iconcave downwards\begin{aligned} f''(x)\ge0&\Rightarrow T_n\ge I && \text{concave upwards}\\ f''(x)\le0&\Rightarrow T_n\le I && \text{concave downwards} \end{aligned}

where I=∫abf(x) dxI=\int_a^b f(x)\,\mathrm dx. The chords lie above the curve in the first case and below in the second.

Not in the formula sheet: learn it.

Use an interval-wide reason, then state “overestimate” or “underestimate” explicitly.

Worked example 5: consistent upward concavity

Question. Decide the trapezium error direction for f(x)=2+1/(1+x)f(x)=2+1/(1+x) on [0,2][0,2].

Plan. Differentiate twice and inspect the sign.

f′′(x)=2(1+x)3>0(0≤x≤2).f''(x)=\frac{2}{(1+x)^3}>0\qquad(0\le x\le2).

Why check the whole interval? One positive second derivative at a single point does not establish consistent concavity.

Answer. The curve is concave upwards, its chords lie above it, and the integral is overestimated.

Check. The denominator is positive everywhere in the interval; there is no excluded point.

Change only the sign of the reciprocal term.

Worked example 6: reversing concavity

Question. Decide the error direction for f(x)=2−1/(1+x)f(x)=2-1/(1+x) on the same interval.

Plan. Reverse the sign of the second derivative.

f′′(x)=−2(1+x)3<0.f''(x)=-\frac{2}{(1+x)^3}<0.

Why does the chord position reverse? Reflecting the varying term reverses the direction of curvature.

Answer. The curve is concave downwards; chords lie below it, giving an underestimate.

Check. Both functions remain positive here, so the signed-integral comparison also describes their geometric areas.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account