Next Article in Journal
Semi-Analytical Analysis of Drug Diffusion through a Thin Membrane Using the Differential Quadrature Method
Previous Article in Journal
Alleviating Long-Tailed Image Classification via Dynamical Classwise Splitting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Numerical Solution and Comparative Study of the Symmetric Rossler Attractor with the Generalized Caputo Fractional Derivative via Two Different Methods

by
Mohamed Elbadri
1,2,
Mohamed A. Abdoon
3,4,
Mohammed Berir
4,5 and
Dalal Khalid Almutairi
6,*
1
Department of Mathematics, Faculty of Sciences and Arts, Jouf University, Tubarjal 74713, Saudi Arabia
2
Department of Mathematic, University of Gezira, Wad Madani 21111, Sudan
3
Department of Basic Sciences (Mathematics), Deanship of Preparatory Year, Shaqra University, Riyadh 15342, Saudi Arabia
4
Department of Mathematics, Faculty of Science, Bakht Al-Ruda University, Duwaym 999129, Sudan
5
Department of Mathematics, Faculty of Science and Arts, Al-Baha University, Baljurashi 1988, Saudi Arabia
6
Department of Mathematics, College of Education (Majmaah), Majmaah University, Al-Majmaah 11952, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(13), 2997; https://doi.org/10.3390/math11132997
Submission received: 9 June 2023 / Revised: 30 June 2023 / Accepted: 3 July 2023 / Published: 5 July 2023
(This article belongs to the Section Computational and Applied Mathematics)

Abstract

:
This study focuses on the solution of the rotationally symmetric Rossler attractor by using the adaptive predictor–corrector algorithm (Apc-ABM-method) and the fractional Laplace decomposition method ( ρ -Laplace DM). Furthermore, a comparison between the proposed methods and Runge–Kutta Fourth Order (RK4) is made. It is discovered that the proposed methods are effective and yield solutions that are identical to the approximate solutions produced by the other methods. Therefore, we can generalize the approach to other systems and obtain more accurate results. In addition to this, it has been shown to be useful for correctly discovering examples via the demonstration of attractor chaos. In the future, the two methods can be used to find the numerical solution to a variety of models that can be used in science and engineering applications.

1. Introduction

In recent decades, fractional differential equations have become prominent due to modeling and chaos [1,2]. Many alternative approaches to solving fractional differential equations have been developed [3,4]. Many fields, from electrical engineering to biology and physics, have found use for the modeling of chaotic and hyper-chaotic systems [5,6,7,8].
Numerous studies cover chaos application: the modeling of electrical circuits. The use of chaotic models is acceptable because of how difficult it is to forecast a wide range of real-world occurrences. Numerous novel methods for assessing chaotic systems have appeared in recent years [9,10,11,12]. Two of these methods, asymptotic stability and Lyapunov exponents, shed light on how the parameters of the model affect the dynamics of the chaotic model. Numerous mathematical and scientific fields make use of fractional calculus. For certain cutting-edge research and applications, scientists, mathematicians, biologists, and those from other fields [13,14,15,16,17,18,19] are increasingly turning to fractional calculus. Given the ambiguity surrounding fractional operators, this finding is noteworthy. Singularity freedom [20] applies to derivatives with exponential and Mittag–Leffler kernels. Since they consider the impact of long-term memory, the fractional derivatives are very useful.
In this study, we use methods to analyze numerical solutions to chaotic systems with generalized fractional derivative orders. By using a rigorous mathematical methodology, we explore cutting-edge and time-saving approaches to the problem of solving the new fractional chaotic model. As a result, there is hope for the investigation of further models using the proposed techniques.
In [21,22,23], Katugampola introduced the generalized fractional integral and the fractional derivative. The generalization of the definition of the fractional derivative with Caputo type is given in [24].
The researchers introduced the fractional Laplace transform to deal with generalized fractional derivatives [25]. It was used to solve some mathematical problems involving generalized fractional derivatives. A coupling of the fractional Laplace transform and some other analytical methods [26,27] was used to solve problems containing a generalized fractional derivative. The results indicated the effectiveness of these methods. One such effective method is the ρ -Laplace decomposition method ( ρ -Laplace DM), which has been used to solve differential equations involving a generalized fraction.
There are many strong reasons to use fractional derivatives in practice. The fractional Rössler system, a set of three nonlinear equations that displays chaotic dynamics, is an extension of the original Rössler system. Some possible physical or chemical phenomena that the fractional Rössler system could describe are a chemical reaction’s fluid flow dynamics and a chaotic oscillator’s dynamics with fractional damping [28], as proven by recent studies [29,30,31]. The literature contains many examples of chaotic systems, such as the Lorenz attractor, the electrical circuit described by Chua, the system described by Chen, the system described by Lu, and the basic chaotic system [32,33,34,35,36].
Systems that have symmetry are more susceptible to multistability because it is guaranteed that each asymmetric attractor will have a twin attractor that is symmetrical with it. Due to the fact that every asymmetric attractor contains a symmetric twin attractor, systems that have symmetry are particularly susceptible to the phenomenon of multistability. However, there is a possibility that there are advantages to multistability, such as the capacity to recreate and research events in the actual world, where they also exist. Consider the rotationally symmetric Rossler attractor [37], which is described by
D 0 α , ρ x ( t ) = y y z , D 0 α , ρ y ( t ) = x + a y , D 0 α , ρ z ( t ) = b + z ( x 2 c ) ,
where D 0 α , ρ ( . ) is the generalized Caputo-type fractional derivative [24], a =0.4, b = 0.4, c = 4.5, and the system is chaotic, with initial conditions x 0 = 1.5 , y 0 = 0.0 , and z 0 = 0.0 .
The prospect of multistability is offered by chaotic dynamical systems that have a symmetric structure thanks to the presence of an independent amplitude control parameter. The development of symmetric Rossler systems results in the production of a symmetric pair of unusual attractors that coexist together. In chaotic systems, symmetry offers a unique amplitude control parameter that may be used independently, which is beneficial for engineering applications.
This study focuses on the solution of the rotationally symmetric Rossler attractor by using two different methods. We are able to apply the methods to various types of systems in order to obtain more accurate findings. We are certain that our approaches will be used in the not-too-distant future to design and simulate a wide range of fractional models. These models have the potential to be utilized in the resolution of increasingly difficult physics, biology, and engineering issues.
The significance of these techniques rests in the fact that they are used to the find a numerical solution in a variety of models, including disease models and chaotic models, and that it can be expanded to incorporate other models in pathology, dynamical models, coding, and hyper-chaos. Both capabilities contribute to the method’s overall usefulness. In addition to this, the method is useful for the discovery of chaos and can be applied to other, more complex models.

2. Basic Definitions

Definition 1.
The generalized fractional integral of the function   f ,   I a + α , ρ f ( t ) , α > 0 , and  ρ > 0 , is given by
I a + α , ρ f t = ρ 1 α Γ α a t   s ρ 1 t ρ s ρ α 1 f s d s ,   α > 0 ,   t > a ,
 for  m 1 < α m , where  m N  (see [21]).
Definition 2.
The generalized Riemann–Liouville fractional derivative (RLFD) of the function  f ,   R D a + α , ρ f ( t ) , of order  α > 0 , is given by [22].
  R D a + α , ρ f ( t ) = ρ α m + 1 Γ m α t 1 ρ d d t m a t   s ρ 1 t ρ s ρ m α 1 f ( s ) d s ,   t > a 0 .
Definition 3.
The generalized Caputo fractional derivative (CFD) operator,   D a + α , ρ ,  α > 0 , is given by
D a + α , ρ f t = ρ α m + 1 Γ m α a t   s ρ 1 t ρ s ρ m α 1 s 1 ρ d d s m f s d s ,   t > a .
where  ρ > 0 , a 0 , and  m 1 < α < m  (see [24]).
Definition 4 ([25]).
The  ρ -Laplace transform of a function  f : [ 0 , ] R  is defined by 
L ρ f ( t ) = 0 e δ t ρ ρ f ( t ) d t t 1 ρ .
The  ρ -Laplace transform of the generalized CFD is defined by
L ρ D 0 α , ρ f ( t ) = δ α L ρ f ( t ) δ α 1   f 0 .

3. Methodology of the Apc-ABM Algorithm

The purpose of this part is to present the algorithm used in the APC-ABM method:
D a + α , ρ y ( t ) = f ( t , y ( t ) ) ,   t [ 0 , T ] , y ( k ) ( a ) = y 0 k ,   k = 0,1 , , α ,
where D a + α , ρ is a generalized Caputo fractional operator.
Then, for n 1 < α n , a > 0 , ρ > 0 , y C n ( [ a , T ] ) , model (8) is equivalent and we obtain
y ( t ) = u ( t ) + ρ 1 α Γ ( α ) a t   s ρ 1 t ρ s ρ α 1 f ( s , y ( s ) ) d s ,
where
u ( t ) = n = 0 m 1     1 ρ n n ! t ρ a ρ n x 1 ρ d d x n y ( x ) x = a .
The first step of our algorithm, under the assumption that the function f is such that a unique solution exists on some interval [ a , T ] , consists of dividing the interval [ a , T ] into N unequal subintervals, t k , t k + 1 , k = 0 , 1 , , N 1 , using the mesh points.
t 0 = a , t k + 1 = t k ρ + h 1 / ρ , k = 0 , 1 , , N 1 ,
where h = T ρ a ρ N . Now, to numerically solve the IVPs, we build approximations y k , k = 0 , 1 , , N . If we have previously assessed the approximations y t j and y j y t j ( j = 1 , 2 , , k ) , we wish to use the integral equation to generate the approximation y k + 1 y t k + 1 .
y t k + 1 = u t k + 1 + ρ 1 α Γ α a t k + 1   s ρ 1 t k + 1 ρ s ρ α 1 f s , y s d s ,
Making the substitution
z = s ρ
we obtain
y t k + 1 = u t k + 1 + ρ α Γ ( α ) a ρ t k + 1 ρ   t k + 1 ρ z α 1 f z 1 / ρ , y z 1 / ρ d z
That is,
y t k + 1 = u t k + 1 + ρ α Γ ( α ) j = 0 k   t j ρ t j + 1 ρ   t k + 1 ρ z α 1 f z 1 / ρ , y z 1 / ρ d z
Subsequently, the trapezoidal quadrature method is employed in consideration of the weight function t k + 1 ρ z α 1 . In order to estimate the integrals on the right-hand side of the equation, a suitable method must be employed (Equation (13)). The corrector formula is derived in the following manner, y t k + 1 ,   k = 0 , 1 , 2 , , N 1 :
y t k + 1 u t k + 1 + ρ α h α Γ ( α + 2 ) j = 0 k   a j , k + 1 f t j , y t j + ρ α h α Γ ( α + 2 ) f t k + 1 , y t k + 1
where
a j , k + 1 = k α + 1 ( k α ) ( k + 1 ) α if   j = 0 , ( k j + 2 ) α + 1 + ( k j ) α + 1 2 ( k j + 1 ) α + 1 if   1 j < k
The ultimate step of our methodology involves substituting the quantity y t k + 1 . The predictor value is located on the right-hand side of Formula (15), y P t k + 1 . The integral equation (Equation (13)) is solved through the utilization of the one-step Adams–Bashforth method. In this scenario, the act of replacing the function with another equivalent one is considered. f z 1 / ρ , y z 1 / ρ at each integral in Equation (16) with the amount f t j , y t j yields
y P t k + 1 u t k + 1 + ρ α Γ ( α ) j = 0 k     t j ρ t j + 1 ρ     t k + 1 ρ z α 1 f t j , y t j d z
y P t k + 1 = u t k + 1 + ρ α h α Γ ( α + 1 ) j = 0 k     ( k + 1 j ) α ( k j ) α f t j , y t j
Thus, the aforementioned formula comprehensively characterizes our adaptive P-C approach for evaluating the approximation, y k + 1 y t k + 1 :
y k + 1 u t k + 1 + ρ α h α Γ ( α + 2 ) j = 0 k   a j , k + 1 f t j , y j + ρ α h α Γ ( α + 2 ) f t k + 1 , y k + 1 p ,
where y j y t j , j = 0 , 1 , , k , and the predicted value y k + 1 P y P t k + 1 can be determined as described in Equation (18) with the weights a j , k + 1 being defined according to (35). The proposed adaptive Apc-ABM-method uses a non-uniform grid t j + 1 = t j ρ + h ρ : j = 0 , 1 , , N 1 with t 0 =   a and h = T ρ a ρ N , where N represents a positive integer. It is posited that the utilization of the Apc-ABM methodology for initial value problems (IVPs) featuring the generalized CFD is rendered infeasible in instances where a uniform grid is employed, as is the case with the aforementioned scenario [38].

4. Applications of the Apc-ABM Algorithm

This section is dedicated to exploring the solutions of Equation (1). With the Apc-ABM method, great results can be achieved when a = b = 0.4 and c = 4.5, with initial conditions x 0 = 1.5 , y 0 = 0 , and z 0 = 0 . By using Equation (18), the approximations x k + 1 , y k + 1 , and z k + 1 , and for N N and T > 0 ,
x k + 1 x 0 + a ρ α h α Γ ( α + 2 ) j = 0 k     a j , k + 1 y j x j + a ρ α h α Γ ( α + 2 ) y k + 1 P x k + 1 P , y k + 1 y 0 + ρ α h α Γ ( α + 2 ) j = 0 k     a j , k + 1 ( c a ) x j x j z j + c y j + ρ α h α Γ ( α + 2 ) ( c a ) x k + 1 P x k + 1 P z k + 1 P + c y k + 1 P , z k + 1 z 0 + ρ α h α Γ ( α + 2 ) j = 0 k     a j , k + 1 x j y j b z j + ρ α h α Γ ( α + 2 ) x k + 1 P y k + 1 P b z k + 1 P ,
where h = T ρ N and
x k + 1 P x 0 + a ρ α h α Γ ( α + 1 ) j = 0 k     ( k + 2 1 j ) α + ( k + j ) α y j x j , y k + 1 P y 0 + ρ α h α Γ ( α + 1 ) j = 0 k     ( k + 2 1 j ) α + ( k + j ) α ( c a ) x j x j z j + c y j , z k + 1 P z 0 + ρ α h α Γ ( α + 1 ) j = 0 k     ( k + 2 1 j ) α + ( k + j ) α x j y j b z j .
Table 1 presents the numerical solution using the Apc-ABM method to Equation (21) when α = 1 , ρ = 1 , ( a , b , c ) = ( 0.5 , 0.5 , 3.5 ) , and x 0 , y 0 , z 0 = ( 1 / 2 , 1 / 2 , 0 ) , and comparing the results with the RK4 method. Table 2 presents the numerical solution for the value of α = 0.95 .

5. ρ-Laplace DM

In this section, we discuss Equation (7) with a generalized CFD to illustrate the algorithm of the ρ -Laplace DM.
The ρ -Laplace DM divides Equation (1) into a linear term A y t , a nonlinear term B y t , and a source function C ( t ) as follows:
D 0 α , ρ y t = A y t + B y t + C ( t ) , ρ > 0 , 0 < α 1 ,
Taking L ρ to Equation (22), we obtain
L ρ y t = y 0 δ + 1 δ α L ρ C ( t ) + 1 δ α   L ρ A y t + B y t
Operating L ρ 1 , we obtain
y t = G ( t ) + L ρ 1 1 δ α   L ρ A y t + B y t
where G t = L ρ 1 y 0 δ + 1 δ α L ρ C t .
The ρ -Laplace DM represents a series solution of y t by
y t = n = 0 y n ( t )
Furthermore, the nonlinear function B y t can be expressed in a series of polynomials:
n = 0 K n
Substituting Equations (25) and (26) into Equation (24) yields
n = 0 y n t = G ( t ) + L ρ 1 1 δ α   L ρ A n = 0 y n t + n = 0 K n .
As a result, the recurrence relation
y 0 τ = G ( t ) ,
y n + 1 τ = L ρ 1 1 δ α   L ρ A y n t + K n ,   n 0
Finally, a series is used to approximate the solution.
ψ M ( t ) = n = 0 M 1 y n ( t )

6. Application of ρ-Laplace DM

According to what is presented in Section 5, after applying the ρ -Laplace DM to each equation of System (1), we obtain
n = 0 x n t = x 0 + L ρ 1 1 δ α   L ρ n = 0 y n t n = 0 U n n = 0 y n t = y 0 + L ρ 1 1 δ α   L ρ n = 0 x n t a n = 0 y n t , n = 0 z n t = b t α ρ ρ α Γ [ 1 + α ] + L ρ 1 1 δ α   L ρ n = 0 V n c n = 0 z n t .
where the nonlinear terms y z and z x 2 are given by
y z = n = 0 U n   and   y z = n = 0 U n
The ρ -Laplace DM provides the recursive relation:
x 0 t = 1.5
y 0 t = 0
z 0 t = b t α ρ ρ α Γ [ 1 + α ]
x n + 1 t = L ρ 1 1 δ α   L ρ y n t U n t , n 0
y n + 1 t = L ρ 1 1 δ α   L ρ x n t a y n t , n 0
z n + 1 t = L ρ 1 1 δ α   L ρ V n t c z n t , n 0
The first few components are
x 0 t = 1.5
y 0 t = 0
z 0 t = b t α ρ ρ α Γ [ 1 + α ]
x 1 t = 0
y 1 t = 1.5 t ρ α ρ α Γ [ 1 + α ]
z 1 t = b ( 2.25   c ) t 2 ρ α ρ 2 α Γ [ 1 + 2 α ]
x 2 t = 1.5 b t 3 α ρ ρ 3 α Γ [ 1 + 2 α ] Γ 2 [ 1 + α ] Γ [ 1 + 3 α ] 1.5 t 2 α ρ ρ 2 α Γ [ 1 + 2 α ]
y 2 t = 1.5 a t 2 α ρ ρ 2 α Γ [ 1 + 2 α ]
z 2 t = b ( 5.0625   4.5 c + c 2 ) t 3 α ρ ρ 3 α Γ [ 1 + 3 α ]
x 3 t = 1.5 a t 3 α ρ ρ 3 α Γ [ 1 + 3 α ] 1.5 b ( 2.25   c ) Γ [ 1 + 3 α ] t 4 α ρ ρ 4 α Γ [ 1 + α ] Γ [ 1 + 2 α ] Γ [ 1 + 4 α ] 1.5 a b Γ [ 1 + 3 α ] t 4 α ρ ρ 4 α Γ [ 1 + α ] Γ [ 1 + 2 α ] Γ [ 1 + 4 α ]
y 3 t = 1.5 t 3 α ρ ρ 3 α Γ [ 1 + 3 α ] + 1.5 a 2 t 3 α ρ ρ 3 α Γ [ 1 + 3 α ] 1.5 b t 4 α ρ ρ 4 α Γ [ 1 + 2 α ] Γ 2 [ 1 + α ] Γ [ 1 + 4 α ]
z 3 t = 4.5 b t 4 α ρ ρ 4 α Γ [ 1 + 3 α ] Γ [ 1 + α ] Γ [ 1 + 2 α ] Γ [ 1 + 4 α ] + 11.390625 b t 4 α ρ ρ 4 α Γ [ 1 + 4 α ] 15.1875 b c t 4 α ρ ρ 4 α Γ [ 1 + 4 α ] + 6.75 b c 2 t 4 α ρ ρ 4 α Γ [ 1 + 4 α ] b c 3 t 4 α ρ ρ 4 α Γ [ 1 + 4 α ] 4.5 b 2 t 5 α ρ ρ 5 α Γ [ 1 + 2 α ] [ 1 + 4 α ] Γ 3 [ 1 + α ] Γ [ 1 + 3 α ] Γ [ 1 + 5 α ]
Therefore, the approximate solution is given as
x t = n = 0 x n t , y t = n = 0 y n t , z t = n = 0 z n t .
Table 3 displays the ρ -Laplace DM solution of Equation (1). When α = 0.90 and ρ = 0.95 , we observed that the ρ -Laplace DM solution given in Table 3 for the fractional order and fractional parameter has the same behavior as the ρ -Laplace DM solution given in Table 4 for integers α = 1 and ρ = 1 .
Table 4 also shows the comparison of outcomes of the solution of system (1) using the ρ -Laplace DM, Apc-ABM method, and RK4 method when α = 1 , ρ = 1 a = b = 0.4, and c = 4.5, with x 0 = 1.5 , y 0 = 0 , and z 0 = 0 . We note that the results obtained when using the ρ -Laplace DM and Apc-ABM methods are very close to those obtained using the RK4 method.
In Figure 1 and Figure 2, solutions were drawn using Equation (1) by using the Apc-ABM-method, with ( a , b , c ) = ( 0.2 , 0.2 , 6.5 ) , x 0 , y 0 , z 0 = ( 1 , 0 , 0 ) , T = 400 , and N = 800 when α = 0.95 and ρ = 0.8 , 1.2 . In Figure 3, solutions were drawn using Equation (1) by using α = 1 and ρ = 1 .

7. Conclusions

This study presents the numerical solution of a fractional system using two different methods and compares the solutions. The APc-ABM method and the ρ -Laplace DM both benefited from the provision of a numerical strategy that was accomplished with the help of software packages. As the step size h fell, the APc-ABM technique produced numerical results that were impressively close to the RK4 solutions. Furthermore, we showed from the comparison that the numerical solutions obtained by the Laplace DM and APc-ABM are identical to the approximate solutions produced by the RK4 method. The obtained numerical results demonstrate that our method carries out its procedures in the context of fractions in a way that satisfies expectations with regard to the degree to which it maintains its numerical stability. We suggest applying this strategy to more complex physics and engineering problems.

Author Contributions

Methodology, M.E., M.A.A. and D.K.A.; Software, M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No data were used to support this study.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research at Majmaah University for supporting this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  2. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  3. Carpinteri, A.; Mainardi, F. (Eds.) Fractals and Fractional Calculus in Continuum Mechanics; Springer: Wien, Austria; New York, NY, USA, 1997; pp. 277–290. [Google Scholar]
  4. Samko, G.; Kilbas, A.; Marichev, O. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: Amsterdam, The Netherlands, 1993. [Google Scholar]
  5. Dudkowski, D.; Jafari, S.; Kapitaniak, T.; Kuznetsov, N.V.; Leonov, G.A.; Prasad, A. Hidden Attractors in Dynamical Systems. Phys. Rep. 2016, 637, 1–50. [Google Scholar] [CrossRef]
  6. Xu, G.; Shekofteh, Y.; Akgül, A.; Li, C.; Panahi, S. A New Chaotic System with a Self-Excited Attractor: Entropy Measurement, Signal Encryption, and Parameter Estimation. Entropy 2018, 20, 86. [Google Scholar] [CrossRef] [Green Version]
  7. Lai, Q.; Akgul, A.; Li, C.; Xu, G.; Çavuşoğlu, Ü. A New Chaotic System with Multiple Attractors: Dynamic Analysis, Circuit Realization and S-Box Design. Entropy 2017, 20, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Petráš, I. A Note on the Fractional-Order Chua’s System. Chaos Solitons Fractals 2006, 38, 140–147. [Google Scholar] [CrossRef]
  9. Abdoon, M.A.; Hasan, F.L.; Taha, N.E. Computational Technique to Study Analytical Solutions to the Fractional Modified KDV-Zakharov-Kuznetsov Equation. Abstr. Appl. Anal. 2022, 2022, 2162356. [Google Scholar] [CrossRef]
  10. Abd El-Maksoud, A.J.; Abd El-Kader, A.A.; Hassan, B.G.; Rihan, N.G.; Tolba, M.F.; Said, L.A.; Radwan, A.G.; Abu-Elyazeed, M.F. FPGA Implementation of Sound Encryption System Based on Fractional-Order Chaotic Systems. Microelectronics 2019, 90, 323–335. [Google Scholar] [CrossRef]
  11. Ahmed, S.A.; Elzaki, T.M.; Elbadri, M.; Mohamed, M.Z. Solution of Partial Differential Equations by New Double Integral Transform (Laplace—Sumudu Transform). Ain Shams Eng. J. 2021, 12, 4045–4049. [Google Scholar] [CrossRef]
  12. Sene, N. Global Asymptotic Stability of the Fractional Differential Equations. J. Nonlinear Sci. Its Appl. 2019, 13, 171–175. [Google Scholar] [CrossRef] [Green Version]
  13. Singh, J.; Kumar, D.; Hammouch, Z.; Atangana, A. A Fractional Epidemiological Model for Computer Viruses Pertaining to a New Fractional Derivative. Appl. Math. Comput. 2017, 316, 504–515. [Google Scholar] [CrossRef]
  14. Kumar, S.; Kumar, R.; Cattani, C.; Samet, B. Chaotic Behaviour of Fractional Predator-Prey Dynamical System. Chaos Solitons Fractals 2020, 135, 109811. [Google Scholar] [CrossRef]
  15. Hammouch, Z.; Mekkaoui, T.; Belgacem, F.B.M. Numerical Simulations for a Variable Order Fractional Schnakenberg Model. AIP Conf. Proc. 2014, 1637, 1450–1455. [Google Scholar]
  16. Abdoon, M.A.; Saadeh, R.; Berir, M.; Guma, F.E. Analysis, Modeling and Simulation of a Fractional-Order Influenza Model. Alex. Eng. J. 2023, 74, 231–240. [Google Scholar] [CrossRef]
  17. Hammouch, Z.; Mekkaoui, T. Traveling-wave solutions of the generalized Zakharov equation with time-space fractional derivatives. Math. Eng. Sci. Aerosp. 2014, 5, 489–498. [Google Scholar]
  18. Belgacem, F.B.M.; Silambarasan, R.; Zakia, H.; Mekkaoui, T. New and Extended Applications of the Natural and Sumudu Transforms: Fractional Diffusion and Stokes Fluid Flow Realms. In Trends in Mathematics; Springer: Singapore, 2017; pp. 107–120. [Google Scholar]
  19. Toufik, M.; Atangana, A. New Numerical Approximation of Fractional Derivative with Non-Local and Non-Singular Kernel: Application to Chaotic Models. Eur. Phys. J. Plus 2017, 132, 444. [Google Scholar] [CrossRef]
  20. Behzad Ghanbari, D. Abdon Atangana, Some New Edge Detecting Techniques Based on Fractional Derivatives with Nonlocal and Non-Singular Kernels. Adv. Differ. Equ. 2020, 2020, 435. [Google Scholar] [CrossRef]
  21. Katugampola, U.N. New Approach to a Generalized Fractional Integral. Appl. Math. Comput. 2011, 218, 860–865. [Google Scholar] [CrossRef] [Green Version]
  22. Katugampola, U.N. A New Approach to Generalized Fractional Derivatives. Bull. Math. Anal. Appl. 2014, 6, 1–15. [Google Scholar]
  23. Katugampola, U.N. Existence and Uniqueness Results for a Class of Generalized Fractional Differential Equations. arXiv 2014, arXiv:1411.5229v2, preprint. [Google Scholar]
  24. Almeida, R.; Malinowska, A.B.; Odzijewicz, T. Fractional Differential Equations with Dependence on the Caputo–Katugampola Derivative. J. Comput. Nonlinear Dyn. 2016, 11, 061017. [Google Scholar] [CrossRef] [Green Version]
  25. Jarad, F.; Abdeljawad, T. A modified Laplace transform for certain generalized fractional operators. Results Nonlinear Anal. 2018, 1, 88–98. [Google Scholar]
  26. Sene, N.; Fall, A.N. Homotopy Perturbation ρ-Laplace Transform Method and Its Application to the Fractional Diffusion Equation and the Fractional Diffusion-Reaction Equation. Fractal Fract. 2019, 3, 14. [Google Scholar] [CrossRef] [Green Version]
  27. Elbadri, M. Initial Value Problems with Generalized Fractional Derivatives and Their Solutions via Generalized Laplace Decomposition Method. Adv. Math. Phys. 2022, 2022, 3586802. [Google Scholar] [CrossRef]
  28. Rysak, A.; Sedlmayr, M.; Gregorczyk, M. Revealing Fractionality in the Rössler System by Recurrence Quantification Analysis. Eur. Phys. J. Spec. Top. 2023, 232, 83–98. [Google Scholar] [CrossRef]
  29. Qazza, A.; Abdoon, M.; Saadeh, R.; Berir, M.A. New Scheme for Solving a Fractional Differential Equation and a Chaotic System. Eur. J. Pure Appl. Math. 2023, 16, 1128–1139. [Google Scholar] [CrossRef]
  30. Saadeh, R.; AAbdoon, M.; Qazza, A.; Berir, M. A Numerical Solution of Generalized Caputo Fractional Initial Value Problems. Fractal Fract. 2023, 7, 332. [Google Scholar] [CrossRef]
  31. Elbadri, M.; Abdoon, M.A.; Berir, M.; Almutairi, D.K. A Symmetry Chaotic Model with Fractional Derivative Order via Two Different Methods. Symmetry 2023, 15, 1151. [Google Scholar] [CrossRef]
  32. Elbadri, M.; Ahmed, S.A.; Abdalla, Y.T.; Hdidi, W. A New Solution of Time-Fractional Coupled KdV Equation by Using Natural Decomposition Method. Abstr. Appl. Anal. 2020, 2020, 3950816. [Google Scholar] [CrossRef]
  33. Sene, N. Mathematical Views of the Fractional Chua’s Electrical Circuit Described by the Caputo-Liouville Derivative. Rev. Mex. Fis. 2021, 67, 91–99. [Google Scholar] [CrossRef]
  34. Morin, O.; Gillis, A.; Chen, J.; Aubin, M.; Bucci, M.K.; Roach, M., 3rd; Pouliot, J. Megavoltage Cone-Beam CT: System Description and Clinical Applications. Med. Dosim. Off. J. Am. Assoc. Med. Dosim. 2005, 31, 51–61. [Google Scholar] [CrossRef]
  35. Zhang, D.; Lu, G. Review of Shape Representation and Description Techniques. Pattern Recognit. 2004, 37, 1–19. [Google Scholar] [CrossRef]
  36. Zhang, X.; Wang, C.; Yao, W.; Lin, H. Chaotic System with Bondorbital Attractors. Nonlinear Dyn. 2019, 97, 2159–2174. [Google Scholar] [CrossRef]
  37. Li, C.; Hu, W.; Sprott, J.C.; Wang, X. Multistability in Symmetric Chaotic Systems. Eur. Phys. J. Spec. Top. 2015, 224, 1493–1506. [Google Scholar] [CrossRef]
  38. Odibat, Z.; Baleanu, D. Numerical Simulation of Initial Value Problems with Generalized Caputo-Type Fractional Derivatives. Appl. Numer. Math. Trans. IMACS 2020, 156, 94–105. [Google Scholar] [CrossRef]
Figure 1. Chaotic attractor of Equation (1), when α = 0.95 and ρ = 0.8 .
Figure 1. Chaotic attractor of Equation (1), when α = 0.95 and ρ = 0.8 .
Mathematics 11 02997 g001
Figure 2. Chaotic attractor of Equation (1), when α = 0.95 and ρ = 1.2 .
Figure 2. Chaotic attractor of Equation (1), when α = 0.95 and ρ = 1.2 .
Mathematics 11 02997 g002
Figure 3. Chaotic attractor of Equation (1), when α = 1 and ρ = 1 .
Figure 3. Chaotic attractor of Equation (1), when α = 1 and ρ = 1 .
Mathematics 11 02997 g003
Table 1. Numerical solutions for a fractional equation (Equation (1)) when α = 1 , ρ = 1 , N = 200 , and t = 2 .
Table 1. Numerical solutions for a fractional equation (Equation (1)) when α = 1 , ρ = 1 , N = 200 , and t = 2 .
h x y z
1/160−1.4392490042450221.9268797538303340.115497632553681
1/320−1.4370278042762521.9197498622741070.115546731072315
1/640−1.4359041008084111.9161919439455590.115570069443809
1/1280−1.4353389885191511.9144147586783360.115581435786859
1/2560−1.4350556191617821.9135266116376550.115587043236967
1/5120−1.4349137314214011.9130826497797140.115589828030596
1/10240−1.4348427368160091.9128606967993130.115591215694437
1/20480−1.4347805887996491.9126665032358240.115592427311921
R K4−1.4347708381598451.9126375980493350.115592614892760
Table 2. Numerical solutions for a fractional equation (Equation (1)) when α = 0.95 , ρ = 1 , and t = 0.5 .
Table 2. Numerical solutions for a fractional equation (Equation (1)) when α = 0.95 , ρ = 1 , and t = 0.5 .
h x y z
1/1601.289215263867680.7940128941943670.113077179858054
1/3201.2878659848438710.7938569548127100.112751179833214
1/6401.2871921290569130.7937766343175180.112588754203015
1/12801.2868553997335130.7937358865105490.112507684690222
1/25601.2866870850436680.7937153657366760.112467185722739
1/51201.2866029402328410.7937050686346850.112446945181642
1/102401.2865608709660920.7936999109052490.112436827146182
1/204801.2865240620756990.7936953928726330.112427975087295
R K41.2865188033328910.7936947470802760.112426709876828
Table 3. ρ -Laplace DM solutions to fractional model equation (Equation (1)) when α = 0.90 and ρ = 0.95 .
Table 3. ρ -Laplace DM solutions to fractional model equation (Equation (1)) when α = 0.90 and ρ = 0.95 .
t x y z
01.50.00.0
0.1 1.4796996703743301 0.2347426924022824 0.05053705282095178
0.2 1.430966701161987 0.4315761993664601 0.07795664830887342
0.3 1.3578018024536942 0.6161357560230636 0.09178543190387903
0.4 1.2621444225163563 0.7907363746739074 0.09111861049966521
0.5 1.145647714788557 0.9552091317148814 0.07230353761707926
Table 4. A comparison of the numerical solutions to fractional model equation (Equation (1)).
Table 4. A comparison of the numerical solutions to fractional model equation (Equation (1)).
t x - ρ -Laplace D x - Apc-ABM x - RK4
01.51.51.5
0.11.4922138751.4922197896425301.492218810160190
0.21.4678221.4678998376211271.467897787108336
0.31.42552387499999991.4259017429142821.425898584622672
0.41.36435199999999981.3655273697753991.365523123300982
0.51.2836718751.2865240620756991.286518803332891
t y -  ρ -Laplace D y - Apc-ABM y - RK4
00.00.00.0
0.10.152785000000000030.1527804185455650.152780702318047
0.20.310240000000000070.3101756374251400.310176020055442
0.30.47092500000000010.4706250992204810.470625370691797
0.40.63328000000000020.6324106598294440.632410592018532
0.50.79562499999999990.7936953928726330.793694747080276
t z -  ρ -Laplace D z - Apc-ABM z - RK4
00.00.00.0
0.10.0357955356250.0357995379763470.035799091697084
0.20.064020890000000010.0641327336306660.064131982855701
0.30.0851356256250.0859616813719140.085960705651558
0.40.098488480.1018931865042050.101892039586943
0.50.1022597656250.1124279750872950.112426709876828
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Elbadri, M.; Abdoon, M.A.; Berir, M.; Almutairi, D.K. A Numerical Solution and Comparative Study of the Symmetric Rossler Attractor with the Generalized Caputo Fractional Derivative via Two Different Methods. Mathematics 2023, 11, 2997. https://doi.org/10.3390/math11132997

AMA Style

Elbadri M, Abdoon MA, Berir M, Almutairi DK. A Numerical Solution and Comparative Study of the Symmetric Rossler Attractor with the Generalized Caputo Fractional Derivative via Two Different Methods. Mathematics. 2023; 11(13):2997. https://doi.org/10.3390/math11132997

Chicago/Turabian Style

Elbadri, Mohamed, Mohamed A. Abdoon, Mohammed Berir, and Dalal Khalid Almutairi. 2023. "A Numerical Solution and Comparative Study of the Symmetric Rossler Attractor with the Generalized Caputo Fractional Derivative via Two Different Methods" Mathematics 11, no. 13: 2997. https://doi.org/10.3390/math11132997

APA Style

Elbadri, M., Abdoon, M. A., Berir, M., & Almutairi, D. K. (2023). A Numerical Solution and Comparative Study of the Symmetric Rossler Attractor with the Generalized Caputo Fractional Derivative via Two Different Methods. Mathematics, 11(13), 2997. https://doi.org/10.3390/math11132997

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