Steffensen's method: Difference between revisions

Revision as of 00:39, 12 July 2008

Revision as of 00:41, 12 July 2008

Line 12:

The function <math>g\ </math> is the average slope of the function <math>f\ </math> between the last sequence point <math>x=x_n,\ y=f(x_n)</math> and the auxilliary point <math>x=x_n + f(x_n),\ y=f(x_n + f(x_n))</math>&nbsp;.

The main advantage of Steffensen's method is that it can find the roots of an equation <math>f\ </math> just as "[[quadratic convergence|quickly]]" as [[Newton's method]] but the formula does not require a separate function for the derivative, so it can be programmed for any generic function. In this case ''[[quadratic convergence|quicly]]'' means that the number of correct digits in the answer doubles with each step. The cost, howeverfor the quick convergence is the double function evaluation: both <math>f(x_n)\ </math> and <math>f(x_n + f(x_n))\ </math> must be evaluatedcalculated, which migh be costly if <math>f\ </math> is a complicated function.

Similar to [[Newton's method]] and most other quadratically convergent methods, the crucial weakness with the method is the choice of the starting value <math>x_0\ </math>&nbsp;. If the value of <math>x_0\ </math> is not "close enough" to the actual solution, the method will fail and the sequence of values <math>x_0,\ x_1,\ x_2,\ x_3 ...</math> will either flip flop between two extremes, or diverge to infinity (possibly both!).

==Generalised definition==