Gauss quadrature for integrals and sums

A. D. Alhaidari

DOI: https://doi.org/10.51483/IJPAMR.3.1.2023.1-10

2023-05-29

Abstract:Gauss quadrature integral approximation is extended to include integrals with a measure consisting of continuous as well as discrete components. That is, we give an approximation for the integral of a function plus its sum over a discrete weighted set.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to extend Gauss quadrature so that it can handle integral measures containing continuous and discrete components. Specifically, the author proposes a Gauss quadrature method for mixed measures to simultaneously approximate the integral of a function and its sum on a discrete weighted set.

### Core of the Problem

1. **Limitations of Traditional Gauss Quadrature**
- Traditional Gauss quadrature is mainly used to handle integral problems with pure continuous measures.
- For integral or summation problems containing discrete components, the traditional method cannot be directly applied.

2. **Requirement for Mixed Measures**
- In many practical application scenarios, the integral measure may contain both continuous and discrete parts. For example, in some physical and engineering problems, it is necessary to handle continuous distributions and discrete point sources simultaneously.
- Therefore, it is very necessary to develop a quadrature method that can handle such mixed measures.

### Solution

The author extends Gauss quadrature through the following steps:

1. **Introducing Mixed Measures**
- Represent the integral measure as a combination of continuous and discrete parts:
\[
\mu(x)=\sigma(x)+\sum_{k = 0}^{M}\xi_k\delta(x - x_k)
\]
where \(\sigma(x)\) is the continuous part, and \(\xi_k\) and \(x_k\) are the weights and positions of the discrete part respectively.

2. **Orthogonal Polynomials and Recurrence Relations**
- Use the orthogonal polynomials related to the mixed measure and their recurrence relations to construct an infinite - dimensional symmetric tridiagonal matrix (Jacobi matrix).
- These polynomials satisfy the following three - term recurrence relation:
\[
x p_n(x)=a_n p_{n + 1}(x)+b_n p_n(x)+a_{n - 1}p_{n - 1}(x)
\]

3. **Numerical Approximation**
- Truncate the above matrix to a finite dimension and calculate its eigenvalues and eigenvectors to obtain Gauss quadrature nodes and weights.
- Finally, transform the integral and summation problems into a finite summation form:
\[
\int_{a}^{b}f(x)d\mu(x)\approx\sum_{n = 0}^{N - 1}\omega_n f(\epsilon_n)+\sum_{k = 0}^{M}\xi_k f(x_k)
\]

### Application Examples

The paper provides several specific examples to verify the effectiveness and accuracy of the proposed method, including:

- **Charlier and Meixner Polynomials**: Used to approximate infinite series.
- **Krawtchouk Polynomials**: Used to approximate finite series.
- **Continuous Double Hahn Polynomials**: Used to handle integral - sum combinations containing continuous and discrete components.

Through these examples, the paper demonstrates the superior performance and wide applicability of the new method in different scenarios.

### Summary

The main contribution of this paper is to provide a new Gauss quadrature method that can effectively handle integral measures containing continuous and discrete components. This provides a powerful tool for solving complex practical problems.