Study of a Fractional-Order Chaotic System Represented by the Caputo Operator

)is paper is presented on the theory and applications of the fractional-order chaotic system described by the Caputo fractional derivative. Considering the new fractional model, it is important to establish the presence or absence of chaotic behaviors. )e Lyapunov exponents in the fractional context will be our fundamental tool to arrive at our conclusions. )e variations of the model’s parameters will generate chaotic behavior, in general, which will be established using the Lyapunov exponents and bifurcation diagrams. For the system’s phase portrait, we will present and apply an interesting fractional numerical discretization. For confirmation of the results provided in this paper, the circuit schematic is drawn and simulated. As it will be observed, the results obtained after the simulation of the numerical scheme and with the Multisim are in good agreement.


Introduction
In the last decade, chaos theory has attracted many researchers. Chaos theory is a field of mathematics that can be applied in many domains: from modeling chaotic financial systems [1], in representing circuit schematics [2,3], and others. e chaotic systems and their circuit schematics have received many investigations in the literature, see, for example, [2][3][4]. With the new development of fractional calculus, some researchers were interested in modeling chaotic systems using the new and old fractional operators. e impact of the fractional derivatives is an interesting question which merits investigations. Many problems are opened in fractional calculus. e first question is how to construct the numerical scheme to obtain the phase portraits of chaotic systems. Many numerical techniques are used in fractional differential equations that can be applied naturally to solve chaotic systems. In [5], Danca finds interesting remarks that the characterization of the hyperchaotic system using two positive Lyapunov exponents is not an adequate definition in a fractional context. erefore, modeling chaos using fractional operators is an open problem that admits some interesting interrogations which are not solved yet.
Fractional calculus finds many recent advancements in fractional operators and their applicability in real-world problems. For recent advances of fractional calculus and its applications, see the following investigations [1,[6][7][8][9][10]. e chaos literature in integer version and fractional context is very long, but we recall the investigations we find essential for our study. In [11], Pacheco et al. study a new analysis of a new chaotic system with different families of hidden and self-excited attractors by using the fractional operator. In [12], Rajagopal et al. analyze using fractionalorder derivative a memristor no equilibrium chaotic system; they proposed synchronizing their introduced model. In other words, they try to find adequate control to synchronize slave and master chaotic systems. In [2], Avalos-Ruiz et al. propose a control implementation on an FPGA for a class of fractional chaotic systems in the context of fractional variable order derivative with Mittag-Leffler kernel. In [13], Perez et al. propose the fractional-order chaotic system in the context of conformable Liouville-Caputo operator. At the same time, they present a novel numerical scheme for the used fractional operator to depict the phase portraits of their proposed chaotic system. In [14], Gomez et al. investigate the chaos in a calcium oscillation model using the Mittag-Leffler fractional derivative operator. In [15], Atangana and Gomez offer a new hyperchaotic system in the context of the Mittag-Leffler fractional operator. ey provided a numerical scheme to depict the considered model's phase portrait in the fractional operators' context. ey mainly apply their numerical procedure to the fractional Chua model. In [16], Li et al. introduce a new four-wing model with integer-order derivative, which will be subject to investigations in the fractional context. In [3], Akgul gives a chaotic system for modeling the memcapacitor model using a fractional-order derivative. In [17], Rajagopal et al. continue their exciting investigation of the fractional-order chaotic system in the fractional calculus field. ey propose a work on the multiscroll chaotic system and give an application to the synchronization of their proposed model. In [18], Akgul et al. address a user interface for the generation of random numbers based on fractional and integer-order chaotic systems. For more investigations in chaos, see in the following investigations [4,[12][13][14]19]. e novelties of the present paper are described in the following lines. is paper focuses on the phase space, the bifurcation diagrams, and the chaos's characterization using the Lyapunov exponents for a fractional version of the Li et al. chaotic system. We analyze the impact of the fractionalorder derivative on the dynamics of a fractional version of the system in [16] with Lyapunov exponents and bifurcation diagram aids. Its possible implementation in terms of electrical modeling has been proposed as well. One of the novelties of the present work is the Lyapunov exponents in the fractional context, which is an open problem. Note that the Lyapunov exponents' properties in the integer context for hyperchaotic and chaotic behaviors are not adequate definitions in the fractional version. is remark was first discussed in [5]. Furthermore, the slight variation of the model's parameters in the dynamics has been investigated using bifurcation diagrams. We provide in this paper that the initial conditions impact the chaotic behaviors in the fractional context. e circuit schematic of the fractional chaotic system and the results after simulation in Multisim confirm the theoretical findings obtained via the numerical scheme. e paper is structured in the following form. In Section 2, we define the fractional operators' necessaries for our work. In Section 3, we present the fractional-order chaotic system described by the Caputo derivative. Section 4 offers the numerical scheme used to obtain the considered fractional system's phase portraits. In Section 5, we illustrate the numerical method by depicting the supposed system's phase portraits with different fractional-order of the Caputo derivative. We also characterize the nature of the chaos using the Lyapunov exponents. In Section 6, we present the bifurcation diagrams and phase space according to the variation of the considered fractional-order chaotic system's parameters. In Section 7, we give the possible implementation of the fractional-order chaotic system in modeling electrical circuits; this section will confirm the paper's theoretical findings. Section 8 provides the stability analysis of the fractional-order chaotic system's equilibrium points with the Matignon criterion used in fractional calculus. In Section 9, we analyze the sensitivity of the initial conditions in the fractional-order chaotic system. We finish with the final remarks and future perspectives of research.

On Fractional Operators in Fractional Calculus
In this paper, we study a fractional differential system; before, it was indispensable to recall the fractional operators used in this work. Many fractional derivatives exist, those with singularities as the Caputo and Riemann-Liouville derivatives, those with nonsingularities as the exponential derivative, and the Mittag-Leffler fractional derivative. Due to space limitations, we provide the Caputo derivative and the Riemann-Liouville derivative and their generalizations.
Definition 1 (see [20,21]). e integral defined by Riemann and Liouville called Riemann-Liouville integral of a given function z: with the function Γ(. . .) symbolizing the Gamma Euler function and under the order α > 0.
ere exists a generalization of the Riemann-Liouville operator recently proposed in the literature; we have the following definition.
Definition 2 (see [22]). e generalized Riemann-Liouville fractional integral of a given function z: where the orders α and ρ satisfy the following relationship α and ρ > 0 and the function gamma is Γ(. . .), for all t > 0.
We can observe that when the order ρ � 1, we get the classical Riemann-Liouville integral operator. Furthermore, when α � ρ � 1, we have the integer version of the integral.
Caputo proposed another derivative due to the inconvenience of the Riemann-Liouville operator. We recall the Caputo derivative in the following definitions and its generalization.
ere exist many fractional operators, and all of them have their advantages. e use of the Caputo derivative in the present work has many motivations. Firstly, all fractional operators take into account the memory effect. Secondly, our main motivation is the initial conditions because the Riemann-Liouville derivative cannot consider them. In the Riemann-Liouville derivative case, the starting condition should be an integral form known in physics to be unrealistic. It cannot be satisfied in many real-world problems. Another basic motivation is in the context of Riemann-Liouville; the derivative of a constant is not null, contrary to the context of the Caputo derivative where the derivative of a constant is null.

Fractional-Order Chaotic System
is section provides a fractional version of the chaotic model based on the system presented by Li et al. [16]. Before the investigation, we recall Li et al. chaotic equations represented with the integer-order derivative by the following equation: x ′ � ax + y + yz, with initial conditions e strange attractor is obtained at a � 1, b � 1, and m � 1 according to Li et al. work. e Lyapunov exponents of the model presented in system (5) are LYE 1 � 0.409, LYE 2 � 0, and LYE 3 � − 1.773. In the original paper, the system presented in equation (5) is called a four-wing chaotic equation and admits three real equilibrium points, which are unstable and have two other complex equilibrium points. For more explanations related to the system in equation (5), see [16]. In our present work, we consider the memory effect generated by the fractional operators. erefore, we investigate the following modified fractional-order chaotic system: with the initial condition defined by the following relation: Our objective is to give a new fractional-order system; therefore, we set the following values for the parameters of model (7); there are a � 1, b � 1, c � 1, and m � 0. In terms of comparison, we can observe that equation (7) is a fractional version of equation (5). Furthermore, we can observe when the order of the Caputo derivative is κ � β � α � 1 in equation (7); we recover the constructive equation of original equation (5) presented in [16]. e difference is in the initial conditions. is paper will be to focus on detecting the chaotic behavior when the fractional-order derivative is utilized. e fractional-order operator impact will be analyzed and proved by the comparison of the attractors' geometries. e variation of the model's parameters such as a and c generate strange actuators and will be confirmed using the bifurcation diagrams. We will also establish via the stability analysis of the equilibrium points of our considered chaotic system; there exists a chaotic region according to the variation of the order of the Caputo derivative. Equation (7) is said the system with commensurate fractional-order when κ � β � α; else, it is said the equation with incommensurate fractional-order.

Numerical Scheme for the Fractional-Order Chaotic System
In this section, we apply the numerical scheme proposed by Garrapa [23] for our proposed fractional-order equation (7) necessary to depict the phase portraits. is numerical discretization uses the numerical approach of the Riemann-Liouville integral. In this section, we consider equation (7) with commensurate fractional-order. e first step will be to determine the analytical solution of fractional differential equation (7) using Riemann-Liouville integral: We set the following function obtained from model (7): For the graphical representations, we introduce the discretization at the point t n ; then, equations (9)-(11) can be written in the following form: Complexity z t n � z(0) + I α w t n , x 1 .

(15)
Introducing the step-size h and t n � nh, the discretization of the fractional integrals can be represented for the functions u, v, and w as the following: where the parameters are given by For the simplification of the numerical scheme of the Riemann-Liouville integral, we introduce the first-order interpolant polynomial of the functions u(τ), v(τ), and w(τ) given by the following relations [23]: Replacing equations (20)-(22) into equations (16)-(18), we get, after recursive summation, more simple formulas for the numerical approach of the fractional integrals: where the parameters are given explicitly as the following form: and furthermore, when n � 1, 2, . . ., then the parameter d can be expressed as the following forms: For the numerical schemes which are used to represent the graphics, it is obtained after plugging equations (23)-(27) into equations (13)-(15); we have the following expressions: For the previous representations, it is necessary to recall the numerical discretization for our functions u, v, and w; we have the following: (29) e convergence and stability are not detailed in this section, see [23]; here, we just give information related to these two properties. We assume that x(t n ), y(t n ), and z(t n ) are the approximate numerical solutions of the fractional system represented in equation (7), and x n , y n , and z n are the exact solutions; the residual function with the Caputo derivative is given by the following expression: We can notice that the convergence of the numerical scheme presented in this section is obtained when the stepsize h converges to 0 [23]. e stability of our numerical scheme is obtained from the Lipchitz criterion of the functions u, v, and w.

Illustration of the Numerical Scheme
is section presents the phase portraits of fractional-order system (7) with commensurate and incommensurate fractional orders. e main objective is to illustrate the numerical scheme presented in previous Section 4. We set the 4 Complexity following values for the parameters of the model a � 1, b � 1, and c � 1 and the initial conditions given by equation (9). For Figure 1, the considered order is α � 0.95. For more illustrations, we consider system (7) with commensurate order α � 0.90. e phase portraits are presented in Figure 2.
An interesting procedure of comparison between the geometry of the attractors can be found in [24]. e method consists of comparing the amplitude of the attractors in different orders of the Caputo derivative. e comparison of the attractors' geometry at the orders α � 0.9 and α � 0.95 can be analyzed for model (7). We adopt the sketch presented in [24]. We consider the phase portraits in plane (x, y) in Figures 1 and 2 and construct polygons, calculate their areas, and compare them. We have following Figures 3 and 4.
In Figure 3, amplitudes are EF � 12.2 cm and AB � CD � 4 cm, and in Figure 4, the amplitudes are KL � 10.88 cm and GH � IJ � 3 cm; referring to the dimension and calculating the areas of the polygons AFBDEC and GLHJKI, we notice that the area(AFBDEC) > area(GLHJKI); thus, there exists a significant difference between the attractors at the orders α � 0.9 and α � 0.95. e method adopted in this section proves the influence of the order of the Caputo derivative. To see more the impact of the order of the Caputo derivative, we consider fractional chaotic equation (7) with incommensurate order: System (31) is obtained when κ � 1 and β � 1 in equation (7). We have the phase portraits in following Figure 5.
e Caputo derivative's impact can be observed by comparing the phase portrait in the commensurate scenario and the phase portrait in the incommensurate scenario, which have different geometries. To detect equation (7) chaotic behavior with the orders α � 0.9 and α � 0.95, we calculate the Lyapunov exponents. e Lyapunov exponents will permit us, in general, to recognize chaotic behaviors.
e Jacobian matrix of fractional differential equation (7) with the same order is given by the following matrix, which is a fundamental tool for finding the Lyapunov exponents: e Lyapunov exponents in the context of fractional calculus in the chaotic system using Matlab code are proposed in Danca's work. Following the same proposed algorithm, the Lyapunov exponents for fractional-order equation (7) with commensurate order when α � 0.95 are given as follows: And, its associated Kaplan-Yorke dimension is given as follows: Following the previous results, we can give the nature of the behaviors. First of all, the fractional system is dissipative because the sum of all Lyapunov exponents gives a negative number. e second property is to detect the chaotic behavior of our presented equation (7) with the same order. We can see with equation (33), we have one positive Lyapunov exponent corresponding to the chaotic behavior of equation (7). We repeat the same analysis with the order α � 0.9; the Lyapunov exponents are represented as the following form: And, its associated Kaplan-Yorke dimension is given as follows: e system is dissipative as well because the addition of all Lyapunov exponents is negative. For the nature of the behaviors, we analyze equation (35); we have the same conclusion as in the previous case, equation (7) has chaotic behavior with the order α � 0.9. It is caused by the existence of one positive Lyapunov exponent.
is section's main conclusion is when the fractional derivative orders vary in the interval [0.9, 1), fractional equation (7) with the same order has chaotic behavior. To confirm this assumption, we calculate in following Table 1 the Lyapunov exponents versus the fractional-order derivative's variation into [0.9, 1).
Checking the results on Table 1, we conclude that fractional-order equation (7) is dissipative at all orders in the interval [0.9, 1). It is because, at all orders, the addition of the Lyapunov exponents is negative. We observe one positive Lyapunov exponent at all orders α, which means the fractional-order system's solutions have chaotic behavior.

Stability of the Equilibrium Points in Fractional Context Calculus
In this part, we discuss the fractional-order system's stability by using the Matignon criterion. It is a famous criterion used in fractional calculus. Stability analysis is one of the fundamental points in the chaotic and hyperchaotic system. In general, all equilibrium points when the system is chaotic fail to be stable. In this section, we consider equation (7) with    6 Complexity commensurate order. It is also interesting to focus on the stability of its equilibrium points. e Jacobian matrix given in the previous section will be used again. We have the matrix described by the following expression: e following points are the real equilibrium points when you exclude the complex eigenvalues: (38) e principle of the method consists of evaluating the Jacobian matrix at all points previously mentioned. For the first point P 1 , the Jacobian matrix is as follows: e previous matrix has eigenvalues as the following numbers λ 1 � 1, λ 2 � 0, and λ 3 � − 1. We notice that the third eigenvalue satisfies |arg(λ 3 )| � π > (απ/2). e first eigenvalue obeys that |arg(λ 1 )| ≈ 0 < (απ/2) for all α. us, the equilibrium point is not stable for all orders α of the Caputo derivative. e matrix in equation (37) at the point P 2 is described as the following form: e previous matrix has eigenvalues as the following numbers λ 1 � 0.1230 + 1.8776i, λ 2 � 0.1230 − 1.8776i, and λ 3 � − 2.246. We remark that the third eigenvalue satisfies |arg(λ 3 )| � π > (απ/2). e first eigenvalue obeys to    Complexity |arg(λ 1 )| ≈ (43π/90) > (απ/2) for all α < 0.9 and the last eigenvalue |arg(λ 2 )| ≈ (43π/90) > (απ/2) for all α. us, the equilibrium point is not stable when α exceeds 0.9.

Bifurcation Diagrams and Lyapunov Exponents
is section analyzes the impact caused by the variation of the parameters a and c on chaotic behavior. Small variations are applied to a and c, and we depict the responses via the bifurcation diagrams plus phase portraits. In this section, for more clarity in the curve, we consider the step-size h � 0.005. All findings are obtained with equation (7) in commensurate order.
We begin our analysis with the variation of the parameter a in the interval (0, 5). We fix b � 1 and c � 1. In Figure 6, we represent the bifurcation diagram versus the parameter a.
We notice that when the parameter a varies in the region [0.6, 5), fractional system (5)-(8) has chaotic behavior. Note that the variation of the parameter a into  (0, 0.5) generates the solutions' convergence to a stable equilibrium point. For an illustration of the chaotic behavior, we represent the phase portraits of our system at a � 1.5 and fixing b � 1 and c � 1 at the order α � 0.95, see Figure 7. We confirm the chaotic behavior with the existence of one positive and large Lyapunov exponent: And, its associated Kaplan-Yorke dimension is given as follows: To validate the convergence of the solutions to a stable equilibrium point, we study the local stability of the nontrivial equilibrium point of equation (7) with commensurate order considering a � 0.4. e equilibrium points of equation (7) e eigenvalues of the previous matrix are given by λ 1 � 0.4, λ 2 � 0, and λ 3 � − 1. We can remark that the first eigenvalue does not obey to Matignon criterion; thus, the point O(0, 0, 0) is not local stable. us, the convergence of the solutions of equation (7) is not at the point O(0, 0, 0). We now try to a point A(− 1.18, − 1.18, − 1.4); thus, the Jacobian matrix is given by We obtain as eigenvalues the following values λ 1 � − 0.0763 + 1.4507i, λ 2 � − 0.0763 − 1.4507i, and λ 3 � − 1.8474. We can observe that all eigenvalues have negative real parts. us, the Matignon criterion is satisfied by all eigenvalues. us, the point A (− 1.18, − 1.18, − 1.4) is a stable equilibrium point. We finish with the point B(1.18, 1.18, − 1.4).
e Jacobian matrix at the point B is described as the following form: We finish this section by the bifurcation diagram generated by the variation of the parameter c ∈ (0, 25), and we illustrate the impact of the parameter c on the phase portraits, see Figure 8, taken at a � 1 and c � 4. e bifurcation diagram shown in Figure 9 indicates that , in our considered interval (0, 25), the system has continuously chaotic behavior. Let us c � 4; the chaotic behavior at this point is confirmed by the existence of one positive Lyapunov exponent; we have the following calculations: And, its associated Kaplan-Yorke dimension is given as follows: For more confirmation of the chaotic region, we depict our fractional system's phase portraits at a � 1 and c � 0.5, see Figure 10.

Electrical Circuit Schematic of the Fractional System
In this section, we give the circuit representation associated to chaotic system (7) with incommensurate order. For more details in the context of the fractional chaotic circuit, we consider equation (7) with κ � β � 1 written using electronic tools: where C i for i � 1, . . . , n denotes the value of the capacitors and R i represents the resistance value. To connect the results in this paper with the experimental results, we draw the circuit schematic and generate the phase portraits using oscilloscopes. We first use our numerical scheme in the previous sections to draw the phase portraits. We consider the order of the Caputo derivative as α � 0.95. e portraits of model (49)-(51) are represented in following Figure 11. For the circuit schematic, the values of the capacitors are given by C 1 � C 2 � C 3 � 1 nF, and the values of the resistors are given R 1 � 250 kΩ, R 2 � 250 kΩ, R 3 � 6.25 kΩ, R 4 � 6.25 kΩ, R 5 � 6.25 kΩ, R 6 � 250 kΩ, and R 7 � 6.25 kΩ.
In equation (52), the third-order approximation is used, but more suitable approximations can also be used. is section's objective is the circuit schematic to validate our theoretical results; therefore, third-or fourth-order approximations can be used for illustration. e method of       (52) is not unique. ere also exists the 4 th -order approximation, which can generate a satisfactory magnitude and phase response over a considerable frequency band [26]. For more related works addressing the approximation of the fractional-order capacitor via fourth-order continued fraction expansion, see Mishra et al.'s investigation [26]. Utilizing the transfer function, the above function can be written as the following formula: Identifying two above equations (45) and (46), we have the following values for the resistors and the capacitors for the fractional integrator, C 3 � 1.27 nF, C 4 � 4.7 nF, C 5 � 3.63 nF, R 9 � 700 MΩ, R 10 � 1.5 MΩ, and R 15 � 15 kΩ. e circuit schematic is done with the aids of 14 resistors, five capacitors, four multipliers for the 4 nonlinear functions, and others, see Figure 12.
e results corresponding to the simulation of the electronic circuit in the Multisim platform are shown in following  We can observe that the results with Matlab after simulating the numerical scheme in Figure 11 and the findings corresponding to the simulation of the electronic circuit in the Multisim platform shown in Figures 13-15 are in good agreement.
is section connects the theoretical results with experimental results.
Before closing this section, we will implement the integer version of model equations (49)-(51). erefore, the order α should be α � 1. For the confirmation of our results obtained after the simulation of the electronic circuit in the Multisim platform, we first simulate the following differential equations:   Under the initial conditions given by system (9) and with a � 1 and c � 1, the phase portraits are represented in following Figure 16. For simulation of system (54)-(56), we rewrite equations (54)-(56) by using resistors and capacitors; we have the following representation: where the values of capacitors are given C 1 � C 2 � C 3 � 1 nF and resistors are given by R 1 � 250 kΩ, R 2 � 250 kΩ, R 3 � 6.25 kΩ, R 4 � 6.25 kΩ, R 5 � 6.25 kΩ, R 6 � 250 kΩ, and R 7 � 6.25 kΩ. e schematic of the above system is represented in following Figure 17. e phase portraits obtained after the simulation of the circuit in Figure 17 are in following  We can observe that the theoretical results 16 and the experimental results 18, 19, and 20 are in good agreement.

Sensitivity to Initial Conditions
In this section, we study the sensitivity of our model when the initial conditions change permanently. By the definition of chaos, the chaotic systems are very sensitive to the starting conditions. In this section, we try to measure the impact of the initial conditions via the bifurcation diagrams. For the rest of the paper, we consider equation (7) with commensurate order α � 0.95. e comparisons are made with the conditions given by We first consider the variation of the parameter a, and we influence the last initial condition z 0 . e bifurcation map21 associated with the initial conditions (0.2, 0.2, 0.2) (blue) and (0.2, 0.2, 0.5) (red) versus the variation of the parameter a is given in following Figure 21.
e present bifurcation diagram, see Figure 21, informs us that the chaotic regions do not change significantly when the parameter a varies. Note that, for the initial condition (0.2, 0.2, 0.2), the chaotic region is (0.6, 5), and for the initial condition (0.2, 0.2, 0.5), the chaotic region is (0.7, 5). e initial conditions' impact can be observed at the beginning of the bifurcation diagram, in the interval (0, 1), where system (7) trajectories converge to a stable equilibrium point. We notice that, with the first initial condition (0.2, 0.2, 0.2), the system enters first in the chaotic region. e small perturbation on the starting condition can generate changes but not significantly when the parameter varies. To see the initial condition's influence according to the variation of a, we illustrate by representing the phase portraits at a � 0.6 and a � 0.7. When a � 0.6, we can observe that, in Figure 22, the phase space in (x, y) plane (blue) corresponding to (0.2, 0.2, 0.2) is chaotic, contrary to the phase space22 in plane (x, y) (red) corresponding to (0.2, 0.2, 0.5), where clearly the solutions of the system converge to an stable equilibrium point. Let us now a � 0.7; in Figure 22, the phase space in the plane (x, y) (blue) corresponding to the condition (0.2, 0.2, 0.2) is chaotic, and the phase space22 in the plane (x, y) (red) corresponding to (0.2, 0.2, 0.5) is also chaotic.
In the second, we suppose y 0 � 0.5, x 0 � 0.2, and z 0 � 0.2 are fixed. e variation of the parameter a is also considered in the bifurcation diagram. e sensitivity to the starting condition can also be observed in the bifurcation diagram. e bifurcation diagram, see Figure 23, associated with the conditions (0.2, 0.2, 0.2) (blue) and (0.2, 0.2, 0.5) (red) versus the variation of the parameter a is given in following Figure 23. e bifurcation diagram, see Figure 23, informs us the chaotic regions (0.6, 5) for the starting condition of equation (7) are given by (0.2, 0.2, 0.2) and (0.7, 5) for the new condition (0.2, 0.5, 0.2). e influence of the initial conditions can be noticed in the interval (0.6, 0.7), after this region equation (7)  We finish this section with the perturbation of the first initial condition. e bifurcation diagram, see Figure 25, 14 Complexity   e bifurcation diagram presented in Figure 25 can be interpreted as with the condition (0.2, 0.2, 0.2), the chaotic region does not change, and with (0.5, 0.2, 0.2), the new chaotic region is also the interval (0.7, 5). In bifurcation diagrams, see Figures 21,23, and 25, the general observation is that the chaotic regions do not significantly change when the parameter a varies. We also observe that the chaotic    16 Complexity   us, the small variation in the starting conditions does not affect the chaotic region of our present model.

Conclusion
In this paper, we have proposed a new system of differential equations that exhibit chaotic behavior and then studied with the Caputo fractional derivative. We have discussed the nature of the chaos for different orders of the fractional derivative using Lyapunov exponents. It is noticed that the fractional-order system considered in this paper has chaotic behavior when the order of the fractional derivative is into (0.9, 1). e small variation of the model's parameters has been studied using bifurcation diagrams, which play an essential role in the chaos characterization. We also find that the initial condition's choice is crucial because our fractional-order system is very sensitive to the initial conditions. e investigations related to the Lyapunov exponents in the context of the fractional derivative are an open problem. e incommensurate fractional-order chaotic system can be investigated for future work. For example, how the phase portraits can be obtained in the context of incommensurate order. How can the Lyapunov exponents be calculated for a fractional-order chaotic system with incommensurate order? What about the variations of the parameters? What will be the impact of the different fractional-order derivatives?

Data Availability
No data were used to support this study.

Conflicts of Interest
e author declares that there are no conflicts of interest.