Boole's rule

The widely propagated typographical error Bode's rule redirects here. For Bode's Law, see Titius–Bode law.

In mathematics, Boole's rule, named after George Boole, is a method of numerical integration.

Formula

Simple Boole's Rule

It approximates an integral: $${\displaystyle \int _{a}^{b}f(x)\,dx}$$ by using the values of f at five equally spaced points: ^[1] $${\displaystyle {\begin{aligned}&x_{0}=a\\&x_{1}=x_{0}+h\\&x_{2}=x_{0}+2h\\&x_{3}=x_{0}+3h\\&x_{4}=x_{0}+4h=b\end{aligned}}}$$

It is expressed thus in Abramowitz and Stegun: ^[2] $${\displaystyle \int _{x_{0}}^{x_{4}}f(x)\,dx={\frac {2h}{45}}{\bigl [}7f(x_{0})+32f(x_{1})+12f(x_{2})+32f(x_{3})+7f(x_{4}){\bigr ]}+{\text{error term}}}$$ where the error term is $${\displaystyle -\,{\frac {8f^{(6)}(\xi )h^{7}}{945}}}$$ for some number ⁠${\displaystyle \xi }$⁠ between ⁠${\displaystyle x_{0}}$⁠ and ⁠${\displaystyle x_{4}}$⁠ where 945 = 1 × 3 × 5 × 7 × 9.

It is often known as Bode's rule, due to a typographical error that propagated from Abramowitz and Stegun. ^[3]

The following constitutes a very simple implementation of the method in Common Lisp which ignores the error term:

Example implementation in Common Lisp

(defunintegrate-booles-rule(fx1x5)"Calculates the Boole's rule numerical integral of the function F in the closed interval extending from inclusive X1 to inclusive X5 without error term inclusion."(declare(type(function(real)real)f))(declare(typerealx1x5))(let((h(/(-x5x1)4)))(declare(typerealh))(let*((x2(+x1h))(x3(+x2h))(x4(+x3h)))(declare(typerealx2x3x4))(*(/(*2h)45)(+(*7(funcallfx1))(*32(funcallfx2))(*12(funcallfx3))(*32(funcallfx4))(*7(funcallfx5)))))))

Composite Boole's Rule

In cases where the integration is permitted to extend over equidistant sections of the interval ${\displaystyle [a,b]}$, the composite Boole's rule might be applied. Given ${\displaystyle N}$ divisions, where ${\displaystyle N}$ mod ${\displaystyle 4=0}$, the integrated value amounts to: ^[4]

$${\displaystyle \int _{x_{0}}^{x_{N}}f(x)\,dx={\frac {2h}{45}}\left(7(f(x_{0})+f(x_{N}))+32\left(\sum _{i\in \{1,3,5,\ldots ,N-1\}}f(x_{i})\right)+12\left(\sum _{i\in \{2,6,10,\ldots ,N-2\}}f(x_{i})\right)+14\left(\sum _{i\in \{4,8,12,\ldots ,N-4\}}f(x_{i})\right)\right)+{\text{error term}}}$$

where the error term is similar to above. The following Common Lisp code implements the aforementioned formula:

Example implementation in Common Lisp

(defunintegrate-composite-booles-rule(fabn)"Calculates the composite Boole's rule numerical integral of the function F in the closed interval extending from inclusive A to inclusive B across N subintervals."(declare(type(function(real)real)f))(declare(typerealab))(declare(type(integer1*)n))(let((h(/(-ba)n)))(declare(typerealh))(flet((f[i](i)(declare(type(integer0*)i))(let((xi(+a(*ih))))(declare(typerealxi))(thereal(funcallfxi)))))(*(/(*2h)45)(+(*7(+(f[i]0)(f[i]n)))(*32(loopforifrom1to(-n1)by2sum(f[i]i)))(*12(loopforifrom2to(-n2)by4sum(f[i]i)))(*14(loopforifrom4to(-n4)by4sum(f[i]i))))))))

Example implementation in R

booleQuad<-function(fx,dx){# Calculates the composite Boole's rule numerical # integral for a function with a vector of precomputed# values fx evaluated at the points in vector dx.n<-length(dx)h<-diff(dx)stopifnot(exprs={length(fx)==nn>8Lh[1L]>=0n>=2Ln%%4L==1LisTRUE(all.equal(h,rep(h[1L],length(h))))})nm2<-n-2Lcf<-double(nm2)cf[seq.int(1,nm2,2L)]<-32cf[seq.int(2,nm2,4L)]<-12cf[seq.int(4,nm2,4L)]<-14cf<-c(7,cf,7)sum(cf*fx)*2*h[1L]/45}

See also

• Newton–Cotes formulas
• Simpson's rule
• Romberg's method

Notes

1. Boole 1880, p. 47, Eq(21).
2. Davis & Polonsky 1983.
3. Weisstein.
4. Sablonnière, Sbibih & Tahrichi 2010, p. 852.

References

• Boole, George (1880) [1860]. A Treatise on the Calculus of Finite Differences (3rd ed.). Macmillan and Company.
• Davis, Philip J.; Polonsky, Ivan (1983) [June 1964]. "Chapter 25, eqn 25.4.14". In Abramowitz, Milton; Stegun, Irene Ann (eds.). Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series. Vol. 55 (Ninth reprint with additional corrections of tenth original printing with corrections (December 1972); first ed.). Washington D.C.; New York: United States Department of Commerce, National Bureau of Standards; Dover Publications. p. 886. ISBN   978-0-486-61272-0. LCCN   64-60036. MR   0167642. LCCN   65-12253.
• Sablonnière, P.; Sbibih, D.; Tahrichi, M. (2010). "Error estimate and extrapolation of a quadrature formula derived from a quartic spline quasi-interpolant". BIT Numerical Mathematics. 50 (4): 843–862. doi:10.1007/s10543-010-0278-0.
• Weisstein, Eric W. "Boole's Rule". MathWorld .