Next Article in Journal
Computer-Aided Laser-Fiber Output Beam 3D Spatial and Angular Design
Previous Article in Journal
RDF 1.1: Knowledge Representation and Data Integration Language for the Web
Previous Article in Special Issue
The Asymmetric Alpha-Power Skew-t Distribution
 
 
Erratum published on 30 October 2020, see Symmetry 2020, 12(11), 1797.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generation of Julia and Mandelbrot Sets via Fixed Points

1
Department of Mathematics, Government College University, Lahore 54000, Pakistan
2
Department of Medical Research, China Medical University No. 91, Hsueh-Shih Road, Taichung 400, Taiwan
3
Department of Sciences and Humanities, National University of Computer and Emerging Sciences, Lahore Campus, Lahore 54000, Pakistan
4
Institute of Research and Development of Processes, University of the Basque Country, Campus of Leioa (Bizkaia), P.O. Box 644, Bilbao, Barrio Sarriena, 48940 Leioa, Spain
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(1), 86; https://doi.org/10.3390/sym12010086
Submission received: 27 November 2019 / Revised: 24 December 2019 / Accepted: 26 December 2019 / Published: 2 January 2020
(This article belongs to the Special Issue Symmetry in Applied Mathematics)

Abstract

:
The aim of this paper is to present an application of a fixed point iterative process in generation of fractals namely Julia and Mandelbrot sets for the complex polynomials of the form T ( x ) = x n + m x + r where m , r C and n 2 . Fractals represent the phenomena of expanding or unfolding symmetries which exhibit similar patterns displayed at every scale. We prove some escape time results for the generation of Julia and Mandelbrot sets using a Picard Ishikawa type iterative process. A visualization of the Julia and Mandelbrot sets for certain complex polynomials is presented and their graphical behaviour is examined. We also discuss the effects of parameters on the color variation and shape of fractals.
MSC:
Primary: 47H10; Secondary:47J25

1. Introduction

Fixed point theory provides a suitable framework to investigate various nonlinear phenomena arising in the applied sciences including complex graphics, geometry, biology and physics [1,2,3,4]. Complex graphical shapes such as fractals, were discovered as fixed points of certain set maps [1]. Informally, fractals can be treated as self similar mathematical structures which have similarity and symmetry such that considerably small parts of the shape are geometrically akin to the whole shape. Fractals are also known as expanding symmetries or unfolding symmetries. Although, fractals do not have a formal definition, however they are identified through their irregular structure that cannot be found in Euclidean geometry. Julia [5] who is considered as one of the pioneers of fractal geometry, studied iterated complex polynomials and introduced Julia set as a classical example of fractals. Let C be the complex space, T : C C be a complex polynomial of degree n 2 with complex coefficients and T i ( x ) be the i t h iterate of x. The behaviour of the iterates T i ( x ) for large i determine the Julia set (see [1,6,7,8]).
Definition 1
([1]). The set of points in C whose orbits do not converge to a point at infinity is known as filled Julia set, K T , that is,
K T = x C : { | T i ( x ) | } i = 0 is bounded .
Julia set of T denoted by J T is the boundary of filled Julia set, that is, J T = K T .
Therefore, we may say that x J T if for every neighborhood of x there exist points w and v such that T i ( w ) and T i ( v ) . The complement of a Julia set is a Fatou set.
Let p C be a fixed point of T and | ( T i ) p | = ρ , where prime denotes the complex differentiation. A point p is called a periodic point if p = T i p for some integer i 0 . Let p , T p , , T i p , be an orbit of p. The point p is called an attracting point if 0 ρ < 1 and a repelling point if ρ > 1 [6,7]. The following result gives a significant connection between repelling points of a polynomial and the Julia set.
Theorem 1
([6]). If T is a complex polynomial, then J T is the closure of the repelling periodic points of T.
Let p be an attracting fixed point of T. Then, the set A ( p ) is called the basin of attraction of p if
A ( p ) = x C : T i x p as i .
The basin of attraction of infinity, A ( ) , is defined in the same way. The following lemma is pivotal in determining Julia sets.
Lemma 1.
[7] Let p be an attracting fixed point of T. Then, J T = A ( p ) .
Thus, the Julia set is the boundary of the basin of attraction of each attracting fixed point of T, including ∞. The existence of the fixed point p for any complex polynomial is guaranteed by Brouwer fixed point theorem [9]. However, the existence of an attracting fixed point depends on the choice of the parameters. Consider the polynomial Q r ( x ) = x 2 + r . Then it has two fixed points excluding infinity. In this case, a fixed point p is attracting if | 2 p | < 1 i.e., | 1 1 4 r | < 1 . Fix v r = 1 4 r , then the set of parameters r such that Q r has an attracting fixed point is given by S = r C : | 1 v r | < 1 . Julia sets, J Q r , on the real axis i.e., r = 0 are reflection symmetric while those with complex parameter values, r C demonstrate rotational symmetry.
Mandelbrot [10] extended the idea of Julia sets and presented the notion of fractals. He investigated the graphical behaviour of connected Julia sets and plotted them for complex function, Q r ( x ) = x 2 + r , where x C is a complex variable and r C is an input parameter. He noted that various geometrical properties involving dimension, symmetry and similarity play consequential role in the study of fractal geometry.
Definition 2
([6]). Let T be any complex polynomial of degree n 2 . A Mandelbrot set M is the set consisting of all parameters r for which the Julia set, J Q r , is connected, that is,
M = r C : J Q r is connected ,
or an equivalent definition is
M = r C : { | Q r n ( 0 ) | } as n .
Mandelbrot [10,11] noted that records of heart beat, irregular coastal structures, variations of traffic flow and many naturally existing textures are examples of fractals.
In order to generate and analyze fractals, various techniques are used such as iterated function systems, random fractals, escape time criterion etc. The escape time algorithm is the stopping criterion that is based on the number of iterations necessary to determine if the orbit sequence tends to infinity or not. This algorithm provides a suitable mechanism used to demonstrate some attributes of dynamic system under iterative process. Generally, the escape criterion for Julia and Mandelbrot sets is given by:
Theorem 2
([6]). For Q r ( x ) = x 2 + r , x , r C , if there exists i 0 such that
| Q r i ( x ) | > max | r | , 2 ,
then Q r i ( x ) as i .
The term max | r | , 2 is also known as escape radius threshold. The escape radius varies in each iteration. The escape radius has a key role in visualizing the fractals.
Historically, Julia and Mandelbrot sets are investigated for the polynomials Q r but the study has been extended to quadratic, cubic, and n t h degree complex polynomials. Lakhtakia et al. [12] explored the Julia sets for general complex function of the form T ( x ) = x n + r where n N . The superior Julia and superior Mandelbrot sets for such complex polynomials in the context of noises arising in the objects were analyzed by Negi et al. [13,14]. Rochon [15] considered a more generalized form of Mandelbrot sets in bi-complex planes, see also [16,17].
Many authors have utilized various iterative processes to generate fractals. Julia and Mandelbrot sets have usually been studied for quadratic, cubic and higher degree polynomials in Picard orbit [8]. Let T : C C and x 0 C . The Picard orbit [6] is a sequence { x i } which is given by
x i + 1 = T ( x i ) ,
where i 0 .
Since the convergence of Picard process is slow, various faster converging iterative processes have been introduced to generate Julia and Mandelbrot sets. Rani and Kumar [18,19] used one-step Mann iterative process to generate superior Julia and Mandelbrot sets for n t h degree complex polynomials of the form T ( x ) = x n + r . The Mann orbit, for any x 0 C , is a sequence { x i } which is given by
x i + 1 = ( 1 α ) ( x i ) + α T ( x i ) ,
where i = 0 , 1 , and α ( 0 , 1 ] .
In 2010, a two-step Ishikawa iteration was used by Rana and Kumar [20] and Chauhan et al. [21] to study relative superior Julia and relative superior Mandelbrot sets, respectively. The dynamics of the nth order complex polynomial for non integer values were investigated in [22]. The authors also obtained new Julia and Mandelbrot sets via Ishikawa orbit. The Ishikawa orbit, for any x 0 C , is a sequence { x i } which is given by
x i + 1 = ( 1 α ) x i + α T y i , y i = ( 1 β ) x i + β T x i ,
where i = 0 , 1 , and α , β ( 0 , 1 ] .
Ashish and Rani [23] investigated the three-step Noor iteration process for Julia and Mandelbrot sets. The Noor orbit, for any x 0 C , is a sequence { x i } which is given by
x i + 1 = ( 1 α ) T x i + α T y i , y i = ( 1 β ) T x i + α T u i , u i = ( 1 γ ) T x i + γ T x i ,
where i = 0 , 1 , and α , β , γ ( 0 , 1 ] .
The modified Ishikawa process, S-iteration, was employed by Kang et al. [24,25] to study relative superior Mandelbrot sets, tricorn and multicorns. The S-orbit, for any x 0 C , is a sequence { x i } given by
x i + 1 = ( 1 α ) x i + α T y i , y i = ( 1 β ) x i + α T x i ,
where i = 0 , 1 , and α , β ( 0 , 1 ] .
Kumari et al. [26] used a four-step iterative process which is faster than of Picard, Mann and S-iteration processes and obtained some generalizations of Julia and Mandelbrot sets for quadratic, cubic and higher degree polynomials.
It is noteworthy that for each iterative process the behaviour and dynamics of the Julia and Mandelbrot sets differ. For some thought-provoking and fascinating comparisons, the reader may refer to [1,24,27,28,29] and references therein.
Complex polynomials of the form T ( x ) = x n + m x + r , where m , r C occur in various engineering problems including digital signal processing. These complex polynomials are used to determine the pole-zero plots for signals and the study of the structure and solutions of linear time invariant (LTI) state-space models, for details see [30]. Thus the study of behaviour of these polynomials and their Julia and Mandelbrot sets has gained immense interest among researchers. Kang et al. [28] introduced Julia and Mandelbrot sets in implicit Jungck Mann and Jungck Ishikawa orbits. Later, several researchers [27,29,31,32,33] employed this implicit iterative process to generate graphs of such complex polynomials. In order to achieve this, they split the polynomial T into two functions T 1 ( x ) = x n + r and T 2 ( x ) = m x . However, the Jungck iterative process and its variants are used to determine the common fixed points of two mappings. Therefore, the question arises whether we can obtain an escape criterion and generate fractals for polynomials of the form T using explicit iterative processes.
The purpose of this paper is to answer this question. In this paper, we discuss the graphical behaviour of the complex polynomial of the form T ( x ) = x n + m x + r where m , r C and n 2 using Picard Ishikawa type fixed point iteration process for the generation of fractals. Note that the Julia and Mandelbrot sets generated have distinctive shapes for the proposed iterative process as compared to already present iterative processes in the literature. Further, we show the effect of change of parameters on color variation and graph of the sets.
The Picard Ishikawa type iteration process was introduced by Piri et al. [34]. They claimed that this iterative process converges faster than Mann and Ishikawa iteration processes. Let D be a subset of a Banach space and f : D D then the three step iteration process is given by
x 1 = x D , x i + 1 = ( 1 α i ) y i + α i f y i , y i = f z i , z i = f ( ( 1 β i ) x i + β i f x i ) , i 0 ,
where α i , β i ( 0 , 1 ] .

2. Main Results

In this section, we use a Picard Ishikawa type iterative process and some prove escape criterions to determine the escape radius for this process. Throughout this paper we assume that for any complex polynomial the parameters are chosen in a way that at the least one attracting fixed point exists.
Let C be a complex space and T C : C C be a complex polynomial with complex coefficients. The Picard Ishikawa type orbit around any x 0 C , is a sequence { x i } given by
x i + 1 = ( 1 α ) y i + α T C y i , y i = T C z i , z i = T C t i , t i = ( 1 β ) x i + β T C x i ,
where i = 0 , 1 , 2 , and α , β ( 0 , 1 ] .
We need the following escape criterions for the quadratics, cubic and higher degree polynomials.

2.1. Escape Criterion for Quadratic Complex Polynomials in a Picard Ishikawa Type Orbit

For the quadratic polynomial T C ( x ) = x 2 + m x + r where m , r C , we have the following result.
Theorem 3.
Suppose that | x | | r | > max 2 ( 1 + | m | ) α , 2 ( 1 + | m | ) β , α , β ( 0 , 1 ] . Define { x i } i N as in (2) where x 0 = x , y 0 = y , z 0 = z and t 0 = t . Then, | x i | as i .
Proof. 
As, T C ( x ) = x 2 + m x + r . From (2), we have
| t | = | ( 1 β ) x + β T C x | = | ( 1 β ) x + β ( x 2 + m x + r ) | | ( 1 β ) x + β ( x 2 + m x ) | β | r | .
The assumption | x | | r | yields
| t | | ( 1 β ) x + β ( x 2 + m x ) | β | x | β | x 2 | ( 1 β + β | m | ) | x | β | x | = β | x 2 | ( 1 + β | m | ) | x | = | x | β | x | ( 1 + β | m | ) .
Since β 1 , we obtain ( 1 + β | m | ) > ( 1 + | m | ) which implies that
| t | | x | β | x | ( 1 + | m | ) .
Thus, we have
| t | | x | ( 1 + | m | ) β | x | 1 + | m | 1 .
Therefore,
| t | | t | ( 1 + | m | ) | x | β | x | 1 + | m | 1 .
From our assumption; | x | > max 2 ( 1 + | m | ) α , 2 ( 1 + | m | ) β , we get
β | x | 1 + | m | 1 > 1 .
Now, (3) gives that
| t | > | x | .
As z = z 0 , (2) gives
| z | = | T C ( t ) | = | t 2 + m t + r | | t 2 + m t | | r | .
Since β 1 , it follows from (5) and assumption | x | | r | that
| z | | t 2 + m t | | x | β | t 2 | | m | | t | | t | = | t | β | t | ( 1 + | m | ) ,
which further implies that
| z | | z | ( 1 + | m | ) | t | β | t | ( 1 + | m | ) 1 .
Using (4) and (5) we have
| t | > | x | β | t | 1 + | m | > β | x | 1 + | m | β | t | 1 + | m | 1 > β | x | 1 + | m | 1 > 1 .
Consequently, (5)–(7) yield
| z | > | x | .
Moreover, let y = y 0 , | y | = | T C ( z ) | = | z 2 + m z + r | . Then, by an assumption | x | | r | , (8) and the fact that β 1 we obtain
| y | | z 2 + m z | | r | β | z | 2 | m | | z | | z | | z | β | z | ( 1 + | m | ) .
This implies
| y | | z | β | z | 1 + | m | 1 .
From (4) and (8) we obtain
| y | | x | β | x | 1 + | m | 1 > | x | .
Finally, we have
| x 1 | = | ( 1 α ) y + α T C ( y ) | = | ( 1 α ) y + α ( y 2 + m y + r ) | .
Furthermore, from | x | | r | and (9) we get that
| x 1 | = | ( 1 α ) y + α ( y 2 + m y + r ) | α | y 2 | ( 1 α + α | m | ) | y | α | r | α | y 2 | ( 1 α + α | m | ) | y | α | y | = α | y 2 | ( 1 + α | m | ) | y | = | y | α | y | ( 1 + α | m | ) .
As α 1 , we obtain
| x 1 | | y | α | y | ( 1 + α | m | ) | y | α | y | ( 1 + | m | ) | y | ( 1 + | m | ) α | y | ( 1 + | m | ) 1 .
By (9), we have
| x 1 | | x | α | x | 1 + | m | 1 .
From our given assumption, we have | x | > 2 ( 1 + | m | ) α and hence α | x | 1 + | m | 1 > 1 . Thus, there exists a real number ρ > 0 such that
α | x | 1 + | m | 1 > 1 + ρ .
It follows that
| x 1 | > ( 1 + ρ ) | x | .
In particular, | x 1 | > | x | . Continuing in the same manner yields
| x i | > ( 1 + ρ ) i | x | .
Therefore, the orbit of x tends to infinity. □
The following corollary is the refinement of the Theorem 3.
Corollary 1.
Suppose that | x i | > max | r | , 2 ( 1 + | m | ) α , 2 ( 1 + | m | ) β where α , β ( 0 , 1 ] then | x i | as i .

2.2. Escape Criterion for Cubic Complex Polynomials in a Picard Ishikawa Type Orbit

For the cubic polynomial T C ( x ) = x 3 + m x + r where m , r C , we have the following result.
Theorem 4.
Suppose | x | | r | > max 2 ( 1 + | m | ) α 1 2 , 2 ( 1 + | m | ) β 1 2 , α , β ( 0 , 1 ] . Define a sequence { x i } i N as in (2) where x 0 = x , y 0 = y , z 0 = z and t 0 = t . Then, | x i | as i .
Proof. 
As T C ( x ) = x 3 + m x + r , from (2) we have
| t | = | ( 1 β ) x + β T C ( x ) | = | ( 1 β ) x + β ( x 3 + m x + r ) | | ( 1 β ) x + β ( x 3 + m x ) | β | r | .
The assumption | x | | r | yields that
| t | | ( 1 β ) x + β ( x 3 + m x ) | β | x | β | x 3 | ( 1 β + β | m | ) | x | β | x | = β | x 3 | ( 1 + β | m | ) | x | = | x | β | x 2 | ( 1 + β | m | ) .
As β 1 ,
| t | | x | β | x 2 | ( 1 + | m | ) .
Therefore,
| t | | t | ( 1 + | m | ) | x | β | x 2 | 1 + | m | 1 .
The assumption, | x | > max 2 ( 1 + | m | ) α 1 2 , 2 ( 1 + | m | ) β 1 2 implies that
β | x 2 | 1 + | m | 1 > 1 .
It follows from (10) that
| t | > | x | .
As z = z 0 , by (2) we have
| z | = | T C ( t ) | | t 3 + m t | | r | .
As β 1 , from (12) and assumption | x | | r | we obtain
| z | | t 3 + m t | | x | = | t | β | t 2 | ( 1 + | m | )
which further implies that
| z | | t | β | t 2 | ( 1 + | m | ) 1 .
Now by (12) and (11), we have
β | t | 2 1 + | m | 1 β | x | 2 1 + | m | 1 > 1 .
Consequently, (5), (20) and (14) imply that
| z | > | x | .
Also, y = y 0 , | y | = | T C ( z ) | = | z 3 + m z + r | . Then, the given assumption | x | | r | , (8) and the fact that β 1 yield
| y | | z 3 + m z | | r | | z | β | z 2 | ( 1 + | m | ) .
Thus
| y | | z | β | z 2 | 1 + | m | 1 .
From (11) and (15), we obtain
| y | | x | β | x 2 | 1 + | m | 1 > | x | .
Lastly, we have
| x 1 | = | ( 1 α ) y + α T C y | = | ( 1 α ) y + α ( y 3 + m y + r ) | .
From | x | | r | , (16) and α 1 , we have
| x 1 | = | ( 1 α ) y + α ( y 3 + m y + r ) | α | y 3 | ( 1 α + α | m | ) | y | α | y | = α | y 2 | ( 1 + α | m | ) | y | | y | α | y 2 | ( 1 + | m | ) | y | ( 1 + | m | ) α | y 2 | ( 1 + | m | ) 1 .
From (16), we have
| x 1 | | x | α | x 2 | 1 + | m | 1 .
By our assumption we have | x | > 2 ( 1 + | m | ) α 1 2 and hence α | x 2 | 1 + | m | 1 > 1 . Thus, there exists a real number ρ > 0 such that
α | x 2 | 1 + | m | 1 > 1 + ρ .
It follows that
| x 1 | > ( 1 + ρ ) | x | .
Continuing in the same manner, we obtain
| x i | > ( 1 + ρ ) i | x | .
Therefore, the orbit of x tends to infinity. □
The following corollary is the refinement of the Theorem 4.
Corollary 2.
Suppose that | x i | > max | r | , 2 ( 1 + | m | ) α 1 2 , 2 ( 1 + | m | ) β 1 2 where α , β ( 0 , 1 ] then | x i | as i .

2.3. Escape Criterion for General Complex Polynomials in a Picard Ishikawa Type Orbit

For the general complex polynomial T C ( x ) = x n + m x + r where m , r C , we have the following result.
Theorem 5.
Suppose | x | | r | > max 2 ( 1 + | m | ) α 1 n 1 , 2 ( 1 + | m | ) β 1 n 1 , with n 2 and α , β ( 0 , 1 ] . Define a sequence { x i } i N as in (2) where x 0 = x , y 0 = y , z 0 = z and t 0 = t . Then, | x i | as i .
Proof. 
Let T C ( x ) = x n + m x + r . Note that (2), assumptions | x | | r | and β 1 give
| t | = | ( 1 β ) x + β T C ( x ) | | ( 1 β ) x + β ( x n + m x ) | β | r | β | x n | ( 1 β + β | m | ) | x | β | x | = | x | β | x n 1 | ( 1 + β | m | ) | x | β | x n 1 | ( 1 + | m | ) .
Therefore,
| t | | x | β | x n 1 | 1 + | m | 1 .
By our assumption, we have | x | > 2 ( 1 + | m | ) β 1 n 1 and hence
β | x n 1 | 1 + | m | 1 > 1 .
It follows from (17) that
| t | > | x | .
Since z = z 0 , so from (2) we obtain
| z | | t n + m t | | r | .
As β 1 , from (19) and assumption | x | | r | , we have
| z | | t | β | t n 1 | ( 1 + | m | ) 1 .
Now by (18) and (19), we have
β | t | n 1 1 + | m | 1 > 1 .
Hence,
| z | > | x | .
As y = y 0 , | y | = | T C ( z ) | = | z n + m z + r | , so using the similar arguments as before we obtain
| y | | x | β | x n 1 | 1 + | m | 1 > | x | .
Also, from | x | | r | , (22), and α 1 we have
| x 1 | = | ( 1 α ) y + α ( y n + m y + r ) | α | y n | ( 1 α + α | m | ) | y | | r | = α | y 2 | ( 1 + α | m | ) | y | = | y | α | y 2 | ( 1 + α | m | ) | x | α | x 2 | 1 + | m | 1 .
Furthermore, from our assumption we have | x | > 2 ( 1 + | m | ) α 1 n 1 and thus α | x n 1 | 1 + | m | 1 > 1 . Thus, there exists a real number ρ > 0 such that
α | x n 1 | 1 + | m | 1 > 1 + ρ .
Finally, we obtain
| x 1 | > ( 1 + ρ ) | x | .
Now, continuing this process
| x i | > ( 1 + ρ ) i | x | .
Therefore, the orbit of x tends to infinity. □
The following corollary is the refinement of the Theorem 5.
Corollary 3.
Suppose that | x i | > max | r | , 2 ( 1 + | m | ) α 1 n 1 , 2 ( 1 + | m | ) β 1 n 1 where n 2 and α , β ( 0 , 1 ] then | x i | as i .
Theorem 6.
Suppose that { x i } i N { 0 } is a sequence in the Picard Ishikawa type orbit for the complex polynomial T C ( x ) = x n + m x + r where m , r C with n 2 such that | x i | as i , then | x | | r | > 2 ( 1 + | m | ) α 1 n 1 and | x | | r | > 2 ( 1 + | m | ) β 1 n 1 , α , β ( 0 , 1 ] .
Proof. 
Let { x i } i N be a sequence in Picard Ishikawa type orbit. First, we prove that | x | | r | . According to hypothesis, | x i | as i , the sequence { | x i | } must be unbounded. Hence, | x i | | r | for all i N { 0 } and therefore | x | | r | . Let T C ( x ) = x n + m x + r , m , r C where t 0 = t , x 0 = x , y 0 = y and z 0 = z , then | x | | r | implies that
| t | = | ( 1 β ) x + β T C x | = | ( 1 β ) x + β ( x n + m x + r ) | | β x n | + ( ( 1 β ) + m β ) x | β | r | β | x n | ( ( 1 β ) + | m | β ) | x | β | x | β | x n | ( 1 + | m | β ) | x | .
Thus,
| t | | x | ( β | x n 1 | ( 1 + | m | ) ) = | x | ( 1 + | m | ) β | x n 1 | 1 + | m | 1
implies
| t | | x | β | x n 1 | 1 + | m | 1 .
Here, we have two possibilities; either β | x n 1 | 1 + | m | 1 1 or β | x n 1 | 1 + | m | 1 > 1 . If β | x n 1 | 1 + | m | 1 1 we have
β | x n 1 | 1 + | m | 2
which implies that
| x n 1 | 2 ( 1 + | m | ) β
and hence
| x | 2 ( 1 + | m | ) β 1 n 1 ,
a contradiction. Indeed, { | x i | } is not bounded where i N { 0 } . Therefore, we must have β | x n 1 | 1 + | m | 1 > 1 . Thus, | x | > 2 ( 1 + | m | ) β 1 n 1 . Now, inequality (23) implies that
| t | | x | β | x n 1 | 1 + | m | 1 > | x | .
Furthermore, β 1 and | x | | r | give
| z | = | T C ( t ) | | t n + m t | | r | β | t n | | m | | t | | x | = | t | ( β | t n 1 | | ( 1 + | m | ) ) .
As β | x n 1 | ( 1 + | m | ) 1 > 1 , so we have
| t | > | x | β | x n 1 | 1 + | m | 1 > | x | .
As a consequence we obtain
| z | | x | β | x n 1 | 1 + | m | 1 ( 1 + | m | ) .
Thus,
| z | | x | β | x n 1 | 1 + | m | 1 > | x | .
Similarly, | y | = | T C ( z ) | = | z n + m z + r | , | x | > | r | and β 1 imply that
| y | | x | β | x n 1 | 1 + | m | 1 .
Consequently,
| y | | x | β | x n 1 | 1 + | m | 1 > | x | .
Finally, we have
| x 1 | = | ( 1 α ) y + α T C ( y ) | = | ( 1 α ) y + α ( y n + m y + r ) | α | y n | ( 1 α + α | m | ) | y | α | r | α | y n | ( 1 α + α | m | ) | y | α | y | α | y n | ( 1 + α | m | ) | y | α | y n | ( 1 + | m | ) | y | = | y | ( α | y n 1 | ( 1 + | m | ) ) | x | ( α | x n 1 | ( 1 + | m | ) ) .
Hence
| x 1 | | x | α | x n 1 | 1 + | m | 1 .
Using arguments similar to those as before, we only have one possibility that α | x n 1 | 1 + | m | 1 > 1 . Therefore, | x | > 2 ( 1 + | m | ) α 1 n 1 . This completes the proof. □

3. Visualization of Fractals

In this section, we present some Julia and Mandelbrot sets for quadratic and higher order polynomials. We found several captivating new fractals having various geometric shapes. However, we have chosen some figures. The color variation occurs due to the change of input parameters. We have also investigated the effect of change of parameters α and β on the shape and the variation of colors. The number of iterations was fixed at 10.

3.1. Generation of Julia sets

Following Algorithm 1 is the pseudocode for the generation of Julia sets. Note that T ( z ) represents the iteration process.
Algorithm 1: Generation of Julia Set.
Symmetry 12 00086 i001
Now, we present quadratic, cubic and septic Julia sets in Picard Ishikawa type orbit for the complex polynomial, T C ( x ) = x n + m x + r .
  • For Figure 1, we consider the polynomial T ( x ) = x 2 + ( 0.5 + 0.7 i ) x + ( 0.01 + 0.18 i ) and A = [ 2.5 , 2.5 ] × [ 2.1 , 2.1 ] . It is easy to see that T has one attracting fixed point, p = 0.1427 + 0.1019 i . Observe that for α = 0.2 , β = 0.097 and α = 0.11 , β = 0.18 we obtain different images due to color variation caused by parameters. It is interesting to note that for α = 1 , β = 1 and α = 10 10 , β = 10 10 we have similar shapes but there is clear variation of colors.
  • For Figure 2, we consider the polynomial T ( x ) = x 3 + ( 0.275 + 0.5 i ) x + ( 0.559 + 0.35 i ) and A = [ 1.5 , 1.5 ] × [ 1.8 , 1.8 ] . The polynomial T has attracting fixed point p 0.6434 + 0.2687 i in A. Note that the cubic Julia sets for α = 0.08 and β = 0.09 have more color variation as compared to the Julia sets for α = 0.1 , and β = 0.2 . Again, for α = 1 , β = 1 and α = 10 10 , β = 10 10 the shapes are same but there is variability in colors.
  • For Figure 3, we input T ( x ) = x 7 + ( 0.23 + 1.2 i ) x + ( 0.5 + 0.7 i ) and A = [ 1.3 , 1.3 ] 2 . The attracting fixed point of the polynomial is p 0.2391 + 0.5835 i . We can see that for α = 0.01 and β = 0.08 the shape is spread and stretched while the shape is dense and neatly packed for α = 0.1 and β = 0.05 . Note the variation of colors in figures (C) and (D) as well.

3.2. Generation of Mandelbrot Sets

Following Algorithm 2 is the pseudocode for the generation of Mandelbrot sets. Note that T ( z ) represents the iteration process.
Algorithm 2: Generation of Mandelbrot set.
Symmetry 12 00086 i002
For Figure 4 we input A = [ 2 , 2 ] × [ 1.2 , 2.5 ] and observe that for α = 0.1 and β = 0.3 , the shape is stretched and the bulb is wider and for α = 0.75 and β = 0.7 the shape is compact with defined bulb. Notice the variation of colors for Mandelbrot sets for α = 1 , β = 1 and α = 0.009 , β = 0.009 . Also, observe that Mandelbot sets generated are symmetric about origin.

4. Conclusions

In this paper, a Picard Ishikawa type orbit was used to study the behaviour of complex poylnomials. We obtained escape criterions for complex quadratic, cubic and higher degree polynomials. Some alluring Julia and Mandelbrot sets have been generated. We also observed that the variation of parameters has shown eminent changes in the Julia and Mandelbrot sets. Our results are different from comparable existing results as we obtain escape criterion and fractals for polynomials of the form T ( x ) = x n + m x + r where m , r C without using the Jungck iterative process. It is also worth mentioning that the behaviour of the polynomial and shape of the fractal generated under the iterative process (2) is different and unique as compared to the iterative process studied before in the literature [1,24,29,32].

Author Contributions

Conceptualization, M.A. and H.I.; methodology, M.A. and H.I.; validation, M.A. and M.D.l.S., formal analysis H.I.; investigation, H.I. and M.A.; writing—original draft preparation, H I.; writing—review and editing, H.I.; visualization, M.D.l.S.; supervision, M.A. and M.D.l.S.; project administration, M.D.l.S.; funding acquisition, M.D.l.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Basque Government through grant number IT 1207-19.

Acknowledgments

All the authors are grateful to the referees for their critical remarks and valuable suggestions which helped to improve the presentation of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barnsley, M. Fractals Everywhere, 2nd ed.; Academic Press: San Diego, CA, USA, 1993. [Google Scholar]
  2. Hundertmark-Zauôkovù, A. On the convergence of fixed point iterations for the moving geometry in a fluid-structure interaction problem. J. Differ. Equ. 2019, 267, 7002–7046. [Google Scholar] [CrossRef]
  3. Rahmani, M.; Koutsopoulos, H.N.; Jenelius, E. Travel time estimation from sparse floating car data with consistent path inference: A fixed point approach. Transp. Res. Part C Emerg. Technol. 2017, 85, 628–643. [Google Scholar] [CrossRef]
  4. Strogatz, S.H. Nonlinear Dynamics and Chaos With Applications to Physics, Biology, Chemistry, and Engineering, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  5. Julia, G. Memoire sur l’iteration des functions rationnelles. J. Math. Pures Appl. 1918, 8, 737–747. [Google Scholar]
  6. Devaney, R.L. A First Course in Chaotic Dynamical Systems: Theory and Experiment, 2nd ed.; Addison-Wesley: Boston, MA, USA, 1992. [Google Scholar]
  7. Falconer, K. Fractal Geometry: Mathematical Foundations and Applications, 2nd ed.; John Wiley & Sons: Chichester, UK, 2004. [Google Scholar]
  8. Frame, M.; Robertson, J. A generalized Mandelbrot set and the role of critical points. Comput. Graph. 1992, 16, 35–40. [Google Scholar] [CrossRef]
  9. Brouwer, L.E.J. Über Abbildungen von Mannigfaltigkeiten. Math. Ann. 1912, 71, 97–115. [Google Scholar] [CrossRef] [Green Version]
  10. Mandelbrot, B.B. The Fractal Geometry of Nature; W.H. Freeman: New York, NY, USA, 1983; Volume 2. [Google Scholar]
  11. Debnath, L. A brief historical introduction to fractals and fractal geometry. Int. J. Math. Educ. Sci. Technol. 2006, 37, 29–50. [Google Scholar] [CrossRef]
  12. Lakhtakia, A.; Varadan, W.; Messier, R.; Varadan, V.K. On the symmetries of the Julia sets for the process zp + c. J. Phys. A Math. Gen. 1987, 20, 3533–3535. [Google Scholar] [CrossRef]
  13. Negi, A.; Rani, M. Midgets of superior Mandelbrot set. Chaos Solitons Fract. 2008, 36, 237–245. [Google Scholar] [CrossRef]
  14. Negi, A.; Rani, M. A new approach to dynamic noise on superior Mandelbrot set. Chaos Solitons Fract. 2008, 36, 1089–1096. [Google Scholar] [CrossRef]
  15. Rochon, D. A generalized Mandelbrot set for bicomplex numbers. Fractals 2000, 8, 355–368. [Google Scholar] [CrossRef]
  16. Wang, X.; Sun, Y. The general quaternionic M-J sets on the mapping zzα + c(α ∈ ℕ). Comput. Math. Appl. 2007, 53, 1718–1732. [Google Scholar] [CrossRef] [Green Version]
  17. Rani, M. Theoretical Framework for Fractal Models under Two-Step Feedback Process. Ph.D. Thesis, Kumaun University, Nainital, India, 2016. [Google Scholar]
  18. Rani, M.; Kumar, V. Superior Julia set. Res. Math. Educ. 2004, 8, 261–277. [Google Scholar]
  19. Rani, M.; Kumar, V. Superior Mandelbrot set. Res. Math. Educ. 2004, 8, 279–291. [Google Scholar]
  20. Rani, M.; Chauhan, Y.S.; Negi, A. Non linear dynamics of Ishikawa iteration. Int. J. Comput. Appl. 2010, 7, 43–49. [Google Scholar] [CrossRef]
  21. Chauhan, Y.S.; Rana, R.; Negi, A. New Julia sets of Ishikawa iterates. Int. J. Comput. Appl. 2010, 7, 34–42. [Google Scholar] [CrossRef]
  22. Chauhan, Y.S.; Rana, R.; Negi, A. Complex dynamics of Ishikawa iterates for non integer values. Int. J. Comput. Appl. 2010, 9, 9–16. [Google Scholar] [CrossRef]
  23. Ashish; Rani, M.; Chugh, R. Julia sets and Mandelbrot sets in Noor orbit. Appl. Math. Comput. 2014, 228, 615–631. [Google Scholar]
  24. Kang, S.M.; Rafiq, A.; Latif, A.; Shahid, A.A.; Ali, F. Fractals through modified iteration scheme. Filomat 2016, 30, 3033–3046. [Google Scholar] [CrossRef] [Green Version]
  25. Kang, S.M.; Rafiq, A.; Latif, A.; Shahid, A.A.; Kwun, Y.C. Tricorns and multicorns of S-iteration scheme. J. Funct. Spaces 2015, 2015, 417167. [Google Scholar]
  26. Kumari, M.; Ashish, R.C. New Julia and Mandelbrot sets for a new faster iterative process. Int. J. Pure Appl. Math. 2016, 107, 161–177. [Google Scholar] [CrossRef] [Green Version]
  27. Kang, S.M.; Nazeer, W.; Tanveer, M.; Shahid, A.A. New fixed point results for fractal generation in Jungck Noor orbit with s-convexity. J. Funct. Spaces 2015, 2015, 963016. [Google Scholar] [CrossRef] [Green Version]
  28. Kang, S.M.; Rafiq, A.; Tanveer, M.; Ali, F.; Kwun, Y.C. Julia and Mandelbrot sets in modified Jungck three-step orbit. Wulfenia J. 2015, 22, 167–185. [Google Scholar]
  29. Kwun, Y.C.; Tanveer, M.; Nazeer, W.; Abbas, M.; Kang, S.M. Fractal generation in modified Jungck-S orbit. IEEE Access 2019, 7, 35060–35071. [Google Scholar] [CrossRef]
  30. Proakis, J.G.; Manolakis, D.G. Digital Signal Processing: Principles, Algorithms and Applications, 4th ed.; Pearson: Bengaluru, India, 2007. [Google Scholar]
  31. Mishra, M.K.; Ojha, D.B.; Sharma, D. Fixed point results in tricorn and multicorns of Ishikawa iteration and s-convexity. IJEST 2011, 2, 157–160. [Google Scholar]
  32. Cho, S.Y.; Shahid, A.A.; Nazeer, W.; Kang, S.M. Fixed point results for fractal generation in noor orbit and s-convexity. SpringerPlus 2016, 5, 1843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Nazeer, W.; Kang, S.M.; Tanveer, M.; Shahid, A.A. Fixed point results in the generation of Julia and Mandelbrot sets. J. Inequal. Appl. 2015, 2015, 298. [Google Scholar] [CrossRef] [Green Version]
  34. Piri, H.; Daraby, B.; Rahrovi, S.; Ghasemi, M. Approximating fixed points of generalized α-nonexpansive mappings Banach spaces by new faster iteration process. Numer. Algorithms 2018, 81, 1129–1148. [Google Scholar] [CrossRef]
Figure 1. Quadratic Julia sets.
Figure 1. Quadratic Julia sets.
Symmetry 12 00086 g001
Figure 2. Cubic Julia sets.
Figure 2. Cubic Julia sets.
Symmetry 12 00086 g002
Figure 3. Septic Julia sets.
Figure 3. Septic Julia sets.
Symmetry 12 00086 g003
Figure 4. Mandelbrot sets.
Figure 4. Mandelbrot sets.
Symmetry 12 00086 g004

Share and Cite

MDPI and ACS Style

Abbas, M.; Iqbal, H.; De la Sen, M. Generation of Julia and Mandelbrot Sets via Fixed Points. Symmetry 2020, 12, 86. https://doi.org/10.3390/sym12010086

AMA Style

Abbas M, Iqbal H, De la Sen M. Generation of Julia and Mandelbrot Sets via Fixed Points. Symmetry. 2020; 12(1):86. https://doi.org/10.3390/sym12010086

Chicago/Turabian Style

Abbas, Mujahid, Hira Iqbal, and Manuel De la Sen. 2020. "Generation of Julia and Mandelbrot Sets via Fixed Points" Symmetry 12, no. 1: 86. https://doi.org/10.3390/sym12010086

APA Style

Abbas, M., Iqbal, H., & De la Sen, M. (2020). Generation of Julia and Mandelbrot Sets via Fixed Points. Symmetry, 12(1), 86. https://doi.org/10.3390/sym12010086

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop