Simpson's Rule
Simpson's rule
[′sim·sənz ‚rül]Simpson’s Rule
a formula for approximating definite integrals. In has the form
where h = (b - a)/2n and fi = f(a + ih), i = 0, 1, 2,…, 2n.
The derivation of Simpson’s rule is based on the replacement of the integrand f(x) on each of the closed intervals [a + 2hk, a + 2h(k + 1)], k = 0, 1,…, n - 1, by a corresponding interpolation polynomial of the second degree (seeINTERPOLATION FORMULAS). In geometric terms, the curve described by the equation y = f(x) is replaced by an approximating curve consisting of segments of parabolas. The error resulting from the use of Simpson’s rule is
where a ≤ ξ ≤ b. If f(x) is a polynomial of degree m ≤ 3, then Simpson’s rule is not approximate but exact, since in this case f(4)(x) = 0.
Simpson’s rule is named after T. Simpson, who devised it in 1743. The rule was, however, known earlier. For example, it was given by J. Gregory in 1668.
Other formulas for the approximate calculation of definite integrals are discussed in the article APPROXIMATE INTEGRATION.