Nyström method
In numerical analysis, the Nyström method[1] or quadrature method seeks the numerical solution of an integral equation by replacing the integral with a representative weighted sum. The continuous problem is broken into discrete intervals, quadrature or numerical integration determines the weights and locations of representative points for the integral. The discrete problem to be solved is now a system of linear equations with n equations and n unknowns.
From the n solved points the function value at other points is interpolated consistent with the chosen quadrature method. Depending on the original problem and the choice of quadrature the problem may be ill-conditioned.
Since the linear equations require operations to solve, hence high-order quadrature rules perform better because low-order quadrature rules require large for a given accuracy. Gaussian quadrature is normally a good choice for smooth, non-singular problems.
Discretization of the integral
- where are the weights of the quadrature rule, and points are the abscissas.
Example
Applying this to the inhomogeneous Fredholm equation of the second kind
- ,
results in
- .
References
- ↑ Nyström, Evert Johannes (1930). "Über die praktische Auflösung von Integralgleichungen mit Anwendungen auf Randwertaufgaben". Acta Mathematica 54 (1): 185–204. doi:10.1007/BF02547521.
- Leonard M. Delves & Joan Walsh (eds): Numerical Solution of Integral Equations, Clarendon, Oxford, 1974.
- Hans-Jürgen Reinhardt: Analysis of Approximation Methods for Differential and Integral Equations, Springer, New York, 1985.