Full linear multistep methods as root-finders

Root-finders based on full linear multistep methods (LMMs) use previous function values, derivatives and root estimates to iteratively find a root of a nonlinear function. As ODE solvers, full LMMs are typically not zero-stable. However, used as root-finders, the interpolation points are convergent so that such stability issues are circumvented. A general analysis is provided based on inverse polynomial interpolation, which is used to prove a fundamental barrier on the convergence rate of any LMM-based method. We show, using numerical examples, that full LMM-based methods perform excellently. Finally, we also provide a robust implementation based on Brents method that is guaranteed to converge.

[1]  D. Anderson,et al.  Algorithms for minimization without derivatives , 1974 .

[2]  Changbum Chun,et al.  An analysis of a new family of eighth-order optimal methods , 2014, Appl. Math. Comput..

[3]  Miquel Grau-Sánchez,et al.  Zero-finder methods derived from Obreshkov's techniques , 2009, Appl. Math. Comput..

[4]  Julio Chaves,et al.  Introduction to Nonimaging Optics , 2008 .

[5]  José Luis Díaz-Barrero,et al.  Zero-finder methods derived using Runge-Kutta techniques , 2011, Appl. Math. Comput..

[6]  José Luis Díaz-Barrero,et al.  Adams-like techniques for zero-finder methods , 2009, Appl. Math. Comput..

[7]  高等学校計算数学学報編輯委員会編 高等学校計算数学学報 = Numerical mathematics , 1979 .

[8]  Young Hee Geum,et al.  A family of optimal sixteenth-order multipoint methods with a linear fraction plus a trivariate polynomial as the fourth-step weighting function , 2011, Comput. Math. Appl..

[9]  Walter Gautschi,et al.  Numerical Analysis , 1978, Mathemagics: A Magical Journey Through Advanced Mathematics.

[10]  Stefan Siegmund,et al.  A new class of three-point methods with optimal convergence order eight and its dynamics , 2014, Numerical Algorithms.

[11]  Young Ik Kim,et al.  An Optimal family of Eighth-order iterative Methods with an inverse interpolatory rational function error corrector for nonlinear equations , 2017, Math. Model. Anal..

[12]  Xiangke Liao,et al.  A new fourth-order iterative method for finding multiple roots of nonlinear equations , 2009, Appl. Math. Comput..

[13]  Arnold R. Miller,et al.  Pathological functions for Newton's method , 1993 .

[14]  Fazlollah Soleymani,et al.  Two new classes of optimal Jarratt-type fourth-order methods , 2012, Appl. Math. Lett..

[15]  J. Butcher The Numerical Analysis of Ordinary Di erential Equa-tions , 1986 .

[16]  J. Butcher The numerical analysis of ordinary differential equations: Runge-Kutta and general linear methods , 1987 .

[17]  W. Marsden I and J , 2012 .

[18]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[19]  Robert A. Adams,et al.  Calculus: A Complete Course , 1994 .

[20]  Andrew S. Glassner,et al.  An introduction to ray tracing , 1989 .

[21]  Alicia Cordero,et al.  Construction of fourth-order optimal families of iterative methods and their dynamics , 2015, Appl. Math. Comput..