next up previous
Next: Convergence Rate of Secant Up: Solution of Non-Linear Equations Previous: Newton's Method
The Secant Method
We can obtain this from Newton's method by replacing the tangent slope f'(x) by the chord or secant slope
. Using a Taylor expansion of f about x, we find
so that the error is O(h) if f'' is bounded. The obvious choice for h is so that if x=xn, the nth iterate, x-h=xn-1, the
previous iterate. Putting this into Newton's method, we obtain
In this method, f has to be evaluated just once on each iteration, while f' does not have to be evaluated at all.
However, it does need two starting values.
The formula can be written with a common denominator as
(2)
open in browser PRO version Are you a developer? Try out the HTML to PDF API pdfcrowd.com
which is recognisable as the formula for inverse interpolation, that is, it linearly interpolates the values of f at xn and
xn-1 to approximate the value of x at which f vanishes. This form should not be used in practical evaluation, since it
computes the ratio of two quantities which are liable to be differences of nearly equal numbers.
This formulation suggests a variant of the method, in which convergence is guaranteed provided starting values can
be found for which the values of f are opposite in sign. This is defined by
where is the greatest index for which . Convergence is guaranteed because there is
always a zero of f between xn and xm. The rate of convergence is better than linear, but not as good as that of the
original secant method.
Example Evaluate .
We use the secant method to solve the equation x2-15=0. The secant formula is
Starting with x0=3, x1=4, we calculate to 9D successively x2=3.857142857,
x4=3.872983871, x5=3.872983347,x6=3.872983347. This looks nearly as good as Newton's method,
but our starting values were fairly close.
next up previous
open in browser PRO version Are you a developer? Try out the HTML to PDF API pdfcrowd.com
Next: Convergence Rate of Secant Up: Solution of Non-Linear Equations Previous: Newton's Method
John Gilbert
5/8/1998
open in browser PRO version Are you a developer? Try out the HTML to PDF API pdfcrowd.com