# 5. The Trapezoidal Rule

by M. Bourne

### Interactive exploration

See an applet where you can explore Simpson's Rule and other numerical techniques:

Problem: Find

int_0^1sqrt(x^2+1)\ dx

We put u = x^2+ 1 then du = 2x\ dx.

But the question does not contain an x\ dx term so we cannot solve it using any of the integration methods we have met so far.

We need to use numerical approaches. (This is usually how software like Mathcad or graphics calculators perform definite integrals).

We can use one of two methods:

Continues below

## The Trapezoidal Rule

We saw the basic idea in our first attempt at solving the area under the arches problem earlier.

Instead of using rectangles as we did in the arches problem, we'll use trapezoids (trapeziums) and we'll find that it gives a better approximation to the area.

y0
Deltax
y1
Deltax
y2
Deltax
y3
Deltax
y4
Deltax
y5
Deltax
y6

Area under curve using trapezoidal Rule

The approximate area under the curve is found by adding the area of all the trapezoids.

(Recall that we write "Δx" to mean "a small change in x".)

### Area of a trapezoid

p
q
h
"Area"=h/2(p+q)

Area of a trapezoid (trapezium)

Now, the area of a trapezoid (trapezium) is given by:

"Area"=h/2(p+q)

We need "right" trapezoids (which means the parallel sides are at right angles to the base), and they are rotated 90° so that their new base is actually h, as follows, and h = Δx.

y0
y1
Deltax

A "typical" trapezoid

So the total area is given by:

"Area"~~1/2(y_0+y_1)Deltax+ 1/2(y_1+y_2)Deltax+ 1/2(y_2+y_3)Deltax+...

We can simplify this to give us the Trapezoidal Rule, for n trapezoids:

"Area"~~ Deltax((y_0)/2+y_1+y_2+y_3+ {:...+(y_n)/2)

To find Δx for the area from x=a to x=b, we use:

Deltax=(b-a)/n

and we also need

y_0= f(a)

y_1= f(a + Δx)

y_2= f(a + 2Δx)

...

y_n= f(b)

### Note

• We get a better approximation if we take more trapezoids [up to a limit!].
• The more trapezoids we take, Delta x will tend to 0, that is, Δx → 0.
• We can write (if the curve is above the x-axis only between x=a and x=b):

"Area"=int_a^b f(x)dx

~~ Deltax((y_0)/2+y_1+{:...+(y_n)/2)

## Don't miss...

There is an interactive applet where you can explore the Trapezoid Rule, here:

Calculus from First Principles applet

## Exercise

Using n= 5, approximate the integral:

int_0^1sqrt(x^2+1)\ dx