Computational Bifurcations Occurring on Red Fixed Components in the λ -Parameter Plane for a Family of Optimal Fourth-Order Multiple-Root Finders under the Möbius Conjugacy Map

: Optimal fourth-order multiple-root ﬁnders with parameter λ were conjugated via the Möbius map applied to a simple polynomial function. The long-term dynamics of these conjugated maps in the λ -parameter plane was analyzed to discover some properties of periodic, bounded and chaotic orbits. The λ -parameters for periodic orbits in the parameter plane are painted in different colors depending on their periods, and the bounded or chaotic ones are colored black to illustrate λ -dependent connected components. When a red ﬁxed component in the parameter plane branches into a q -periodic component, we encounter geometric bifurcation phenomena whose characteristics determine the desired boundary equation and bifurcation point. Computational results along with illustrated components support the bifurcation phenomena underlying this paper.


Introduction
The governing equations of natural phenomena are usually represented as nonlinear equations in the form of f (x) = 0, when solving the relevant algebraic or differential equations obtained through the scientific observations, experiments, experiences as well as mathematical modeling. Since the exact solutions of such nonlinear governing equations are rarely feasible, it is considered best to rely on numerical solution techniques such as multi-point iterative methods [1][2][3][4], splitting methods [5,6], semi-explicit methods [7,8] as well as analytical-numerical methods [9]. Among such techniques, the most widely used multi-point iterative method under general circumstances is Newton's method which converges quadratically to a simple root and linearly to a multiple root. In order to more accurately solve the nonlinear governing equations, many researchers have been developing new iteration methods with high convergence orders, including Newton's method as a particular case. Motivation for writing this paper is to analyze the convergence behavior of these iterative numerical methods from the viewpoint of the limit behavior under the iteration dynamics with illustrative computational results.
Although the limit behavior of the aforementioned numerical solution techniques draws great attention to the readers, our study will focus on the limit behavior of a family of Newton-like iterative methods with a simple structure of a weight function possessing a λ-control parameter. The convergence behavior of these iterative numerical methods equipped with λ-control parameters can be analyzed from a viewpoint of iteration dynamics. Let us begin with a family of nonlinear governing equations with parameter λ as represented by the following nonlinear recurrence relation: x n+1 = W f (λ; x n ), n = 0, 1, 2, · · · , (1) where W f is called the iteration function or fixed point operator [10]. If iteration index n is regarded as an indicator of time, then the nonlinear recurrence relation (1) would represent itself the relevant discrete dynamical system. Consequently, pursuing the convergence of the root of the nonlinear recurrence relation (1) is initiated by dealing with its limit behavior under the action of the nonlinear discrete dynamical system (1). The governing equation to be dealt with in this paper is a family of optimal [11] two-point Newton-like fourth-order multiple-root finders with parameter λ studied in Reference [3]: where Q f : C 2 → C is analytic [12] in a neighborhood of (λ, 0) with λ ∈ C fixed and t is taken on the principal analytic branch of the m-th root. In the current study, we will employ a specific form of Q f with a controllable parameter λ ∈ C as follows: . (3) The advantage of choosing such a rational function Q f is to allow the conjugated W f to contain the low-order polynomial factor T(λ; z) with exact λ-dependent roots, that is, the exact fixed points of F from fixed-point relation (6) under the Möbius conjugacy to be introduced in Section 4. If a different type of Q f , say, a polynomial were chosen, then the higher order of T(λ; z) would hinder us from locating the exact fixed points.
In Reference [13], as a part of long-term dynamics for the conjugated map of (2), a study was conducted on bifurcation [14,15] phenomena occurring on the boundary curve of the yellow fixed component related to λ-free strange fixed point z = 1. The number of the component was turned out to be only one. Moreover, the equation of the boundary curve turned out to be a circle as one of the known geometric shapes. This simple circle allowed us to analyze the global bifurcations occurring on it without difficulties.
The aim of this paper is to study of the fixed points under the Möbius conjugacy map and to investigate in depth the boundary curves and other geometric properties of the hyperbolic components [16] in the λ-parameter plane (to be defined in Section 4). The current study will cover a challenging task not handled in the previous study in Reference [13]. We will extensively investigate the bifurcation phenomena of the red fixed component R f (whose notation and analysis will be described in Sections 3 and 4) related to the strange fixed point z = z(λ), which is λ-dependent, that is, not λ-free as in the case of the previous study for the yellow fixed component. In particular, it is our curiosity and concern to find out whether the number of R f is more than one, and whether the equation of the boundary curve represents a known geometric figure. From the perspective of plane geometric curves, our main interest lies in the computational, theoretical and visual exploration of bifurcation phenomena occurring on the boundary curve of R f .
A theoretical formulation of the boundary equation of R f developed here will serve as a basis for describing bifurcation phenomena occurring in another discrete dynamical system. For example, it is helpful to establish the boundary equations of hyperbolic components with small periods in the Mandelbrot set [17] for a better analysis of bifurcation phenomena occurring there.
The remainder of this paper further describes the important points for the analysis of bifurcation phenomena. Described in Section 2 are preliminary results on the long-term behavior of a dynamical system and local bifurcations. In Section 3, we describe various properties of the fixed and critical points under the Möbius conjugacy map. Section 4 investigates a long-term dynamical behavior and the symmetry of the parameter plane. The location of the bifurcation points is extensively explored based upon the existence of a common osculating point between the fixed component (maximal connected set) [12] and q-periodic component. Section 5 describes the boundary curves of the fixed components along with their geometric properties. Finally in Section 6, we summarize overall studies and state the possible future work on the long-term dynamics of period-2 components.

Preliminary Results
Although iteration function W f in (1) is in general meromorphic, it suffices to treat the dynamics of rational functions only, in view of Lemma of 2.1 of Reference [13]. Indeed, we apply f (z) = ((z − a)(z − b)) m to obtain the conjugation of rational W f , namely, F = h • W f • h −1 in the form of a rational function to be shown in (4) via the first-order rational Möbius conjugacy map h(z) = z−b z−a , with a = b [18].
The following definition of conjugacy and other assertions are useful for continuing discussion of the relevant dynamics.
Such h is called a conjugacy [19].
According to Theorem of 2.1 in Reference [13], we find: Theorem 1. Let F , K ∈ C 1 and be conjugate to each other via the diffeomorphic conjugacy h. In addition, let ξ be a fixed point of K. Then the following hold: (a) The fixed point property remains invariant under a topological conjugacy h, that is, (b) The Poincaré characteristic multiplier [15] of ξ by K, denoted by m(K, ξ), is invariant under a diffeomorphic conjugacy h, that is, Hence, the topological conjugacy h gives an isomorphism between two orbits of F and K.
Then the orbit behavior of F exhibits a behavior similar to that of the iterative function W f . According to Theorem 1, the invariant properties of the fixed point and the Poincaré characteristic multiplier after the conjugacy transformation address that F and W f have the advantage that their orbit dynamics is invariant. Such an advantage may give the favorable geometric characteristics induced by the orbit behavior of F . Lemma 1. Suppose that ψ : Ω ⊂ C → Ω is meromorphic and has a fixed point ξ ∈ Ω with |ψ (ξ)| < 1. Then ψ has a unique fixed point ξ such that sequence {z n+1 = ψ(z n )} ∞ 0 converges to ξ for any given z 0 ∈ Ω.
Proof. (i) For ξ ∈ Ω ⊂ C, the proof is done by Lemma 1 of Reference [20]. (ii) For ψ(∞) = ∞, it is possible to make |ψ (∞)| < 1 under some restrictions on the form of ψ, which can be justified based on the following description. By definition from Section 3.1 of Reference [21], ψ is said to be analytic at ∞ if H • ψ • H −1 is analytic at 0 with H(z) = 1/z. Due to Lemma 2.1 of Reference [13], we can express ψ(z) = a 0 +a 1 z+···+a n z n b 0 +b 1 z+···+b m z m as a rational function with a n b m = 0, n > m for n, m ∈ N ∪ {0}. We further define a function if n > m + 1 by direct computation or from Section 2.6 of Reference [19]. Hence if restriction |ψ (∞)| = | b m a n | < 1 is made, then ∞ becomes an attracting fixed point of ψ.
If we select z 0 as a critical point z * of ψ satisfying ψ (z * ) = 0, then Lemma 1 with z 0 = z * yields the following corollary. Corollary 1. Let q ∈ N and ξ be an attracting q-periodic point of ψ in Lemma 1. Then every critical orbit of ψ tends to ξ. Corollary 1 will theoretically support the formation of hyperbolic components serving as important parts of the structure of the parameter space to be introduced later in Section 4.
We now define a function G(λ; z) = ψ λ (z) for a given parameter λ ∈ C. Let ξ and ω respectively denote a fixed point and a critical point of G(λ; z). Then, in view of Lemma 1, | d dz G(λ; ξ)| = |G (λ; ξ)| = |ψ λ (ξ)| = 1 would influence the long-term orbit behavior of G(λ; ω) as the parameter λ varies across the boundary of region Ω ⊂ C. The unit circle characterizes the stability G (λ; ξ) and is better designated as the stability unit circle S for further analysis.
An abrupt change of such a long-term orbit behavior of G(λ; ω) often arises at λ = λ * , where the value of G (λ; ξ) crosses the boundary of S. This kind of abrupt change in its dynamical behavior is said to be a bifurcation. Consulting the notions of such bifurcations in Reference [22], the following three types of bifurcation may occur: (1) If G (λ * ; ξ) = 1, then the bifurcation is of (cyclic) fold(saddle-node).
(3) If G (λ * ; ξ) = e iθ , then the bifurcation is of Neimark-Sacker(secondary Hopf), Such a control parameter λ * in the parameter plane is called a bifurcation point whose location will be treated in Section 4.2.
According to Theorem 1-(a), the fixed points 0 and They are super-attractive, free of λ and also imply the critical points of F (λ; z) in view of (5).
Nevertheless, because their orbits approach themselves, the dynamical impact is expected to be small. Fixed points excluding such fixed points {0, ∞} are termed as strange fixed points [28] which are not the roots of f (z) = (z − a) m (z − b) m . Notice that z = 1 is a λ-free strange fixed point whose dynamics was explored in Reference [13]. The roots of T(λ; z) = 0 in (6) are indeed λ-dependent strange fixed points. Simple computations directly prove the following lemma.
hold for any λ ∈ C and any z ∈ C.
We seek λ-dependent strange fixed points ξ λ by solving T(λ; ξ λ ) = 0 for a fixed λ ∈ C. According to Proposition 1-(a), T(λ; z) can be written in the form: Hence, we are able to express all the desired λ-dependent strange fixed points of F as ξ (k) λ and 1/ξ (k) λ with the notation ξ Figure 4 of Reference [13] illustrates the stability surfaces of λ-dependent strange fixed points ξ λ The super-attractors of F (λ; z) can be found by solving T(λ; z) = Q(λ; z) = 0 for z and λ. Eliminating z in T(λ; z) = Q(λ; z) = 0 yields a 4-degree polynomial from which we find λ ∈ {−4.7034, 3.86866, −2, −1.83192} for the desired super-attractors. Note that the cubic factor in (8) with positive discriminant value of 1, 594, 800 has three distinct real roots. On the other hand, eliminating λ in T(λ; z) = Q(λ; z) = 0 yields an 8-degree polynomial which induces 8 supper-attractors denoted by 4 pairs of (z j , 1/z j ) for 1 ≤ j ≤ 4 with The critical points of F are given by the roots of F (λ; z) = 0. A quick glance of (5) reveals that z = −1 is a λ-free critical point. Evidently z = 0 and z = ∞ are critical points that respectively correspond to the roots a and b of the polynomial (z − a) m (z − b) m . The critical points excluding the both roots a and b are termed as free critical points [28]. Since the orbit behavior of λ-dependent critical point of F (λ; z) draws our attention, we will treat free critical points that are λ-dependent.

Parameter Plane and Long-Term Orbit Behavior
One way to efficiently handle the long-term orbit dynamics of the iterative map F (λ; z) in (4) is to visually illustrate its dynamical orbit behavior in the λ-parameter plane whose definition will be made in this section. The long-term critical orbit under the action of F (λ; z) will induce attracting periodic, finitely bounded or chaotic orbits for each given λ selected in the λ-parameter plane. It is certainly interesting to paint the point λ for which the long-term orbit approaches an attracting periodic orbit, with different color depending on its period. Corollary 2 and Remark 2 well state the orbit behavior between two critical points z and 1/z of F .

Corollary 2.
Let q ∈ N be given. If z ∈ C is a q-periodic point of F (λ; z), then so is 1 z .
Proof. Lemma 2 directly shows that F (λ; z) = z and F (λ; 1/z) = 1/z for a fixed point z of F with a fixed λ ∈ C. Hence there exists a topological conjugacy h(z) = 1/z such that z , completing the proof.
In light of Remark 2, it suffices to consider only one branch of the critical points ω 1 for its typical orbit behavior. We define the parameter plane(space) P for iterative map F (λ; z) as follows: P = {λ ∈ C : an orbit of a critical point ω tends to a number γ ∈ C under the action of F (λ; ω)}.
From now on, let ω denote the critical point ω 1 for convenience. Consider the limit behavior of an orbit of F (λ; ω) approaching a number γ ∈ C in the long run. If γ = ∞, then it draws very little attention to the limit behavior since F (λ; ∞) = ∞. On the contrary, if γ is finite, then it will create a diversity of the limit behavior as λ ∈ C varies. More importantly, if γ is finite as well as bounded, then Bolzano-Weierstass Theorem guarantees the existence of a subsequence of the orbit of F (λ; ω) converging to the fixed point ξ(λ). In other words, we can find a q-periodic point ξ(λ) = γ for some q ∈ N in such a way that by expressing = p j + q for any , p, j ∈ N with p < and q ∈ {0, 1, 2, · · · , p − 1}; if q = 0, then it suggests that a non-periodic bounded orbit could exist. As a result, the q-periodic orbit of ξ(λ) = γ forms a q-cycle: ¶ ξ(λ), F (λ; ξ(λ)), F 2 (λ; ξ(λ)), · · · , F q−1 (λ; ξ(λ)) © .
Taking the discussion in Section 3 into consideration, we find the possible fixed points We now define a function G(λ; z) = F q (λ; z) − z for any given q ∈ N. As a consequence, G(λ; z) establishes a bivariate rational function of z and λ on C. Since the rational function G(λ; z) has a finite number of zeros z = ξ(λ), that is, the fixed points of F q (λ; z) for any given q ∈ N, we can find a λ for any periodic-point z = ξ(λ) being approached by the long-term critical orbit of F (λ; ω).
To debate the stability of the q-periodic point ξ(λ), we first consider the case for q = 1. Due to the fact that F (λ; 0) = 0 and F (λ; ∞) = ∞, we find that 0 and ∞ are the fixed points as well as super-attractors in view of Proposition 1-(d). In this case, regions of corresponding λ-values are respectively colored in cyan and magenta in the parameter space P in Figure 1. For the fixed point ξ(λ) = 1, we have F (λ; 1) = 1 and |F (λ; 1)| < 1 for λ in the yellow region of P in Figure 1.  -3. -2. -1.

4.
5. Let us zoom in on the three regions indicated by the dotted rectangles in Figure 1 so that the red fixed components denoted by R  λ for 1 ≤ k ≤ 2. We next consider the case for q ≥ 2. Observe that the attracting fixed point ξ(λ) maintains its stability for λ in the red region R f i . As λ crosses the finite boundary of R f i , that is, the fixed point ξ(λ) loses its stability, but other attracting q-periodic points with begin to emerge along the boundaries of R f i . Hyperbolic components of λ-values for these q-periodic points are painted according to the color chart in Figure 5 of Reference [13]. Some of q-periodic components in P are indicated by arrow numbers. If λ is chosen in such a q-periodic component, then the orbit of critical point ω attracts to a q-periodic point of F . To generate the λ-parameter plane P, Mathematica software has been used with 16-digits of precision for real and complex numbers. The desired convergence criterion for a q-periodic orbit was 10 −6 after maximum 3000 iterations.
It is evident that proper components of P can be found to select best members of a family (2) locating the desired roots z = 0 or z = ∞ with acceptable numerical stability, provided λ-parameters belong to such components.
According to Lemma 3.5 in Reference [29] and Theorem 5.2 in Reference [13], the parameter space P is symmetric about its horizontal axis.

Bifurcation Points along a Red Fixed Component R f
Suppose that a point λ(t) on a piecewise analytic curve in the complex plane satisfies relation: An exemplary curve consisting of such points λ is sketched in Figure 3 along with attracting q-periodic satellite [17] components denoted by C q and identified by numbers q ∈ N. Observe that λ (t) = dλ dt represents a tangential vector at a parametric location t ∈ [t i , t f ]. We extensively pursue the case when λ(t) traces a boundary curve of a typical red fixed component R f as t ∈ [t i , t f ] varies. For convenience of notation, we let ∂R f denote the finite boundary of R f , and write F λ (z) ≡ F (λ; z) for simple treatment of derivatives. Let λ be an osculating point ∈ ∂R f where the λ-dependent critical orbit of F approaches an attracting (hyperbolic) q-periodic point ξ. We find that the q-periodic point ξ is dependent upon parameter t, since the boundary ∂R f depends on t ∈ [t i , t f ]. Hence, we express ξ = |ξ|e i t with t ∈ [0, 2π).
Given q ∈ N, the q-periodic point ξ satisfies the following set of two equations at the common boundary point λ ∈ ∂R f ∂C q as shown in Figure 3.
Solving for λ from the fixed point relation F λ (ξ(t)) = ξ(t) via (4), we are able to obtain: where H(z) is a rational function in z in view of (4) and differentiable at points except at a finite number of poles. This λ will trace a curve in the complex plane as parameter t varies. Hence, by the chain rules, we find the desired derivative at the fixed point ξ: from which we have α q = 1 in view of the fact that F q λ (ξ) = ξ. Accordingly, Equation (14) reduces to: Since dλ dt represents a tangent line at a location t, Equation (16) states the common tangent line of ∂R f and ∂C q . The discussion thus far leads us to the notion of /q-bifurcation point below.
By solving α(λ) = e i2π /q for λ, we are able to find the /q-bifurcation point λ = λ ,q in terms of and q, provided that a set of coprime values of ( , q) are given. If q = 1 then = 0 must occur and gives λ ,q = 1, which is the fold(saddle-node) bifurcation point in the sense of local bifurcations discussed in Section 2. At this kind of fold bifurcation point, another fixed component different from R f might arise or the same fixed component as R f itself merge together due to the connectedness of the component. In the latter case, such a fold bifurcation point λ ,q = 1 might be a cusp-like point where R f itself becomes a primitive component. If ( , q) = (1, 2), then λ ,q = −1 becomes the flip (period-doubling) bifurcation point. If q > 2 and = 0, then λ ,q becomes the Neimark-Sacker (secondary Hopf) bifurcation point.
We now conveniently denote the closed boundary curves of R Moreover, coprime pairs of ( , q) give us the number of bifurcation points λ ,q as the length |F n | of the Farey [30] sequence F n of order n. Theorem 4 in Reference [31] directly gives |F n |: with φ(j) as Euler's totient function.   Table 3. Typical λ ,q along B 3 at t = 2π q , with λ 3 (t) for t ∈ [0, π] and λ 2 (t) for t ∈ (π, 2π).

Establishing Parametric Equations of ∂R f
We establish a boundary equation for such an attracting red fixed component R f i associated with a typical λ-dependent fixed point ξ given by a root of T(λ; ξ) = 0. During the course of development, we will show the existence of three red fixed components R f , each boundary of which consists of two analytic branches forming a piecewise smooth closed curve. Details of such a boundary curve along with the location of bifurcation points will be discussed in the next section.
Consulting (5) and (6), we find that the fixed point ξ satisfies a set of equations with λ ∈ ∂R f as follows: where c k are the two roots being dependent on λ from (7); moreover, for preferable algebraic treatments, expressions Q(λ; z), q(λ; z) and F (λ; z) are written as Q λ (z), q λ (z) and F λ (z) in order. We reconstruct the second equation of (19) with a real parameter t ∈ [0, 2π): To make notations short, let us write ξ = z, c k = c, and compute λ from (7) in terms of c: which defines a branch curve of ∂R f . It is our interest to express λ in terms of the parameter t. By use of the symbolic Mathematica commands PolynomialQuotient and PolynomialRemainder described in Section 3.3.4 of Reference [32], we find: where Ω q (z) and Θ q (z) are 9-and 6-degree polynomials in order; Ω r (z) = 2(18 − 7c − 4c 2 + c 3 )δ(z) and Θ r (z) = −4(c − 2)δ(z) with δ(z) = ( c−2 c−3 ) 2 c 3 − 2c + (1 − 3c 2 + c 4 )z . Due to the fact that (z 2 + cz + 1) = 0, relation (20) directly yields: where c = c k for k ∈ {1, 2} can be expressed in terms of λ from (7). Since F λ (ξ) is not explicitly dependent on fixed points ξ themselves but only on coefficient c, the stabilities of all four fixed points ξ found from T(λ; ξ) = 0 remain equal to one another. This kind of favorable dependence on c clearly reduces the possible lengthy algebraic work required for each fixed point. Evidently, three roots C 1 (t), C 2 (t), C 3 (t) of c = c k (t) can be found for each k ∈ {1, 2} in terms of t ∈ [0, 2π) from the cubic equation of (23): with σ = 53 − 18e i t + 3 √ 6 √ −886 + 421e i t − 68e 2i t + 4e 3i t 1/3 . Further refined computations readily offer: where η reveals the complex conjugate of η. We say the points t at which a complex function f (z) is not analytic on C are called singular points or singularities of f (z). Preferably, the continuity relations hold at the singular points t = 0 and t = π: where a + and b − imply the respective limits from the above of a and the below of b.
In view of (21) and (24), we can express three λ-values at C 1 (t), C 2 (t) and C 3 (t). Denoting λ j (t) = λ(C j (t)) for t ∈ [0, 2π) leads us to the parametric expression for each j ∈ {1, 2, 3} as follows:  The result of Equation (27) suggests that the number of the red fixed components is three and enables us to construct the three boundaries of R f parametrically as t ∈ [0, 2π) varies. Indeed, C j (t) and λ(C j (t)) have also singular points at t = 0 and t = π.
Since λ k (t) in (21) is a function of c which is piecewise analytic with such singularities, we shall have same properties given by (25) and (26) with letter C j replaced by λ j . In light of this observation along with the symmetry of P, we characterize the three closed boundary curves B j of R f as B j = {β j (t), t ∈ [0, 2π)} for (1 ≤ j ≤ 3), with Any of β j (t) is generally comprised of piecewise analytic arcs with singularities at t = 0 and t = π. The functions β j (t), (1 ≤ j ≤ 3) are plotted for t ∈ [0, 2π) in Figure 5. Let B j be a closed plane curve traced by β j (t) for t ∈ [0, 2π), which represents the boundary of j-th red fixed component R

Geometric Features of Boundaries B j
In this section, various typical geometric features, including curvatures, perimeters and bounding areas of the closed boundary curves B j , will be treated from a perspective of Euclidean planar curves. Let x j (t) = (λ j (t)) and y j (t) = (λ j (t)) respectively denote the real and imaginary part of λ j (t). By elaborate computing with approximately 10 digits of accuracy, we obtain the arc lengths of B j : as well as the bounding area A j enclosed by B j : The curvatures of B j calculated when t = π 4 , π 3 , π 2 , 2π 3 , 3π 4 , and the average curvatures are shown in Table 4. In addition, the bounding areas and perimeters are listed together with them in Table 4. The following lemma immediately follows from the geometric properties listed in Table 4.

Remark 3.
Consider a circle T = {(a c + r c cos t, r c sin t), t ∈ [0, 2π)} osculating at the topmost point (a c , r c ) of B 1 , with a c = −4.6900936467, and r c = 0.28847354384. Let ρ 1 (t) = |β 1 (t) − a c | be the distance from the center (a c , 0) of T to a point β 1 (t) at a location t. Then by elaborate computations, we have ρ 1 (t) ≥ r c for all t ∈ [0, 2π) and find T is completely inscribed in B 1 .

Conclusions
In the previous study in Reference [13], the only yellow fixed component related to λ-free strange fixed point z = 1 had a circle as its boundary curve, which described branching phenomena without difficulties. On the other hand, in the present study, the formation of the red fixed component R f is more complicated because the strange fixed point z = z(λ) (given by the roots of degree-4 polynomial T(λ; z) in (6)) depends on λ, and thus the analysis of the desired bifurcation phenomena requires a higher level of tools. One such tool is to introduce a suitable parameter, say, t, to represent a curve on the complex plane. The expression λ(t) = x(t) + iy(t), given at the beginning of Section 4.2, will be very helpful in describing the branching phenomena occurring in the λ-parameter plane P. In particular, it has enabled us to draw some geometric properties of R f observed in P.
Through the Möbius-conjugacy map applied to f (z) = (z − a) m (z − b) m , we have studied the long-term orbit behavior of a λ-family of fourth-order multiple-root finders in P. We have theoretically formulated the boundary equation of the red fixed hyperbolic component R f , by representing a point λ(t) with parameter t ∈ [0, 2π). The boundaries of R f have not turned out to be known geometric figures. In fact, they have proven to be neither circles nor cardioids. The number of R f is 3 higher than that of the yellow fixed component studied in Reference [13]. Equating the derivative of a q-periodic point to that of a λ-dependent fixed point at the common osculating point λ between R f and q-periodic components have induced the relation F λ (ξ) q = 1, which has enabled us to successfully locate the bifurcation points λ ,q directly from F λ (ξ) = e i 2π /q with and q defined in Definition 2. With λ selected in the magenta or cyan region of P, iterative scheme W f in (2) stably converges to the desired root. However, members of the family with λ selected in other regions would exhibit unfavorable numerical behavior. According to Proposition 1-a, T(λ; z) is composed of factors of the form (z 2 + c k z + 1) for k = 1, 2, · · · . Therefore, if the order of the governing equation (for instance, (7)) for the coefficients c k is not greater than 4, that is, the order of T(λ; z) is not greater than 8, then it is possible to find the exact λ-dependent strange fixed point z = z(λ). Otherwise, the bifurcation phenomena will be resolved neither theoretically nor computationally. As our future study, we will challenge to investigate the bifurcation phenomena of the period-2 orange components shown in Figure 2 and locate the desired global bifurcation points by a theoretical formulation of the desired boundary equation, provided the order of T(λ; z) is hopefully not greater than 8.