Search This Blog

Monday, November 14, 2011

The Trapezoidal Rule

Problem: Find
math expression
We put u = x2 + 1 and as x = 0 → 1,u = 1 → 2
If u = x2 + 1 then du = 2x dx.
But the question does not contain an x dxterm so we cannot solve it using any of the normal integration methods.
We need to use numerical approaches. When software like Mathcad or graphics calculators perform definite integrals, they use numerical methods given earlier.
We can use one of two methods:

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.
math expression
Recall that we write "Δx" to mean "a small change in x".
Let's see this in LiveMath. Note that our approximation is much better than using rectangles.
Now, the area of a trapezoid (trapezium) is given by:
math expression
math expression
So the approximate area under the curve is found by adding the area of the trapezoids. (Our trapezoids are rotated 90° so that their new base is actually the height. So h = Δx.)
Area ≈
math expression
We can simplify this to give us the Trapezoidal Rule, for n trapezoids:
math expression
To find Δx for the area from x = a to x = b, we use:
math expression
and we also need
y0 = f(a)
y1 = f(a + Δx)
y2 = f(a + 2Δx)
yn = f(b)

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

Exercise: Using n = 5, approximate the integral:
math expression
Here, a = 0 and b = 1.
math expression
y0 = f(a) = f(0) = math expression = 1
y1 = f(a + Δx) = f(0.2) = math expression
y2 = f(a + 2Δx) = f(0.4) = math expression
y3 = f(a + 3Δx) = f(0.6) = math expression
y4 = f(a + 4Δx) = f(0.8) = math expression
y5 = f(b) = f(1) = math expression
So the area ≈
math expression
So math expression 1.150

No comments:

Post a Comment