Abstract and Applied Analysis

Volume 2014, Article ID 276279, 10 pages

http://dx.doi.org/10.1155/2014/276279

## Analytical Study of Fractional-Order Multiple Chaotic FitzHugh-Nagumo Neurons Model Using Multistep Generalized Differential Transform Method

^{1}Department of Mathematics, Faculty of Science, University of Jordan, Amman 11942, Jordan^{2}Nonlinear Analysis and Applied Mathematics (NAAM) Research Group, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia^{3}Pioneer Center for Gifted Students, Ministry of Education, Jerash 26110, Jordan^{4}Applied Science Department, Ajloun College, Al-Balqa Applied University, Ajloun 26816, Jordan

Received 8 March 2014; Accepted 13 May 2014; Published 12 June 2014

Academic Editor: Dumitru Baleanu

Copyright © 2014 Shaher Momani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The multistep generalized differential transform method is applied to solve the fractional-order multiple chaotic FitzHugh-Nagumo (FHN) neurons model. The algorithm is illustrated by studying the dynamics of three coupled chaotic FHN neurons equations with different gap junctions under external electrical stimulation. The fractional derivatives are described in the Caputo sense. Furthermore, we present figurative comparisons between the proposed scheme and the classical fourth-order Runge-Kutta method to demonstrate the accuracy and applicability of this method. The graphical results reveal that only few terms are required to deduce the approximate solutions which are found to be accurate and efficient.

#### 1. Introduction

Mathematical modeling method of real-life phenomena is widely applied in medicine and biology. Specifically, the understanding of neural system model plays an important role in several branches of medical science and technology such as neuroscience, brain activity, chemical reaction kinetics, and behavior of cardiac tissue [1–6]. So it has attracted many medical researchers over the last two decades in order to understand the biogenesis, mechanism, and function. Despite that the formulation of such systems is considerably simple, the lacking understanding of their complex behaviors remains to be a very challenging task, especially when the results are expected in a very short time. However, the advent of computers significant progress has been made recently to reduce this gap.

Furthermore, the FHN neural system is one of the best mathematical models describing the electrical activity in the field of electrocardiology, which is a simplified model for the qualitative characteristics and dynamical and neuronal investigations of electrical propagation in the myocardium. For a comprehensive introduction in this field, we refer to [7–16].

In this paper, the chaotic FHN neurons model under external electrical stimulation is given by the following three coupled equations with different gap junctions: where and represent the state variables of a neuron representing the activation potential and the recovery voltage, respectively; , , and represent the states of the master, the first slave, and the second slave FHN neuron, respectively; , , and represent the strengths of gap junctions between the master and the first slave neurons, between the master and the second slave neurons, and between the two slave neurons, respectively. Disturbances at the master, the first slave, and the second slave neurons are represented by , , and , respectively. The term represents the external stimulation current with time and angular frequency . Here, we use the angular frequency and the amplitude as dimensionless quantities as specified for FHN neurons model.

The literature on this subject is quite vast, for example, the full FitzHugh model on an infinite domain has been studied in [17]. In [18], the Hopf bifurcations have analyzed FHN model for nerve conduction. The dynamics of uncertain coupled chaotic delayed FHN neurons with various parametric variations under external electrical stimulation have been investigated in [19], where separate conditions for single-input and multiple-input control schemes for synchronization of a wide class of FHN systems were provided. In [20], the authors have discussed the synchronization of three coupled chaotic FHN neurons under external electrical stimulation with different gap junctions. Moreover, numerical simulation of the FHN equations has been presented using the variational iteration method and Adomian decomposition method [21]. Whilst, the analytical solutions for the FHN model in the case where a collection of unstable cells is surrounded by a collection of stable cells have been generated in [22].

Nowadays, fractional calculus has been used to model physical and engineering processes, which are found to be best described by fractional differential equations [23–26]. It is worth noting that the standard mathematical models of integer-order derivatives, including nonlinear models, do not work adequately in many cases. More recently, fractional calculus has become a powerful tool to describe the dynamics of chaotic neurons system, which appear frequently in many branches of medical science. Chaotic neurons systems have a profound effect on its approximate solutions and are highly sensitive to time step sizes. Thus, it will be beneficial to find a reliable analytical tool to test its long-term accuracy and efficiency. The multistep generalized differential transform method (MSGDTM) is powerful in investigating approximate solutions of various kinds of these systems.

In this paper, the attention is given to obtain the approximate solution of the fractional-order multiple chaotic FHN neurons model under external electrical stimulation with different gap junctions using the MSGDTM. This method is only a simple modification of the generalized differential transform method (GDTM), in which it is treated as an algorithm in a sequence of small intervals (i.e., time step) for finding accurate approximate solutions to the corresponding systems. The approximate solutions obtained by using the GDTM are valid only for a short time. The ones obtained by using the MSGDTM are more valid and accurate during a long time and are in good agreement with the classical Runge-Kutta method numerical solution when the order of the derivative is one.

The remainder of this paper is organized as follows. In next section, we present basic facts, definitions, and notations related to the fractional calculus and MSGDTM. In Section 3, the MSGDTM is applied to the fractional-order multiple chaotic FHN neurons model. In Section 4, numerical simulation is shown graphically to illustrate the feasibility and effectiveness of the proposed method. Finally, the conclusions are drawn in Section 5.

#### 2. The Multistep Generalized Differential Transform Method (MSGDTM)

To describe the MSGDTM [26–29], we consider the following initial value problem for systems of fractional differential equations: subject to the initial conditions where is the Caputo fractional derivative of order , and , for , and the Caputo fractional derivative of of order with is defined as

for , , , . For more details about the fractional calculus theory, see [30–32].

Let be the interval over which we want to find the solution of the initial value problem (2)-(3). The differential transform of the -order derivative of a function on a subinterval is defined as follows: Using (5), one can easily prove the following corollary.

Corollary 1. *If , then , while if , then .*

In actual applications of GDTM, the -order approximate solution of the initial value problem (2)-(3) can be expressed by the finite series where satisfied the recurrence relation and and are the differential transform of function for .

Assume that the interval is divided into subintervals , , of equal step size by using the nodes . The main ideas of the MSGDTM are as follows. Firstly, we will apply the GDTM to the initial value problem (2)-(3) over the interval . Then, we will obtain the approximate solution , , using the initial condition , for . For and at each subinterval , we will use the initial condition and apply the GDTM to the initial value problem (2)-(3) over the interval . The process is repeated and generates a sequence of approximate solutions , , for . Finally, the MSGDTM assumes the following solution: where

The new algorithm of the MSGDTM is simple for computational performance for all values of . As we will see in the next section, the main advantage of the new algorithm is that the obtained solution converges for wide time regions.

#### 3. Applications of the MSGDTM for the Fractional-Order Multiple Chaotic FHN Neurons Model

To demonstrate the applicability, accuracy, and efficiency of the MSGDTM for solving linear and nonlinear fractional-order equations, we applied this scheme to the fractional-order model of three coupled chaotic FHN neurons with different gap junctions [20], which is the lowest-order chaotic system among all the chaotic systems. Where the integer-order derivatives are replaced by the fractional-order derivatives as follows: where , , and are the strengths of gap junctions between the master and the first slave neurons, between the master and the second slave neurons, and between the two slave neurons, respectively; , , and are the system parameters and and are the state variables of a neuron representing the activation potential and the recovery voltage, respectively; , , and are the states of the master, the first slave, and the second slave FHN neuron, respectively; and , are parameters describing the order of the fractional time-derivatives in the Caputo sense.

By applying the MSGDT algorithm to obtain the numerical solution for the fractional-order multiple chaotic FHN neurons model, the system (10) gives where , , and and are the differential transformation of and , , respectively. The differential transform of the initial conditions are given by , , , , , and . In view of the differential inverse transform, the differential transform series solution for the system (10) can be obtained as

According to the MSGDTM, the series solution for the system (10) is suggested by where , , , , , and , for , satisfy the following recurrence relations: such that , , , , , and .

Finally, starting with , , , , and and using the recurrence relation given in (14), the multistep solution can be obtained as in (13).

#### 4. A Test Problem for the Fractional-Order Chaotic FHN Neurons Model

In this work, we carefully propose the MSGDTM, a reliable modification of the GDTM that improves the convergence of the series solution. The method provides immediate and visible symbolic terms of analytic solutions as well as numerical approximate solutions to both linear and nonlinear differential equations. Moreover, we shall demonstrate the accuracy of the MSGDT scheme against the Mathematica built-in fourth-order Runge-Kutta (RK4) procedure for the solutions of multiple chaotic FHN neurons model in the case of integer order derivatives. The MSGDT scheme is coded in the computer algebra package Mathematica. The Mathematica environment variable digits controlling the number of significant digits are set to 20 in all the calculations done in this paper. The time range studied in this work is and the step size . In this regard, we take the initial condition for chaotic FHN neurons model such as , , , , , and with parameters , and , whilst , , , , , , , , and .

Figure 1 shows the phase portrait for the classical multiple chaotic FHN neurons model, when , using the MSGDT and RK4 methods. However, it can be seen that the results obtained using the MSGDTM match the results of the RK4 method very well, which implies that the MSGDTM can predict the behavior of these variables accurately for the region under consideration. Additionally, Figures 2, 3, and 4 show the phase portrait for the fractional multiple chaotic FHN neurons using the MSGDTM. From the numerical results in Figures 2, 3, and 4, it is clear that the approximate solutions depend continuously on the time-fractional derivative , . The effective dimension of (10) is defined as the sum of orders In the meantime, we can see that the chaos exists in the fractional-order multiple chaotic FHN neurons model with order as low as .

#### 5. Conclusions

In this paper, a multistep generalized differential transform method has been successfully applied to find the numerical solutions of the fractional-order multiple chaotic FitzHugh-Nagumo neurons model. This method has the advantage of giving an analytical form of the solution within each time interval which is not possible using purely numerical techniques like the fourth-order Runge-Kutta method (RK4). We conclude that MSGDT method is a highly accurate method in solving a broad array of dynamical problems in fractional calculus due to its consistency used in a longer time frame.

The reliability of the method and the reduction in the size of computational domain give this method a wider applicability. Many of the results obtained in this paper can be extended to significantly more general classes of linear and nonlinear differential equations of fractional order.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The authors would like to express their thanks to the unknown referees for their careful reading and helpful comments.

#### References

- D. Noble, “Modeling the heart,”
*Physiology*, vol. 19, no. 4, pp. 191–197, 2004. View at Google Scholar · View at Scopus - R. FitzHugh, “Impulses and physiological states in theoretical models of nerve membrane,”
*Biophysical Journal*, vol. 1, no. 6, pp. 445–466, 1961. View at Google Scholar - B. Jobst, “Brain stimulation for surgical epilepsy,”
*Epilepsy Research*, vol. 89, no. 1, pp. 154–161, 2010. View at Publisher · View at Google Scholar · View at Scopus - M. Aqil, K.-S. Hong, M.-Y. Jeong, and S. S. Ge, “Detection of event-related hemodynamic response to neuroactivation by dynamic modeling of brain activity,”
*NeuroImage*, vol. 63, no. 1, pp. 553–568, 2012. View at Publisher · View at Google Scholar · View at Scopus - Y. Rudy and J. R. Silva, “Computational biology in the study of cardiac ion channels and cell electrophysiology,”
*Quarterly Reviews of Biophysics*, vol. 39, no. 1, pp. 57–116, 2006. View at Publisher · View at Google Scholar · View at Scopus - O. Abu Arqub, A. El-Ajou, S. Momani, and N. Shawagfeh, “Analytical solutions of fuzzy initial value problems by HAM,”
*Applied Mathematics & Information Sciences*, vol. 7, no. 5, pp. 1903–1919, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - M. Al-Smadi, O. Abu Arqub, and S. Momani, “A computational method for two-point boundary value problems of fourth-order mixed integrodifferential equations,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 832074, 10 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - J. Wang, T. Zhang, and B. Deng, “Synchronization of FitzHugh-Nagumo neurons in external electrical stimulation via nonlinear control,”
*Chaos, Solitons & Fractals*, vol. 31, no. 1, pp. 30–38, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G.-C. Wu, D. Baleanu, and S.-D. Zeng, “Discrete chaos in fractional sine and standard maps,”
*Physics Letters A*, vol. 378, no. 5-6, pp. 484–487, 2014. View at Publisher · View at Google Scholar · View at MathSciNet - C. Li and G. Chen, “Chaos in the fractional order Chen system and its control,”
*Chaos, Solitons & Fractals*, vol. 22, no. 3, pp. 549–554, 2004. View at Publisher · View at Google Scholar · View at Scopus - G. C. Wu and D. Baleanu, “Chaos synchronization of the discrete fractional logistic map,”
*Signal Processing*, vol. 102, pp. 96–99, 2014. View at Google Scholar - C. Li and G. Chen, “Chaos and hyperchaos in the fractional-order Rössler equations,”
*Physica A*, vol. 341, no. 1–4, pp. 55–61, 2004. View at Publisher · View at Google Scholar · View at MathSciNet - G.-C. Wu and D. Baleanu, “Discrete fractional logistic map and its chaos,”
*Nonlinear Dynamics*, vol. 75, no. 1-2, pp. 283–287, 2014. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. Li and G. Peng, “Chaos in Chen's system with a fractional order,”
*Chaos, Solitons & Fractals*, vol. 22, no. 2, pp. 443–450, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. H. Deng and C. P. Li, “Chaos synchronization of the fractional Lü system,”
*Physica A*, vol. 353, no. 1–4, pp. 61–72, 2005. View at Publisher · View at Google Scholar · View at Scopus - O. Abu Arqub, Z. Abo-Hammour, and S. Momani, “Application of continuous genetic algorithm for nonlinear system of second-order boundary value problems,”
*Applied Mathematics and Information Sciences*, vol. 8, no. 1, pp. 253–248, 2014. View at Google Scholar - J. Rinzel and J. P. Keaner, “Hopf bifurcation to repetitive activity in nerve,”
*SIAM Journal on Applied Mathematics*, vol. 43, no. 4, pp. 907–922, 1983. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Rinzel and J. B. Keller, “Traveling wave solutions of a nerve conduction equation,”
*Biophysical Journal*, vol. 13, no. 12, pp. 1313–1337, 1973. View at Google Scholar - M. Rehan and K.-S. Hong, “Robust synchronization of delayed chaotic FitzHugh-Nagumo neurons under external electrical stimulation,”
*Computational and Mathematical Methods in Medicine*, vol. 2012, Article ID 230980, 11 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Rehan, K.-S. Hong, and M. Aqil, “Synchronization of multiple chaotic FitzHugh-Nagumo neurons with gap junctions under external electrical stimulation,”
*Neurocomputing*, vol. 74, no. 17, pp. 3296–3304, 2011. View at Publisher · View at Google Scholar · View at Scopus - A. A. Soliman, “Numerical simulation of the FitzHugh-Nagumo equations,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 762516, 13 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Artebrant, A. Tveito, and G. T. Lines, “A method for analyzing the stability of the resting state for a model of pacemaker cells surrounded by stable cells,”
*Mathematical Biosciences and Engineering*, vol. 7, no. 3, pp. 505–526, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Freihat and S. Momani, “Application of multistep generalized differential transform method for the solutions of the fractional-order Chua's system,”
*Discrete Dynamics in Nature and Society*, vol. 2012, Article ID 427393, 12 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Freihat and M. Al-Smadi, “A new reliable algorithm using the generalized differential transform method for the numeric-analytic solution of fractional-order Liu chaotic and hyperchaotic systems,”
*Pensee Journal*, vol. 75, no. 9, pp. 263–276, 2013. View at Google Scholar - Z. Odibat and S. Momani, “A generalized differential transform method for linear partial differential equations of fractional order,”
*Applied Mathematics Letters*, vol. 21, no. 2, pp. 194–199, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Momani and Z. Odibat, “A novel method for nonlinear fractional partial differential equations: combination of DTM and generalized Taylor's formula,”
*Journal of Computational and Applied Mathematics*, vol. 220, no. 1-2, pp. 85–95, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. S. Erturk, S. Momani, and Z. Odibat, “Application of generalized differential transform method to multi-order fractional differential equations,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 13, no. 8, pp. 1642–1654, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Freihat and S. Momani, “Adaptation of differential transform method for the numeric-analytic solution of fractional-order Rössler chaotic and hyperchaotic systems,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 934219, 13 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Odibat, S. Momani, and V. S. Ertürk, “Generalized differential transform method: application to differential equations of fractional order,”
*Applied Mathematics and Computation*, vol. 197, no. 2, pp. 467–477, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. S. Miller and B. Ross,
*An Introduction to the Fractional Calculus and Fractional Differential Equations*, A Wiley-Interscience Publication, John Wiley & Sons, New York, NY, USA, 1993. View at MathSciNet - C. Li and W. Deng, “Remarks on fractional derivatives,”
*Applied Mathematics and Computation*, vol. 187, no. 2, pp. 777–784, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Sabatier, O. P. Agrawal, and J. A. Tenreiro Machado,
*Advances in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering*, Springer, Dordrecht, The Netherlands, 2007. View at Publisher · View at Google Scholar · View at MathSciNet