Skip to main content
\(\renewcommand{\chaptername}{Unit} \newcommand{\derivativehomeworklink}[1]{\href{http://db.tt/cSeKG8XO}{#1}} \newcommand{\chpname}{unit} \newcommand{\sageurlforcurvature}{http://bmw.byuimath.com/dokuwiki/doku.php?id=curvature_calculator} \newcommand{\uday}{ \LARGE Day \theunitday \normalsize \flushleft \stepcounter{unitday} } \newcommand{\sageworkurl}{http://bmw.byuimath.com/dokuwiki/doku.php?id=work_calculator} \newcommand{\sagefluxurl}{http://bmw.byuimath.com/dokuwiki/doku.php?id=flux_calculator} \newcommand{\sageworkfluxurl}{http://bmw.byuimath.com/dokuwiki/doku.php?id=both_flux_and_work} \newcommand{\sagelineintegral}{http://bmw.byuimath.com/dokuwiki/doku.php?id=line_integral_calculator} \newcommand{\sagephysicalpropertiestwod}{http://bmw.byuimath.com/dokuwiki/doku.php?id=physical_properties_in_2d} \newcommand{\sagephysicalpropertiesthreed}{http://bmw.byuimath.com/dokuwiki/doku.php?id=physical_properties_in_3d} \newcommand{\sageDoubleIntegralCheckerURL}{http://bmw.byuimath.com/dokuwiki/doku.php?id=double_integral_calculator} \newcommand{\myscale}{1} \newcommand{\ds}{\displaystyle} \newcommand{\dfdx}[1]{\frac{d#1}{dx}} \newcommand{\ddx}{\frac{d}{dx}} \newcommand{\ii}{\vec \imath} \newcommand{\jj}{\vec \jmath} \newcommand{\kk}{\vec k} \newcommand{\vv}{\mathbf{v}} \newcommand{\RR}{\mathbb{R}} \newcommand{\R}{ \mathbb{R}} \newcommand{\inv}{^{-1}} \newcommand{\im}{\text{im }} \newcommand{\colvec}[1]{\begin{bmatrix}#1\end{bmatrix} } \newcommand{\cl}[1]{ \begin{matrix} #1 \end{matrix} } \newcommand{\bm}[1]{ \begin{bmatrix} #1 \end{bmatrix} } \DeclareMathOperator{\rank}{rank} \DeclareMathOperator{\rref}{rref} \DeclareMathOperator{\vspan}{span} \DeclareMathOperator{\trace}{tr} \DeclareMathOperator{\proj}{proj} \DeclareMathOperator{\curl}{curl} \newcommand{\blank}[1]{[14pt]{\rule{#1}{1pt}}} \newcommand{\vp}{^{\,\prime}} \newcommand{\lt}{<} \newcommand{\gt}{>} \newcommand{\amp}{&} \)

Section1.4Optional: Higher Order Approximations

When you ask a calculator to tell you what \(e^{.1}\) means, your calculator uses an extension of differentials to give you an approximation. The calculator only uses polynomials (multiplication and addition) to give you an answer. This same process is used to evaluate any function that is not a polynomial (so trig functions, square roots, inverse trig functions, logarithms, etc.) The key idea needed to approximate functions is illustrated by the next problem.

Exercise1.4.1

Let \(f(x)=e^x\text{.}\) You should find that your work on each step can be reused to do the next step.

(a)

Find a first degree polynomial \(P_1(x)=a+bx\) so that \(P_1(0)=f(0)\) and \(P'_1(0)=f'(0)\text{.}\) In other words, give me a line that passes through the same point and has the same slope as \(f(x)=e^x\) does at \(x=0\text{.}\) Set up a system of equations and then find the unknowns \(a\) and \(b\text{.}\) The next two are very similar.

(b)

Find a second degree polynomial \(P_2(x)=a+bx+cx^2\) so that \(P_2(0)=f(0)\text{,}\) \(P'_2(0)=f'(0)\text{,}\) and \(P''_2(0)=f''(0)\text{.}\) In other words, give me a parabola that passes through the same point, has the same slope, and has the same concavity as \(f(x)=e^x\) does at \(x=0\text{.}\)

(c)

Find a third degree polynomial \(P_3(x)=a+bx+cx^2+dx^3\) so that \(P_3(0)=f(0)\text{,}\) \(P'_3(0)=f'(0)\text{,}\) \(P''_3(0)=f''(0)\text{,}\) and \(P'''_3(0)=f'''(0)\text{.}\) In other words, give me a cubic that passes through the same point, has the same slope, the same concavity, and the same third derivative as \(f(x)=e^x\) does at \(x=0\text{.}\)

(d)

Now compute \(e^{.1}\) with a calculator. Then compute \(P_1(.1)\text{,}\) \(P_2(.1)\text{,}\) and \(P_3(.1)\text{.}\) How accurate are the line, parabola, and cubic in approximating \(e^{.1}\text{?}\)

Exercise1.4.2

Now let \(f(x)=\sin x\text{.}\) Find a 7th degree polynomial so that the function and the polynomial have the same value and same first seven derivatives when evaluated at \(x=0\text{.}\) Evaluate the polynomial at \(x=0.3\text{.}\) How close is this value to your calculator's estimate of \(\sin(0.3)\text{?}\) You may find it valuable to use the notation

\begin{equation*} P(x) = a_0+a_1x+a_2x^2+a_3x^3 +\cdots+a_7 x^7. \end{equation*}

The polynomial you are creating is often called a Taylor polynomial. (I'm giving you the name so that you can search online for more information if you are interested.)

The previous two problems involved finding polynomial approximations to the function at \(x=0\text{.}\) The next problem shows how to move this to any other point, such as \(x=1\text{.}\)

Exercise1.4.3

Let \(f(x)=e^x\text{.}\)

(a)

Find a second degree polynomial

\begin{equation*} T(x)=a+bx+cx^2 \end{equation*}

so that \(T(1)=f(1)\text{,}\) \(T'(1)=f'(1)\text{,}\) and \(T''(1)=f''(1)\text{.}\) In other words, give me a parabola that passes through the same point, has the same slope, and the same concavity as \(f(x)=e^x\) does at \(x=1\text{.}\)

(b)

Notice that we just replaced \(x\) with \(x-1\text{.}\) This centers, or shifts, the approximation to be at \(x=1\text{.}\) The first part will be much simpler now when you let \(x=1\text{.}\)

Find a second degree polynomial written in the form

\begin{equation*} S(x)=a+b(x-1)+c(x-1)^2 \end{equation*}

so that \(S(1)=f(1)\text{,}\) \(S'(1)=f'(1)\text{,}\) and \(S''(1)=f''(1)\text{.}\) In other words, find a quadratric that passes through the same point, has the same slope, and the same concavity as \(f(x)=e^x\) does at \(x=1\text{.}\)

(c)

Find a third degree polynomial written in the form

\begin{equation*} P(x)=a+b(x-1)+c(x-1)^2+d(x-1)^3 \end{equation*}

so that \(P(1)=f(1)\text{,}\) \(P'(1)=f'(1)\text{,}\) \(P''(1)=f''(1)\text{,}\) and \(P'''(1)=f'''(1)\text{.}\) In other words, give me a cubic that passes through the same point, has the same slope, the same concavity, and the same third derivative as \(f(x)=e^x\) does at \(x=1\text{.}\)

This example refers back to Exercise 1.2.3. We wanted a spherical tank of radius 3m, but due to manufacturing error the radius was slightly off. Let's now illustrate how we can use polynomials to give a first, second, and third order approximation of the volume if the radius is 3.02m instead of 3m.

We start with \(V=\frac{4}{3} \pi r^3\) and then compute the derivatives

\begin{equation*} V'=4\pi r^2, V''=8\pi r, \text{ and } V'''=8\pi. \end{equation*}

Because we are approximating the increase in volume from \(r=3\) to something new, we'll create our polynomial approximations centered at \(r=3\text{.}\) We'll consider the polynomial

\begin{equation*} P(r)=a_0+a_1(r-3)+a_2(r-3)^2+a_3(r-3)^3, \end{equation*}

whose derivatives are

\begin{equation*} P'=a_1+2a_2(r-3)+3a_3(r-3)^2, P''=2a_2+6a_3(r-3), P'''=6a_3. \end{equation*}

So that the derivatives of the volume function match the derivatives of the polynomial (at \(r=3\)), we need to satisfy the equations in the table below.

\(k\) Value of \(V\) at the \(k\)th derivative Value of \(P\) at the the \(k\)th derivative Equation
\(0\) \(V(3) = \frac{4}{3}\pi (3)^3 = 36\pi\) \(P(3) = a_0\) \(a_0=36\pi\)
\(1\) \(V'(3) = 4\pi (3)^2=36\pi\) \(P'(3) = a_1\) \(a_1=36\pi\)
\(2\) \(V''(3) = 8\pi (3)=24\pi\) \(P''(3) = 2a_2\) \(2a_2=24\pi\)
\(3\) \(V'''(3) = 8\pi\) \(P'''(3) = 6a_3\) \(6a_3=8\pi\)
Table1.4.2

This tells us that the third order polynomial is

\begin{equation*} P(r)=a_0+a_1(r-3)+a_2(r-3)^2+a_3(r-3)^3 =36\pi+36\pi(r-3)+12\pi(r-3)^2+\frac{4}{3}\pi(r-3)^3 . \end{equation*}

We wanted to approximate the volume if \(r=3.2\text{,}\) so our change in \(r\) is \(dr=3.2-3=0.2\text{.}\) We can rewrite our polynomial as

\begin{equation*} P(r)=36\pi+36\pi(dr)+12\pi(dr)^2+\frac{4}{3}\pi(dr)^3. \end{equation*}

We are now prepared to approximate the volume using a first, second, and third order approximation.

  1. A first order approximation yields \(P=36\pi+ 36\pi\cdot 0.02 =36.72\pi.\) The volume increased by \(0.72\pi\) m\(^3\text{.}\)

  2. A second order approximation yields

    \begin{equation*} P=36\pi+ 36\pi\cdot 0.02 +12\pi (0.02)^2 =36.7248\pi. \end{equation*}
  3. A third order approximation yields

    \begin{equation*} P=36\pi+36\pi\cdot 0.02 +12\pi (0.02)^2+\frac{4}{3}\pi(0.02)^3 =36.724810\bar6\pi. \end{equation*}

With each approximation, we add on a little more volume to get closer to the actual volume of a sphere with radius \(r=3.02\text{.}\) The actual volume of a sphere involves a cubic function, so when we approximate the volume with a cubic, we should get an exact approximation (and \(V(3.02) = \frac 43 \pi (3.02)^3 =(36.724810\bar6)\pi\text{.}\))

We'll end this section with a problem to practice the example above.

Exercise1.4.4

As a challenge, try to draw a 3D graph which illustrates the volume added on by each successive approximation. If you have it before we get here as a class, let me know and I'll let you share what you have discovered with the class.

Suppose you are constructing a cube whose side length should be \(s=2\) units. The manufacturing process is not exact, but instead creates a cube with side lengths \(s=2+ds\) units. (You should assume that all sides are still the same, so any error on one side is replicated on all. We have to assume this for now, but before the semester ends we'll be able to do this with high dimensional calculus.)

Suppose that the machine creates a cube with side length \(2.3\) units instead of 2 units. Note that the volume of the cube is \(V=s^3\text{.}\) Use a first, second, and third order approximation to estimate the increase in volume caused by the .3 increase in side length. Then compute the actual increase in volume \(V(2.3)-V(2)\text{.}\)