The Effect of a Non-Local Fractional Operator in an Asymmetrical Glucose-Insulin Regulatory System: Analysis, Synchronization and Electronic Implementation

: For studying biological conditions with higher precision, the memory characteristics deﬁned by the fractional-order versions of living dynamical systems have been pointed out as a meaningful approach. Therefore, we analyze the dynamics of a glucose-insulin regulatory system by applying a non-local fractional operator in order to represent the memory of the underlying system, and whose state-variables deﬁne the population densities of insulin, glucose, and β -cells, respectively. We focus mainly on four parameters that are associated with different disorders (type 1 and type 2 diabetes mellitus, hypoglycemia, and hyperinsulinemia) to determine their observation ranges as a relation to the fractional-order. Like many preceding works in biosystems, the resulting analysis showed chaotic behaviors related to the fractional-order and system parameters. Subsequently, we propose an active control scheme for forcing the chaotic regime (an illness) to follow a periodic oscillatory state, i.e., a disorder-free equilibrium. Finally, we also present the electronic realization of the fractional glucose-insulin regulatory model to prove the conceptual ﬁndings.


Introduction
Homeostasis is the tendency of organisms to auto-regular and maintain their internal environment in a stable state [1], For instance, an excellent model to describe the homeostatic process in the organism is the glucose-insulin system [1,2]. On one side, when the glucose level is low, arises diverse pathologies (anxiety, tremors, obfuscation, coma, etc.). On the other side, microvascular damages in the retina, kidney, and neuronal injuries, which lead to chronic renal insufficiency and blindness, are originated by high glucose concentrations. The principal pathology of glucose homeostasis is diabetes. Diabetes can be stated as a chronic disorder provoked when the pancreas does not produce insulin in enough quantities, and also if the body cannot successfully process it, resulting in atypical high blood sugar levels. Because of the autoimmune annihilation of β-cells, the insulin released by the pancreas to the human organism is not enough to maintain a certain healthy level, and therefore, the Type 1 Diabetes Mellitus (T1DM) may be manifested. On the contrary, when the beta-cells can produce supranormal, or average concentrations of insulin, but cannot be adequately used in reducing glycemia, the Type 2 Diabetes Mellitus (T2DM) appears [3].
Health agencies [4] demonstrate people, mainly adults, living with diabetes has almost quadrupled from 1980 to 422 million. This rise is primarily due to the rise in T2DM and factors driving it, including overweight and obesity. Diverse mathematical models using ODEs and OdEs have provided a common path to understand multiple complex systems [5,6]. In the last years, because of the long-memory of fractional-order operators, fractional-order systems have gained extensive attention for describing and understanding physical and biological systems [5,[7][8][9][10]. For instance, the analysis of diverse biological systems, such as bio-impedance, drug diffusion, respiratory tissue, and so forth, was reported in [11] using fractional calculus. In [12], the authors characterize the trade-offs between HIV infection and the tumor-immune system employing fractional-order biological systems. Kheiri and Jafari [13] formulated a fractional-order theory that was focused on the multi-patch HIV/AIDS model to analyze whether human migration has effects on the propagation of the HIV/AIDS outbreak. In [14], the authors reported a fractional-order Izhikevich system to obtain insights on distinct neuronal spike responses, including bursting, fast-, regular-spiking, and chattering, as the no-integer order varies. The functioning at electrical level of a fractional-order system of an isolated β-cell is presented in [15,16]. In [17], was reported the comparison among diverse scenarios, such as integer, constant, no-constant and fractional-order derivatives, in order to explain the memory index. In [18], the transmission issues of a susceptible-infected-recovered model were analyzed. They found a proper yield of memory of the fractional-order systems to forecast the pandemic spread. Finally, [19] discovered that the neuron's firing rate could be emulated with a fractional derivative and a slowly varying of the parameters. Therefore, the unique rat neocortical pyramidal neurons have several time-scales.
Therefore, the study of fractional-order biological models continues been critical for an accurate analysis of several health conditions, as well as being essential to understanding this significant open-topic.
Additionally, the synchronization and collective dynamics play an essential role both in physical [20][21][22] and biological systems [15,[23][24][25][26]. By one side, the synchronization is necessary for systems with stable behavior, for instance, the human heartbeat and respiration [27,28]; nevertheless, on the other hand, the synchronized state could provoke severe pathologies such as the Parkinson's disease where the excessive synchronization correlates with a motor deficit. For the glucose-insulin regulatory system, a synchronization state is mandatory, since this condition is closely related to a disorder-free state where the glucose-insulin concentrations synchronize. In [29], the β-cell synchronization is fundamental to effectuate a pulsatile insulin liberation, which carries more substantial hypoglycemic results compared to constant secretion. In [30], the synchronization between gap-junctions and adjacent cells is essential to limit the heterogeneity and biological noise, thence obtaining a robust activation of the β-cells population within the islet. Pecora and Carroll in [31] showed that two nonlinear dynamical systems could be synchronized by introducing appropriate coupling, since then, a variety of approaches were proposed to deal with synchronized states of fractional-order chaos generators, and these include backstepping, adaptive, and active controls [32,33]. Regarding the chaos synchronization applying the active control method [34], it has been demonstrated that an active controller can be straightforwardly designed to achieve synchronization globally if the nonlinearity of the system is known. Thereupon, it is considered to be a promising control strategy due to its straightforwardness. [33].
In addition, the experimental verification of fractional-order models is a topic that has been attracting the attention of researchers [35][36][37][38][39]. In this scenario, ARM-based embedded systems have become in a central block integrated with non-embedded technologies, such as FPGAs, DSPs, and microcontrollers. Subsequently, now it is possible to implement complex software and hardware functionality on a single chip [40]. Several digital hardware have been reported and verified contemplating integer-order chaotic systems [41][42][43]. Notwithstanding, just a short-list of papers have studied the implementation of fractional chaos oscillators on ARM platforms [44,45].
Motivated by the discussion mentioned above, we analyze the effect of a non-local fractional operator in an asymmetrical glucose-insulin regulatory system. More specifically, the system parameters that were related to hypoglycemia, hyperinsulinemia, T1DM, and T2DM, were studied by using both analytical (stability of equilibrium) and numerical (bifurcation diagrams and basins of attraction) techniques. Additionally, we found that the system presents a chaotic attractor, as demonstrated by its phase portraits and Lyapunov exponent spectrum. Besides, the theoretical insights of the glucose-insulin interaction can be validated by using the proposed electronic implementation based on the embedded SoC ARM Broadcom BCM2837B0. Latter, the synchronization between a chaotic behavior (disorder) and periodic behaviors (average condition) were achieved by applying a simple control strategy.
The sections of the manuscript contain the following. The fractional-order glucose-insulin system, in the sense of Caputo, is presented in Section 2. The stability of the equilibrium points is analyzed in Section 3. Several dynamical analysis for the hypoglycemia, hyperinsulinemia, T1D, and T2D are presented in Section 4. Two synchronization schemes using the active control approach are introduced in Section 5. The ARM implementation is carried out in Section 6, which is followed by the conclusions in Section 7.

Fractional-Order Glucose-Insulin Regulatory System
Because diverse biological systems present memory effects, direct generalization maybe by applying fractional-order differential equations, i.e., arbitrary (non-integer) orders. In this manner, we can get a closer insight into the real phenomena. The benefits of fractional-order are mainly to capture the entire time evolution for physical processes, being a kind of memory index, as well as more degrees of freedom for the resulting models. The analysis conducted in this work is inspired by the glucose-insulin regulatory system in [46], as given in Equation (1). This model was derived from the Ackerman, Bajaj, and Rao, and predator-prey Volterra models. As a consequence, the resulting system is formulated as the trade-off relating the glucose and insulin, but including the beta-cells interaction. Apart from normal metabolic conditions described by the classical relation between prey (glucose) and predator (insulin), the model (1) also characterizes the metabolic disorders as chaotic oscillations. Abnormal biological conditions, including glucose-insulin interactions, can be stated as chaotic evolutions of the dynamical systems, as given in [47][48][49][50], to mention a few. From experimental data, Ref. [51] also discussed that the T1DM may have chaotic behaviors. In this manner, the proposed analysis of the glucose-insulin regulatory system is in the same line of previous studies, but instead applying fractional-calculus theory.
x = −a 1 x + 0.1xy + 1.09y 2 − 1.08y 3 + 0.03z − 0.06z 2 + a 7 z 3 − 0.19, z = a 15 y − 1.35y 2 + 0.5y 3 + 0.42z + 0.15yz, we now introduce a fractional-order version [52], but instead applying the Caputo's definition. Let us consider the fractional differential operator D q f (t) = J n−γ D n f (t), with q > 0 and n ∈ N, where D n describes the n-order derivative, while J γ defines the γ-order integral operator in the sense of Riemann-Liouville as For a function f (t), the fractional-order (γ) integral, with γ ∈ R + can be determined by where t ≥ t 0 and Γ(·) is the Gamma function, defined as Γ(s) = ∞ 0 t s−1 e −t dt.
In this manner, the Caputo operator of fractional-order can be established by Definition 2. The fractional derivative of Caputo with order q for a function f (t) ∈ C n ([t 0 , ∞), R) is given as and when 0 < q < 1 is given, as follows By substituting Equation (4) in Equation (1), we get the fractional-order glucose-insulin regulatory system that is given by where D q i (x, y, z) with i = 1, 2, 3 is the the fractional derivative operator in Caputo's sense [53], q i is the fractional-order satisfying 0 < q i ≤ 1, (x, y, z) describe the population density of insulin, the population size of glucose, and the population size of β-cells, respectively. The representative reduction of insulin levels for glucose deficiency is described by a 1 , whereas the parameter a 8 stands for the impact of insulin on glucose. Additionally, the increment rate of insulin concentrations released by b-cells is provided by a 7 . Finally, the increasing rate of the β-cells provoked, when the glucose levels also grow, is given by a 15 [46]. From the nature of how this system was conceived, the parameters must be nonnegative [46]. The discrete-time version for the proposed system (5) is set as: where with and note that x 0 , y 0 , z 0 are the initial values for, 56 and f 3 (x, y, z) = a 15 y − 1.35y 2 + 0.5y 3 + 0.42z + 0.15yz, whether q i = q, i = 1, 2, 3 the system is said to be commensurate then the convergence order is described as

Stability Analysis of Fractional-Order Glucose-Insulin Model
A non-local differential equation with fractional-order q ∈ (0, 1), normally has a stability region larger than that of the integer-order version with q = 1 [57]. We introduce the following theorems and definitions to discuss the stability of the proposed fractional-order system.

Definition 3. A system denoted as
with trajectory x(t) = 0, where c i are the starting conditions is t −q asymptotically stable if there is a nonnegative real q, so that: Subsequently, fractional-order systems are known to have long memory, since they obey a behavior like t −q , which slowly tends to 0 for the solutions x(t). A particular case of Mittag-Leffler stability is stated as power-law stability t −q . [58,59].

Theorem 3.
The equilibria E * with a instability measure defined as is asymptotically stable if and only if (12) is defined nonpositive. Where λ i 's are the roots of characteristic equation:

Remark 1.
If ρ is nonnegative, the underlying system can show a chaotic evolution since the equilibrium E * is unstable [64,65].
The first step consists of finding the equilibrium points of (5), which are calculated via f(x) = 0, as follows Because of the biological interpretation of state variables [46], the stability analysis for E * = (x * , y * , z * ) is only performed for nonnegative fixed points. Equation (5) has the following Jacobian By setting a 1 = 1.3, a 7 = 2.01, a 8 = 0.22, a 15 = 0.3, we compute the equilibrium points and eigenvalues, as shown in Table 1. As can be noted, the system (5) has two nonnegative fixed points characterized for saddle points of index-1 and index-2, respectively. According to Theorem 1, chaos behavior may arise in the fractional-order glucose-insulin system when | arg(λ)| > qπ/2, holds. Therefore, the minimum commensurate fractional-order leading to chaotic oscillations is q > 0.7638, which is remarkably lower than that reported in [52]. Table 1. Nonnegative equilibrium points and eigenvalues of the system (5).

E i Equilibrium Point Eigenvalues
The chaotic attractor existence and local stability of the equilibria are presented below. The Jacobian matrix considering the equilibria E 1 and parameter a 1 is which generates the characteristic equation given, as follows .395, and b 5 = 90.7232, when considering the Routh-Hurwitz theory, obtaining the following conditions there is a root positive which is unstable and a pair of complex conjugate with the negative real part. Thus E 1 is unstable being a saddle point of index 1.
The presence of a Hopf bifurcation in the system (5) at the equilibrium point E 1 is analyzed by replacing λ = iω in (16) obtaining taking into account the real part of (18), therefore, because a 1 + b 1 > 0, and b 5 , b 4 , a 1 > 0, hence ω 2 < 0, so it is not possible, then, the equilibrium point E 1 does not suffer a Hopf bifurcation.
In order to study the stability of E 1 we consider the discriminant D(P) of characteristic equation P(λ), given as follows where Following the Routh-Hurwitz stability conditions for fractional-order differential equations [66], we establish the following terms: then the equilibrium point E 1 is locally asymptotically stable ∀q ∈ (0, 1). (ii) d 3 > 0 is the necessary condition for the equilibrium point E 1 to be locally asymptotically stable.
When q 1 = q 2 = q 3 ≡ q = 0.9 = 9/10, with w = 10, the characteristic equation of (5) in the equilibrium point E 1 = (0.802, 1.853, 1.286) is given by Theorem 3, as follows if we set a 1 = 1.3 the characteristic equation has an unstable root λ = 1.0434 and the | arg(λ)| < π/20, moreover the E 1 is saddle point. While the characteristic equation at the equilibrium point by considering a 1 = 1.3 we obtain unstable roots λ 1,2 = 1.0771 ± 0.1444i. thence, the instability measure of the system is ρ = (π/2w) − 0.1333 > 0. Therefore, the fractional-order system (5) satisfies the necessary condition for exhibiting a chaotic attractor according to Theorem 3. This can be understood by locating the respective eigenvalues in the complex plane. Figure 1, displays the 27 eigenvalues of the system in the complex plane, the unstable region is delimited by the red lines which is denoted by π/2w. The eigenvalues for E 1 are given in Figure 1a, we can observe that the system has one eigenvalue in the unstable region, it is a saddle point of index 1. Meanwhile, Figure 1b shows the eigenvalues of equilibrium point E 2 ; it has two unstable eigenvalues (saddle point of index 2).   (5) for fractional-orders q 1 = q 2 = q 3 = 0.9. The results were computed by applying the predictor-corrector scheme Adams-Bashforth-Moulton (ABM) [54][55][56]. Due to Caputo's fractional differential operator (4) permits to select both homogeneous and inhomogeneous initial conditions, the ABM algorithm can be executed without particular constraints [56].
As well known, the differential equations that model dynamical systems regularly present symmetries and, therefore, it is rare the solutions do not evolve in a symmetrical orientation. The most common effect is asymmetric attractors for a proper range of parameters, which almost all converges to a singlesymmetric attractor [67]. For dynamical systems with three state variables x, y, and z, there are three types of involuntary symmetries: inversion, rotation, and reflection. We consider the following transformations (x, y, z) → (−x, −y, −z), (x, y, z) → (−x, −y, z), and (x, y, z) → (−x, y, z) corresponding to the invariance of the equations with changes of sign in three, two, and one variable, respectively. As a result, the system (5) does not have a symmetry under the proposed transformations, and its trajectories (x(t), y(t), z(t)) cannot cross the (0, 0, 0) coordinate.

Numerical Analysis of the Non-Local Fractional Operators on Hypoglycemia, Hyperinsulinemia, T1DM, and T2DM
As reported in [46][47][48], when the chaotic behavior appears, it could mean the existence of some disorders in the inherent dynamics of a biological system. In this manner, we analyze four health disorders related to fractional-order glucose-insulin model (5) through one-and two-dimensional bifurcation diagrams and Lyapunov exponents. In particular, we construct a map that relates the fractional-order derivative with a specific parameter, such as a 1 (hypoglycemia), a 7 (hyperinsulinemia), a 8 (T1DM), and a 15 (T2DM).

Hypoglycemia: Parameter a 1 as a Function of Fractional-Order q
For patients with diabetes, hypoglycemia emerges when the reduction of blood glucose concentration reduces below 185 3.9 mmol/L (70 mg/dL) [68]. This is a critical condition, since hypoglycemia may lead to a life-alarming state. In Equation (5), this complication is analyzed when considering the parameter a 1 and the fractional-order q. It means that, if the rate of insulin decrease, which is represented by a 1 in the system (5), gets low, then the hypoglycemia phenomenon emerge. Therefore, we suppose that the underlying system converges into a chaotic behavior as shown in Figure 2. Figure 3a exhibits the bifurcation diagram of system (5) with a fixed fractional-order and considering a 1 as a critical parameter. The bifurcation diagram was made when considering the following: when the state-variable x intersects the Poincaré plane provided by x − p x = 0 with p x = 0.5, the measure r = y 2 + z 2 is delineated. It can be observed that system is stable for values of parameter a 1 > 1.5 but as the parameter diminishes the behavior turns chaotic. Additionally, we observed that the fractional-order q produces a shift concerning the bifurcation diagram showed in Ref. [46]. This consideration exemplifies the importance of considering a fractional-order derivative in the dynamical system, i.e., when values lesser than a 1 = 2.3 are set in the integer-order system [46], chaotic behavior was observed; however, this limit is different for the fractional-order model (a 1 ≤ 1.45). It is at this moment when we could mention that a disorder appears. The numerical results of Lyapunov exponents denoted by λ i with i = 1, 2, 3 are shown in Figure 3b for a 1 = 1.3 and q = 0.9 by applying Wolf's algorithm [69]. The fractional-order glucose-insulin system is chaotic because of the exponents are λ 1 > 0, λ 2 = 0 and λ 3 < 0 with |λ 1 | < |λ 2 + λ 3 |. Those results imply that the system is sensitive to tiny variations of its initial conditions [70,71].
Besides, a two-dimensional bifurcation diagram between the hypoglycemia parameter a 1 and q is presented in Figure 4. The unbounded behavior is represented by green regions, whereas chaos regions are denoted by red color. The black regions lead to healthy behavior (free of hypoglycemia). We found the lower the fractional-order, the lower the effect of a 1 . The basin of attraction in the plane x(0) − y(0) for z(0) = 1, q = 0.9 and a 1 = 1.3, is plotted in Figure 4b, the yellow region stand for a chaotic attractor shown in Figure 2, whereas initial conditions from blue region converge into a unbounded behavior. Finally, Figure 5a-c and Figure 6, presents the phase portraits and Lyapunov exponents, respectively, of healthy behavior for system (5) obtaining a Lyapunov exponent with magnitude zero and two negatives.

Hyperinsulinemia: Parameter a 7 as a Function of Fractional-Order q
Hyperinsulinemia means the quantity of insulin in the blood is higher than normal levels. Hyperinsulinemia is most often caused by insulin resistance, both humans and animals [72]. A condition in which the body is not capable of acts in the right form to the effects of insulin. Consequently, in order to compensate the high blood glucose levels, the pancreatic β-cells irrigate more insulin [73][74][75][76]. Hyperinsulinemia condition is analyzed in the fractional-order glucose-insulin regulatory system (5) by the parameter a 7 . Figure 7a shows the bifurcation diagram for different values of a 7 as a function of fractional order q. From Figure 7a, we can observe that the system is stable when the values for a 7 are small, which describes the increased rate of insulin. If a 7 increases, the system becomes in a chaotic manner, which can be proved by Proposition 1 as follows Proposition 1. When q 1 = q 2 = q 3 ≡ q = 0.95 and a 1 = 2.04, a 7 = 2.4, a 8 = 0.22, a 15 = 0.3, the system (5) exhibits a chaotic attractor.
Besides, the phenomenon antimonotonicity is stated in Figure 7a, which refers to the creation of period orbits followed by their nullification with reverse bifurcation sequences [77]. This phenomenon is one of the most common paths to chaos [78,79]. Antimonotonicity was found in Equation (5) by sweeping a 7 in the interval 2.6 ≤ a 7 ≤ 3.2 with q = 0.95. Additionally, we obtain the Lyapunov exponents λ 1 = 1.6617, λ 2 =0, λ 3 =-24.7646) indicating chaos.
On the other hand, Figure 7b gives the two-dimensional bifurcation diagram between the hyperinsulinemia parameter a 7 and the fractional-order q. The unbounded behavior is represented by the green regions; chaotic behavior is denoted by red regions and the healthy behavior (free of hyperinsulinemia) is given by the blue regions. We found that hyperinsulinemia disorder depends on the value of fractional-order. For values q < 0.92, the hyperinsulinemia tends to periodic oscillations.

Type-2 Diabetes Mellitus: Parameter a 7 as a Function of Fractional-Order q
The abnormal insulin secretion of the pancreatic β-cells is commonly related to T2DM or non-insulin-dependent diabetes mellitus, which is known to be a disorder with insulin resistance [80][81][82]. The interconnection among T2DM, insulin resistance, and obesity relies on the β-cell dysfunction [80,83]. T2DM condition is characterized by the parameter a 8 in (5). Figure 8a shows the bifurcation diagram for the parameter a 8 with a fractional-order q = 0.98. The parameter a 8 is appropriate to understand the insulin resistance of the human body since it describes the effect of emitted insulin on glucose level [46]. In the bifurcation diagram, that phenomenon is detected when a 8 < 0.37, which is associated with chaotic behavior, as demonstrated by Proposition 2. Two-dimensional bifurcation diagram for a 8 and fractional-order q, the chaotic behavior is denoted by red regions, and the periodic behavior (healthy behavior) is given by the blue regions.
Additionally, we compute the Lyapunov exponents λ 1 = 0.56, λ 2 = 0, λ 3 = −21.03. Figure 8b sketches the two-dimensional bifurcation diagram for the T2DM parameter a 8 and the fractional-order q. Analogous previous cases, the red areas evolve to chaos, whereas the blue regions converge to a stable behavior (healthy condition). There is a linear fit between the fractional-order and T2DB. The lower the fractional-order, the lower the value for a 8 , where the T2DM disorder is observed. Besides, for q < 0.97 the T2DM disappear for 0.25 < a 8 < 0.7. These results suggest that the T2DM is not presented in the glucose-insulin system (5) for lowers fractional-orders.

Type-1 Diabetes Mellitus: Parameter a 15 as a Function of Fractional-Order q
T1DM is a common autoimmune disease that originates when the pancreatic β-cells cannot produce insulin at normal levels, and patients will require hormone dosage for their entire life. [84]. The fractional-order system (5) exhibits this condition when the density of β-cells distinguished by a 15 reduces and, therefore, the pancreas may not secrete sufficient insulin to stabilize the glucose concentration.
The bifurcation diagram of Equation (5) is shown in Figure 9a when considering a 15 as critical parameter with q = 0.95. As can be seen, the system (5) exhibits different types of steady behaviors for specific values of a 15 . However, whether this parameter decreases, the system behaves chaotically, as is demonstrated in Proposition 3.
The Lyapunov exponents when q = 0.95 and a 1 = 2.04, a 7 = 2.01, a 8 = 0.22, a 15 = 0.26 are λ 1 = 1.7733, λ 2 = 0, and λ 3 = −24.3966. Similarly previous case, Figure 9b presents the two-dimensional bifurcation diagram relating a 15 and the fractional-order q. The green, red, and blue colors denote unbounded, chaotic, and steady-state behaviors, respectively. A healthy behavior, free of T1DM, is found for fractional-orders lowers than q < 0.925. Those results may imply that lower fractional-orders mitigate the effect of the reduction of population density of β-cells for T1DM.

Synchronization between Fractional-Order Glucose Insulin Systems
Synchronization is a nonlinear phenomenon that was observed in biological systems; it is seen on isolated cells [15], clusters of cells as in organisms, and even in collective dynamics of populations [25]. Regarding the glucose-insulin system, it has been shown pancreatic β-cells also present a collective behavior whose synchronization underlies a small-world functional organization [24][25][26]. Thus, the synchronization is crucial to effectuate a pulsatile insulin liberation in cells, which guarantees more substantial hypoglycemic effects. Hence, we study the synchronization between fractional-order glucose-insulin regulatory systems. We expect that the synchronization state converges into a periodic behavior, because it is the typical response in a suband blood glucose concentrations. We define thject with normal metabolic conditions, allowing with this, the synchronization between the insulin e drive and response system, as follows D q 2 y 1 = −a 8 x 1 y 1 + 3.84x 2 1 + 1.2x 3 1 + 0.3y 1 (1 − y 1 ) − 1.37z 1 + 0.3z 2 1 − 0.22z 3 1 − 0.56, D q 3 z 1 = a 15 y 1 − 1.35y 2 1 + 0.5y 3 1 + 0.42z 1 + 0.15y 1 z 1 , and where u 1 , u 2 , u 3 in (30) represents the unknown control terms, and the error can be defined by To achieve the synchronization, it is essential that the errors e i → 0 as t → ∞ with i = 1, 2, 3. Equation (31), together with (29) and (30), yield the error system D q 1 e 1 = −â 1 x 2 + 0.1x 2 y 2 + 1.09y 2 2 − 1.08y 3 2 + 0.03z 2 − 0.06z 2 2 +â 7 z 3 2 + +a 1 x 1 − 0.1x 1 y 1 − 1.09y 2 1 + 1.08y 3 Let us define the active control functions u i with i = 1, 2, 3 where the linear functions V 1 , V 2 , V 3 are given by By using (33) and (34), the error system (32) becomes The synchronization error vanishes eventually because of the eigenvalues are −1, −1, −1 in Equation (35).
The synchronization scenario is as follows. The drive system has a periodic behavior, while the response system is in a chaotic state. We study the Type-1 Diabetes Mellitus (parameter a 8 ), since it is the most common disorder, and affects most world population as well as it is correlated with obesity. For this case, a 8 = 0.5 andâ 8 = 0.27 for drive and response systems, respectively, while a 1 =â 1 , a 7 =â 7 , and a 15 =â 15 . The fractional-order are q 1 = q 2 = q 3 = q = 0.95 in both systems with x 1 (0) = 0.53, y 1 (0) = 1.31, z 1 (0) = 1.03 and x 2 (0) = 0.5, y 2 (0) = 1.1, z 2 (0) = 1.3, for drive and response systems, respectively. Figure 10a-c show the phase planes between the periodic (free of T1DM) and chaotic (with T1DM) systems. Additionally, the synchronization error by considering (36) is given in Figure 11. Due to the error tends to zero as time evolves, we infer that the proposed control strategy is suitable for forcing the system with the disorder to a state free of T1DM. It is worth noting that the control strategy can be extended to incorporate uncertainties and improve the robustness of the synchronization using other approaches, as shown in [22,85]. From a practical biological point of view, for instance, recent works have employed optical-based control using a light-activated Na+ channel, to attain insulin from β-cells both in-vitro and vivo [86,87]. Therefore, our results could be useful for future works where the glucose-insulin system could be controlled with an artificial control signal.

Physical Realization of the Fractional-Order Glucose-Insulin System Based on an ARM Processor
As well known, the experimental realization of fractional-order dynamical systems is a hot topic that has been attracting the attention of researchers since it is a path for demonstrating the complex dynamics, including chaos [35][36][37][38][39]44,45]. For fractional-order systems, there three typical approaches for getting electronic circuits: frequency-domain approximation, numerical algorithms, and the Adomian decomposition method [35][36][37][38][39]44,45]. The first-mentioned is not recommended for chaos detection, since it may induce incorrect results [62,63]. On the other hand, the second and latter approaches are good options for physically implementing fractional-order systems in re-programmable digital hardware [44,45]. Therefore, we chose the numerical algorithm approach for programming the ABM method. Subsequently, we select herein the ARM SoC Broadcom BCM2837B0 for the experimental verification of the fractional-order glucose-insulin regulatory system. The SoC contains an ARM core with 64-bit. An SDRAM LPDDR2 with 1GB. The ARM cores are capable of running at up to 1.4 GHz. It's possible to create an interface by using the GPIO port with a 16-bit monotonic voltage output D/A converter AD569. Figure 12a,b present the block diagram of the working principle and the main instructions of the pseudo-code, respectively.
After initializing the ARM processor, we set the initial values, h, q, x 0 , y 0 , z 0 . Because of the negative values of the system (1), a positive integer φ is needed to offset the time-series of the state-variables to avoid losing data. In this manner, all computed data are now positive. Next, we multiply the data by a positive integer γ to fit them to the DAC resolution of 2 14 bits. Finally, the obtained results visualize in an oscilloscope, as shown in Figures 13 and 14. We analyze the scenario related to hypoglycemia. First, we implement the case where system (1) presents the hypoglycemia condition, as given in Figure 13. As can be seen, the experimental phase portraits are pretty similar to those that are shown in Figure 1. Finally, the case when the fractional-order system (1) is free of hypoglycemia, i.e., an steady-state and, therefore, convergers to a periodic attractor, is given in Figure 14. Similar to the previous case, the experimental results are in good agreement with Figure 4. Subsequently, it indicates that the fractional-order glucose-insulin regulatory system was successfully realized on an ARM digital platform.

Conclusions
The dynamical analysis, synchronization, and physical realization of a glucose-insulin regulatory system has been presented by using Caputo's non-local fractional-order operator. In particular, we studied four common disorders, such as T1DM, T2DM, Hypoglycemia, and Hyperinsulinemia. We found that the fractional-order system switches between a chaotic behavior (a health disorder) and a disorder-free state, not only for the values of systems parameters, but also as a function of the fractional-order, due it adds more degrees of freedom in the model.To understand that insight, we computed two-dimensional bifurcations diagrams, which demonstrated the importance of considering the fractional-order (memory index) for getting a higher approximation of the observed behavior because fractional-order systems describe the whole-time domain in the solution, while the integer-order model is related to the local properties. Additionally, a phenomenon of antimonotonicity was observed in the parameter related to the hyperinsulinemia case. Besides, by applying the straightforward active control method, we showed that stable behavior in the fractional glucose-insulin system under the T1DM condition could be attained when it synchronizes with a disorder-free system. We remark that the synchronization can be extended to the remaining conditions. Finally, the electronics approach-based validation of chaotic and periodic behaviors was shown using an ARM digital platform. The experimental observations were in good agreement with the theoretical findings.
In this manner, the system-of-a-chip circuit designs are the right candidate for exploiting the advantages of the fractional-order models due to their simplicity, programmable characteristics, and portability, therefore increasing the fractional-order-based oncoming applications. As future work, an analysis related to robustness of the synchronization scheme will be developed.