Numerical Integration

Definition

Numerical integration is a fundamental technique in numerical analysis used to calculate the approximate numerical value of a definite integral. When an integrand is difficult or impossible to integrate analytically (using standard calculus rules), or when the function is known only through a set of discrete data points, numerical integration methods allow us to estimate the area under the curve.

Main Content

1. The Trapezoidal Rule

This method approximates the region under the graph of the function as a trapezoid and calculates its area.
It is based on a first-order polynomial approximation, where the function is assumed to be linear between two consecutive points.

2. Simpson’s 1/3 Rule

This method provides a more accurate approximation by using a second-order polynomial (a parabola) to connect three consecutive points.
It requires an even number of sub-intervals (or an odd number of data points) to function correctly.

3. Simpson’s 3/8 Rule

This is a variation of Simpson’s rule that uses a third-order polynomial to connect four consecutive points.
It is often used when the number of sub-intervals is a multiple of three, providing higher accuracy than the 1/3 rule for specific curve types.

Working / Process

1. Partitioning the Interval

Divide the interval $[a, b]$ into $n$ equal sub-intervals of width $h$, where $h = \frac{b - a}{n}$.
Define the coordinates as $x_0, x_1, x_2, \dots, x_n$ such that $x_i = a + i \cdot h$.

2. Evaluating Function Values

Calculate the height of the function at each partition point: $y_i = f(x_i)$.
Create a table of values to organize these outputs for easier calculation.

3. Applying the Integration Formula

Insert the values into the chosen integration formula. For the Trapezoidal Rule: $I \approx \frac{h}{2} [y_0 + 2(y_1 + y_2 + \dots + y_{n-1}) + y_n]$.
Sum the weighted values to obtain the final approximate area.

Visual representation of the Trapezoidal Rule:
y |    /|
  |   / |
  |  /  |
  | /   |
  |/____|__ x
    a   b
(Area of the trapezoid approximates the integral)

Advantages / Applications

Allows for solving integrals of complex functions where no antiderivative exists (e.g., $e^{-x^2}$).
Enables processing of experimental data where the function formula is unknown but data points are available.
Essential in engineering simulations, physics modeling, and financial data analysis where high-speed computation is required.

Summary

Numerical integration is the process of approximating the definite integral of a function using discrete points rather than symbolic integration. It simplifies complex mathematical problems into manageable arithmetic operations.

Key terms:

Interval: The range $[a, b]$ over which integration occurs.
Sub-interval: Smaller segments of width $h$ used to approximate the area.
Nodes: The discrete $x$-coordinates where the function value is evaluated.