Open Journal of Mathematical Sciences
ISSN: 2523-0212 (Online) 2616-4906 (Print)
DOI: 10.30538/oms2018.0017
Fixed point results for the complex fractal generation in the S-iteration orbit with s-convexity
Krzysztof Gdawiec\(^{1}\), Abdul Aziz Shahid
Institute of Computer Science, University of Silesia, B¸edzi´nska 39, 41-200 Sosnowiec, Poland. (K.G)
Department of Mathematics and Statistics, University of Lahore, Lahore, Pakistan.(A.A.S)
\(^{1}\)Corresponding Author: kgdawiec@ux2.math.us.edu.pl
Abstract
Keywords:
1. Introduction
Mandelbrot and Julia sets are some of the best known illustrations of a highly complicated chaotic systems generated by a very simple mathematical process. They were introduced by Benoit Mandelbrot in the late 1970's [1], but Julia sets were studied much earlier, namely in the early 20th century by French mathematicians Pierre Fatou and Gaston Julia. Mandelbrot working at IBM has studied their works and plotted the Julia sets for \(z^2 + c\) and corresponding to them the Mandelbrot set. He was surprised by the result that he obtained. Since then many mathematicians have studied different properties of Mandelbrot and Julia sets and proposed various generalizations of those sets. The first and the most obvious generalization was the use of \(z^p + c\) function instead of the quadratic one used by Mandelbrot [2, 3]. Then some other types of functions were studied: rational [4], transcendental [5], elliptic [6], anti-polynomials [7] etc. Another step in the studies on Mandelbrot and Julia sets was the extension from complex numbers to other algebras, e.g., quaternions [8], octonions [9], bicomplex numbers [10] etc.Another interesting generalization of Mandelbrot and Julia sets is the use of the results from fixed point theory. In the fixed point theory there exist many approximate methods of finding fixed points of a given mapping, that are based on the use of different feedback iteration processes. These methods can be used in the generalization of Mandelbrot and Julia sets. In 2004, Rani and Kumar [11, 12] introduced superior Julia and Mandelbrot sets using Mann iteration scheme. Chauhan et al. in [13] introduced the relative superior Julia sets using Ishikawa iteration scheme. Also, relative superior Julia sets, Mandelbrot sets and tricorn, multicorns by using the \(S\)-iteration scheme were presented in [14, 15]. Recently, Ashish et al. in [16] introduced Julia and Mandelbrot sets using the Noor iteration scheme, which is a three-step iterative procedure. The junction of a \(s\)-convex combination [17] and various iteration schemes was studied in many papers. Mishra et al. [18, 19] developed fixed point results in relative superior Julia sets, tricorn and multicorns by using the Ishikawa iteration with \(s\)-convexity. In [20] Kang et al. introduced new fixed point results for fractal generation using the implicit Jungck-Noor orbit with \(s\)-convexity, whereas Nazeer et al. in [21]NKTS15} used the Jungck-Mann and Jungck-Ishikawa iterations with \(s\)-convexity. The use of Noor iteration and \(s\)-convexity was shown in [22].
In this paper, we present some fixed point results for Julia and Mandelbrot sets by using the \(S\)-iteration scheme with \(s\)-convexity. We derive the escape criteria for quadratic, cubic and the \((k+1)\)th degree complex polynomial.
The remainder of this paper is outlined as follows. In Sec. 2 we present some basic definitions used in the paper. Next, in Sec.3, we present the main results of this paper, namely we prove the escape criteria for the quadratic, cubic and the \((k+1)\)th degree complex polynomial in the \(S\)-iteration with \(s\)-convexity. Some graphical examples of complex fractals generated using the escape criteria are presented in Sec.4. The paper ends with some concluding remarks and future work (Sec.5).
2. Preliminaries
Definition2.1. (see [23], Julia set) Let \(f: \mathbb{C} \rightarrow \mathbb{C}\) be a polynomial of degree \(\geq 2\). Let \(F_{f}\) be the set of points in \(\mathbb{C}\) whose orbits do not converge to the point at infinity, i.e., \(F_{f} = \{z \in \mathbb{C}: \{\left\vert f^{n}(z)\right\vert \}_{n=0}^{\infty} \text{ is bounded} \}.\) \(F_{f}\)\ is called as filled Julia set of the polynomial \(f\). The boundary points of \(F_{f}\) are called the points of Julia set of the polynomial \(f\) or simply the Julia set.
Definition 2.2.(see [24], Mandelbrot set) The Mandelbrot set \(M\) consists of all parameters \(c\) for which the filled Julia set of \(Q_{c}(z)=z^{2}+c\) is connected, i.e.,
Definition 2.3. Let \(C \subset \mathbb{C}\) be a non-empty set and \(f:C\rightarrow C\). For any point \(z_{0} \in C\) the Picard orbit is defined as the set of iterates of the point \(z_{0}\), i.e.,
Definition 2.4.(see [25], \(S\)-iteration orbit) Consider a sequence \(\{z_{n}\}\) of iterates for initial point \(z_{0}\in C\) such that
Theorem 2.5.[Escape Criterion for \(S\)-iteration] Let \(Q_c(z) = z^{k+1} + c\), where \(k = 1, 2, 3, \ldots\) and \(c \in \mathbb{C}\). Iterate \(Q_c\) using \eqref{eq:s} with \(\gamma_n = \gamma, \delta_n = \delta\), where \(\gamma, \delta \in (0, 1]\). Suppose that
3. Main results
In each of the two steps of \(S\)-iteration we use a convex combination of two elements. In the literature we can find some generalizations of the convex combination. One of such generalizations is the \(s\)-convex combination.Definition 3.1.[\(s\)-convex combination [17]] Let \(z_1, z_2, \ldots, z_n \in \mathbb{C}\) and \(s \in (0, 1]\). The \(s\)-convex combination is defined in the following way:
3.1. Escape criterion for quadratic function
Theorem 3.2. Assume that \(|z|\geq |c|>\frac{2}{s\gamma }\) and \(|z|\geq |c|>\frac{2}{s\delta }\), where \(\gamma, \delta, s \in (0, 1]\) and \(c \in \mathbb{C}\). Let \(z_0 = z\). Then for (7) with \(Q_c(z) = z^2 + c\) we have \(|z_{n}| \rightarrow \infty\) as \(n \rightarrow \infty\).
Proof. Consider \[ \left\vert w\right\vert =\left\vert (1-\delta )^{s}z+\delta^{s}Q_{c}(z)\right\vert. \] For \(Q_{c}(z)=z^{2}+c\), \begin{eqnarray*} \left\vert w\right\vert &=&\left\vert (1-\delta )^{s}z+\delta^{s}(z^{2}+c)\right\vert \\ &=&\left\vert (1-\delta )^{s}z+(1-(1-\delta ))^{s}(z^{2}+c)\right\vert. \end{eqnarray*} By binomial expansion upto linear terms of \(\delta \) and \((1-\delta ),\) we obtain
Corollary 3.3. Suppose that
Corollary 3.4.[Escape Criterion] Let \(\gamma, \delta, s \in (0, 1]\). Suppose that
3.2. Escape criterion for cubic function
Theorem 3.5. Assume that \(|z| \geq |c|> (\frac{2}{s\gamma })^{\frac{1}{2}}\) and \(|z|\geq |c| > (\frac{2}{s\delta })^{\frac{1}{2}}\), where \(c \in \mathbb{C}\) and \(\gamma, \delta, s \in (0, 1]\). Let \(z_0 = z\). Then for (7) with \(Q_c(z) = z^3 + c\) we have \(\left\vert z_{n}\right\vert \rightarrow \infty\) as \(n \rightarrow \infty\).
Proof. Consider \[ \left\vert w\right\vert =\left\vert (1-\delta )^{s}z+\delta^{s}Q_{c}(z)\right\vert. \] For \(Q_{c}(z)=z^{3}+c\), we have \begin{eqnarray*} \left\vert w\right\vert &=&\left\vert (1-\delta )^{s}z+\delta^{s}(z^{3}+c)\right\vert \\ &=&\left\vert (1-\delta )^{s}z+(1-(1-\delta ))^{s}(z^{3}+c)\right\vert. \end{eqnarray*} By binomial expansion upto linear terms of \(\delta \) and \((1-\delta ),\) we obtain
Corollary 3.6. Suppose that
Corollary 3.7.[Escape Criterion] Let \(\gamma, \delta, s \in (0, 1]\). Suppose that
3.3. A general escape criterion
Theorem 3.8. Assume that \(|z|\geq |c|>(\frac{2}{s\gamma })^{\frac{1}{k+1}}\) and \(|z|\geq |c|>(\frac{2}{s\delta })^{\frac{1}{k+1}}\), where \(k = 1, 2, \ldots\), \(c \in \mathbb{C}\) and \(\gamma, \delta, s \in (0, 1]\). Let \(z_0 = z\). Then for (7) with \(Q_c(z) = z^{k+1} + c\) we have \(\left\vert z_{n}\right\vert \rightarrow \infty \) as \(n \rightarrow \infty\).
Proof. Consider \[ \left\vert w\right\vert =\left\vert (1-\delta )^{s}z+\delta^{s}Q_{c}(z)\right\vert. \] For \(Q_{c}(z)=z^{k+1}+c\), we have \begin{eqnarray*} \left\vert w\right\vert &=&\left\vert (1-\delta )^{s}z+\delta^{s}(z^{k+1}+c)\right\vert \\ &=&\left\vert (1-\delta )^{s}z+(1-(1-\delta ))^{s}(z^{k+1}+c)\right\vert. \end{eqnarray*} By binomial expansion upto linear terms of \(\delta \) and \((1-\delta ),\) we obtain
Corollary3.9.[Escape Criterion] Let \(\gamma, \delta, s \in (0, 1]\). Suppose that
4. Generation of Mandelbrot and Julia sets
In this section we present some graphical examples of Mandelbrot and Julia sets. To generate the images we used the escape time algorithms which were implemented in Mathematica 9.0. Pseudocode of the Mandelbrot set generation algorithm is presented in Algorithm 1, whereas Algorithm 2 presents the pseudocode for the Julia set generation algorithm.4.1. Mandelbrot sets for the quadratic polynomial \(Q_{c}(z)=z^{2}+c\)
Examples of a quadratic Mandelbrot set, i.e., the set for \(Q_c(z) = z^2 + c\), in the \(S\) orbit with \(s\)-convexity are presented in Fig. 1. The common parameters used to generate the images were the following: \(K = 50\), \(\gamma = 0.7\), \(\delta = 0.2\). Whereas, the varying parameters were the following: (a) \(A = [-2, 0.3] \times [-1.2, 1.2]\), \(s = 0.1\), (b) \(A = [-2.4, 0.5] \times [-1.5, 1.5]\), \(s = 0.2\), (c) \(A = [-2.6, 0.5] \times [-1.8, 1.8]\), \(s = 0.3\), (d) \(A = [-3, 0.9] \times [-2, 2]\), \(s = 0.4\), (e) \(A = [-3.5, 0.9] \times [-2, 2]\), \(s = 0.5\), (f) \(A = [-4, 1] \times [-2.5, 2.5]\), \(s = 0.6\), (g) \(A = [-4.5, 1] \times [-2.7, 2.7]\), \(s = 0.7\), (h) \(A = [-6, 1] \times [-3, 3]\), \(s = 0.8\), (i) \(A = [-6, 1] \times [-3, 3]\), \(s = 0.9\), (j) \(A = [-6, 1] \times [-3, 3]\), \(s = 1.0\).Figure 1. Quadratic Mandelbrot set for \(\gamma = 0.7, \delta = 0.2\) and varying \(s\)
4.2. Mandelbrot sets for the cubic polynomial \(Q_{c}(z)=z^{3}+c\)
Examples of a cubic Mandelbrot set, i.e., the set for \(Q_c(z) = z^3 + c\), in the \(S\) orbit with \(s\)-convexity are presented in Fig. 2. The common parameters used to generate the images were the following: \(K = 50\), \(\gamma = 0.6\), \(\delta = 0.4\). Whereas, the varying parameters were the following: (a) \(A = [-0.7, 0.7] \times [-1, 1]\), \(s = 0.1\), (b) \(A = [-0.8, 0.8] \times [-1.3, 1.3]\), \(s = 0.2\), (c) \(A = [-1, 1] \times [-1.5, 1.5]\), \(s = 0.3\), (d) \(A = [-1, 1] \times [-2, 2]\), \(s = 0.4\), (e) \(A = [-1, 1] \times [-2, 2]\), \(s = 0.5\), (f) \(A = [-1, 1] \times [-2, 2]\), \(s = 0.6 ,(g), A = [-1, 1] \times [-2.2, 2.2]\), \(s = 0.7\), (h) \(A = [-1, 1] \times [-2.5, 2.5]\), \(s = 0.8\), (i) \(A = [-1, 1] \times [-2.5, 2.5]\), \(s = 0.9\), (j) \(A = [-1.5, 1.5] \times [-2.5, 2.5]\), \(s = 1.0\).Figure 2. Cubic Mandelbrot set for \(\gamma = 0.6, \delta = 0.4\) and varying \(s\).
4.3. Mandelbrot sets for the polynomial \(Q_{c}(z)=z^{k+1}+c\), where \(k=3\)
Examples of the fourth order Mandelbrot set, i.e., the set for \(Q_c(z) = z^4 + c\), in the \(S\) orbit with \(s\)-convexity are presented in Fig. 3 . The common parameters used to generate the images were the following: \(K = 50\), \(\gamma = 0.5\), \(\delta = 0.3\). Whereas, the varying parameters were the following: (a) \(A = [-0.9, 0.7] \times [-0.7, 0.7]\), \(s = 0.1\), (b) \(A = [-1.1, 0.8] \times [-0.8, 0.8]\), \(s = 0.2\), (c) \(A = [-1.3, 1] \times [-0.9, 0.9]\), \(s = 0.3\), (d) \(A = [-1.4, 1] \times [-1.1, 1.1]\), \(s = 0.4\), (e) \(A = [-1.5, 1.1] \times [-1.3, 1.3]\), \(s = 0.5\), (f) \(A = [-1.7, 1.2] \times [-1.4, 1.4]\), \(s = 0.6\), (g) \(A = [-1.8, 1.3] \times [-1.5, 1.5]\), \(s = 0.7\), (h) \(A = [-1.8, 1.4] \times [-1.6, 1.6]\), \(s = 0.8\), (i) \(A = [-1.8, 1.5] \times [-1.6, 1.6]\), \(s = 0.9\), (j) \(A = [-1.8, 1.5] \times [-1.6, 1.6]\), \(s = 1.0\).Figure 3. Mandelbrot set for \(k = 3\), \(\gamma = 0.5, \delta = 0.3\) and varying \(s\).
4.4. Julia sets for the quadratic polynomial \(Q_{c}(z)=z^{2}+c\)
Examples of a quadratic Julia set, i.e., the set for \(Q_c(z) = z^2 + c\), in the \(S\) orbit with \(s\)-convexity are presented in Fig. 4. The common parameters used to generate the images were the following: \(c = -1.45 + 0.3\mathbf{i}\), \(K = 50\), \(\gamma = 0.3\), \(\delta = 0.5\). Whereas, the varying parameters were the following: (a) \(A = [-2.1, 1.3] \times [-1, 1]\), \(s = 0.1\), (b) \(A = [-2.5, 1.5] \times [-1.4, 1.4]\), \(s = 0.2\), (c) \(A = [-2.3, 1.5] \times [-1.4, 1.4]\), \(s = 0.3\), (d) \(A = [-2.5, 1.7] \times [-1.6, 1.6]\), \(s = 0.4\), (e) \(A = [-2.5, 1.7] \times [-1.8, 1.8]\), \(s = 0.5\), (f) \(A = [-2.5, 1.7] \times [-2, 2]\), \(s = 0.6\), (g) \(A = [-2.5, 1.8] \times [-2.3, 2.3]\), \(s = 0.7\), (h) \(A = [-2.5, 1.8] \times [-2.6, 2.6]\), \(s = 0.8\), (i) \(A = [-2.5, 1.8] \times [-2.9, 2.9]\), \(s = 0.9\), (j) \(A = [-2.8, 2.3] \times [-3.4, 3.4]\), \(s = 1.0\).Figure 4. Quadratic Julia set for \(c = -1.45 + 0.3\mathbf{i}\), \(\gamma = 0.3, \delta = 0.5\) and varying \(s\).
5. Conclusions
In this paper we proposed the use instead of a convex combination the \(s\)-convex combination in the \(S\)-iteration. Using this modified \(S\)-iteration we derived escape criteria for the polynomial function of the form \(z^{k+1} + c\), where \(c \in \mathbb{C}\) and \(k = 1, 2, \ldots\). Moreover, we presented some graphical examples of Mandelbrot and Julia sets generated using the escape time algorithm with the derived escape criteria. Very interesting changes in the Mandelbrot and Julia sets can be observed when \(s\) varies from low to high values. Because the \(s\)-convex combination is a generalization of the convex combination, thus the results presented in this paper are generalization of the results obtained by Kang et al. in [14].
In our further work we will try to derive the escape criteria in the \(S\)-iteration with \(s\)-convexity for functions of other classes than the polynomial one, e.g., trigonometric. Moreover, in the fixed point literature we can find many different iteration methods that can be used in the study of Julia and Mandelbrot sets. A review of explicit iterations and their dependencies can be found in [26].
Competing Interest
The authors declare no competing interest.References
- Mandelbrot, B. B. (1983). The fractal geometry of nature (Vol. 173). New York: WH freeman.[Google Scholor]
- Dhurandhar, S. V., Bhavsar, V. C., & Gujar, U. G. (1993). Analysis of \(z\)-plane fractal images from \(z\leftarrow z^{\alpha}+ c\) for \(a< 0\). Computers & Graphics, 17(1), 89-94.
- Lakhtakia, A., Varadan, V. V., Messier, R., & Varadan, V. K. (1987). On the symmetries of the Julia sets for the process \(z\leftarrow z^p+ c\). Journal of Physics A: Mathematical and General, 20(11), 3533.[Google Scholor]
- Peherstorfer, F., & Stroh, C. (2001). Connectedness of Julia sets of rational functions. Computational Methods and Function Theory, 1(1), 61-79. [Google Scholor]
- Domínguez, P., & Fagella, N. (2008, January). Residual Julia sets of rational and transcendental functions. In Transcendental dynamics and complex analysis (pp. 138-164). [Google Scholor]
- Koss, L. (2016). Elliptic Functions with Disconnected Julia Sets. Int. J. Bifurcation Chaos, 26(6), 1650095 (2016) [10 pages].[Google Scholor]
- Crowe, W. D., Hasson, R., Rippon, P. J., & Strain-Clark, P. E. D. (1989). On the structure of the Mandelbar set. Nonlinearity, 2(4), 541.[Google Scholor]
- Dang, Y., Kauffman, L. H., & Sandin, D. J. (2002). Hypercomplex Iterations: Distance Estimation and Higher Dimensional Fractals (Vol. 1). World Scientific. ISBN: 978-981-02-3296-2 [Google Scholor]
- Griffin, C. J., & Joshi, G. C. (1992). Octonionic Julia sets. Chaos, Solitons & Fractals, 2(1), 11-24.[Google Scholor]
- Katunin, A. (2016, September). Analysis of 4D Hypercomplex Generalizations of Julia Sets. In International Conference on Computer Vision and Graphics (pp. 627-635). Springer, Cham. [Google Scholor]
- Rani, M., & Kumar, V. (2004). Superior Julia set. Research in Mathematical Education, 8(4), 261-277.[Google Scholor]
- Rani, M., & Kumar, V. (2004). Superior Mandelbrot set. Research in Mathematical Education, 8(4), 279-291.[Google Scholor]
- Chauhan, Y. S., Rana, R., & Negi, A. (2010). New Julia sets of Ishikawa iterates. International Journal of Computer Applications, 7(13), 34-42.[Google Scholor]
- Kang, S. M., Rafiq, A., Latif, A., Shahid, A. A., & Ali, F. (2016). Fractals through Modified Iteration Scheme. Filomat, 30(11), 3033-3046. [Google Scholor]
- Kang, S. M., Rafiq, A., Latif, A., Shahid, A. A., & Kwun, Y. C. (2015). Tricorns and Multicorns of-Iteration Scheme. Journal of Function Spaces, Volume 2015, Article ID 417167, 7 pages. [Google Scholor]
- Rani, M., & Chugh, R. (2014). Julia sets and Mandelbrot sets in Noor orbit. Applied Mathematics and Computation, 228, 615-631. [Google Scholor]
- Pinheiro, M. R. (2008). S-convexity (Foundations for Analysis). Differ. Geom. Dyn. Syst, 10, 257-262.[Google Scholor]
- Mishra, M. K., Ojha, D. B., & Sharma, D. (2011). Fixed point results in tricorn and multicorns of Ishikawa iteration and s-convexity. IJEST, 2(2), 157-160.[Google Scholor]
- Mishra, M. K., & Ojha, D. B. Some Common Fixed Point Results In Relative Superior Julia Sets With Ishikawa Iteration And S-Convexity. IJAEST 2(2), 174 - 179 [Google Scholor]
- Kang, S. M., Nazeer, W., Tanveer, M., & Shahid, A. A. (2015). New Fixed Point Results for Fractal Generation in Jungck Noor Orbit with-Convexity. Journal of Function Spaces, Volume 2015 (2015), Article ID 963016, 7 pages.[Google Scholor]
- Nazeer, W., Kang, S. M., Tanveer, M., & Shahid, A. A. (2015). Fixed point results in the generation of Julia and Mandelbrot sets. Journal of Inequalities and Applications, 2015(1), 298.[Google Scholor]
- Cho, S. Y., Shahid, A. A., Nazeer, W., & Kang, S. M. (2016). Fixed point results for fractal generation in Noor orbit and s-convexity. SpringerPlus, 5(1), 1843.[Google Scholor]
- Barnsley, M. F. (1988). Fractals Everywhere, Acad. Press, New York.[Google Scholor]
- Devaney, R. L. (1993). A First Course in Chaotic Dynamical Systems: Theory and Experiment., Avalon Publishing,CA 94608-1077 USA.[Google Scholor]
- Agarwal, R. P., O Regan, D., & Sahu, D. R. (2007). Iterative construction of fixed points of nearly asymptotically nonexpansive mappings. Journal of Nonlinear and convex Analysis, 8(1), 61-79.[Google Scholor]
- Gdawiec, K., & Kotarski, W. (2017). Polynomiography for the polynomial infinity norm via Kalantari’s formula and nonstandard iterations. Applied Mathematics and Computation, 307, 17-30.[Google Scholor]