# Study of a Biparametric Family of Iterative Methods

Authors: B. Campos; A. Cordero; Á. A. Magreñán; J. R. Torregrosa; P. Vindel
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

The dynamics of a biparametric family for solving nonlinear equations is studied on quadratic polynomials. This biparametric family includes thec-iterative methods and the well-known Chebyshev-Halley family. We find the analytical expressions for the fixed and critical points by solving 6-degree polynomials. We use the free critical points to get the parameter planes and, by observing them, we specify some values of (
α
,
c
) with clear stable and unstable behaviors.


---

## Body

## 1. Introduction

Iterative methods are needed for solving the most of the nonlinear equationsf
x
=
0, because they are difficult or impossible to solve exactly by means of analytical methods. When they are applied on polynomials, they give rise to rational functions whose behavior is not well known except in small areas. A way to extend these regions is by studying the dynamical behavior of the rational functions associated to the iterative methods. In some previous papers, we focus on this option and we start with the dynamical study of Chebyshev-Halley family ([1–3]), the King’s family [4], and the c-family [5].The(
α
,
c
)-family is a biparametric family of iterative methods that includes the Chebyshev-Halley ([6, 7]) and c-families when one of the parameters is fixed.In [8] the authors establish the conjugacy classes for (
α
,
c
)-family and they find explicit formulations for fixed and critical points of the Chebyshev, Halley, and super Halley methods and also for the c-family, applied on an arbitrary polynomial.In this paper, we consider the(
α
,
c
)-family and we calculate the fixed and critical points of this family applied on quadratic polynomials. By solving 6-degree polynomials we obtain the analytical expressions for the fixed and critical points of this family in terms of both parameters.From these analytical expressions we make a study of the points in the(
α
,
c
)-plane for α and c real. The analysis of the critical points helps us to obtain the parameter planes. From these parameter planes we specify some values of (
α
,
c
) with clear stable and unstable behavior and show the dynamical planes of the corresponding iterative methods.The(
α
,
c
)-family (see [9]) is a two-parametric third-order class of iterative methods defined by
(1)
z
n
+
1
=
z
n
-
1
+
1
2
L
f
z
n
1
-
α
L
f
z
n
+
c
L
f
z
n
2
f
z
n
f
′
z
n
,
where
(2)
L
f
z
=
f
z
f
′
′
z
f
′
z
2
,
and parameters α and c are complex. This family includes Chebyshev-Halley family for c
=
0 and c-family when α
=
0. As it is seen throughout this paper, the dynamical behavior of this family is much more complicated than previously studied because it includes two parameters α and c.For beginning the study of the dynamics of this family, we apply the operator associated to (1) on quadratic polynomials. Taking into account that this family satisfies the scaling theorem (see [9]), the rational operators obtained when applying the method on every quadratic polynomial are conjugated to the one obtained on p
z
=
z
2
+
a. For this polynomial, the operator M
p
z
,
α
,
c
,
a associated to (1) is a rational function depending on three complex parameters, a, α, and c:
(3)
M
z
,
α
,
c
,
a
=
-
α
c
a
4
+
2
-
1
+
2
α
c
a
3
z
2
+
2
-
1
+
2
α
-
3
c
+
3
α
c
a
2
z
4
+
2
-
6
-
3
c
+
2
α
c
a
z
6
+
6
-
4
α
-
2
c
+
α
c
z
8
×
8
z
5
α
a
-
2
z
2
+
α
z
2
-
1
.Rational functionsR and g
∘
R
∘
g
-
1 are conjugated, g being a Möebius map. Considering g
(
z
)
=
(
z
+
i
a
)
/
(
z
-
i
a
) the parameter a can be obviated and the operator M
z
,
α
,
c
,
a is conjugated to
(4)
O
p
(
z
,
α
,
c
)
=
z
3
1
+
z
4
-
2
+
2
α
-
z
+
4
c
1
+
z
2
-
2
α
+
z
1
+
z
4
2
α
z
-
1
-
2
z
+
4
c
z
3
1
+
z
2
-
8
α
c
z
4
.
By using this Möebius map, the roots of p
z become 0 and ∞; meanwhile g
(
∞
)
=
1.Iterative methods are used for finding the roots of a nonlinear equation and, from a dynamical point of view, these roots are included in the set of the fixed points of the operator associated with the method. So, from this point of view, our greatest interest is to find the fixed points of operator (4) and to study their stability.The natural space to study the dynamics of a rational mapR is the Riemann sphere C
^, R
:
C
^
→
C
^ (see [10, 11], e.g.). The dynamics of this rational map induces a subdivision of the complex sphere into two sets, namely, the Fatou set F
R and the Julia set J
R of R. The Fatou set of a nonconstant rational function R is the maximal open subset of C
^ on which the family of iterates {
R
n
} is equicontinuous (roughly speaking, they preserve the proximity of points) and the Julia set is its complement in C
^. From this definition, we know that F
R is open and J
R is compact. Moreover, if S
=
g
∘
R
∘
g
-
1, where g is a Möebius map, then F
S
=
g
F
R and J
S
=
g
J
R ([11, Theorem 3.1.4]).Let us recall some dynamical concepts of complex dynamics that we use in this paper (see [10]). Given a rational function R
:
C
^
→
C
^, theorbit of a point
z
0
∈
C
^ is defined as
(5)
{
z
0
,
R
z
0
,
R
2
z
0
,
…
,
R
n
z
0
,
…
}
.
We analyze the phase plane of the map R by classifying the starting points from the asymptotic behavior of their orbits. A point z
0
∈
C
^ is called afixed point if R
z
0
=
z
0. Aperiodic point
z
0 of period p
>
1 is a point such that R
p
z
0
=
z
0 and R
k
z
0
≠
z
0, for k
<
p. Apreperiodic point is a point z
0 that is not periodic, but there exists a k
>
0 such that R
k
z
0 is periodic. A point z
0 is acritical point of the rational map R if R fails to be injective in any neighborhood of z
0. Moreover, a fixed point z
0 is calledattractor if R
′
z
0
<
1,superattractor if R
′
z
0
=
0,repulsor if R
′
z
0
>
1, andparabolic if R
′
z
0
=
1. The fixed points different from the roots of the polynomial p
(
z
) are calledstrange fixed points.The basin of attraction of an attractor α is defined as
(6)
A
α
=
{
z
0
∈
C
^
:
R
n
z
0
⟶
α
,
n
⟶
∞
}
.The Fatou set of the rational functionR is the set of points z
∈
C
^ whose orbits tend to an attractor (fixed point, periodic orbit, or infinity). Since its complement in C
^ is the Julia set the basin of attraction of any fixed point belongs to the Fatou set and the boundaries of these basins of attraction belong to the Julia set.The rest of the paper is organized as follows. In Section2 we make a preliminary study of the fixed points, obtaining their analytic expressions and presenting the regions of the (
α
,
c
)-plane where they are complex or real. Section 3 is devoted to the analysis of the critical points, which will be useful in Section 5 to obtain the parameter planes. In Section 4 we check that the obtained results are consistent with the known ones of Chebyshev-Halley class and c-family. Section 6 shows some dynamical planes where interesting qualitative behaviors appear. Finally, in Section 7, some conclusions are stated.


## 2. Study of the Fixed Points

The fixed points of the operator (4) are obtained by solving the equation O
p
z
,
α
,
c
=
z. Evaluating O
p
z
,
α
,
c
-
z we have to solve the equation
(7)
-
z
(
z
-
1
)
P
z
,
α
,
c
1
+
z
4
(
2
α
z
-
1
-
2
z
)
+
4
c
z
3
1
+
z
2
-
8
α
c
z
4
=
0
,
where P
z
,
a
,
c is the symmetric 6-degree polynomial:
(8)
P
z
,
a
,
c
=
1
+
7
-
2
α
z
+
19
-
8
α
+
4
c
z
2
+
26
-
12
α
+
8
c
-
8
α
c
z
3
+
19
-
8
α
+
4
c
z
4
+
7
-
2
α
z
5
+
z
6
.The solutionsz
=
0 and z
=
∞ of (7) are associated to the roots of p
(
z
). We observe that z
=
1 is a strange fixed point for every value of the parameters. The other fixed points come from the roots of P
z
,
a
,
c. These roots must include the strange fixed points of Chebyshev-Halley family for c
=
0 and the respective ones of c-family for α
=
0, as we will see in Section 4. In the following result we present the explicit expressions of the strange fixed points. An exhaustive proof can be found in [12].Theorem 1.
The fixed points of the operator (4), associated to the biparametric family of iterative methods (1) on quadratic polynomials, are:(i)
z
=
0 and z
=
∞, corresponding to the roots of the polynomial;
(ii)
seven strange fixed points:z
=
1 and the following six points:
(9)
z
1
α
,
c
=
x
1
α
,
c
+
x
1
α
,
c
2
-
4
2
,
z
2
α
,
c
=
x
1
α
,
c
-
x
1
α
,
c
2
-
4
2
,
z
3
α
,
c
=
x
2
α
,
c
+
x
2
α
,
c
2
-
4
2
,
z
4
α
,
c
=
x
2
α
,
c
-
x
2
α
,
c
2
-
4
2
,
z
5
α
,
c
=
x
3
α
,
c
+
x
3
α
,
c
2
-
4
2
,
z
6
α
,
c
=
x
3
α
,
c
-
x
3
α
,
c
2
-
4
2
,where
(10)
x
1
α
,
c
=
1
3
2
α
-
7
+
1
3
f
α
,
c
3
+
g
α
,
c
3
,
x
2
α
,
c
=
1
3
2
α
-
7
-
1
6
f
α
,
c
3
+
g
α
,
c
3
-
i
3
6
f
α
,
c
3
-
g
α
,
c
3
,
x
3
α
,
c
=
1
3
2
α
-
7
-
1
6
f
α
,
c
3
+
g
α
,
c
3
+
i
3
6
f
α
,
c
3
-
g
α
,
c
3
,
f
α
,
c
=
-
1
+
2
α
3
+
18
1
+
4
α
c
i
+
6
3
c
2
α
-
1
+
2
α
3
+
c
-
1
+
40
α
+
32
α
2
+
16
c
,
g
α
,
c
=
-
1
+
2
α
3
+
18
1
+
4
α
c
i
-
6
3
c
2
α
-
1
+
2
α
3
+
c
-
1
+
40
α
+
32
α
2
+
16
c
.Proof.
The fixed points satisfyO
p
z
,
α
,
c
=
z. As we have seen, the numerator of O
p
z
,
α
,
c
-
z can be factorized as
(11)
-
z
z
-
1
P
z
,
a
,
c
,
where P
z
,
a
,
c is the polynomial given in (8).
The fixed pointsz
=
0 and z
=
1 are trivially obtained. Now, we look for the roots of the 6-degree polynomial (8). The change of variable x
=
z
+
1
/
z transforms polynomial (8) into the cubic one:
(12)
x
3
+
7
-
2
α
x
2
+
16
-
8
α
+
4
c
x
+
12
-
8
α
+
8
c
-
8
α
c
,
whose roots are
(13)
x
1
α
,
c
=
1
3
-
7
+
2
α
+
1
3
f
α
,
c
3
+
g
α
,
c
3
,
x
2
α
,
c
=
1
3
-
7
+
2
α
-
1
6
f
α
,
c
3
+
g
α
,
c
3
-
i
3
6
f
α
,
c
3
-
g
α
,
c
3
,
x
3
α
,
c
=
1
3
-
7
+
2
α
-
1
6
f
α
,
c
3
+
g
α
,
c
3
+
i
3
6
f
α
,
c
3
-
g
α
,
c
3
,
where
(14)
f
α
,
c
=
-
1
+
2
α
3
+
18
1
+
4
α
c
-
6
3
Δ
,
g
α
,
c
=
-
1
+
2
α
3
+
18
1
+
4
α
c
+
6
3
Δ
,
Δ
=
c
2
α
-
1
+
2
α
3
+
c
-
1
+
40
α
+
32
α
2
+
16
c
.
Undoing the change of variable, the six roots of the polynomialP
(
z
,
α
,
c
) are obtained by solving the quadratic equation:
(15)
z
2
-
z
x
+
1
=
0
,
whose solutions are
(16)
z
±
=
x
±
x
2
-
4
2
.
That is, the six roots of (8) are obtained by substituting in (16) the values of x given in (13).From Cardano’s method we know that,α and c being real:(i)
ifΔ
>
0 there is a unique real solution and a pair of complex conjugated roots;
(ii)
ifΔ
=
0 there are two different real solutions, one of them double; and there is one real solution with multiplicity three (x
=
(
-
7
+
2
α
)
/
3) when the following system of equations is fulfilled:
(17)
Δ
=
0
,
-
1
+
2
α
3
+
18
1
+
4
α
c
=
0
;
the only values of the parameters that satisfy this system are(18)
α
=
1
2
,
c
=
0
,
a
=
-
1
16
,
c
=
27
256
;
(iii)
ifΔ
<
0 there are three different real solutions.So,(
α
,
c
)-plane is divided into two regions, as we see in Figure 1. In this figure, the orange region corresponds to Δ
<
0 and the blue region corresponds to Δ
>
0. They are delimited by the axis c
=
0 and the curve Δ
=
0. This curve is formed by the graphics of the two functions:
(19)
C
±
=
1
-
40
a
-
32
a
2
±
16
a
+
1
3
32
.Figure 1
The curvesC
- and C
+.Moreover, from (16), if a point x
i
(
α
,
c
), i
=
1,2
,
3, is real, it must be satisfied x
i
(
α
,
c
)
≥
2 or x
i
(
α
,
c
)
≤
-
2 in order to be real the corresponding fixed points z
i
(
α
,
c
). So, we look for the separating curves x
1
(
α
,
c
)
=
2 and x
1
(
α
,
c
)
=
-
2:
(20)
x
1
α
,
c
=
2
⟹
1
3
-
7
+
2
α
+
1
3
f
α
,
c
3
+
g
α
,
c
3
=
2
.
Raising to the cube both sides of the equality and taking into account that
(21)
f
α
,
c
3
+
g
α
,
c
3
3
=
f
α
,
c
+
g
α
,
c
+
3
f
α
,
c
g
α
,
c
3
f
α
,
c
3
+
g
α
,
c
3
=
f
α
,
c
+
g
α
,
c
+
3
1
-
4
α
+
4
α
2
-
12
c
f
α
,
c
3
+
g
α
,
c
3
,
we obtain
(22)
216
(
α
(
4
+
c
)
-
2
(
5
+
c
)
)
=
0
⟹
c
=
2
-
5
+
2
α
2
-
α
.
Similarly,
(23)
x
1
(
α
,
c
)
=
-
2
⟹
216
α
c
=
0
⟹
c
=
0
or
α
=
0
.It can be checked that the same separating curves appear when we considerx
2
α
,
c
=
±
2 and x
3
α
,
c
=
±
2. In Figure 2, we observe the different regions separated by these curves and the number of real and complex fixed points in each region. Let us note that the fixed points change their nature, from complex to real or vice versa, when they cross these curves. For example, for α
=
-
1, there are six real roots for negative values of c below the curve c
=
2
(
-
5
+
2
α
)
/
(
2
-
α
); two of them change to complex when they cross this curve and other two roots change from real to complex when crossing the axis c
=
0 (see Figure 3). An exhaustive study of these bifurcations has been made in [12].Figure 2
Bifurcation curves of fixed points.Figure 3
Bifurcation diagram of fixed points forα
=
-
1.


## 3. Study of Critical Points

Critical points satisfyO
p
′
z
,
α
,
c
=
0, where the prime means the derivative of O
p
z
,
α
,
c with respect to z. The expression of this derivative is
(24)
O
p
′
(
z
,
α
,
c
)
=
-
2
z
2
1
+
z
4
Q
z
,
α
,
c
1
+
z
4
2
α
z
-
1
-
2
z
+
4
c
z
3
1
+
z
2
-
8
α
c
z
4
2
,
where
(25)
Q
z
,
a
,
c
=
b
0
+
b
1
z
+
b
2
z
2
+
b
3
z
3
+
b
2
z
4
+
b
1
z
5
+
b
0
z
6
,
(26)
b
0
=
-
3
+
3
α
+
6
c
,
b
1
=
-
18
+
20
α
-
4
α
2
+
16
c
-
24
α
c
,
b
2
=
-
45
+
53
α
-
16
α
2
+
10
c
-
16
α
c
+
24
α
2
c
,
b
3
=
-
60
+
72
α
-
24
α
2
+
16
α
c
-
32
α
2
c
.The solutionsz
=
0 and z
=
∞ are associated to the roots of p
(
z
), which are always super attractors because methods of family (1) have order of convergence three. The expressions of the other critical points come from the roots of the 6-degree symmetric polynomial Q
z
,
a
,
c given in the following result.Theorem 2.
The critical points of the operatorO
p
(
z
,
α
,
c
), associated to the biparametric family of iterative methods (1) on quadratic polynomials, are(i)
z
=
0 and z
=
∞, which are associated to the nonstrange fixed points;
(ii)
seven free critical points:z
=
-
1 and
(27)
z
1
α
,
c
=
x
1
α
,
c
+
x
1
α
,
c
2
-
4
2
,
z
2
α
,
c
=
x
1
α
,
c
-
x
1
α
,
c
2
-
4
2
,
z
3
α
,
c
=
x
2
α
,
c
+
x
2
α
,
c
2
-
4
2
,
z
4
α
,
c
=
x
2
α
,
c
-
x
2
α
,
c
2
-
4
2
,
z
5
α
,
c
=
x
3
α
,
c
+
x
3
α
,
c
2
-
4
2
,
z
6
α
,
c
=
x
3
α
,
c
-
x
3
α
,
c
2
-
4
2
,where
(28)
x
1
α
,
c
=
-
b
1
3
b
0
+
s
1
+
s
2
,
x
2
α
,
c
=
-
b
1
3
b
0
-
1
2
s
1
+
s
2
+
i
3
2
s
1
-
s
2
,
x
3
α
,
c
=
-
b
1
3
b
0
-
1
2
s
1
+
s
2
-
i
3
2
s
1
-
s
2
,
s
1
=
-
q
+
D
2
3
,
p
=
-
3
+
b
2
b
0
-
b
1
2
3
b
0
2
,
s
2
=
-
q
-
D
2
3
,
q
=
b
3
-
b
1
b
0
-
b
1
b
2
3
b
0
2
+
2
b
1
3
27
b
0
3
,
D
=
4
27
p
3
+
q
2
.Proof.
The critical pointsz
=
0
,
-
1
,
∞ are obtained immediately. So, we search the roots of the polynomial Q
(
z
,
α
,
c
) in operator (24). For b
0
≠
0, the solutions of Q
z
,
a
,
c
=
0 are the same of
(29)
1
+
a
1
z
+
a
2
z
2
+
a
3
z
3
+
a
2
z
4
+
a
1
z
5
+
z
6
=
0
,
where a
i
=
b
i
/
b
0, i
=
1,2
,
3, and α
≠
1
-
2
c.
Asz
=
0 is not a root of the polynomial (25), we can divide the previous equation by z
3 and use the change of variable x
=
z
+
1
/
z, which leads to the cubic polynomial equation:
(30)
x
3
+
a
¯
x
2
+
b
¯
x
+
c
¯
=
0
,
where
(31)
a
¯
=
a
1
=
b
1
b
0
,
b
¯
=
-
3
+
b
2
b
0
,
c
¯
=
b
3
-
2
b
1
b
0
.
The quadratic term is eliminated by means of the change(32)
x
=
y
-
a
¯
3
,
and the final polynomial equation is
(33)
y
3
+
p
y
+
q
=
0
,
where
(34)
p
=
b
¯
-
a
¯
2
3
=
-
3
+
b
2
b
0
-
b
1
2
3
b
0
2
,
q
=
2
a
¯
3
27
-
a
¯
b
¯
3
+
c
¯
=
b
3
-
b
1
b
0
-
b
1
b
2
3
b
0
2
+
2
b
1
3
27
b
0
3
.
By using the changey
=
s
1
+
s
2, y
3
=
s
1
3
+
s
2
3
+
3
s
1
s
2
y is obtained. By identifying the coefficients of this equation with the coefficients of (33) we obtain
(35)
s
1
3
+
s
2
3
=
-
q
,
s
1
3
s
2
3
=
-
p
3
27
.
So, s
1
3 and s
2
3 are solutions of the quadratic equation W
2
+
q
W
-
p
3
/
27
=
0; that is,
(36)
s
1
3
,
s
2
3
=
-
q
±
q
2
+
4
p
3
/
27
2
.
Therefore, the three solutions of (33) are
(37)
y
1
=
s
1
+
s
2
,
y
2
=
-
1
2
s
1
+
s
2
+
i
3
2
s
1
-
s
2
,
y
3
=
-
1
2
s
1
+
s
2
-
i
3
2
s
1
-
s
2
,
where
(38)
s
1
=
-
q
+
D
2
3
,
s
2
=
-
q
-
D
2
3
,
D
=
4
27
p
3
+
q
2
.
Undoing all the changes, the six roots of (25) are obtained from
(39)
z
=
x
±
x
2
-
4
2
,
for the different values of variable x.The number of real and complex roots of (33) depends on the sign of D. If D
>
0, there are one real and two complex roots, meanwhile if D
≤
0, the three roots are real.As(40)
D
=
-
64
α
3
c
-
10
α
2
-
1
+
2
α
3
+
α
-
3575
+
5360
α
-
2696
α
2
+
384
α
3
+
36
α
4
c
-
8
-
2
+
α
-
5
+
3
α
3
c
2
×
2187
-
1
+
α
+
2
c
4
-
1
,
the curve D
=
0 implies α
=
0, c
=
0, or
(41)
-
10
α
2
-
1
+
2
α
3
+
α
-
3575
+
5360
α
-
2696
α
2
+
384
α
3
+
36
α
4
c
-
8
-
2
+
α
-
5
+
3
α
3
c
2
=
0
.This last equation defines two functions:(42)
C
±
=
-
α

±
α
3
5
-
6
α
+
2
α
2
95
-
58
α
+
6
α
2
3
-
3575
+
5360
α
-
2696
α
2
+
384
α
3
+
36
α
4
w
±
α
3
5
-
6
α
+
2
α
2
95
-
58
α
+
6
α
2
3

-
α
±
α
3
5
-
6
α
+
2
α
2
95
-
58
α
+
6
α
2
3
-
3575
+
5360
α
-
2696
α
2
+
384
α
3
+
36
α
4
×
16
-
2
+
α
-
5
+
3
α
3
-
1
,
whose domain is {
α
:
α
≤
(
1
/
6
)
29
-
271
∪
α
≥
(
1
/
6
)
29
+
271
}. We can check that α
=
2 is an asymptote of c
+ and α
=
5
/
3 is an asymptote of C
-; moreover, C
- has the horizontal asymptote c
=
20
/
9 and C
+ has not any horizontal asymptote.Taking into account the change (39), the number of real and complex roots also depends on whether x
≥
2 or x
≤
-
2. Taking x
=
2, we get
(43)
s
1
+
s
2
=
6
+
a
1
3
⟹
s
1
+
s
2
3
=
6
+
a
1
3
3
.
As s
1
s
2
=
-
p
/
3 and s
1
3
+
s
2
3
=
-
q, we obtain the equation
(44)
16
(
α
-
2
)
(
6
+
α
(
-
4
+
c
)
-
2
c
)
3
(
-
1
+
α
+
2
c
)
=
0
.
So, if -
1
+
α
+
2
c
≠
0, we obtain the bifurcation curves α
=
2 and c
=
2
2
α
-
3
/
(
α
-
2
), α
≠
2.Takingx
=
-
2, then
(45)
-
80
α
2
c
3
(
-
1
+
α
+
2
c
)
=
0
.So, if-
1
+
α
+
2
c
≠
0, we obtain the bifurcation curves α
=
0 and c
=
0. All these bifurcation curves can be seen in Figure 4.Figure 4
Bifurcation curves of critical points.These curves separate the(
α
,
c
)-plane into different regions with different number of real and complex critical points. These regions are shown in Figures 5, 6, and 7.Figure 5
Critical points in the(
α
,
c
)-plane for α
<
0.Critical points in the(
α
,
c
)-plane for 0
<
α
<
2.
(a)
0
<
α
<
2
(b)
Detail of the zone for0
<
α
<
0.5Critical points in the(
α
,
c
)-plane for 1.9
<
α and c
>
0.
(a)
1.9
<
α
<
2.2
(b)
α
>
2.2Let us remark that there are only three free critical points because they satisfy the relationships:(46)
z
1
z
2
=
1
,
z
3
z
4
=
1
,
z
5
z
6
=
1
,
for every values (
α
,
c
).If-
1
+
α
+
2
c
=
0, then polynomial (25) becomes a product of z and a symmetric 4-degree polynomial. Following the same steps as before, on the line α
=
1
-
2
c and α
≠
±
5
/
2, we have for α
<
(
1
/
2
)
(
-
5
-
65
) two real and two complex roots, for (
1
/
2
)
(
-
5
-
65
)
<
α
<
0 four real roots, for 0
<
α
<
1 four complex roots, for 1
<
α
<
(
1
/
2
)
(
-
5
+
65
) two real and two complex roots, for (
1
/
2
)
(
-
5
+
65
)
<
α
<
2 four complex roots, and for α
>
2 two real and two complex roots. Finally, for α
=
±
5
/
2, polynomial (25) yields a symmetric 2-degree polynomial. For α
=
-
5
/
2 both roots are complex and for α
=
5
/
2 both roots are real.


## 4. Chebyshev-Halley Class andc-Family of Iterative Methods

We know that Chebyshev-Halley class andc-family of iterative methods are included in the biparametric family (1). In this section, we prove that the 6-degree polynomial (8) can be factorized as a product of a polynomial of degree two and a polynomial of degree four, whose roots correspond to the fixed points of Chebyshev family for c
=
0 and c-family for α
=
0, respectively.So,(47)
P
z
,
α
,
c
=
h
1
z
,
α
,
c
h
2
z
,
α
,
c
,
where
(48)
h
1
z
,
α
,
c
=
1
+
d
z
+
z
2
,
h
2
z
,
α
,
c
=
1
+
d
1
z
+
d
2
z
2
+
d
3
z
3
+
z
4
.To fulfill this relationship, the coefficients must satisfy the relations(49)
d
1
+
d
=
7
-
2
α
,
d
2
+
d
d
1
+
1
=
19
-
8
α
+
4
c
,
d
3
+
d
d
2
+
d
1
=
26
-
12
α
+
8
c
-
8
α
c
,
1
+
d
d
3
+
b
2
=
19
-
8
α
+
4
c
,
d
+
d
3
=
7
-
2
α
.
From the first and the last relations, it follows that d
1
=
d
3. Clearing d
1 and d
2 in terms of d in the first two equations, we get
(50)
d
1
=
7
-
2
α
-
d
,
d
2
=
18
-
8
α
+
4
c
-
7
d
+
2
d
α
+
d
2
,
and substituting in the third one, we obtain again the cubic equation:
(51)
d
3
+
-
7
+
2
α
d
2
+
16
-
8
α
+
4
c
d
-
12
+
8
α
-
8
c
-
8
α
c
=
0
.For(
α
,
c
) real, the real solution of this polynomial is
(52)
d
α
,
c
=
1
3
7
-
2
α
-
1
3
f
α
,
c
3
+
g
α
,
c
3
,
with f and g previously defined.Forc
=
0, then f
(
α
,
0
)
=
g
(
α
,
0
)
=
-
1
+
2
α
3; then,
(53)
h
1
z
,
α
,
0
=
1
+
b
z
+
z
2
=
1
+
3
-
2
α
z
+
z
2
,
whose zeros coincide with the strange fixed points of the Chebyshev-Halley family (see [1]).Moreover, forα
=
0, the cubic equation (51) is
(54)
d
3
-
7
d
2
+
16
+
4
c
d
-
8
c
-
12
=
0
,
whose solutions are
(55)
d
=
2
,
d
=
5
2
+
1
2
1
-
16
c
,
d
=
5
2
-
1
2
1
-
16
c
,
and if d
=
2, then d
1
=
5 and d
2
=
8
+
4
c. So, the zeros of the function,
(56)
h
2
z
,
0
,
c
=
1
+
5
z
+
8
z
2
+
4
z
2
c
+
5
z
3
+
z
4
,
coincide with the strange fixed points z
i, i
=
1,2
,
3,4, of the c-family (see [5]).On the other hand, forc
=
0, the critical points are
(57)
z
1
(
α
,
0
)
=
3
-
4
α
+
2
α
2
+
α
(
-
2
+
α
)
(
-
3
+
2
α
)
(
-
1
+
2
α
)
3
(
-
1
+
α
)
,
z
2
(
α
,
0
)
=
3
-
4
α
+
2
α
2
-
α
(
-
2
+
α
)
(
-
3
+
2
α
)
(
-
1
+
2
α
)
3
(
-
1
+
α
)
,
z
3
(
α
,
0
)
=
z
4
=
z
5
=
z
6
=
-
1
,
which are the critical points of Chebyshev-Halley family (see [1]).Finally, forα
=
0, the critical points are
(58)
z
1
(
0
,
c
)
=
3
+
4
c
+
2
5
c
(
3
-
c
)
3
(
-
1
+
2
c
)
,
z
2
(
0
,
c
)
=
3
+
4
c
-
2
5
c
(
3
-
c
)
3
(
-
1
+
2
c
)
,
z
3
(
0
,
c
)
=
z
4
=
z
5
=
z
6
=
-
1
,
which coincide with the critical points of c-family (see [5]).


## 5. Study of the Parameter Space

As we have said, the critical points whose analytical expression has been previously obtained satisfy thatz
2
k
-
1
(
α
,
c
)
=
1
/
z
2
k
(
α
,
c
), for k
=
1,2
,
3. So, only three of them are free critical points.It is a known fact that any basin of attraction must hold, at least, one critical point (see, e.g., [10] or [11]). Then, the knowledge of the critical points and their behavior as a starting point of the iterative process is of vital importance. In this section, the parameter plane associated to each free critical point is obtained, by using them as starting points of the (
α
,
c
)-family and taking note of their orbits depending on the real values of the parameters involved. So, in Figure 8 each point (
a
,
b
) of the real plane [
-
10,10
]
×
[
-
10,10
] represents a scheme of the (
α
,
c
)-family, for α
=
a and c
=
b. Then, the color of this point is drawn in red if this specific method converges to one of the roots (0 and ∞, after Möbius transformation), or in white in other case. Besides, the color will be brighter if the number of iterations needed to reach one of the roots is small and darker if it is high. The routines used to get these pictures are slightly modified from those appearing in [13]. Similar features can be found in [14].Parameter planes corresponding to different critical points.
(a)
Critical pointsz
1 and z
2
(b)
Critical pointsz
3 and z
4
(c)
Critical pointsz
5 and z
6So, we will observe the different parameter planes in order to find the stable behavior (in red) and the “dangerous” regions (in white). In the last ones, divergence, attracting strange fixed points or even attracting periodic orbits, can appear.In Figure8(a), the parameter plane associated to critical points z
1 and z
2 can be observed. Wide white areas of no convergence to the roots are found, but also the number of elements of the (
α
,
c
)-family that can be used to converge to them is high.A similar behavior can be found in Figures8(b) and 8(c). Stable behavior appears in wide areas of the region where α
>
0 and c
>
0 and also for α
<
0 and 0
<
c
<
3. In the following, we select different points in stable and unstable regions of the parameter planes (i.e., we choose some iterative methods of the class) and we see their specific dynamical behavior on quadratic polynomials.


## 6. Dynamical Planes

The analysis of dynamical planes of a parametric family or a method is a usual technique for improving the knowledge of a method; see, for example, [15–18] and references therein.The generation of dynamical planes is very similar to the one of parameter spaces. In case of dynamical planes, the value of each parameter,α and c, is constant. So, each point of the complex plane is considered as a starting point of concrete element of the family of iterative methods and it is painted in different colors depending on the point which it has converged to. If it converges to 0, then the point is painted in orange, when it converges to ∞, it appears in clear blue. If the iterates converge to a strange fixed point it is painted in different colors, as green or dark blue. Finally, the basins of attraction of periodic orbits are painted in black. As in the case of the parameter planes, routines appearing in [13] have been used.Firstly, let us take into account the information obtained in the parameter planes (Figure8). The most stable behavior corresponds to those (
α
,
c
) values that remain in red at every parameter plane: they correspond to methods of the family that, using any free critical point as an initial estimation, it only converges to zero or infinite, that is, to any of the roots of the original quadratic polynomial. However, it can be seen at Figure 9 that all these dynamical planes are not equally useful, as the wideness of the immediate basin of attraction is clearly bigger in some cases.Stable behavior.
(a)
α
=
9, c
=
3
(b)
α
=
1
/
2, c
=
-
1
/
5
(c)
α
=
-
10, c
=
7
/
2
(d)
α
=
-
5, c
=
-
1
/
5With respect to the unstable behavior, in Figures10 to 12, different situations can be observed. Specifically, in Figures 10(a) and 10(b), two periodic orbits of period 4 are found for values of the parameters α
=
5 and c
=
5. Indeed, in Figure 11(b), a two-periodic orbit is observed for α
=
c
=
10, with a wide basin of attraction in black. In Figure 11(a), a two-periodic orbit is found, but also an attracting strange fixed point whose basin of attraction is much bigger than the one of z
=
0. Finally, in Figure 12, some dynamical planes with four different basins of attraction are shown, two of them corresponding to attracting strange fixed points.Unstable behavior: two 4-periodic orbits.
(a)
α
=
5, c
=
5
(b)
α
=
5, c
=
5Unstable behavior: one 2-periodic orbit.
(a)
α
=
5, c
=
-
5
(b)
α
=
10, c
=
10Unstable behavior: attracting strange fixed points.
(a)
α
=
5, c
=
-
1
(b)
α
=
2.5, c
=
-
1Let us remark that in all cases, a cubic symmetry is observed. This fact is due to the symmetry of the polynomials (with respect to their coefficients) involved in the rational operator associated to the family.


## 7. Conclusions

In this paper, the dynamical behavior of the(
α
,
c
)-family has been studied, confirming that it is much more than the union of cases α
=
0 (c-family) and c
=
0 (Chebyshev-Halley family), as there exist hybrid iterative methods with specially stable behavior. The fixed and critical points of the rational operator associated to the family on quadratic polynomials were obtained, and the regions in which these points are real or complex have been shown. From critical points, the parameter planes of the family were drawn and some regions of very stable and unstable behavior have been detected. Dynamical planes for specific values of the parameters in these regions were shown.



---

*Source: 141643-2014-12-24.xml*