The ApproximateInt(f(x), x = a..b, method = simpson, opts) command approximates the integral of f(x) from a to b by using Simpson's rule. The first two arguments (function expression and range) can be replaced by a definite integral.

If the independent variable can be uniquely determined from the expression, the parameter x need not be included in the calling sequence.

Given a partition $P=\left(a=x_0,x_1,...,x_N=b\right)$ of the interval $\left(a\,b\right)$, Simpson's rule approximates the integral on each subinterval $\left(x_{i-1}\,x_i\right)$ by integrating the quadratic function that interpolates the three points $\left(x_{i-1}\,f\left(x_{i-1}\right)\right)$, $\left(\frac{x_{i-1}}{2}+\frac{x_i}{2}\,f\left(\frac{x_{i-1}}{2}+\frac{x_i}{2}\right)\right)$, and $\left(x_i\,f\left(x_i\right)\right)$.  This value is

In the case that the widths of the subintervals are equal, the approximation can be written as

Traditionally, Simpson's rule is written as: given N where N is an even integer and given equally spaced points $a=x_0,x_1,x_2,...,x_N=b$, an approximation to the integral ${\int }_a^{b}f\left(x\right)\ⅆx$ is

By default, the interval is divided into $10$ equal-sized subintervals.

For the options opts, see the ApproximateInt help page.

This rule can be applied interactively, through the ApproximateInt Tutor.

$\mathrm{polynomial}≔\mathrm{CurveFitting}\left[\mathrm{PolynomialInterpolation}\right]\left(\left[x\left[0\right]\,\frac{x\left[0\right]+x\left[1\right]}{2}\,x\left[1\right]\right]\,\left[f\left(0\right)\,f\left(\frac{1}{2}\right)\,f\left(1\right)\right]\,z\right)\:$

$\mathrm{integrated}≔\mathrm{int}\left(\mathrm{polynomial}\,z=x\left[0\right]..x\left[1\right]\right)\:$

$\mathrm{factor}\left(\mathrm{integrated}\right)$

$-\frac{\left(x_0-x_1\right)\left(f\left(0\right)+f\left(1\right)+4f\left(\frac{1}{2}\right)\right)}{6}$

$\mathrm{with}\left(\mathrm{Student}\left[\mathrm{Calculus1}\right]\right)\:$

$\mathrm{ApproximateInt}\left(\sin \left(x\right)\,x=0..5\,\mathrm{method}=\mathrm{simpson}\right)$

$\frac{\sin \left(\frac{19}{4}\right)}{3}+\frac{\sin \left(5\right)}{12}+\frac{\sin \left(3\right)}{6}+\frac{\sin \left(\frac{13}{4}\right)}{3}+\frac{\sin \left(\frac{7}{2}\right)}{6}+\frac{\sin \left(\frac{15}{4}\right)}{3}+\frac{\sin \left(4\right)}{6}+\frac{\sin \left(\frac{17}{4}\right)}{3}+\frac{\sin \left(\frac{9}{2}\right)}{6}+\frac{\sin \left(\frac{5}{4}\right)}{3}+\frac{\sin \left(\frac{3}{2}\right)}{6}+\frac{\sin \left(\frac{7}{4}\right)}{3}+\frac{\sin \left(2\right)}{6}+\frac{\sin \left(\frac{9}{4}\right)}{3}+\frac{\sin \left(\frac{5}{2}\right)}{6}+\frac{\sin \left(\frac{11}{4}\right)}{3}+\frac{\sin \left(\frac{1}{4}\right)}{3}+\frac{\sin \left(\frac{1}{2}\right)}{6}+\frac{\sin \left(\frac{3}{4}\right)}{3}+\frac{\sin \left(1\right)}{6}$

$\mathrm{ApproximateInt}\left(x\left(x-2\right)\left(x-3\right)\,x=0..5\,\mathrm{method}=\mathrm{simpson}\,\mathrm{output}=\mathrm{plot}\right)$

$\mathrm{ApproximateInt}\left(\tan \left(x\right)-2x\,x=-1..1\,\mathrm{method}=\mathrm{simpson}\,\mathrm{output}=\mathrm{plot}\,\mathrm{partition}=50\right)$

$\mathrm{ApproximateInt}\left(\ln \left(x\right)\,1..100\,\mathrm{method}=\mathrm{simpson}\,\mathrm{output}=\mathrm{animation}\right)$