Last modified: January 27, 2025

This article is written in: 🇺🇸

Trapezoidal Rule

The Trapezoidal Rule is a fundamental numerical integration technique employed to approximate definite integrals, especially when an exact antiderivative of the function is difficult or impossible to determine analytically. This method is widely used in various fields such as engineering, physics, and economics due to its simplicity and effectiveness. The core idea behind the Trapezoidal Rule is to estimate the area under a curve by dividing it into a series of trapezoids, rather than relying solely on simpler geometric shapes like rectangles. By approximating the region under the graph of a function with trapezoids, this method provides a better fit to the actual curve, thereby enhancing the accuracy of the integral approximation compared to basic methods like the Rectangle (Midpoint) Rule.

The Trapezoidal Rule is particularly advantageous for functions that are relatively smooth and do not exhibit extreme fluctuations or sharp bends within the interval of integration. By constructing trapezoids that align with the function's endpoints, the method captures the linear trend between successive points, effectively averaging the function's behavior over each subinterval. This approach reduces the potential for significant overestimation or underestimation of the area, leading to more reliable integral approximations. Additionally, the Trapezoidal Rule is easy to implement computationally, making it a practical choice for both manual calculations and computer-based algorithms.

Mathematical Formulation

The mathematical foundation of the Trapezoidal Rule is based on the assumption that the region under the graph of a function $f\left(x\right)$ over a given interval $[a,b]$ can be approximated by a trapezoid. A trapezoid is a quadrilateral with at least one pair of parallel sides, which in this context are the lines connecting the points $(a,f(a\left)\right)$ and $(b,f(b\left)\right)$. The area of this trapezoid provides an estimate of the definite integral of the function over the interval $[a,b]$.

The basic formula for the Trapezoidal Rule when applied to a single interval $[a,b]$ is:

$$\frac{b-a}{2}\left[f\right(a)+f(b\left)\right]$$

This formula calculates the area of the trapezoid by taking the average of the function values at the endpoints $f\left(a\right)$ and $f\left(b\right)$, and then multiplying by the width of the interval $(b-a)$. This approximation assumes that the function behaves linearly between the two endpoints, which is a reasonable assumption for small intervals or functions that do not vary drastically.

To improve the accuracy of the approximation, the interval $[a,b]$ can be divided into $n$ smaller subintervals of equal width. Let $h$ represent the width of each subinterval, calculated as:

$$h=\frac{b-a}{n}$$

Within each subinterval $[x_i,x_{i+1}]$, the function $f\left(x\right)$ is evaluated at the endpoints $x_i$ and $x_{i+1}$. The area of each trapezoid within these subintervals is then computed using the formula:

$$\frac{f\left(x_{i+1}\right)+f\left(x_i\right)}{2}h$$

The integral $I$ over the entire interval $[a,b]$ is then approximated by summing the areas of all the trapezoids:

$$I=\sum _{i=0}^{n-1}\frac{f\left(x_{i+1}\right)+f\left(x_i\right)}{2}h$$

This summation effectively aggregates the contributions of each trapezoid to estimate the total area under the curve $f\left(x\right)$ across the interval. The graphical representation below illustrates how the Trapezoidal Rule approximates the area under a curve by stacking trapezoids side by side:

In this diagram, each trapezoid corresponds to a subinterval $[x_i,x_{i+1}]$, with the parallel sides representing the function values $f\left(x_i\right)$ and $f\left(x_{i+1}\right)$. The height of each trapezoid is the width $h$ of the subinterval, and the area of each trapezoid contributes to the overall approximation of the integral.

Algorithm Steps

Implementing the Trapezoidal Rule involves a systematic process to ensure an accurate approximation of the definite integral. The following steps outline the procedure:

I. Partition the Interval:

• Begin by dividing the interval $[a,b]$ into $n$ equal subintervals. The number of subintervals $n$ determines the granularity of the approximation; a larger $n$ generally leads to a more accurate result.
• Calculate the width of each subinterval using the formula:

$$h=\frac{b-a}{n}$$

This step ensures that the interval is evenly partitioned, facilitating consistent application of the Trapezoidal Rule across all subintervals.

II. Evaluate Function at Endpoints:

• For each subinterval $[x_i,x_{i+1}]$, compute the function values at the endpoints $f\left(x_i\right)$ and $f\left(x_{i+1}\right)$.
• These evaluations are crucial as they determine the heights of the trapezoids that approximate the area under the curve within each subinterval.
• Organizing these function values systematically allows for efficient computation of the trapezoidal areas.

III. Apply the Trapezoidal Rule Formula:

• Use the Trapezoidal Rule formula to calculate the area of each trapezoid within the subintervals:

$$\frac{f\left(x_{i+1}\right)+f\left(x_i\right)}{2}h$$

• This formula takes the average of the function values at the endpoints and multiplies it by the width $h$ to determine the area of the trapezoid.
• Applying this formula uniformly across all subintervals ensures consistency and accuracy in the approximation process.

IV. Aggregate the Results:

Sum the areas of all the trapezoids computed in the previous step to obtain the total integral approximation:

$$I=\sum _{i=0}^{n-1}\frac{f\left(x_{i+1}\right)+f\left(x_i\right)}{2}h$$

• This aggregation provides the final estimate of the definite integral over the interval $[a,b]$.
• The sum effectively captures the cumulative area under the curve as approximated by the series of trapezoids.

By meticulously following these steps, the Trapezoidal Rule systematically approximates the integral by leveraging the linear behavior of the function within each subinterval. This method balances simplicity and accuracy, making it a practical choice for a wide range of numerical integration problems.

Example

To illustrate the application of the Trapezoidal Rule, consider the function $f\left(x\right)=x^2$. This example demonstrates how the method approximates the integral over a specific interval using a manageable number of subintervals.

I. Choose Interval and Subdivisions:

• Let $a=0$ and $b=2$, and partition this interval into $n=2$ equal subintervals. This choice ensures that the interval is evenly divided, facilitating the application of the Trapezoidal Rule.
• Calculate the width of each subinterval:

$$h=\frac{2-0}{2}=1$$

This results in two subintervals: $[0,1]$ and $[1,2]$, each with a width of 1.

II. Evaluate Function at Endpoints:

Compute the function values at the endpoints of each subinterval:

For the first subinterval $[0,1]$:

• $f\left(0\right)=0^2=0$
• $f\left(1\right)=1^2=1$

For the second subinterval $[1,2]$:

• $f\left(1\right)=1^2=1$
• $f\left(2\right)=2^2=4$

These evaluations provide the necessary data points for calculating the areas of the trapezoids.

III. Apply the Trapezoidal Rule Formula:

Calculate the area of each trapezoid using the Trapezoidal Rule formula:

For $[0,1]$:

$$\frac{1-0}{2}\left[f\right(0)+f(1\left)\right]=\frac{1}{2}(0+1)=0.5$$

For $[1,2]$:

$$\frac{2-1}{2}\left[f\right(1)+f(2\left)\right]=\frac{1}{2}(1+4)=2.5$$

These calculations yield the areas of the trapezoids corresponding to each subinterval.

IV. Sum the Areas:

Add the areas of the individual trapezoids to obtain the total integral approximation:

$$I=0.5+2.5=3$$

This sum represents the estimated value of the definite integral of $f\left(x\right)=x^2$ over the interval $[0,2]$ using the Trapezoidal Rule.

V. Compare with Exact Integral:

For validation, compute the exact value of the integral analytically:

$${\int }_0^2{x}^{2}dx={\left[\frac{x^3}{3}\right]}_0^2=\frac{8}{3}\approx 2.6667$$

The Trapezoidal Rule approximation of 3 is reasonably close to the exact value of approximately 2.6667, demonstrating the method's effectiveness even with a small number of subintervals.

This example highlights the practical application of the Trapezoidal Rule and illustrates how it provides a balance between computational simplicity and accuracy. While the approximation is not exact in this case, it is sufficiently close for many practical purposes, especially when dealing with functions that do not exhibit extreme behavior within the interval of integration.

Advantages

• The Trapezoidal Rule is straightforward to understand and implement, making it accessible to students and practitioners alike. The method relies on basic arithmetic operations and function evaluations at specific points, without requiring complex algorithms or iterative procedures.
• Compared to the Rectangle Rule, which approximates the area under a curve using rectangles based on function values at specific points (such as midpoints), the Trapezoidal Rule often provides a more accurate approximation. This increased accuracy stems from the method's ability to account for the linear trend between two consecutive points, effectively averaging the function's behavior over each subinterval.
• The Trapezoidal Rule is versatile and can be applied to a broad spectrum of functions, including those that are linear, polynomial, or exhibit smooth curvature. It is particularly effective for functions that do not have extreme fluctuations or sharp bends within the interval of integration.
• The method is computationally efficient, especially for functions that are computationally inexpensive to evaluate. Since the Trapezoidal Rule requires only two function evaluations per subinterval, it scales well with increasing numbers of subintervals, maintaining reasonable computational demands even for large-scale problems.
• The Trapezoidal Rule serves as a foundational technique in numerical integration, providing a basis for more advanced methods such as Simpson's Rule and Gaussian Quadrature. Understanding the principles of the Trapezoidal Rule facilitates the comprehension and implementation of these more sophisticated techniques.

Limitations

• The Trapezoidal Rule assumes that the function behaves linearly within each subinterval, which can lead to significant approximation errors for functions that exhibit high curvature or rapid changes within the interval of integration. In such cases, the linear approximation may not capture the function's true behavior, resulting in overestimation or underestimation of the area under the curve.
• Achieving a high level of precision with the Trapezoidal Rule often requires dividing the interval $[a,b]$ into a large number of subintervals. While increasing the number of subintervals can enhance the accuracy of the approximation, it also increases the computational effort, as more function evaluations and calculations are necessary.
• The accuracy of the Trapezoidal Rule is highly dependent on the behavior of the function being integrated. For functions that are not smooth or that have irregular behavior within the interval, the linear approximations used by the Trapezoidal Rule may not adequately capture the true nature of the function.
• When dealing with a large number of subintervals, especially in computational implementations, the Trapezoidal Rule can be susceptible to the accumulation of numerical errors. These errors may arise from factors such as round-off errors in floating-point arithmetic or from inaccuracies in function evaluations.
• The Trapezoidal Rule is primarily designed for definite integrals over finite intervals. It is not inherently suited for improper integrals, where the interval of integration is infinite or where the integrand has singularities within the interval.
• The method assumes that the function $f\left(x\right)$ is continuous over the interval $[a,b]$. For functions with discontinuities or points of non-differentiability within the interval, the Trapezoidal Rule may fail to provide a meaningful approximation.