In this work we show a general procedure to obtain optimal derivative free iterative methods (Kung and Traub (1974))  for nonlinear equations f(x) = 0, applying polynomial interpolation to a generic optimal derivative free iterative method of lower order. Let us consider an optimal method of order q=2n-1, v=φn(x), that uses n functional evaluations. Performing a Newton step w=v-f(v)f'(v) one obtains a method of order 2n, that is not optimal because it introduces two new functional evaluations. Instead, we approximate the derivative by using a polynomial of degree n that interpolates n+1 already known functional values and keeps the order 2n. We have applied this idea to Steffensen's method, (Ortega and Rheinboldt (1970)) , obtaining a family of optimal derivative free iterative methods of arbitrary high order. We provide different numerical tests, that confirm the theoretical results and compare the new family with other well known family of similar characteristics. © 2012 Elsevier Ltd.
Cordero, A., Hueso, J. L., Martínez, E., & Torregrosa, J. R. (2013). Generating optimal derivative free iterative methods for nonlinear equations by using polynomial interpolation. Mathematical and Computer Modelling, 57(7–8), 1950–1956. https://doi.org/10.1016/j.mcm.2012.01.012