Skip to main content\(\require{physics}\require{upgreek}\everymath{\displaystyle}
\newcommand{\N}{\mathbb N}
\newcommand{\Z}{\mathbb Z}
\newcommand{\Q}{\mathbb Q}
\newcommand{\R}{\mathbb R}
\newcommand{\inv}{^{-1}}
\newcommand{\DS}{\displaystyle}
\newcommand{\eps}{\varepsilon}
\newcommand{\tsub}[1]{_{\mathrm{#1}}}
\newcommand{\ee}{\mathrm{e}}
\newcommand{\ii}{\mathrm{i}}
\newcommand{\limit}[1]{\lim\limits_{#1}}
\newcommand{\resid}[1]{\underset{#1}{\Res}}
\DeclareMathOperator{\sinc}{sinc}
\DeclareMathOperator{\sgn}{sgn}
\newcommand{\pii}{\pi}
\DeclareMathOperator{\Prob}{P}
\DeclareMathOperator{\EV}{E}
\DeclareMathOperator{\Var}{Var}
\newcommand{\bv}[1]{\boldsymbol{#1}}
\newcommand{\uv}[1]{\hat{\bv{#1}}}
\newcommand{\cl}[1]{\mathcal{#1}}
\newcommand{\bb}[1]{\mathbb{#1}}
\DeclareMathOperator{\Cis}{cis}
\DeclareMathOperator{\RE}{Re}
\DeclareMathOperator{\IM}{Im}
\newcommand{\xd}{\mathbf{d}}
\newcommand{\seq}[3]{{#1}_{#2},\cdots,{#1}_{#3}}
\newcommand{\psup}[1]{^{(#1)}}
\newcommand{\hypext}{{}^*}
\DeclareMathOperator{\st}{st}
\newcommand{\set}[1]{\left\{#1\right\}}
\DeclareMathOperator{\Sin}{Sin}
\DeclareMathOperator{\Cos}{Cos}
\DeclareMathOperator{\Tan}{Tan}
\DeclareMathOperator{\Sec}{Sec}
\DeclareMathOperator{\Csc}{Csc}
\DeclareMathOperator{\Cot}{Cot}
\DeclareMathOperator{\Log}{Log}
\DeclareMathOperator{\Arg}{Arg}
\DeclareMathOperator{\Ln}{Ln}
\DeclareMathOperator{\Grad}{grad}
\DeclareMathOperator{\Div}{div}
\DeclareMathOperator{\Curl}{curl}
\newcommand{\rd}{\textstyle\mathop{}\!\mathrm{d}^{\!\!\!-}\hspace{-0.0555 em}}
\newcommand{\rpd}{\textstyle\mathop{}\!\partial^{\hspace{-0.5 em}-}\hspace{-0.1666 em}}
\newcommand{\rdv}[2]{\frac{\rd{#1}}{\rd{#2}}}
\newcommand{\rpdv}[2]{\frac{\rpd{#1}}{\rpd{#2}}}
\newcommand{\lt}{<}
\newcommand{\gt}{>}
\newcommand{\amp}{&}
\definecolor{fillinmathshade}{gray}{0.9}
\newcommand{\fillinmath}[1]{\mathchoice{\colorbox{fillinmathshade}{$\displaystyle \phantom{\,#1\,}$}}{\colorbox{fillinmathshade}{$\textstyle \phantom{\,#1\,}$}}{\colorbox{fillinmathshade}{$\scriptstyle \phantom{\,#1\,}$}}{\colorbox{fillinmathshade}{$\scriptscriptstyle\phantom{\,#1\,}$}}}
\)
Section 9.2 Numerical integration
Objectives
-
Approximate integrals numerically when no closed-form antiderivative is available.
-
Compare rectangular, trapezoidal, and parabolic approximations of integrals.
-
Use concavity and derivative bounds to identify when each method over- or underestimates the true value.
-
Frame this as how we can deal with functions with no closed-form antiderivative, like the ones we just saw.
-
Left Rectangular Approximation:
\begin{equation*}
L_{n}=\sum_{i=1}^{n} f(x_{i-1})\,\Delta x_{i}
\end{equation*}
If \(f\) is increasing/decreasing, this will be an under/overestimate.
-
Right Rectangular Approximation:
\begin{equation*}
R_{n}=\sum_{i=1}^{n} f(x_{i})\,\Delta x_{i}
\end{equation*}
If \(f\) is increasing/decreasing, this will be an over/underestimate.
-
Midpoint Rectangular Approximation:
\begin{equation*}
M_{n}=\sum_{i=1}^{n} f(\bar x_{i})\,\Delta x_{i}
\end{equation*}
(Here, \(\bar x_{i}=\dfrac{x_{i-1}+x_{i}}{2}\) is the midpoint of \(x_{i-1}\) and \(x_{i}\text{.}\))
If
\(f\) is concave up/down, this will be an under/overestimate.
Error if
\(|f''(x)|\le K\) for
\(a\le x\le b\text{:}\) \(|E_{M}|\le\dfrac{K(b-a)^{3}}{24n^{2}}\)
-
Trapezoidal Approximation:
\begin{equation*}
T_{n}=\sum_{i=1}^{n} \dfrac{f(x_{i-1})+f(x_{i})}{2}\,\Delta x_{i}=\dfrac{L_{n}+R_{n}}{2}
\end{equation*}
If \(f\) is concave up/down, this will be an over/underestimate.
Error if
\(|f''(x)|\le K\) for
\(a\le x\le b\text{:}\) \(|E_{T}|\le\dfrac{K(b-a)^{3}}{12n^{2}}\)
-
Simpson’s rule:
\begin{equation*}
\int_{a}^{b} f(x)\dd{x}\approx\dfrac{b-a}{6}\left(f(a)+4f\left(\tfrac{a+b}2\right)+f(b)\right)
\end{equation*}
Relation to trapezoid and midpoint rules: \(S_{2n}=\dfrac{T_{n}+2M_{n}}{3}\)
Error if
\(|f^{(4)}(x)|\le K\) for
\(a\le x\le b\text{:}\) \(|E_{T}|\le\dfrac{K(b-a)^{5}}{180n^{4}}\)