Newton’s Method

Learning Objectives

  • Describe the steps of Newton’s method.
  • Explain what an iterative process means.
  • Recognize when Newton’s method does not work.
  • Apply iterative processes to various situations.

In many areas of pure and applied mathematics, we are interested in finding solutions to an equation of the form [latex]f\left(x\right)=0.[/latex] For most functions, however, it is difficult—if not impossible—to calculate their zeroes explicitly. In this section, we take a look at a technique that provides a very efficient way of approximating the zeroes of functions. This technique makes use of tangent line approximations and is behind the method used often by calculators and computers to find zeroes.

Describing Newton’s Method

Consider the task of finding the solutions of [latex]f\left(x\right)=0.[/latex] If [latex]f[/latex] is the first-degree polynomial [latex]f\left(x\right)=ax+b,[/latex] then the solution of [latex]f\left(x\right)=0[/latex] is given by the formula [latex]x=-\frac{b}{a}.[/latex] If [latex]f[/latex] is the second-degree polynomial [latex]f\left(x\right)=a{x}^{2}+bx+c,[/latex] the solutions of [latex]f\left(x\right)=0[/latex] can be found by using the quadratic formula. However, for polynomials of degree [latex]3[/latex] or more, finding roots of [latex]f[/latex] becomes more complicated. Although formulas exist for third- and fourth-degree polynomials, they are quite complicated. Also, if [latex]f[/latex] is a polynomial of degree [latex]5[/latex] or greater, it is known that no such formulas exist. For example, consider the function

[latex]f\left(x\right)={x}^{5}+8{x}^{4}+4{x}^{3}-2x-7.[/latex]

No formula exists that allows us to find the solutions of [latex]f\left(x\right)=0.[/latex] Similar difficulties exist for nonpolynomial functions. For example, consider the task of finding solutions of [latex]\text{tan}\left(x\right)-x=0.[/latex] No simple formula exists for the solutions of this equation. In cases such as these, we can use Newton’s method to approximate the roots.

Newton’s method makes use of the following idea to approximate the solutions of [latex]f\left(x\right)=0.[/latex] By sketching a graph of [latex]f,[/latex] we can estimate a root of [latex]f\left(x\right)=0.[/latex] Let’s call this estimate [latex]{x}_{0}.[/latex] We then draw the tangent line to [latex]f[/latex] at [latex]{x}_{0}.[/latex] If [latex]{f}^{\prime }\left({x}_{0}\right)\ne 0,[/latex] this tangent line intersects the [latex]x[/latex]-axis at some point [latex]\left({x}_{1},0\right).[/latex] Now let [latex]{x}_{1}[/latex] be the next approximation to the actual root. Typically, [latex]{x}_{1}[/latex] is closer than [latex]{x}_{0}[/latex] to an actual root. Next we draw the tangent line to [latex]f[/latex] at [latex]{x}_{1}.[/latex] If [latex]{f}^{\prime }\left({x}_{1}\right)\ne 0,[/latex] this tangent line also intersects the [latex]x[/latex]-axis, producing another approximation, [latex]{x}_{2}.[/latex] We continue in this way, deriving a list of approximations: [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…}.[/latex] Typically, the numbers [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…}[/latex] quickly approach an actual root [latex]x*,[/latex] as shown in the following figure.

The approximations [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…}[/latex] approach the actual root [latex]x*.[/latex] The approximations are derived by looking at tangent lines to the graph of [latex]f.[/latex]
This function f(x) is drawn with points (x0, f(x0)), (x1, f(x1)), and (x2, f(x2)) marked on the function. From (x0, f(x0)), a tangent line is drawn, and it strikes the x axis at x1. From (x0, f(x0)), a tangent line is drawn, and it strikes the x axis at x2. If a tangent line were drawn from (x2, f(x2)), it appears that it would come very close to x*, which is the actual root. Each tangent line drawn in this order appears to get closer and closer to x*.

Now let’s look at how to calculate the approximations [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…}.[/latex] If [latex]{x}_{0}[/latex] is our first approximation, the approximation [latex]{x}_{1}[/latex] is defined by letting [latex]\left({x}_{1},0\right)[/latex] be the [latex]x[/latex]-intercept of the tangent line to [latex]f[/latex] at [latex]{x}_{0}.[/latex] The equation of this tangent line is given by

[latex]y=f\left({x}_{0}\right)+{f}^{\prime }\left({x}_{0}\right)\left(x-{x}_{0}\right).[/latex]

Therefore, [latex]{x}_{1}[/latex] must satisfy

[latex]f\left({x}_{0}\right)+{f}^{\prime }\left({x}_{0}\right)\left({x}_{1}-{x}_{0}\right)=0.[/latex]

Solving this equation for [latex]{x}_{1},[/latex] we conclude that

[latex]{x}_{1}={x}_{0}-\frac{f\left({x}_{0}\right)}{f\prime \left({x}_{0}\right)}.[/latex]

Similarly, the point [latex]\left({x}_{2},0\right)[/latex] is the [latex]x[/latex]-intercept of the tangent line to [latex]f[/latex] at [latex]{x}_{1}.[/latex] Therefore, [latex]{x}_{2}[/latex] satisfies the equation

[latex]{x}_{2}={x}_{1}-\frac{f\left({x}_{1}\right)}{f\prime \left({x}_{1}\right)}.[/latex]

In general, for [latex]n>0,{x}_{n}[/latex] satisfies

[latex]{x}_{n}={x}_{n-1}-\frac{f\left({x}_{n-1}\right)}{f\prime \left({x}_{n-1}\right)}.[/latex]

Next we see how to make use of this technique to approximate the root of the polynomial [latex]f\left(x\right)={x}^{3}-3x+1.[/latex]

Finding a Root of a Polynomial

Use Newton’s method to approximate a root of [latex]f\left(x\right)={x}^{3}-3x+1[/latex] in the interval [latex]\left[1,2\right].[/latex] Let [latex]{x}_{0}=2[/latex] and find [latex]{x}_{1},{x}_{2},{x}_{3},{x}_{4},[/latex] and [latex]{x}_{5}.[/latex]

From [link], we see that [latex]f[/latex] has one root over the interval [latex]\left(1,2\right).[/latex] Therefore [latex]{x}_{0}=2[/latex] seems like a reasonable first approximation. To find the next approximation, we use [link]. Since [latex]f\left(x\right)={x}^{3}-3x+1,[/latex] the derivative is [latex]{f}^{\prime }\left(x\right)=3{x}^{2}-3.[/latex] Using [link] with [latex]n=1[/latex] (and a calculator that displays [latex]10[/latex] digits), we obtain

[latex]{x}_{1}={x}_{0}-\frac{f\left({x}_{0}\right)}{f\prime \left({x}_{0}\right)}=2-\frac{f\left(2\right)}{f\prime \left(2\right)}=2-\frac{3}{9}\approx 1.666666667.[/latex]

To find the next approximation, [latex]{x}_{2},[/latex] we use [link] with [latex]n=2[/latex] and the value of [latex]{x}_{1}[/latex] stored on the calculator. We find that

[latex]{x}_{2}={x}_{1}=\frac{f\left({x}_{1}\right)}{f\prime \left({x}_{1}\right)}\approx 1.548611111.[/latex]

Continuing in this way, we obtain the following results:

[latex]\begin{array}{c}{x}_{1}\approx 1.666666667\hfill \\ {x}_{2}\approx 1.548611111\hfill \\ {x}_{3}\approx 1.532390162\hfill \\ {x}_{4}\approx 1.532088989\hfill \\ {x}_{5}\approx 1.532088886\hfill \\ {x}_{6}\approx 1.532088886.\hfill \end{array}[/latex]

We note that we obtained the same value for [latex]{x}_{5}[/latex] and [latex]{x}_{6}.[/latex] Therefore, any subsequent application of Newton’s method will most likely give the same value for [latex]{x}_{n}.[/latex]

The function [latex]f\left(x\right)={x}^{3}-3x+1[/latex] has one root over the interval [latex]\left[1,2\right].[/latex]
The function f(x) = x3 – 3x + 1 is drawn. It has roots between −2 and −1, 0 and 1, and 1 and 2.

Letting [latex]{x}_{0}=0,[/latex] let’s use Newton’s method to approximate the root of [latex]f\left(x\right)={x}^{3}-3x+1[/latex] over the interval [latex]\left[0,1\right][/latex] by calculating [latex]{x}_{1}[/latex] and [latex]{x}_{2}.[/latex]

[latex]{x}_{1}\approx 0.33333333,{x}_{2}\approx 0.347222222[/latex]

Hint

Use [link].

Newton’s method can also be used to approximate square roots. Here we show how to approximate [latex]\sqrt{2}.[/latex] This method can be modified to approximate the square root of any positive number.

Finding a Square Root

Use Newton’s method to approximate [latex]\sqrt{2}[/latex] ([link]). Let [latex]f\left(x\right)={x}^{2}-2,[/latex] let [latex]{x}_{0}=2,[/latex] and calculate [latex]{x}_{1},{x}_{2},{x}_{3},{x}_{4},{x}_{5}.[/latex] (We note that since [latex]f\left(x\right)={x}^{2}-2[/latex] has a zero at [latex]\sqrt{2},[/latex] the initial value [latex]{x}_{0}=2[/latex] is a reasonable choice to approximate [latex]\sqrt{2}.\right)[/latex]

For [latex]f\left(x\right)={x}^{2}-2,{f}^{\prime }\left(x\right)=2x.[/latex] From [link], we know that

[latex]\begin{array}{cc}\hfill {x}_{n}& ={x}_{n-1}-\frac{f\left({x}_{n-1}\right)}{f\prime \left({x}_{n-1}\right)}\hfill \\ & ={x}_{n-1}-\frac{{x}^{2}{}_{n-1}-2}{2{x}_{n-1}}\hfill \\ & =\frac{1}{2}{x}_{n-1}+\frac{1}{{x}_{n-1}}\hfill \\ & =\frac{1}{2}\left({x}_{n-1}+\frac{2}{{x}_{n-1}}\right).\hfill \end{array}[/latex]

Therefore,

[latex]\begin{array}{}\\ \\ {x}_{1}=\frac{1}{2}\left({x}_{0}+\frac{2}{{x}_{0}}\right)=\frac{1}{2}\left(2+\frac{2}{2}\right)=1.5\hfill \\ {x}_{2}=\frac{1}{2}\left({x}_{1}+\frac{2}{{x}_{1}}\right)=\frac{1}{2}\left(1.5+\frac{2}{1.5}\right)\approx 1.416666667.\hfill \end{array}[/latex]

Continuing in this way, we find that

[latex]\begin{array}{c}{x}_{1}=1.5\hfill \\ {x}_{2}\approx 1.416666667\hfill \\ {x}_{3}\approx 1.414215686\hfill \\ {x}_{4}\approx 1.414213562\hfill \\ {x}_{5}\approx 1.414213562.\hfill \end{array}[/latex]

Since we obtained the same value for [latex]{x}_{4}[/latex] and [latex]{x}_{5},[/latex] it is unlikely that the value [latex]{x}_{n}[/latex] will change on any subsequent application of Newton’s method. We conclude that [latex]\sqrt{2}\approx 1.414213562.[/latex]

We can use Newton’s method to find [latex]\sqrt{2}.[/latex]
The function y = x2 – 2 is drawn. A dashed line comes up from x0 = 2, and a tangent line is drawn down from there. It touches x1 = 1.5, which is near x* = the square root of 2.

Use Newton’s method to approximate [latex]\sqrt{3}[/latex] by letting [latex]f\left(x\right)={x}^{2}-3[/latex] and [latex]{x}_{0}=3.[/latex] Find [latex]{x}_{1}[/latex] and [latex]{x}_{2}.[/latex]

[latex]{x}_{1}=2,{x}_{2}=1.75[/latex]

Hint

For [latex]f\left(x\right)={x}^{2}-3,[/latex] [link] reduces to [latex]{x}_{n}=\frac{{x}_{n-1}}{2}+\frac{3}{2{x}_{n-1}}.[/latex]

When using Newton’s method, each approximation after the initial guess is defined in terms of the previous approximation by using the same formula. In particular, by defining the function [latex]F\left(x\right)=x-\left[\frac{f\left(x\right)}{{f}^{\prime }\left(x\right)}\right],[/latex] we can rewrite [link] as [latex]{x}_{n}=F\left({x}_{n-1}\right).[/latex] This type of process, where each [latex]{x}_{n}[/latex] is defined in terms of [latex]{x}_{n-1}[/latex] by repeating the same function, is an example of an iterative process. Shortly, we examine other iterative processes. First, let’s look at the reasons why Newton’s method could fail to find a root.

Failures of Newton’s Method

Typically, Newton’s method is used to find roots fairly quickly. However, things can go wrong. Some reasons why Newton’s method might fail include the following:

  1. At one of the approximations [latex]{x}_{n},[/latex] the derivative [latex]{f}^{\prime }[/latex] is zero at [latex]{x}_{n},[/latex] but [latex]f\left({x}_{n}\right)\ne 0.[/latex] As a result, the tangent line of [latex]f[/latex] at [latex]{x}_{n}[/latex] does not intersect the [latex]x[/latex]-axis. Therefore, we cannot continue the iterative process.
  2. The approximations [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…}[/latex] may approach a different root. If the function [latex]f[/latex] has more than one root, it is possible that our approximations do not approach the one for which we are looking, but approach a different root (see [link]). This event most often occurs when we do not choose the approximation [latex]{x}_{0}[/latex] close enough to the desired root.
  3. The approximations may fail to approach a root entirely. In [link], we provide an example of a function and an initial guess [latex]{x}_{0}[/latex] such that the successive approximations never approach a root because the successive approximations continue to alternate back and forth between two values.
If the initial guess [latex]{x}_{0}[/latex] is too far from the root sought, it may lead to approximations that approach a different root.
A function is drawn with two roots, labeled root sought and root found. A point x0 is chosen such that when the tangent of x0 is taken, even though it is nearer to the root sought, the tangent points to the root found.
When Newton’s Method Fails

Consider the function [latex]f\left(x\right)={x}^{3}-2x+2.[/latex] Let [latex]{x}_{0}=0.[/latex] Show that the sequence [latex]{x}_{1},{x}_{2}\text{,…}[/latex] fails to approach a root of [latex]f.[/latex]

For [latex]f\left(x\right)={x}^{3}-2x+2,[/latex] the derivative is [latex]{f}^{\prime }\left(x\right)=3{x}^{2}-2.[/latex] Therefore,

[latex]{x}_{1}={x}_{0}-\frac{f\left({x}_{0}\right)}{{f}^{\prime }\left({x}_{0}\right)}=0-\frac{f\left(0\right)}{{f}^{\prime }\left(0\right)}=-\frac{2}{-2}=1.[/latex]

In the next step,

[latex]{x}_{2}={x}_{1}-\frac{f\left({x}_{1}\right)}{f\prime \left({x}_{1}\right)}=1-\frac{f\left(1\right)}{{f}^{\prime }\left(1\right)}=1-\frac{1}{1}=0.[/latex]

Consequently, the numbers [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…}[/latex] continue to bounce back and forth between [latex]0[/latex] and [latex]1[/latex] and never get closer to the root of [latex]f[/latex] which is over the interval [latex]\left[-2,-1\right][/latex] (see [link]). Fortunately, if we choose an initial approximation [latex]{x}_{0}[/latex] closer to the actual root, we can avoid this situation.

The approximations continue to alternate between [latex]0\phantom{\rule{0.2em}{0ex}}\text{and}\phantom{\rule{0.2em}{0ex}}1[/latex] and never approach the root of [latex]f.[/latex]
The function f(x) = x3 – 2x + 2 is drawn, which has a root between −2 and −1. The tangent from x = 0 goes to x = 1, and the tangent from x = 1 goes to x = 0.

For [latex]f\left(x\right)={x}^{3}-2x+2,[/latex] let [latex]{x}_{0}=-1.5[/latex] and find [latex]{x}_{1}[/latex] and [latex]{x}_{2}.[/latex]

[latex]{x}_{1}\approx -1.842105263,{x}_{2}\approx -1.772826920[/latex]

Hint

Use [link].

From [link], we see that Newton’s method does not always work. However, when it does work, the sequence of approximations approaches the root very quickly. Discussions of how quickly the sequence of approximations approach a root found using Newton’s method are included in texts on numerical analysis.

Other Iterative Processes

As mentioned earlier, Newton’s method is a type of iterative process. We now look at an example of a different type of iterative process.

Consider a function [latex]F[/latex] and an initial number [latex]{x}_{0}.[/latex] Define the subsequent numbers [latex]{x}_{n}[/latex] by the formula [latex]{x}_{n}=F\left({x}_{n-1}\right).[/latex] This process is an iterative process that creates a list of numbers [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…},{x}_{n}\text{,…}.[/latex] This list of numbers may approach a finite number [latex]x*[/latex] as [latex]n[/latex] gets larger, or it may not. In [link], we see an example of a function [latex]F[/latex] and an initial guess [latex]{x}_{0}[/latex] such that the resulting list of numbers approaches a finite value.

Finding a Limit for an Iterative Process

Let [latex]F\left(x\right)=\frac{1}{2}x+4[/latex] and let [latex]{x}_{0}=0.[/latex] For all [latex]n\ge 1,[/latex] let [latex]{x}_{n}=F\left({x}_{n-1}\right).[/latex] Find the values [latex]{x}_{1},{x}_{2},{x}_{3},{x}_{4},{x}_{5}.[/latex] Make a conjecture about what happens to this list of numbers [latex]{x}_{1},{x}_{2},{x}_{3}\text{…},{x}_{n}\text{,…}[/latex] as [latex]n\to \infty .[/latex] If the list of numbers [latex]{x}_{1},{x}_{2},{x}_{3}\text{,…}[/latex] approaches a finite number [latex]x*,[/latex] then [latex]x*[/latex] satisfies [latex]x*=F\left(x*\right),[/latex] and [latex]x*[/latex] is called a fixed point of [latex]F.[/latex]

If [latex]{x}_{0}=0,[/latex] then

[latex]\begin{array}{}\\ \\ {x}_{1}=\frac{1}{2}\left(0\right)+4=4\hfill \\ {x}_{2}=\frac{1}{2}\left(4\right)+4=6\hfill \\ {x}_{3}=\frac{1}{2}\left(6\right)+4=7\hfill \\ {x}_{4}=\frac{1}{2}\left(7\right)+4=7.5\hfill \\ {x}_{5}=\frac{1}{2}\left(7.5\right)+4=7.75\hfill \\ {x}_{6}=\frac{1}{2}\left(7.75\right)+4=7.875\hfill \\ {x}_{7}=\frac{1}{2}\left(7.875\right)+4=7.9375\hfill \\ {x}_{8}=\frac{1}{2}\left(7.9375\right)+4=7.96875\hfill \\ {x}_{9}=\frac{1}{2}\left(7.96875\right)+4=7.984375.\hfill \end{array}[/latex]

From this list, we conjecture that the values [latex]{x}_{n}[/latex] approach [latex]8.[/latex]

[link] provides a graphical argument that the values approach [latex]8[/latex] as [latex]n\to \infty .[/latex] Starting at the point [latex]\left({x}_{0},{x}_{0}\right),[/latex] we draw a vertical line to the point [latex]\left({x}_{0},F\left({x}_{0}\right)\right).[/latex] The next number in our list is [latex]{x}_{1}=F\left({x}_{0}\right).[/latex] We use [latex]{x}_{1}[/latex] to calculate [latex]{x}_{2}.[/latex] Therefore, we draw a horizontal line connecting [latex]\left({x}_{0},{x}_{1}\right)[/latex] to the point [latex]\left({x}_{1},{x}_{1}\right)[/latex] on the line [latex]y=x,[/latex] and then draw a vertical line connecting [latex]\left({x}_{1},{x}_{1}\right)[/latex] to the point [latex]\left({x}_{1},F\left({x}_{1}\right)\right).[/latex] The output [latex]F\left({x}_{1}\right)[/latex] becomes [latex]{x}_{2}.[/latex] Continuing in this way, we could create an infinite number of line segments. These line segments are trapped between the lines [latex]F\left(x\right)=\frac{x}{2}+4[/latex] and [latex]y=x.[/latex] The line segments get closer to the intersection point of these two lines, which occurs when [latex]x=F\left(x\right).[/latex] Solving the equation [latex]x=\frac{x}{2}+4,[/latex] we conclude they intersect at [latex]x=8.[/latex] Therefore, our graphical evidence agrees with our numerical evidence that the list of numbers [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…}[/latex] approaches [latex]x*=8[/latex] as [latex]n\to \infty .[/latex]

This iterative process approaches the value [latex]x*=8.[/latex]
The function F(x) = (1/2)x + 4 is graphed along with y = x. From x0, which appears to be at the origin, a line is drawn to the function F(x) at x1 = F(x0). Then a line is drawn to the right from here to the line y = x, at which point a line is drawn up to x2 = F(x1). Then a line is drawn to the right from here to the line y = x, at which point a line is drawn up to x3 = F(x2). Continuing this process would converge on the two lines’ intersection point at x* = 8.

Consider the function [latex]F\left(x\right)=\frac{1}{3}x+6.[/latex] Let [latex]{x}_{0}=0[/latex] and let [latex]{x}_{n}=F\left({x}_{n-1}\right)[/latex] for [latex]n\ge 2.[/latex] Find [latex]{x}_{1},{x}_{2},{x}_{3},{x}_{4},{x}_{5}.[/latex] Make a conjecture about what happens to the list of numbers [latex]{x}_{1},{x}_{2},{x}_{3}\text{,…}{x}_{n}\text{,…}[/latex] as [latex]n\to \infty .[/latex]

[latex]{x}_{1}=6,{x}_{2}=8,{x}_{3}=\frac{26}{3},{x}_{4}=\frac{80}{9},{x}_{5}=\frac{242}{27};x*=9[/latex]

Hint

Consider the point where the lines [latex]y=x[/latex] and [latex]y=F\left(x\right)[/latex] intersect.

Iterative Processes and Chaos

Iterative processes can yield some very interesting behavior. In this section, we have seen several examples of iterative processes that converge to a fixed point. We also saw in [link] that the iterative process bounced back and forth between two values. We call this kind of behavior a [latex]2[/latex]-cycle. Iterative processes can converge to cycles with various periodicities, such as [latex]2-\text{cycles},\phantom{\rule{0.2em}{0ex}}4-\text{cycles}[/latex] (where the iterative process repeats a sequence of four values), 8-cycles, and so on.

Some iterative processes yield what mathematicians call chaos. In this case, the iterative process jumps from value to value in a seemingly random fashion and never converges or settles into a cycle. Although a complete exploration of chaos is beyond the scope of this text, in this project we look at one of the key properties of a chaotic iterative process: sensitive dependence on initial conditions. This property refers to the concept that small changes in initial conditions can generate drastically different behavior in the iterative process.

Probably the best-known example of chaos is the Mandelbrot set (see [link]), named after Benoit Mandelbrot (1924–2010), who investigated its properties and helped popularize the field of chaos theory. The Mandelbrot set is usually generated by computer and shows fascinating details on enlargement, including self-replication of the set. Several colorized versions of the set have been shown in museums and can be found online and in popular books on the subject.

The Mandelbrot set is a well-known example of a set of points generated by the iterative chaotic behavior of a relatively simple function.
A very complicated and organic looking fractal.

In this project we use the logistic map

[latex]f\left(x\right)=rx\left(1-x\right),\phantom{\rule{0.2em}{0ex}}\text{where}\phantom{\rule{0.2em}{0ex}}x\in \left[0,1\right]\phantom{\rule{0.2em}{0ex}}\text{and}\phantom{\rule{0.2em}{0ex}}r>0[/latex]

as the function in our iterative process. The logistic map is a deceptively simple function; but, depending on the value of [latex]r,[/latex] the resulting iterative process displays some very interesting behavior. It can lead to fixed points, cycles, and even chaos.

To visualize the long-term behavior of the iterative process associated with the logistic map, we will use a tool called a cobweb diagram. As we did with the iterative process we examined earlier in this section, we first draw a vertical line from the point [latex]\left({x}_{0},0\right)[/latex] to the point [latex]\left({x}_{0},f\left({x}_{0}\right)\right)=\left({x}_{0},{x}_{1}\right).[/latex] We then draw a horizontal line from that point to the point [latex]\left({x}_{1},{x}_{1}\right),[/latex] then draw a vertical line to [latex]\left({x}_{1},f\left({x}_{1}\right)\right)=\left({x}_{1},{x}_{2}\right),[/latex] and continue the process until the long-term behavior of the system becomes apparent. [link] shows the long-term behavior of the logistic map when [latex]r=3.55[/latex] and [latex]{x}_{0}=0.2.[/latex] (The first [latex]100[/latex] iterations are not plotted.) The long-term behavior of this iterative process is an [latex]8[/latex]-cycle.

A cobweb diagram for [latex]f\left(x\right)=3.55x\left(1-x\right)[/latex] is presented here. The sequence of values results in an [latex]8[/latex]-cycle.
In the first quadrant, f(x) = 3.55x(1 – x) is graphed as is y = x. From some point on the x axis, a line is drawn up to the line y = x, at which point it turns to be horizontal and continues until it touches the outside edge of f(x), at which point it turns again to be vertical until it each the line y = x. This process continues a number of times and creates an interesting series of boxes.
  1. Let [latex]r=0.5[/latex] and choose [latex]{x}_{0}=0.2.[/latex] Either by hand or by using a computer, calculate the first [latex]10[/latex] values in the sequence. Does the sequence appear to converge? If so, to what value? Does it result in a cycle? If so, what kind of cycle (for example, [latex]2-\text{cycle},\phantom{\rule{0.2em}{0ex}}4-\text{cycle}.\right)?[/latex]
  2. What happens when [latex]r=2?[/latex]
  3. For [latex]r=3.2[/latex] and [latex]r=3.5,[/latex] calculate the first [latex]100[/latex] sequence values. Generate a cobweb diagram for each iterative process. (Several free applets are available online that generate cobweb diagrams for the logistic map.) What is the long-term behavior in each of these cases?
  4. Now let [latex]r=4.[/latex] Calculate the first [latex]100[/latex] sequence values and generate a cobweb diagram. What is the long-term behavior in this case?
  5. Repeat the process for [latex]r=4,[/latex] but let [latex]{x}_{0}=0.201.[/latex] How does this behavior compare with the behavior for [latex]{x}_{0}=0.2?[/latex]

Key Concepts

  • Newton’s method approximates roots of [latex]f\left(x\right)=0[/latex] by starting with an initial approximation [latex]{x}_{0},[/latex] then uses tangent lines to the graph of [latex]f[/latex] to create a sequence of approximations [latex]{x}_{1},{x}_{2},{x}_{3}\text{,…}.[/latex]
  • Typically, Newton’s method is an efficient method for finding a particular root. In certain cases, Newton’s method fails to work because the list of numbers [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…}[/latex] does not approach a finite value or it approaches a value other than the root sought.
  • Any process in which a list of numbers [latex]{x}_{0},{x}_{1},{x}_{2}\text{,…}[/latex] is generated by defining an initial number [latex]{x}_{0}[/latex] and defining the subsequent numbers by the equation [latex]{x}_{n}=F\left({x}_{n-1}\right)[/latex] for some function [latex]F[/latex] is an iterative process. Newton’s method is an example of an iterative process, where the function [latex]F\left(x\right)=x-\left[\frac{f\left(x\right)}{{f}^{\prime }\left(x\right)}\right][/latex] for a given function [latex]f.[/latex]

For the following exercises, write Newton’s formula as [latex]{x}_{n+1}=F\left({x}_{n}\right)[/latex] for solving [latex]f\left(x\right)=0.[/latex]

[latex]f\left(x\right)={x}^{2}+1[/latex]

[latex]f\left(x\right)={x}^{3}+2x+1[/latex]

[latex]F\left({x}_{n}\right)={x}_{n}-\frac{{x}_{n}{}^{3}+2{x}_{n}+1}{3{x}_{n}{}^{2}+2}[/latex]

[latex]f\left(x\right)=\text{sin}\phantom{\rule{0.1em}{0ex}}x[/latex]

[latex]f\left(x\right)={e}^{x}[/latex]

[latex]F\left({x}_{n}\right)={x}_{n}-\frac{{e}^{{x}_{n}}}{{e}^{{x}_{n}}}[/latex]

[latex]f\left(x\right)={x}^{3}+3x{e}^{x}[/latex]

For the following exercises, solve [latex]f\left(x\right)=0[/latex] using the iteration [latex]{x}_{n+1}={x}_{n}-cf\left({x}_{n}\right),[/latex] which differs slightly from Newton’s method. Find a [latex]c[/latex] that works and a [latex]c[/latex] that fails to converge, with the exception of [latex]c=0.[/latex]

[latex]f\left(x\right)={x}^{2}-4,[/latex] with [latex]{x}_{0}=0[/latex]

[latex]|c|>0.5[/latex] fails, [latex]|c|\le 0.5[/latex] works

[latex]f\left(x\right)={x}^{2}-4x+3,[/latex] with [latex]{x}_{0}=2[/latex]

What is the value of [latex]\text{“}c\text{”}[/latex] for Newton’s method?

[latex]c=\frac{1}{{f}^{\prime }\left({x}_{n}\right)}[/latex]

For the following exercises, start at

a. [latex]{x}_{0}=0.6[/latex] and

b. [latex]{x}_{0}=2.[/latex]

Compute [latex]{x}_{1}[/latex] and [latex]{x}_{2}[/latex] using the specified iterative method.

[latex]{x}_{n+1}={x}_{n}{}^{2}-\frac{1}{2}[/latex]

[latex]{x}_{n+1}=2{x}_{n}\left(1-{x}_{n}\right)[/latex]

a. [latex]{x}_{1}=\frac{12}{25},{x}_{2}=\frac{312}{625};[/latex] b. [latex]{x}_{1}=-4,{x}_{2}=-40[/latex]

[latex]{x}_{n+1}=\sqrt{{x}_{n}}[/latex]

[latex]{x}_{n+1}=\frac{1}{\sqrt{{x}_{n}}}[/latex]

a. [latex]{x}_{1}=1.291,{x}_{2}=0.8801;[/latex] b. [latex]{x}_{1}=0.7071,{x}_{2}=1.189[/latex]

[latex]{x}_{n+1}=3{x}_{n}\left(1-{x}_{n}\right)[/latex]

[latex]{x}_{n+1}={x}_{n}{}^{2}+{x}_{n}-2[/latex]

a. [latex]{x}_{1}=-\frac{26}{25},{x}_{2}=-\frac{1224}{625};[/latex] b. [latex]{x}_{1}=4,{x}_{2}=18[/latex]

[latex]{x}_{n+1}=\frac{1}{2}{x}_{n}-1[/latex]

[latex]{x}_{n+1}=|{x}_{n}|[/latex]

a. [latex]{x}_{1}=\frac{6}{10},{x}_{2}=\frac{6}{10};[/latex] b. [latex]{x}_{1}=2,{x}_{2}=2[/latex]

For the following exercises, solve to four decimal places using Newton’s method and a computer or calculator. Choose any initial guess [latex]{x}_{0}[/latex] that is not the exact root.

[latex]{x}^{2}-10=0[/latex]

[latex]{x}^{4}-100=0[/latex]

[latex]3.1623\phantom{\rule{0.2em}{0ex}}\text{or}\phantom{\rule{0.2em}{0ex}}-3.1623[/latex]

[latex]{x}^{2}-x=0[/latex]

[latex]{x}^{3}-x=0[/latex]

[latex]0,-1\phantom{\rule{0.2em}{0ex}}\text{or}\phantom{\rule{0.2em}{0ex}}1[/latex]

[latex]x+5\phantom{\rule{0.1em}{0ex}}\text{cos}\left(x\right)=0[/latex]

[latex]x+\text{tan}\left(x\right)=0,[/latex] choose [latex]{x}_{0}\in \left(-\frac{\pi }{2},\frac{\pi }{2}\right)[/latex]

[latex]0[/latex]

[latex]\frac{1}{1-x}=2[/latex]

[latex]1+x+{x}^{2}+{x}^{3}+{x}^{4}=2[/latex]

[latex]0.5188\phantom{\rule{0.2em}{0ex}}\text{or}\phantom{\rule{0.2em}{0ex}}-1.2906[/latex]

[latex]{x}^{3}+{\left(x+1\right)}^{3}={10}^{3}[/latex]

[latex]x={\text{sin}}^{2}\left(x\right)[/latex]

[latex]0[/latex]

For the following exercises, use Newton’s method to find the fixed points of the function where [latex]f\left(x\right)=x;[/latex] round to three decimals.

[latex]\text{sin}\phantom{\rule{0.1em}{0ex}}x[/latex]

[latex]\text{tan}\left(x\right)[/latex] on [latex]x=\left(\frac{\pi }{2},\frac{3\pi }{2}\right)[/latex]

[latex]4.493[/latex]

[latex]{e}^{x}-2[/latex]

[latex]\text{ln}\left(x\right)+2[/latex]

[latex]0.159,3.146[/latex]

Newton’s method can be used to find maxima and minima of functions in addition to the roots. In this case apply Newton’s method to the derivative function [latex]{f}^{\prime }\left(x\right)[/latex] to find its roots, instead of the original function. For the following exercises, consider the formulation of the method.

To find candidates for maxima and minima, we need to find the critical points [latex]{f}^{\prime }\left(x\right)=0.[/latex] Show that to solve for the critical points of a function [latex]f\left(x\right),[/latex] Newton’s method is given by [latex]{x}_{n+1}={x}_{n}-\frac{{f}^{\prime }\left({x}_{n}\right)}{f\text{″}\left({x}_{n}\right)}.[/latex]

What additional restrictions are necessary on the function [latex]f?[/latex]

We need [latex]f[/latex] to be twice continuously differentiable.

For the following exercises, use Newton’s method to find the location of the local minima and/or maxima of the following functions; round to three decimals.

Minimum of [latex]f\left(x\right)={x}^{2}+2x+4[/latex]

Minimum of [latex]f\left(x\right)=3{x}^{3}+2{x}^{2}-16[/latex]

[latex]x=0[/latex]

Minimum of [latex]f\left(x\right)={x}^{2}{e}^{x}[/latex]

Maximum of [latex]f\left(x\right)=x+\frac{1}{x}[/latex]

[latex]x=-1[/latex]

Maximum of [latex]f\left(x\right)={x}^{3}+10{x}^{2}+15x-2[/latex]

Maximum of [latex]f\left(x\right)=\frac{\sqrt{x}-\sqrt[3]{x}}{x}[/latex]

[latex]x=5.619[/latex]

Minimum of [latex]f\left(x\right)={x}^{2}\text{sin}\phantom{\rule{0.1em}{0ex}}x,[/latex] closest non-zero minimum to [latex]x=0[/latex]

Minimum of [latex]f\left(x\right)={x}^{4}+{x}^{3}+3{x}^{2}+12x+6[/latex]

[latex]x=-1.326[/latex]

For the following exercises, use the specified method to solve the equation. If it does not work, explain why it does not work.

Newton’s method, [latex]{x}^{2}+2=0[/latex]

Newton’s method, [latex]0={e}^{x}[/latex]

There is no solution to the equation.

Newton’s method, [latex]0=1+{x}^{2}[/latex] starting at [latex]{x}_{0}=0[/latex]

Solving [latex]{x}_{n+1}=\text{−}{x}_{n}{}^{3}[/latex] starting at [latex]{x}_{0}=-1[/latex]

It enters a cycle.

For the following exercises, use the secant method, an alternative iterative method to Newton’s method. The formula is given by

[latex]{x}_{n}={x}_{n-1}-f\left({x}_{n-1}\right)\frac{{x}_{n-1}-{x}_{n-2}}{f\left({x}_{n-1}\right)-f\left({x}_{n-2}\right)}.[/latex]

Find a root to [latex]0={x}^{2}-x-3[/latex] accurate to three decimal places.

Find a root to [latex]0=\text{sin}\phantom{\rule{0.1em}{0ex}}x+3x[/latex] accurate to four decimal places.

[latex]0[/latex]

Find a root to [latex]0={e}^{x}-2[/latex] accurate to four decimal places.

Find a root to [latex]\text{ln}\left(x+2\right)=\frac{1}{2}[/latex] accurate to four decimal places.

[latex]-0.3513[/latex]

Why would you use the secant method over Newton’s method? What are the necessary restrictions on [latex]f?[/latex]

For the following exercises, use both Newton’s method and the secant method to calculate a root for the following equations. Use a calculator or computer to calculate how many iterations of each are needed to reach within three decimal places of the exact answer. For the secant method, use the first guess from Newton’s method.

[latex]f\left(x\right)={x}^{2}+2x+1,{x}_{0}=1[/latex]

Newton: [latex]11[/latex] iterations, secant: [latex]16[/latex] iterations

[latex]f\left(x\right)={x}^{2},{x}_{0}=1[/latex]

[latex]f\left(x\right)=\text{sin}\phantom{\rule{0.1em}{0ex}}x,{x}_{0}=1[/latex]

Newton: three iterations, secant: six iterations

[latex]f\left(x\right)={e}^{x}-1,{x}_{0}=2[/latex]

[latex]f\left(x\right)={x}^{3}+2x+4,{x}_{0}=0[/latex]

Newton: five iterations, secant: eight iterations

In the following exercises, consider Kepler’s equation regarding planetary orbits, [latex]M=E-\epsilon \phantom{\rule{0.1em}{0ex}}\text{sin}\left(E\right),[/latex] where [latex]M[/latex] is the mean anomaly, [latex]E[/latex] is eccentric anomaly, and [latex]\epsilon [/latex] measures eccentricity.

Use Newton’s method to solve for the eccentric anomaly [latex]E[/latex] when the mean anomaly [latex]M=\frac{\pi }{3}[/latex] and the eccentricity of the orbit [latex]\epsilon =0.25;[/latex] round to three decimals.

Use Newton’s method to solve for the eccentric anomaly [latex]E[/latex] when the mean anomaly [latex]M=\frac{3\pi }{2}[/latex] and the eccentricity of the orbit [latex]\epsilon =0.8;[/latex] round to three decimals.

[latex]E=4.071[/latex]

The following two exercises consider a bank investment. The initial investment is [latex]\text{💲}10,000.[/latex] After [latex]25[/latex] years, the investment has tripled to [latex]\text{💲}30,000.[/latex]

Use Newton’s method to determine the interest rate if the interest was compounded annually.

Use Newton’s method to determine the interest rate if the interest was compounded continuously.

[latex]4.394\text{%}[/latex]

The cost for printing a book can be given by the equation [latex]C\left(x\right)=1000+12x+\left(\frac{1}{2}\right){x}^{2\text{/}3}.[/latex] Use Newton’s method to find the break-even point if the printer sells each book for [latex]\text{💲}20.[/latex]

Glossary

iterative process
process in which a list of numbers [latex]{x}_{0},{x}_{1},{x}_{2},{x}_{3}\text{…}[/latex] is generated by starting with a number [latex]{x}_{0}[/latex] and defining [latex]{x}_{n}=F\left({x}_{n-1}\right)[/latex] for [latex]n\ge 1[/latex]
Newton’s method
method for approximating roots of [latex]f\left(x\right)=0;[/latex] using an initial guess [latex]{x}_{0};[/latex] each subsequent approximation is defined by the equation [latex]{x}_{n}={x}_{n-1}-\frac{f\left({x}_{n-1}\right)}{f\prime \left({x}_{n-1}\right)}[/latex]