# Simpson’s Rule

## Solved Problems

Click or tap a problem to see the solution.

### Example 7

Approximate the integral $\int\limits_0^2 {{x^3}dx}$ using Simpson's Rule with n = 4 subintervals.

### Example 8

Using Simpson's Rule with n = 4 subintervals, evaluate the integral $\int\limits_0^1 {{e^x}dx}.$ Round the answer to 3 decimal places.

### Example 9

Using Simpson's Rule with $$n = 2$$ subintervals, evaluate the integral $\int\limits_1^3 {\ln xdx}.$ Round the answer to $$3$$ decimal places.

### Example 10

Using Simpson's Rule with $$n = 4$$ subintervals approximate the area under the sine curve $f\left( x \right) = \sin x$ between $$x = 0$$ and $$x = \pi.$$

### Example 11

Using Simpson's Rule with $$n = 4$$ subintervals approximate the area under the inverse cosine curve $f\left( x \right) = \arccos x$ between $$x = -1$$ and $$x = 1.$$

### Example 7.

Approximate the integral $\int\limits_0^2 {{x^3}dx}$ using Simpson's Rule with $$n = 4$$ subintervals.

Solution.

The Simpson's Rule formula for $$n = 4$$ subintervals has the form

${S_4} = \frac{{\Delta x}}{3}\left[ {f\left( {{x_0}} \right) + 4f\left( {{x_1}} \right) + 2f\left( {{x_2}} \right) + 4f\left( {{x_3}} \right) + 2\left( {{x_4}} \right)} \right].$

Determine the width of the subinterval:

$\Delta x = \frac{{b - a}}{n} = \frac{{2 - 0}}{4} = \frac{1}{2}.$

The values of the cubic function at the points $${x_i} = \left\{ {0,\frac{1}{2},1,\frac{3}{2},2} \right\}$$ are

$f\left( {{x_0}} \right) = f\left( 0 \right) = {0^3} = 0;$
$f\left( {{x_1}} \right) = f\left( {\frac{1}{2}} \right) = {\left( {\frac{1}{2}} \right)^3} = \frac{1}{8};$
$f\left( {{x_2}} \right) = f\left( 1 \right) = {1^3} = 1;$
$f\left( {{x_3}} \right) = f\left( {\frac{3}{2}} \right) = {\left( {\frac{3}{2}} \right)^3} = \frac{{27}}{8};$
$f\left( {{x_4}} \right) = f\left( 2 \right) = {2^3} = 8.$

Hence

$\int\limits_0^2 {{x^3}dx} \approx {S_4} = \frac{1}{6}\left[ {0 + 4 \cdot \frac{1}{8} + 2 \cdot 1 + 4 \cdot \frac{{27}}{8} + 8} \right] = \frac{1}{6}\left[ {\frac{1}{2} + 2 + \frac{{27}}{2} + 8} \right] = \frac{{24}}{6} = 4.$

### Example 8.

Using Simpson's Rule with $$n = 4$$ subintervals, evaluate the integral $\int\limits_0^1 {{e^x}dx}.$ Round the answer to $$3$$ decimal places.

Solution.

We evaluate the given integral by the formula

${S_4} = \frac{{\Delta x}}{3}\left[ {f\left( {{x_0}} \right) + 4f\left( {{x_1}} \right) + 2f\left( {{x_2}} \right) + 4f\left( {{x_3}} \right) + 2\left( {{x_4}} \right)} \right].$

Determine the width of the subinterval:

$\Delta x = \frac{{b - a}}{n} = \frac{{1 - 0}}{{4}} = \frac{1}{4}.$

Compute the function values at the endpoints of the subintervals:

$f\left( {{x_0}} \right) = f\left( 0 \right) = {e^0} = 1;$
$f\left( {{x_1}} \right) = f\left( {\frac{1}{4}} \right) = {e^{\frac{1}{4}}} = \sqrt[4]{e} \approx 1.2840;$
$f\left( {{x_2}} \right) = f\left( {\frac{1}{2}} \right) = {e^{\frac{1}{2}}} = \sqrt e \approx 1.6487;$
$f\left( {{x_3}} \right) = f\left( {\frac{3}{4}} \right) = {e^{\frac{3}{4}}} = \sqrt[4]{{{e^3}}} \approx 2.1170;$
$f\left( {{x_4}} \right) = f\left( 1 \right) = {e^1} = e \approx 2.7183;$

Plugging in the function values into our equation, we get:

$\int\limits_0^1 {{e^x}dx} \approx {S_4} = \frac{1}{{12}}\left[ {1 + 4 \times 1.2840 + 2 \times 1.6487 + 4 \times 2.1170 + 2.7183} \right] = \frac{1}{{12}} \times 20.6197 = 1.7183 \approx 1.718$

### Example 9.

Using Simpson's Rule with $$n = 2$$ subintervals, evaluate the integral $\int\limits_1^3 {\ln xdx}.$ Round the answer to $$3$$ decimal places.

Solution.

We evaluate the integral by the following approximate formula:

${S_2} = \frac{{\Delta x}}{3}\left[ {f\left( {{x_0}} \right) + 4f\left( {{x_1}} \right) + f\left( {{x_2}} \right)} \right].$

The width of the subinterval is

$\Delta x = \frac{{b - a}}{n} = \frac{{2 - 0}}{{2}} = 1.$

Calculate the function values at the points $${x_i} = \left\{ {1,2,3} \right\}:$$

$f\left( {{x_0}} \right) = f\left( 1 \right) = \ln 1 = 0;$
$f\left( {{x_1}} \right) = f\left( 2 \right) = \ln 2 \approx 0.6931;$
$f\left( {{x_2}} \right) = f\left( 3 \right) = \ln 3 \approx 1.0986;$

Hence

$\int\limits_1^3 {\ln xdx} \approx {S_2} = \frac{1}{3}\left[ {0 + 4\ln 2 + \ln 3} \right] \approx \frac{1}{3}\left[ {4 \times 0.6931 + 1.0986} \right] = \frac{1}{3} \times 3.8710 = 1.290$

### Example 10.

Using Simpson's Rule with $$n = 4$$ subintervals approximate the area under the sine curve $f\left( x \right) = \sin x$ between $$x = 0$$ and $$x = \pi.$$

Solution.

The Simpson's Rule formula with $$n = 4$$ segments is written in the form

${S_4} = \frac{{\Delta x}}{3}\left[ {f\left( {{x_0}} \right) + 4f\left( {{x_1}} \right) + 2f\left( {{x_2}} \right) + 4f\left( {{x_3}} \right) + \left( {{x_4}} \right)} \right].$

The width of the subinterval is

$\Delta x = \frac{{b - a}}{n} = \frac{{\pi - 0}}{4} = \frac{\pi }{4}.$

Calculate the values of the sine function at the points $${x_i} = \left\{ {0,\frac{\pi }{4}, \frac{\pi }{2}, \frac{{3\pi }}{4}, \pi } \right\}:$$

$f\left( {{x_0}} \right) = f\left( 0 \right) = \sin 0 = 0;$
$f\left( {{x_1}} \right) = f\left( {\frac{\pi }{4}} \right) = \sin \frac{\pi }{4} = \frac{{\sqrt 2 }}{2};$
$f\left( {{x_2}} \right) = f\left( {\frac{\pi }{2}} \right) = \sin \frac{\pi }{2} = 1;$
$f\left( {{x_3}} \right) = f\left( {\frac{{3\pi }}{4}} \right) = \sin \frac{{3\pi }}{4} = \frac{{\sqrt 2 }}{2};$
$f\left( {{x_4}} \right) = f\left( \pi \right) = \sin \pi = 0.$

Now we can substitute these values into our equation and find the approximate value of the area:

$A = \int\limits_0^\pi {\sin xdx} \approx {S_4} = \frac{\pi }{{12}}\left[ {0 + 4 \cdot \frac{{\sqrt 2 }}{2} + 2 \cdot 1 + 4 \cdot \frac{{\sqrt 2 }}{2} + 0} \right] = \frac{\pi }{{12}}\left( {4\sqrt 2 + 2} \right) = \frac{\pi }{6}\left( {2\sqrt 2 + 1} \right).$

### Example 11.

Using Simpson's Rule with $$n = 4$$ subintervals approximate the area under the inverse cosine curve $f\left( x \right) = \arccos x$ between $$x = -1$$ and $$x = 1.$$

Solution.

The Simpson's Rule formula with $$n = 4$$ segments has the form

${S_4} = \frac{{\Delta x}}{3}\left[ {f\left( {{x_0}} \right) + 4f\left( {{x_1}} \right) + 2f\left( {{x_2}} \right) + 4f\left( {{x_3}} \right) + \left( {{x_4}} \right)} \right].$

Determine the width of the subinterval:

$\Delta x = \frac{{b - a}}{n} = \frac{{1 - \left({-1}\right)}}{4} = \frac{1}{2}.$

Calculate the values of the inverse cosine at the endpoints of the subintervals:

$f\left( {{x_0}} \right) = f\left( { - 1} \right) = \arccos \left( { - 1} \right) = \pi ;$
$f\left( {{x_1}} \right) = f\left( { - \frac{1}{2}} \right) = \arccos \left( { - \frac{1}{2}} \right) = \frac{{2\pi }}{3};$
$f\left( {{x_2}} \right) = f\left( 0 \right) = \arccos 0 = \frac{\pi }{2};$
$f\left( {{x_3}} \right) = f\left( {\frac{1}{2}} \right) = \arccos \frac{1}{2} = \frac{\pi }{3};$
$f\left( {{x_4}} \right) = f\left( 1 \right) = \arccos 1 = 0.$

Substituting these values into our equation, we find the area under the curve:

$A = \int\limits_{ - 1}^1 {\arccos xdx} \approx {S_4} = \frac{1}{6}\left[ {\pi + 4 \cdot \frac{{2\pi }}{3} + 2 \cdot \frac{\pi }{2} + 4 \cdot \frac{\pi }{3} + 0} \right] = \frac{1}{6} \cdot 6\pi = \pi$