Numerical Integration: Trapezoidal Rule, Simpson’s 1/3rd, and 3/8th Rules

Definition

Numerical integration (often called quadrature) is the set of algorithms used to calculate the definite integral of a function when an analytical solution is difficult or impossible to obtain, or when the data is provided as a discrete set of points. The Trapezoidal Rule, Simpson’s 1/3rd Rule, and Simpson’s 3/8th Rule are closed Newton-Cotes formulas that approximate the area under a curve by dividing the interval into sub-intervals and approximating the function within those intervals using polynomials.

Main Content

1. The Trapezoidal Rule

This rule approximates the area under a curve by dividing the total area into a series of trapezoids. It connects each pair of consecutive points with a straight line.
It is a first-order method, meaning it assumes the function behaves linearly between points.

2. Simpson’s 1/3rd Rule

This rule approximates the area by using second-order (parabolic) polynomials to connect three consecutive points at a time.
It requires an even number of sub-intervals (odd number of data points) to function correctly.

3. Simpson’s 3/8th Rule

This rule uses third-order (cubic) polynomials to connect four consecutive points.
It is generally more accurate than the 1/3rd rule for higher-order curves and requires the number of sub-intervals to be a multiple of three.

Visual representation of approximation methods:

Trapezoidal (Linear)     Simpson's 1/3rd (Parabolic)
      |  /                   |  .-'|
      | /                    | /   |
      |/                     |/    |
      +----                  +-----+

Working / Process

1. Pre-processing and Subdivision

Divide the interval $[a, b]$ into $n$ equal sub-intervals of width $h$, where $h = (b - a) / n$.
Identify the functional values $y_0, y_1, y_2, \dots, y_n$ at each equidistant point $x_0, x_1, \dots, x_n$.

2. Applying the Formula

For Trapezoidal: $\int_a^b f(x)dx \approx \frac{h}{2} [ (y_0 + y_n) + 2(y_1 + y_2 + \dots + y_{n-1}) ]$.
For Simpson's 1/3rd: $\int_a^b f(x)dx \approx \frac{h}{3} [ (y_0 + y_n) + 4(\text{sum of odd } y) + 2(\text{sum of even } y) ]$.
For Simpson's 3/8th: $\int_a^b f(x)dx \approx \frac{3h}{8} [ (y_0 + y_n) + 3(y_1 + y_2 + y_4 + y_5 + \dots) + 2(y_3 + y_6 + \dots) ]$.

3. Calculation and Error Estimation

Compute the sum using the selected rule.
Calculate the error term; note that the error in Trapezoidal rule is $O(h^2)$, while Simpson's rules offer higher-order accuracy.

Advantages / Applications

Engineering: Used to calculate work done by variable forces or cross-sectional areas of civil structures.
Data Analysis: Essential for processing experimental data where the underlying mathematical function $f(x)$ is unknown.
Efficiency: These methods are computationally inexpensive compared to symbolic integration, making them ideal for computer-based numerical solvers.

Summary

Numerical integration techniques are computational methods used to approximate definite integrals by partitioning an area into geometric shapes. The Trapezoidal rule uses linear segments, while Simpson’s rules use quadratic and cubic splines for higher precision. Essential terms include the step size ($h$), nodes (data points), and the order of accuracy (error term). These methods are fundamental tools in computational mathematics and engineering for solving integrals that lack closed-form solutions.