Secant method: Difference between revisions

From Calculus
(Created page with "{{root-finding algorithm}} ==Definition== The '''secant method''' is a root-finding algorithm.")
 
No edit summary
 
(7 intermediate revisions by the same user not shown)
Line 2: Line 2:
==Definition==
==Definition==


The '''secant method''' is a root-finding algorithm.
The '''secant method''' is a root-finding algorithm that makes successive point estimates for the value of a root of a [[continuous function]]. In general, the secant method is not guaranteed to converge towards a root, but under some conditions, it does. A slight variant of this method, called the [[false position method]], functions very similarly to the [[bisection method]].
 
==Summary==
 
{| class="sortable" border="1"
! Item !! Value
|-
| Initial condition || we are given a function, typically a [[continuous function]], and two initial guesses <math>x_0</math> and <math>x_1</math> for roots of the function.
|-
| Iterative step || At stage <math>n</math> (<math>n \ge 2</math>) we compute <math>x_n</math> in terms of <math>x_{n-1}, x_{n-2}, f(x_{n-1}), f(x_{n-2})</math>.
|-
| Convergence rate || The [[order of convergence]] is the [[golden ratio]]
|-
| Computational tools needed || Function evaluation at particular points, multiplication, subtraction, and division
|-
| Termination || We may terminate based on known bounds on the size of the derivative and the function value coming in sufficient proximity to the actual value.
|}
 
==Initial condition==
 
The secant method requires ''two'' initial guesses for the root, say <math>x_0</math> and <math>x_1</math>.
 
==Iterative step==
 
For <math>n \ge 2</math>, we define <math>x_n</math> as the following affine combination of the previous two guesses <math>x_{n-1}</math> and <math>x_{n-2}</math>:
 
<math>x_n := \frac{x_{n-2}f(x_{n-1}) - x_{n-1}f(x_{n-2})}{f(x_{n-1}) - f(x_{n-2})}</math>
 
We can therefore think of <math>x_n</math> as an ''affine linear combination'' of <math>x_{n-2}</math> and <math>x_{n-1}</math> with the following respective coefficients:
 
<math>\frac{f(x_{n-1})}{f(x_{n-1}) - f(x_{n-2})}, \frac{-f(x_{n-2})}{f(x_{n-1}) - f(x_{n-2})}</math>
 
Geometrically, this can be interpreted as follows: we make a line through the points <math>(x_{n-2},f(x_{n-2}))</math> and <math>(x_{n-1},f(x_{n-1}))</math> in the <math>(x,f(x))</math>-plane, and define <math>x_n</math> as the <math>x</math>-coordinate of the intersection of this line with the <math>x</math>-axis.
 
==Convergence rate==

Latest revision as of 16:47, 26 April 2014

This article is about a root-finding algorithm. See all root-finding algorithms

Definition

The secant method is a root-finding algorithm that makes successive point estimates for the value of a root of a continuous function. In general, the secant method is not guaranteed to converge towards a root, but under some conditions, it does. A slight variant of this method, called the false position method, functions very similarly to the bisection method.

Summary

Item Value
Initial condition we are given a function, typically a continuous function, and two initial guesses x0 and x1 for roots of the function.
Iterative step At stage n (n2) we compute xn in terms of xn1,xn2,f(xn1),f(xn2).
Convergence rate The order of convergence is the golden ratio
Computational tools needed Function evaluation at particular points, multiplication, subtraction, and division
Termination We may terminate based on known bounds on the size of the derivative and the function value coming in sufficient proximity to the actual value.

Initial condition

The secant method requires two initial guesses for the root, say x0 and x1.

Iterative step

For n2, we define xn as the following affine combination of the previous two guesses xn1 and xn2:

xn:=xn2f(xn1)xn1f(xn2)f(xn1)f(xn2)

We can therefore think of xn as an affine linear combination of xn2 and xn1 with the following respective coefficients:

f(xn1)f(xn1)f(xn2),f(xn2)f(xn1)f(xn2)

Geometrically, this can be interpreted as follows: we make a line through the points (xn2,f(xn2)) and (xn1,f(xn1)) in the (x,f(x))-plane, and define xn as the x-coordinate of the intersection of this line with the x-axis.

Convergence rate