3
<title>Numerical Integration Introduction - GNU Scientific Library -- Reference Manual</title>
4
<meta http-equiv="Content-Type" content="text/html">
5
<meta name="description" content="GNU Scientific Library -- Reference Manual">
6
<meta name="generator" content="makeinfo 4.8">
7
<link title="Top" rel="start" href="index.html#Top">
8
<link rel="up" href="Numerical-Integration.html#Numerical-Integration" title="Numerical Integration">
9
<link rel="next" href="QNG-non_002dadaptive-Gauss_002dKronrod-integration.html#QNG-non_002dadaptive-Gauss_002dKronrod-integration" title="QNG non-adaptive Gauss-Kronrod integration">
10
<link href="http://www.gnu.org/software/texinfo/" rel="generator-home" title="Texinfo Homepage">
12
Copyright (C) 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006 The GSL Team.
14
Permission is granted to copy, distribute and/or modify this document
15
under the terms of the GNU Free Documentation License, Version 1.2 or
16
any later version published by the Free Software Foundation; with the
17
Invariant Sections being ``GNU General Public License'' and ``Free Software
18
Needs Free Documentation'', the Front-Cover text being ``A GNU Manual'',
19
and with the Back-Cover Text being (a) (see below). A copy of the
20
license is included in the section entitled ``GNU Free Documentation
23
(a) The Back-Cover Text is: ``You have freedom to copy and modify this
24
GNU Manual, like GNU software.''-->
25
<meta http-equiv="Content-Style-Type" content="text/css">
26
<style type="text/css"><!--
27
pre.display { font-family:inherit }
28
pre.format { font-family:inherit }
29
pre.smalldisplay { font-family:inherit; font-size:smaller }
30
pre.smallformat { font-family:inherit; font-size:smaller }
31
pre.smallexample { font-size:smaller }
32
pre.smalllisp { font-size:smaller }
33
span.sc { font-variant:small-caps }
34
span.roman { font-family:serif; font-weight:normal; }
35
span.sansserif { font-family:sans-serif; font-weight:normal; }
41
<a name="Numerical-Integration-Introduction"></a>
42
Next: <a rel="next" accesskey="n" href="QNG-non_002dadaptive-Gauss_002dKronrod-integration.html#QNG-non_002dadaptive-Gauss_002dKronrod-integration">QNG non-adaptive Gauss-Kronrod integration</a>,
43
Up: <a rel="up" accesskey="u" href="Numerical-Integration.html#Numerical-Integration">Numerical Integration</a>
47
<h3 class="section">16.1 Introduction</h3>
49
<p>Each algorithm computes an approximation to a definite integral of the
51
where w(x) is a weight function (for general integrands w(x)=1).
52
The user provides absolute and relative error bounds
53
<!-- {$(\hbox{\it epsabs}, \hbox{\it epsrel}\,)$} -->
54
(epsabs, epsrel) which specify the following accuracy requirement,
56
<!-- {$\hbox{\it RESULT}$} -->
57
RESULT is the numerical approximation obtained by the
58
algorithm. The algorithms attempt to estimate the absolute error
59
<!-- {$\hbox{\it ABSERR} = |\hbox{\it RESULT} - I|$} -->
60
ABSERR = |RESULT - I| in such a way that the following inequality
62
The routines will fail to converge if the error bounds are too
63
stringent, but always return the best approximation obtained up to that
66
<p>The algorithms in <span class="sc">quadpack</span> use a naming convention based on the
69
<pre class="display"> <code>Q</code> - quadrature routine
71
<code>N</code> - non-adaptive integrator
72
<code>A</code> - adaptive integrator
74
<code>G</code> - general integrand (user-defined)
75
<code>W</code> - weight function with integrand
77
<code>S</code> - singularities can be more readily integrated
78
<code>P</code> - points of special difficulty can be supplied
79
<code>I</code> - infinite range of integration
80
<code>O</code> - oscillatory weight function, cos or sin
81
<code>F</code> - Fourier integral
82
<code>C</code> - Cauchy principal value
84
<p class="noindent">The algorithms are built on pairs of quadrature rules, a higher order
85
rule and a lower order rule. The higher order rule is used to compute
86
the best approximation to an integral over a small range. The
87
difference between the results of the higher order rule and the lower
88
order rule gives an estimate of the error in the approximation.
91
<li><a accesskey="1" href="Integrands-without-weight-functions.html#Integrands-without-weight-functions">Integrands without weight functions</a>
92
<li><a accesskey="2" href="Integrands-with-weight-functions.html#Integrands-with-weight-functions">Integrands with weight functions</a>
93
<li><a accesskey="3" href="Integrands-with-singular-weight-functions.html#Integrands-with-singular-weight-functions">Integrands with singular weight functions</a>