The dynamical behavior of different Steffensen-type methods is analyzed. We check the chaotic behaviors alongside the convergence radii (understood as the wideness of the basin of attraction) needed by Steffensen-type methods, that is, derivative-free iteration functions, to converge to a root and compare the results using different numerical tests. We will conclude that the convergence radii (and the stability) of Steffensen-type methods are improved by increasing the convergence order. The computer programming package MATHEMATICA provides a powerful but easy environment for all aspects of numerics. This paper puts on show one of the application of this computer algebra system in finding fixed points of iteration functions.
1. Introduction
Suppose that we wish to find a solution ξ of a nonlinear equation f(z)=0 numerically, where z∈ℂ and the function f:ψ⊆ℂ→ℂ is an analytical complex function. Starting from some z0, the Steffensen’s method [1] uses the iteration:
(1)zn+1=zn-f(zn)2f(zn+f(zn))-f(zn),n=0,1,2,….
If the initial value z0 is close enough to a simple root of f(z), this iteration converges quadratically (see [2, 3] for more information about Steffensen-type methods). If the root is not simple, the convergence becomes linear. Since then, a tremendous amount of effort has been made in the direction of improving the convergence and/or the simplicity of the method resulting in modified Steffensen method in the divided difference form [4] as follows (m is the multiplicity):
(2)zn+1=zn-mf(zn)f[zn,wn],n=0,1,2,…,
wherein wn=zn+f(zn) and f[zn,wn] is the two-point divided difference. Notice that the modification for multiple roots given by (2) is the same as that given by Schröder for Newton’s method. We remark that the notation of divided differences will be used throughout this paper.
Note that the Steffensen method could be written in a more generalized form with one free nonzero parameter [5] as follows:
(3)zn+1=zn-f(zn)f[zn,wn],n=0,1,2,…,
wherein wn=zn+βf(zn) and β∈ℂ∖{0}. It is known that this method has second order of convergence, for every nonzero value of the parameter β. However, as we will see in this paper, the parameter β plays an important role for choosing the initial estimation in the stability of the method, and so forth.
Usually, efficiency indices are used to compare the behavior of iterative methods. Traub in [5] used the operational efficiency index p1/op, where p is the convergence rate and op is the number of products/quotients per iteration. Ostrowski in [6] defined the classical efficiency index p1/d, wherein d is the number of functional evaluations per iteration. Similarly, one may use the informational index defined as p/d. In this work, we are going to introduce another criterion for comparing the iterative methods used: the basins of attraction, (see, e.g., [7–13]).
The dynamical analysis of the rational function associated with an iterative method gives us important information about its stability and reliability. There exists a huge number of publications related to the numerical and dynamical properties of iteration functions (see, e.g., [14–18]).
It is known that the basin of attraction of an iterative method on a particular function is an interesting tool to visually get to know how an iteration function behaves as a function of different initial points [19].
The remaining sections of this paper are organized as follows. Section 2 provides some basic concepts in order to deal with dynamics associated with iterative functions. It is followed by Section 3, wherein a discussion over the basins of attraction for the uniparametric Steffensen scheme (3) on higher degree polynomials, using the software Mathematica 8, is given. Section 4 gives the list of methods that we are going to compare and their dynamics on quadratic polynomials is analyzed. Next in Section 5, we present the basins of attraction for various high-order Steffensen-type methods. Some conclusions and discussions are illustrated in Section 6 to end the paper.
2. Basic Concepts
We need some basic definitions and notions before going to Section 3. Thus, now we remind them briefly.
Let R:ℂ^→ℂ^ be a rational map on the Riemann sphere. Then, a periodic point z0 of period m is such that Rm(z0)=z0, where m is the smallest such integer. If m=1, z0 is called fixed point of R and also point z0 is called attracting if |R′(z0)|<1, repelling if |R′(z0)|>1, and neutral if |R′(z0)|=1. If the derivative is zero then the point is called super-attracting [20]. If R is the rational function associated with an iterative method on a function f, the fixed points of R different from the roots of f(z)=0, are called strange fixed points.
Note that the Julia set of a nonlinear map R(z), denoted 𝒥(R), is the closure of the set of its repelling periodic points. The complementary of 𝒥(R) in the Riemann sphere is the Fatou set ℱ(R).
A point z0 is in the Julia set if and only if dynamics in a neighborhood of z0 shows strong dependence on the starting conditions, so that nearby initial conditions yield to wildly different behavior after a number of full iterations.
Definition 1 (basin of attraction [21]).
If a fixed point p of R is attracting, then all nearby points of p are attracted toward p under the action of R, in the sense that their iterates under R converge to p. The collection of all points whose iterates under R converge to p is called the basin of attraction of p, denoted by Bp={x∈ℂ:limk→∞Rk(x)=p} with Rk(x)=R(R(⋯R(x)⋯)) as the k-fold composite map of x under R.
Lemma 2 (see [20]).
Every attracting periodic orbit is contained in the Fatou set of R. In fact, the entire basin of attraction Bp, which is an open set, for an attracting periodic orbit is contained in the Fatou set. However, every repelling periodic orbit is contained in the Julia set.
Kalantari in [22] coined the term “polynomiography” to be the art and science of visualization in the approximation of roots of polynomial using iteration functions. Note that a polynomiograph may or may not result in a fractal image. Even when a polynomiograph is a fractal image it does not diminish its uniqueness.
As we endeavor to solve increasingly complex problems, computer algebra systems (CAS) are becoming more and more important to our work. We use CAS to help with calculations too time consuming. Mathematica 8 is one of the most popular CASs available today [23]. Due to its wide applicability and power along with easiness [24], we apply this programming package in finding the chaotic behaviors of Steffensen-type methods.
3. Basins of Attractions for the Uniparametric Steffensen Scheme
A wide dynamical study of Steffensen’s method (1) on quadratic and cubic polynomials has been developed in [12]. In this paper, it was showed that no Scaling Theorem is possible for this derivative-free scheme and, therefore, the analysis was made in particular polynomials. In general, it can be concluded that the stability of derivative-free schemes is worse than the one of methods with derivatives. In fact, in the dynamical planes associated with some of these methods, frequently the basin of attraction of the infinity appears. In fact, the infinity can be a strange fixed point: if R(z) is the rational function associated to the method, then we prove that the infinity is a fixed point if Q(0)=0, where Q(z)=1/R(1/z). If this happens, the infinity (as any other strange fixed point) can be attractive, repulsive, or neutral. This character is obtained analyzing the value of Q′(0).
Before providing the chaotic behaviors of multipoint Steffensen-type methods, we conclude some points on the uniparametric Steffensen scheme (3).
Lemma 3.
The associated operator of the uniparametric Steffensen scheme on z2-1 is OSt(x)=(1-βx+x2+βx3)/(-β+2x+βx2). The only strange fixed point of OSt(z) is the infinity, and its character is neutral.
This result can be understood as the behavior of the infinity as fixed point can be attractive or repulsive, and it can be observed also when this method is applied on other functions, as we will see in Figures 5 to 10.
Toward this end and heretofore, we take a rectangle D=[-3,3]×[-3,3]∈ℂ and we assign a color to each point z0∈D according to the simple root at which the corresponding iterative method starting from z0 converges, and we mark the point as black if the method does not converge after the maximum number of iterations. In this way, we distinguish the attraction basins by their colors for different methods.
The criteria we have used in our Mathematica 8 codes [25–27] are that the maximum number of iterations is 100. That is to say, if the method does not reach the considered accuracy after 100 of its full iteration steps, we allocate the black color and also considered it as Not Convergent. The considered accuracy is the obtained residual of the function to be less than 10-4. Note that the small white points will be shown in the exact location of the simple zeros in our fractal patterns.
In what follows, we define the test problems in this paper. For the first test, we have taken the following function with roots {-1.,-i,i,1.}:(4)p1(z)=z4-1.
The second test problem is a polynomial as follows when the simple zeros are {-2.20663,-0.45318-0.78493i,-0.45318+0.78493i,0.906359,1.10332-1.911i,1.10332+1.911i}:
(5)p2(z)=z6+10z3-8.
The third test problem is chosen as follows:
(6)p3(z)=z4-1z.
The roots are {0.309017+0.951057i,0.309017-0.951057I,1., -0.809017+0.587785i, −0.809017 − 0.587785i}.
The fourth test problem is taken into account as
(7)p4(z)=z4-z+i.
The roots are {-0.759845+0.592595i,-0.532605-1.08829i, 0.181924+0.732098I, 1.11052 − 0.236405i}.
The last test problem we have chosen is as follows when the roots are {-1.09112+0.629961i, 0.-1.25992i,1.09112 + 0.629961i}:
(8)p5(z)=z3-2i.
We can observe on Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10 that scheme (3) for the choice β=0.001 is very robust in contrast to the other cases. Method (3) has different radii of convergence according to the free nonzero parameter β. This is also true for all the test functions for Steffensen’s scheme (β=1).
Scheme (3) for the test problem 1.
β=1
β=0.5
β=0.1
Scheme (3) for the test problem 1.
β=-0.1
β=0.01
β=0.001
Scheme (3) for the test problem 2.
β=1
β=0.5
β=0.1
Scheme (3) for the test problem 2.
β=-0.1
β=0.01
β=0.001
Scheme (3) for the test problem 3.
β=1
β=0.5
β=0.1
Scheme (3) for the test problem 3.
β=-0.1
β=0.01
β=0.001
Scheme (3) for the test problem 4.
β=1
β=0.5
β=0.1
Scheme (3) for the test problem 4.
β=-0.1
β=0.01
β=0.001
Scheme (3) for the test problem 5.
β=1
β=0.5
β=0.1
Scheme (3) for the test problem 5.
β=-0.1
β=0.01
β=0.001
From the graphical comparisons in this section, it is obvious that in the Steffensen uniparametric scheme (3), the basins of attractions can be widely improved by choosing a small value for the nonzero free parameter β. In fact, from the Taylor series expansion of the divided difference and the derivative around the solution ξ of the nonlinear equation, we have
(9)f[xn,wn]=f(xn+βf(xn))-f(xn)βf(xn)=f′(ξ)[1+(2+βf′(ξ))c2ennnnnl..+(βf′(ξ)c22nnnnnnn+(3+3βf′(ξ)+β2f′(ξ)2)c3βf′(ξ)c22)ek2]+O(en3),f′(xn)=f′(ξ)[1+2f′(ξ)c2en+3f′(ξ)c3en2]+O(ek3).
Let us note that, when β→0, both Taylor expansions coincide. Note that it is experimentally observable that (almost always) when the forward finite difference is used at the denominator of Steffensen’s scheme, that is, f[xn,wn]=(f(xn+βf(xn))-f(xn))/βf(xn), then very small negative values for β results in larger basins of attractions and higher speed, while when the backward finite difference is used at the denominator, that is, f[xn,wn]=(f(xn)-f(xn-βf(xn)))/βf(xn), then very small positive values for β yield in larger basins of attractions and higher speed.
This finding is very useful when we extend Steffensen’s scheme for solving systems of nonlinear equations by defining the first order divided difference operator properly. In fact, by choosing very small values for the free nonzero parameter the convergence radii will be similar to Newton’s method for solving systems of nonlinear equations, while there is no need to compute the Jacobian matrix which is costly.
Even by applying an iteration on β and make the process with memory, we can obtain the cubical order of convergence. Another interesting point which should be included is the computational time required for obtaining the fractal behaviors of such schemes. In fact, in our numerical experiments, the computational CPU time has a dramatically fall when β is chosen as a very small value while choosing β=1 takes more time to find the fixed points of the iteration functions.
Remark 4.
The chaotic behavior of Steffensen’s method can be simply improved by choosing very small entries for the free nonzero parameter (see Figures 6(c), 8(c), and 10(c)). In this case, the fractal behavior of the scheme tends to Newton’s fractal. In fact, we want the methods to have very few black points on all examples not just one.
4. Behavior of Several Steffensen-Like Methods on Quadratic Polynomials
Many iterative methods have been improved by using various techniques [28, 29]. A drawback of Newton’s method is that for many particular choices of the function f, especially in hard problems, the calculations of the derivatives take some deal of time. That is why higher order derivative-free methods are better root solvers and are in focus recently. The most important merit of Steffensen’s method is that it has quadratic convergence like Newton’s method. For this reason, in this section, we list some of the multipoint derivative-free methods we consider for comparisons and a brief dynamical analysis is made for these methods applied on quadratic polynomials.
Kung and Traub in the pioneer paper [30] provided the following two- and three-step derivative-free families (β∈ℂ∖{0}) of methods with orders four and eight, respectively,
(10)yn=xn+βf(xn),zn=yn-βf(xn)f(yn)f(yn)-f(xn),xn+1=zn-f(xn)f(yn)f(zn)-f(xn)[1f[yn,xn]-1f[zn,yn]],(11)yn=xn+βf(xn),zn=yn-βf(xn)f(yn)f(yn)-f(xn),wn=zn-f(xn)f(yn)f(zn)-f(xn)[1f[yn,xn]-1f[zn,yn]],xn+1=wn-f(xn)f(yn)f(zn)f(wn)-f(xn)×[1f(wn)-f(yn){1f[wn,zn]-1f[zn,yn]}nnn-1f(zn)-f(xn){1f[zn,yn]-1f[yn,xn]}].
We would like to study the general convergence of methods (10) and (11) for quadratic polynomials. To be more precise (see [31, 32]), a given method is generally convergent if the scheme converges to a root for almost every starting point and for almost every polynomial of a given degree. The main problem is that no Scaling Theorem can be stated for derivative-free methods (see, e.g., [33]). So, only the behavior on specific polynomials can be analyzed. In this paper, the polynomial p(z)=z2-1 will be used for this purpose. Therefore, let us denote by OKT4(z) (OKT8(z)) the operator associated with the fourth-order (resp., eighth-order) scheme by Kung and Traub on p(z).
The fixed points of OKT4(z) are the roots of the equation OKT4(z)=z, that is, z=-1, z=1, and the strange fixed points that are the zeros of the polynomial 1+β2+β4-4β(1+2β2)z + (2+21β2-4β4)z2-24β(1-β2)z3+(13-45β2+6β4)z4+(28β-24β3)z5+(23β2-4β4)z6 + 8β3z7 + β4z8.
The stability of the strange fixed points can be deduced from graphical analysis of the respective stability function, that is, representing graphically the regions of the complex plane in which the absolute value of the derivative of the operator evaluated at the strange fixed point is lower than one.
Lemma 5.
The number of simple strange fixed points of operator OKT4(z) is eight, and their stability is described in the following cases.
Four of them are always repulsive, so they remain in the Julia set.
Two of them can be attractive (but not super attractive) in a small complex neighborhood of the origin.
Finally, two of the strange fixed points are attractive in complex region around the origin and are super attractive for the following values of the parameter: β≈-0.406175-1.63893i, β≈-0.406175+1.63893i, β≈0.276582-0.409103i, and β≈0.276582+0.409103i.
The stability region of the complex plane where these four fixed points are attractive is represented in Figure 11.
Stability function of attractive strange fixed points of KT4.
In order to determine the critical points, we calculate the first derivative of OKT4(z). A classical result establishes that there is at least one critical point associated with each invariant Fatou component. It is clear that z=-1 and z=1 are critical points and give rise to their respective Fatou components, but there exist in the family, some free critical points, that is, critical points different from the roots, some of them depending on the value of the parameter. If any of these free critical points is near a stable strange fixed point, then the last one would have its own basin of attraction.
Lemma 6.
Analyzing the equation OKT4′(z)=0, we obtain fourteen free critical points, roots of the polynomial:
(12)2+7β2-3β4-β6-β8+(-28β+26β3+10β5+16β7)z+(26-96β2-25β4-102β6+7β8)z2+(140β-24β3+350β5-96β7)z3+(-58+166β2-750β4+525β6-21β8)z4+(-164β+1068β3-1500β5+240β7)z5+(30-904β2+2430β4-1060β6+35β8)z6+(308β-2168β3+2300β5-320β7)z7+(827β2-2495β4+1065β6-35β8)z8+(1098β3-1550β5+240β7)z9+(843β4-534β6+21β8)z10+(390β5-96β7)z11+(107β6-7β8)z12+16β7z13+β8z14.
Moreover, we can state the following facts.
There is no value of β that makes coincide a free critical point with a strange fixed point.
For β≈0.05 one free critical point coincides with the root z=-1, so in this case the number of free critical points is reduced to 13.
The real values of the critical points for a real range of values of the parameter are showed in Figure 12.
Real values of free critical points of OKT4(z).
Let us state that the strange fixed points of OKT8(z) are the zeros of a polynomial of degree 39, whose coefficients depend on β.
The stability of the strange fixed points is described in the following lemma.
Lemma 7.
The number of simple strange fixed points of operator OKT8(z) is thirty-nine, and their stability is described in the following cases.
Three of them are always repulsive, so they remain in the Julia set.
The other strange fixed points can be attractive or super attractive, in different complex regions.
The stability region of some of these fixed points is represented in Figure 13.
Stability regions of strange fixed points from KT8.
By analyzing these stability functions, it is deduced that there are wide regions of the complex plane where it is easy to find two or more attractive strange fixed points (as, e.g., β=-5i and β=5i). Moreover, there exist also wide regions of stability near the origin. In fact, no attractive strange fixed points can be found for |β|<3/2.
Analyzing the equation OKT8′(z)=0, we obtain sixty free critical points that will coincide with some of the strange fixed points in case these are super attracting.
An optimal fourth-order derivative-free method [34] was introduced by Liu et al. in the following form:
(13)yn=xn-f(xn)2f(xn+f(xn))-f(xn),xn+1=yn-f[xn,yn]-f[yn,zn]+f[xn,zn]f[xn,yn]2f(yn),
where zn=xn+f(xn). It is denoted by L4. The rational function associated with this method on p(z) is the operator OL4(z),
(14)OL4(z)=(1-5z+19z2-27z3-z4+23z5nnl.+5z6+7z7+8z8+2z9)×((-1+2z+z2)(1-2z+3z2+2z3)2)-1.
Lemma 8.
The quantity of simple strange fixed points of operator OL4(z) is seven, {-3.4251, -1.98276-0.105743i, -1.98276+0.105743i, 0.109953-0.282276i, 0.109953+0.282276i, 0.585354-0.246666i, 0.585354+0.246666i}, and their stability is described in the following cases.
All the complex strange fixed points are repulsive, so they remain in the Julia set.
The real strange fixed point is an attractor, but not super attractor.
Lemma 9.
The equation OL4′(z)=0 yields the polynomial 1+5x-31x2+23x3+64x4+30x5+4x6 whose roots are the free critical points, {-3.94434,-2.25291, -1.85186,-0.114261, 0.331686-0.151469i, 0.331686+0.151469i}.
Let us remark that the existence of a free critical point near the real attracting fixed point will derive its own basin of attraction, as can be seen in Figure 14.
Dynamical plane of OL4(z) with three basins of attraction.
Cordero et al. in [35] proposed a seventh-order method, which is free from derivative:
(15)yn=xn-f(xn)2f(wn)-f(xn),zn=yn-f(yn)×((f(yn)-f(wn))(yn-wn)-f(yn)(yn-wn))-1,xn+1=zn-f(zn)×((f(zn)-f(yn))(zn-yn)-f(zn)(yn-zn)nnnn-(f(yn)-f(wn))(yn-wn))-1
wherein wn=xn+f(xn). We will denote it by C7.
The rational function associated C7 on p(z) is OC7(z),
(16)OC7(z)=(1-9z+13z2+65z3-112z4-64z5+119z6n-22z7-98z8+71z9+152z10+85z11+21z12n+2z13)×((-1+2z+z2)3+23z4+12z5+6z6)nnn×(1-4z+14z2+8z3(-1+2z+z2)3(1-4z+14z2+8z3+23z4+12z5+6z6))-1.
Lemma 10.
The simple strange fixed points of operator OC7(z) are
{-2.30882-0.732518i, -2.30882+0.732518i, -2.05354-0.276071i, -2.05354+0.276071i, -1.87555, -0.411291, 0.200847-0.0392285i, 0.200847+0.0392285i, 0.638264-0.0474315i, 0.638264+0.0474315i}, and all of them are finite and repulsive as the derivable operator evaluated in each of them is greater than one, in absolute value;
the infinity is also a strange fixed point and its character is neutral.
Lemma 11.
The free critical points of OC7(z) are the roots of the polynomial 1-9z+32z2+96z3-399z4-2323z5+1676z6+9274z7+2233z8-13025z9-12426z10+1194z11+8849z12+6817z13+2694z14+612z15+76z16+4z17.
As all the strange fixed points are repulsive and, therefore, they lay on the Julia set, the position of the free critical points in the complex plane has no interest.
Zheng et al. in [36] presented the following eighth-order derivative-free family without memory:
(17)yn=xn-f(xn)f[xn,wn],wn=xn+βf(xn),nnnnnnnnnnnnnnnnnnnnnnllllllβ∈ℂ∖{0},zn=yn-f(yn)f[xn,yn]+f[yn,wn]-f[xn,wn],xn+1=zn-f(zn)×(f[zn,yn]+f[zn,yn,xn](zn-yn)nnl+f[zn,yn,xn,wn](zn-yn)(zn-xn))-1,
which will be denoted by Z8.
Lemma 12.
Operator OZ8(z) has ten simple strange fixed points, roots of the equation 1 + 6β2 + β4 + (-24β-16β3)z + (21+48β2-β4)z2 + (-56β+32β3)z3 + (35-84β2-6β4)z4 + 56βz5 + (7+14β4)z6 + (24β-32β3)z7 + (30β2-11β4)z8 + 16β3z9 + 3β4z10. Nine of these fixed points are repulsive for all values of parameter β; however, one of them can be attractive (even super attractive) in a small complex region around the origin. In Figure 15, the stability function of these specific fixed point is shown.
Stability function of the only attractive strange fixed point of Z8.
Again, the existence of values of the parameter that yields attracting strange fixed points forces us to analyze the possibility of free critical points. As we know, if both elements coexist, basins of attraction of fixed points different from the roots appear.
Lemma 13.
Analyzing the equation OZ8′(z)=0, we obtain the free critical points:
(18)cr1(β)=-1-ββ,cr2(β)=-1+ββ,cr3(β)=-2β-2β2+β4β2,cr4(β)=-2β+2β2+β4β2.
Moreover,
if β=1/2, then cr2(1/2)=cr4(1/2)=-1 and there are only two free critical points;
also when β=-1/2, then cr1(1/2)=cr3(1/2)=1 and there are only two free critical points;
finally, if β=-2i or β=2i, then cr3(±2i)=cr4(±2i)=±2i and the number of free critical points is three.
Soleymani et al. in [37] proposed the following optimal three-step iteration family, including four function evaluations, just like (11) and (17):
(19)yn=xn-f(xn)f[xn,wn],wn=xn+βf(xn),nnnnnnnnnnnnnnnnnnnnnllllllllβ∈ℂ∖{0},zn=yn-f(yn)f[xn,yn]+f[yn,wn]-f[xn,wn],xn+1=zn-(1+a3(zn-xn))2f(zn)×(+a2a3(zn-xn)2a1-f(xn)a3+2a2(zn-xn)nnn+a2a3(zn-xn)2)-1,
where a1=f[wn,xn]+f(wn)a3-(wn-xn)a2, a2=f[wn,xn,yn]+a3f[wn,yn], and
(20)a3=(wn(f[zn,xn]-f[yn,xn])-f[zn,xn]ynn+f[wn,xn](yn-zn)+f[yn,xn]zn)×((zn-yn)f(wn)+(wn-zn)f(yn)nnl+(yn-wn)f(zn))-1.
Let us denote this method by S8; the rational function associated with S8 when it is applied on the quadratic polynomial p(z), OS8(z), is the same as the one of Z8. Then, their dynamics is the same (for this polynomial). Nevertheless, it is very different for other functions, as we will see in the following sections.
5. Attraction Basins for Various Steffensen-Type Methods
The aim herein is to use the basin of attraction as another tool for comparing the iteration algorithms given in Section 4.
We have used methods (10) with β=1, (10) with β=-0.001, (11) with β=1, (11) with β=-0.001, (13), (15), and (17) with β=1, (17) with β=-0.001, (19) with β=1, and (19) with β=-0.001, for the test problems of Section 4. The fractal behavior of these comparisons is furnished in Figures 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, and 27. Let us remark that the dynamics of Test 3 is similar to the other test functions, so it is not included.
Basins of attraction for the test problem 1.
Method (10), β=1
Method (10), β=-0.001
Method (11), β=1
Basins of attraction of (11), (13), and (15) for the test problem 1.
Method (11), β=-0.001
Method (13), β=-0.001
Method (17), β=1
Basins of attraction for the test problem 1.
Method (17), β=-0.001
Method (19), β=1
Method (19), β=-0.001
Basins of attraction for the test problem 2.
Method (10), β=1
Method (10), β=-0.001
Method (11), β=1
Basins of attraction of (11), (13), and (15) for the test problem 2.
Method (11), β=-0.001
Method (13), β=-0.001
Method (17), β=1
Basins of attraction for the test problem 2.
Method (17), β=-0.001
Method (19), β=1
Method (19), β=-0.001
Basins of attraction for the test problem 4.
Method (10), β=1
Method (10), β=-0.001
Method (11), β=1
Basins of attraction of (11), (13), and (15) for the test problem 4.
Method (11), β=-0.001
Method (13), β=-0.001
Method (17), β=1
Basins of attraction for the test problem 4.
Method (17), β=-0.001
Method (19), β=1
Method (19), β=-0.001
Basins of attraction for the test problem 5.
Method (10), β=1
Method (10), β=-0.001
Method (11), β=1
Basins of attraction of (11), (13), and (15) for the test problem 5.
Method (11), β=-0.001
Method (13), β=-0.001
Method (17), β=1
Basins of attraction for the test problem 5.
Method (17), β=-0.001
Method (19), β=1
Method (19), β=-0.001
Remark 14.
According to the discussion at the end of Section 4, the Steffensen-type methods, in which there is no nonzero free parameter in their structures, are not competitive and we expect to have small basins of attraction for them.
It is clear for the compared tests that even the higher order Steffensen-type methods (13) and (15) without the free nonzero parameter have improved basins in contrast to Steffensen’s scheme (β=1). Furthermore, it should be noted that among all the methods compared in this section, Soleymani et al. optimal eighth-order method (19) with β=-0.001 has the best performance, followed by (17) with β=-0.001. In this work, the computer specifications are Microsoft Windows XP Intel(R), Pentium(R) 4 CPU, 3.20 GHz with 4 GB of RAM.
Remark 15.
Unlike the Newton-type methods in which whatever the order is higher, the convergence radius is smaller, in the multipoint high-order efficient Steffensen-type methods the increase of convergence order will automatically improve the convergence radius, though the chaotic behavior of the schemes for unappropriate values of the free nonzero parameter is too much. To avoid this chaotic behavior, one may follow Remark 4.
In order to summarize these results, we have attached a weight to the quality of the results obtained by each method. The weight of 1 is for the smallest Julia set and a weight of 4 for scheme with chaotic behavior alongside the convergence behavior. We then averaged those results to come up with the smallest value for the best method overall and the highest for the worst. These data are presented in Table 1.
Results of chaotic comparisons for different derivative-free methods.
Method
Test 1
Test 2
Test 4
Test 5
Average
(10) with β=1
3
4
3
3
13/4
(10) with β=-0.001
2
3
2
2
9/4
(11) with β=1
3
4
3
3
13/4
(11) with β=-0.001
2
3
2
2
9/4
(13)
4
4
4
3
15/4
(15)
4
4
4
3
15/4
(17) with β=1
3
3
3
3
12/4
(17) with β=-0.001
2
2
1
2
7/4
(19) with β=1
2
2
3
3
10/4
(19) with β=-0.001
1
1
1
1
4/4
Notice again that the figures show how fast the method converges to a root based on shading to indicate speed of convergence.
6. Conclusions
In this paper, we have analyzed the dynamics of different Steffensen-type methods, firstly on quadratic polynomials and afterwards on other functions. We have concluded that in Steffensen-type methods whatever the order is higher, the convergence radius will be bigger. We have used Mathematica 8 for finding the fixed and critical points of the rational functions associated with the iterations functions and for drawing the basins of attraction. Besides, if the free nonzero parameter for the families analyzed tends to 0, then its fractal tends to be the same as Newton’s fractal. Choosing very small magnitudes for the free nonzero parameter gives us the ability to avoid computation of the Jacobian matrix when dealing with systems of nonlinear equations and have an acceptable convergence radius. Although we have discussed simple zeros of nonlinear functions, such remarks are valid for Steffensen-type methods when finding multiple roots as well.
Conflict of Interests
The authors declare that they do not have any conflict of interests in their submitted paper.
Acknowledgments
The authors are indebted to the referees for some interesting comments and suggestions. This research was supported by Ministerio de Ciencia y Tecnología MTM2011-28636-C02-02.
SteffensenJ. F.Remarks on iteration1933166472SoleymaniF.On a bi-parametric class of optimal eighth-order derivative-free methods20117212737MR2908187ZBL1248.65050SoleymaniF.Optimal fourth-order iterative methods free from derivatives2011122255264MR2880423ZBL1265.65101ZhengQ.ZhaoP.ZhangL.MaW.Variants of Steffensen-secant method and applications2010216123486349610.1016/j.amc.2010.04.058MR2661705ZBL1200.65036TraubJ. F.1964Englewood Cliffs, NJ, USAPrentice-Hallxviii+310Prentice-Hall Series in Automatic ComputationMR0169356OstrowskiA. M.1964Englewood Cliffs, NJ, USAPrentice-HallNetaB.ScottM.ChunC.Basins of attraction for several methods to find simple roots of nonlinear equations201221821105481055610.1016/j.amc.2012.04.017MR2927071ZBL1250.65067NetaB.ScottM.On a family of Halley-like methods to find simple roots of nonlinear equations2013219157940794410.1016/j.amc.2013.02.035MR3037505NetaB.ChunC.On a family of Laguerre methods to find multiple roots of nonlinear equations201321923109871100410.1016/j.amc.2013.05.002MR3069009NetaB.ChunC.ScottM.Basins of attraction for optimal eighth order methods to find simple roots of nonlinear equations201422756759210.1016/j.amc.2013.11.017MR3146342AmatS.BusquierS.PlazaS.Dynamics of the King and Jarratt iterations200569321222310.1007/s00010-004-2733-yMR2139283ZBL1068.30019ChicharroF.CorderoA.GutiérrezJ. M.TorregrosaJ. R.Complex dynamics of derivative-free methods for nonlinear equations2013219127023703510.1016/j.amc.2012.12.075MR3027866CorderoA.García-MaimóJ.TorregrosaJ. R.VassilevaM. P.VindelP.Chaos in King's iterative family201326884284810.1016/j.aml.2013.03.012MR3066701ChunC.LeeM. Y.NetaB.DžunićJ.On optimal fourth-order iterative methods free from second derivative and their dynamics2012218116427643810.1016/j.amc.2011.12.013MR2879123ZBL1277.65031CorderoA.TorregrosaJ. R.VindelP.Dynamics of a family of Chebyshev-Halley type methods2013219168568858310.1016/j.amc.2013.02.042MR3044065SoleimaniF.SoleymaniF.ShateyiS.Some iterative methods free from derivatives and their basins of attraction for nonlinear equations201320131010.1155/2013/301718301718MR3044881ZBL1264.65076SusantoH.KarjantoN.Newton's method's basins of attraction revisited200921531084109010.1016/j.amc.2009.06.041MR2568964ZBL1264.65076VrscayE. R.GilbertW. J.Extraneous fixed points, basin boundaries and chaotic dynamics for Schröder and König rational iteration functions198852111610.1007/BF01401018MR918313ZBL1264.65076OttE.Basin of attraction2006171701BlanchardP.Complex analytic dynamics on the Riemann sphere19841118514110.1090/S0273-0979-1984-15240-6MR741725ZBL1264.65076GulickD.1992Prentice-HallKalantariB.2009Hackensack, NJ, USAWorld Scientific Publishingxxiv+467MR2483211WagonS.20103rdSpringer
http://www.wolfram.com/mathematica/new-in-8/new-and-improved-core-algorithms/index.htmlGetzC.HelmstedtJ.2004Amsterdam, The NetherlandsElsevier B. V.xii+332MR2162550McClureM.Newton's method for complex polynomialspreprint version of a “Mathematical graphics” column from Mathematica in Education and Research, 1–15VaronaJ. L.Graphic and numerical comparison between iterative methods2002241374610.1007/BF03025310MR1889920ZBL1003.65046IgnatovaB.KyurkchievN.IlievA.Multipoint algorithms arising from optimal in the sense of Kung-Traub iterative procedures for numerical solution of nonlinear equations201164579SoleymaniF.SoleimaniF.Novel computational derivative-free methods for simple roots2012131247258MR2953688ZBL1272.65043KungH. T.TraubJ. F.Optimal order of one-point and multipoint iteration197421643651MR0353657ZBL0289.65023McMullenC.Families of rational maps and iterative root-finding algorithms1987125346749310.2307/1971408MR890160ZBL0634.30028SmaleS.On the efficiency of algorithms of analysis19851328712110.1090/S0273-0979-1985-15391-1MR799791ZBL0592.65032AmatS.BusquierS.PlazaS.Review of some iterative root-finding methods from a dynamical point of view200410335MR2127479ZBL1137.37316LiuZ.ZhengQ.ZhaoP.A variant of Steffensen's method of fourth-order convergence and its applications201021671978198310.1016/j.amc.2010.03.028MR2647066ZBL1208.65064CorderoA.HuesoJ. L.MartínezE.TorregrosaJ. R.A family of derivative-free methods with high order of convergence and its application to nonsmooth equations201220121583690110.1155/2012/836901MR2922950ZBL1246.65079ZhengQ.LiJ.HuangF.An optimal Steffensen-type family for solving nonlinear equations2011217239592959710.1016/j.amc.2011.04.035MR2811234ZBL1227.65044SoleymaniF.Karimi VananiS.PaghalehM. J.A class of three-step derivative-free root solvers with optimal convergence order201220121556874010.1155/2012/568740MR2889100ZBL1235.65047