International Journal of Engineering Mathematics
Volume 2015 (2015), Article ID 107978, 35 pages
http://dx.doi.org/10.1155/2015/107978

Research Article

New Improvement of the Expansion Methods for Solving the Generalized Fitzhugh-Nagumo Equation with Time-Dependent Coefficients

Jalil Manafian and Mehrdad Lakestani

Department of Applied Mathematics, Faculty of Mathematics Science, University of Tabriz, Tabriz 51569 84687, Iran

Received 30 May 2015; Revised 10 September 2015; Accepted 10 September 2015

Academic Editor: Josè A. Tenereiro Machado

Copyright © 2015 Jalil Manafian and Mehrdad Lakestani. 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

An improvement of the expansion methods, namely, the improved -expansion method, for solving nonlinear second-order partial differential equation, is proposed. The implementation of the new approach is demonstrated by solving the generalized Fitzhugh-Nagumo equation with time-dependent coefficients. As a result, many new and more general exact travelling wave solutions are obtained including periodic function solutions, soliton-like solutions, and trigonometric function solutions. The exact particular solutions contain four types: hyperbolic function solution, trigonometric function solution, exponential solution, and rational solution. We obtained further solutions comparing this method with other methods. The results demonstrate that the new -expansion method is more efficient than the Ansatz method and Tanh method applied by Triki and Wazwaz (2013). Recently, this method is developed for searching exact travelling wave solutions of nonlinear partial differential equations. Abundant exact travelling wave solutions including solitons, kink, and periodic and rational solutions have been found. These solutions might play an important role in engineering fields. It is shown that this method, with the help of symbolic computation, provides a straightforward and powerful mathematical tool for solving the nonlinear physics.

1. Introduction

Nonlinear evolution equations (NLEEs) are very important model equations in mathematical physics, engineering, and applied mathematics for describing diverse types of physical mechanisms of natural phenomena in the field of applied sciences, biochemistry, and engineering. Nonlinear wave equations play a major role in various fields such as plasma physics, fluid mechanics, optical fibers, solid state physics, chemical kinetics, geochemistry, and nonlinear optics [1, 2]. Much work has been done over the years on the subject of obtaining the analytical solutions to the nonlinear partial differential equations (NPDEs). One of the most exciting advances of nonlinear science and theoretical physics has been the development of methods to look for exact solutions for NPDEs. The advent of symbolic computation also enables performing some complicated and tedious algebraic processes coupled with differential calculations. With the rapid development of nonlinear sciences based on computer algebraic system, many effective methods have been presented, such as the homotopy analysis method [3, 4], the variational iteration method [5–7], the homotopy perturbation method [8, 9], the sine-cosine method [10], the tanh-coth method [11–13], the modified extended tanh-function method [14, 15], the Exp-function method [16, 17], the exp()-expansion method [18, 19], the -expansion method [20, 21], the modified simple equation method [22], the novel -expansion method [23], the new approach of the generalized -expansion method [21], the Jacobi elliptic function method [24], and the homogeneous balance method [25]. Recently, Naher and Abdullah [21] presented an effective and straightforward method, called the new approach of generalized -expansion method, to obtain exact travelling wave solutions of nonlinear evolution equations. Here, we use an effective method, namely, the improved -expansion method, for constructing a range of exact solutions for the following ordinary partial differential equations that in this paper we developed solutions as well. In this paper, we put forth the new approach of improved -expansion method to construct exact travelling wave solutions including solitons, kink, and periodic and rational solutions to the generalized Fitzhugh-Nagumo equation with time-dependent coefficients. In [22, 26], the standard form of the Fitzhugh-Nagumo equation is given:

where and is of the unknown function depending on the temporal variable and the spatial variable . When , (1) reduces to the real Newell-Whitehead equation which describes the dynamical behavior near the bifurcation point for the Rayleigh-Bénard convection of binary fluid mixtures [27]. The FN equation has various applications in the fields of flame propagation, logistic population growth, neurophysiology, branching Brownian motion process, autocatalytic chemical reaction, and nuclear reactor theory [28]. This equation combines diffusion and nonlinearity which is controlled by the term . Many physicists and mathematicians have paid much attention to the Fitzhugh-Nagumo equation in recent years due to its importance in mathematical physics. Remarkably, this nonlinear evolution equation is an important nonlinear reaction-diffusion equation and is usually used to model the transmission of nerve impulses [29–31]. Shih et al. utilized the approximate conditional symmetry method to determine approximate solutions admitted by a perturbation of (1). Equation (1) has been presented by Hariharan and Kannan via Haar wavelet method [32]. Kawahara and Tanaka [33] have found new exact solutions of (1), by applying the nonclassical symmetry reductions approach by using Hirota method. In [34], Nucci and Clarkson have obtained some new exact solutions with Jacobi elliptic function. Li and Guo [35] have obtained the new exact solutions of the Fitzhugh-Nagumo equation by using first integral method. In this paper, we consider that the generalized Fitzhugh-Nagumo equation with time-dependent coefficients is given as

where , , and are arbitrary functions of . The nonlinear models with variable coefficients are needed to describe the propagation of pulses. Triki and Wazwaz [36] obtained a new variety of soliton solutions by means of specific solitary wave Ansatz and the tanh method for (2). Also, Jiwari et al. [37] applied polynomial differential quadrature method to find the numerical solution of the generalized Fitzhugh-Nagumo equation with time-dependent coefficients in one-dimensional space. In [38], the effect of diffusion on pattern formation in FitzHugh-Nagumo model has been surveyed. Van Gorder applied the method of homotopy analysis to study the Fitzhugh-Nagumo equation [39]. In [40], a FitzHugh-Nagumo monodomain model has been used to describe the propagation of the electrical potential in heterogeneous cardiac tissue. Abbasian et al. have searched symmetric bursting behaviors in the generalized FitzHugh-Nagumo model [41]. Ray and Sahooz used the fractional subequation method for solving the space-time fractional Zakharov-Kuznetsov and modified Zakharov-Kuznetsov equations [42]. Also, the authors of [43] applied the fractional subequation method to the time fractional KdV-Zakharov-Kuznets and space-time fractional modified KdV-Zakharov-Kuznetsov equations. The purpose of this paper is to obtain exact solutions of the generalized Fitzhugh-Nagumo equation and to determine the accuracy of the improved -expansion method in solving these kinds of problems. The paper is organized as follows. In Section 2, we describe the improved -expansion method. In Section 3, we examine the generalized Fitzhugh-Nagumo equation with method introduced in Section 2. Also, conclusion and advantages are given in Section 4. Finally, some references are given at the end of this paper.

2. Description of Improved -Expansion Technique

Step 1. We suppose the given nonlinear fractional partial differential equation for to be in the form

which can be converted to an ODE:

the transformation, , is wave variable. Also, is constant to be determined later.

Step 2. Suppose the traveling wave solution of (4) can be expressed as follows:

where and are constants to be determined, such that , and satisfies the following ordinary differential equation:

We will consider the following special solutions of (6):

 Family  1: When and , then  Family  2: When and , then  Family  3: When and , then  Family  4: When and , then  Family  5: When and , then  Family  6: When and , then  Family  7: When and , then  Family  8: When , then  Family  9: When , then  Family  10: When and , then  Family  11: When , then  Family  12: When , then  Family  13: When , then  Family  14: When , then  Family  15: When and , then   Family  16: When and , then  Family  17: When and , then ,

where , and are constants to be determined later. But the positive integer can be determined by considering the homogeneous balance between the highest order derivatives and nonlinear terms appearing in (6). If is not an integer, then a transformation formula should be used to overcome this difficulty.

Step 3. Substituting (5) into (4) with the value of obtained in Step 2, collecting the coefficients of , and then setting each coefficient to zero, we can get a set of overdetermined partial differential equations for , and with the aid of symbolic computation Maple.

Step 4. Solve the algebraic equations in Step 3, and then substitute in (5).

3. The Generalized Fitzhugh-Nagumo Equation

We consider the generalized Fitzhugh-Nagumo equation with time-dependent coefficients as follows:

where , and are arbitrary functions of . By using the wave variable , reduce it to an ODE as follows:

where and . Balancing the the linear term of the highest order with the highest order nonlinear term by using homogenous principal, we have

Then, the trail solutions are

Substituting (10) and (6) into (8) and collecting all terms with the same order of together, the left-hand side of (10) is converted into polynomial in . Setting each coefficient of each polynomial to zero, we derive a set of algebraic equations for , and as follows:

Solving the set of algebraic equations using Maple, we get the following results.

Set 1. We have the following

where , and are arbitrary constants. By using (13), Families  7, 9, and 13 can be written, respectively, as

where .

Set 2. We have the following:

where , and are arbitrary constants. By using (17), Families 7, 9, and 13 can be written, respectively, as

where .

Set 3. We have the following:

where , and are arbitrary constants. By using (20), Families 7, 9, and 13 can be written, respectively, as

where .

Set 4. We have the following:

where , and are arbitrary constants. By using (23), Families 7, 9, and 13 can be written, respectively, as

where .

Set 5. We have the following:

where , and are arbitrary constants. By using (26), Families 7, 9, and 13 can be written, respectively, as

where .

Set 6. We have the following:

where , and are arbitrary constants. By using (29), Families 7, 9, and 13 can be written, respectively, as

where .

Set 7. We have the following:

where , and are arbitrary constants. By using (32), Family 2 can be written as

where .

Set 8. We have the following:

where , and are arbitrary constants. By using (35), Family 2 can be written as

where .

Set 9. We have the following:

where , and are arbitrary constants. By using (38), Family 2 can be written as

where + + .

Set 10. We have the following:

where , and are arbitrary constants. By using (41), Family 2 can be written as

where .

Set 11. We have the following:

where , and are arbitrary constants. By using (44), Family 2 can be written as

where .

Set 12. We have the following:

where , and are arbitrary constants. By using (47), Family 2 can be written as

where .

Set 13. We have the following:

By using (50), Families 1, 2, 7, and 14 can be written, respectively, as

where .

Set 14. We have the following:

By using (50), Families 1, 2, 7, and 14 can be written, respectively, as

where .

Set 15. We have the following:

By using (50), Families 1, 2, 7, and 14 can be written, respectively, as

where .

Set 16. We have the following:

By using (50), Families 1, 2, 7, and 14 can be written, respectively, as

where .

Set 17. We have the following:

By using (50), Families 1, 2, 7, and 14 can be written, respectively, as

where .

Set 18. We have the following:

By using (50), Families 1, 2, 7, and 14 can be written, respectively, as

where .

Set 19. We have the following:

By using (50), Families 1, 2, and 14 can be written, respectively, as

where .

Set 20. We have the following:

By using (50), Families 1, 2, and 14 can be written, respectively, as

where .

Set 21. We have the following:

By using (50), Families 1, 2, and 14 can be written, respectively, as

where .

Set 22. We have the following:

By using (50), Families 1, 2, and 14 can be written, respectively, as

where .

Set 23. We have the following:

By using (50), Families 1, 2, and 14 can be written, respectively, as

where .

Set 24. We have the following:

By using (50), Families 1, 2, and 14 can be written, respectively, as

where .

Set 25. We have the following:

By using (50), Families 1, 2, and 14 can be written, respectively, as

where .

Set 26. We have the following:

where , and are arbitrary constants. By using (77), Families 7 and 9 can be written, respectively, as

where .

Set 27. We have the following:

where , and are arbitrary constants. By using (80), Families 7 and 9 can be written, respectively, as

where .

Set 28. We have the following:

where , and are arbitrary constants. By using (83), Families 7 and 9, respectively, can be written as

where .

Set 29. We have the following:

where , and are arbitrary constants. By using (86), Families 7 and 9 can be written, respectively, as

where .

Set 30. We have the following:

where , and are arbitrary constants. By using (89), Families 7 and 9 can be written, respectively, as

where .

Set 31. We have the following:

where , and are arbitrary constants. By using (92), Families 7 and 9 can be written, respectively, as

where .

Set 32. We have the following:

where is an arbitrary constant. By using (95), Families 1, 2, and 14 can be written, respectively, as

where .

Set 33. We have the following:

where is an arbitrary constant. By using (99), Families 1, 2, and 14 can be written, respectively, as

where .

Set 34. We have the following:

where is an arbitrary constant. By using (102), Families 1, 2, and 14 can be written, respectively, as

where .

Set 35. We have the following:

where is an arbitrary constant. By using (105), Families 1, 2, and 14 can be written, respectively, as

where .

Set 36. We have the following:

where is an arbitrary constant. By using (108), Families 1, 2, and 14 can be written, respectively, as

where .

Set 37. We have the following:

where is an arbitrary constant. By using (111), Families 1, 2, and 14 can be written, respectively, as

where .

Set 38. We have the following:

where is an arbitrary constant. By using (114), Families 1, 2, and 14 can be written, respectively, as

where .

Set 39. We have the following:

where is an arbitrary constant. By using (117), Families 1, 2, and 14 can be written, respectively, as

where .

Set 40. We have the following:

where is an arbitrary constant. By using (120), Families 1, 2, and 14 can be written, respectively, as

where .

Remark 1. In Figure 1, we plot two- and three-dimensional graphics of imaginary and real values of (14), which denote the dynamics of solutions with appropriate parametric selections. Also, in Figure 2, we draw two- and three-dimensional graphics of imaginary and real values of (15), which demonstrate the dynamics of solutions with convenient parametric choices. In Figure 3, we plot two- and three-dimensional graphics of (96), which represent the dynamics of solutions with proper parametric values. Moreover, in Figure 4, we draw two- and three-dimensional graphics of imaginary and real values of (97), which display the dynamics of solutions with suitable parametric selections. We plot two- and three-dimensional graphics of Figures 1–4, when and for Figures 5–8 when . To the best of our knowledge, these complex exponential function solutions have not been submitted to the literature in advance. The analytical solutions and figures obtained in this paper give us a different physical interpretation for the nonlinear generalized Fitzhugh-Nagumo equation.

(a)

(b)

(c)

(d)

Figure 1: Graphs of ((a) and (b)) real values and ((c) and (d)) imaginary values of (14) are demonstrated at , and and by considering the values ((a) and (c)) and ((b) and (d)) .

(a)

(b)

(c)

(d)

Figure 2: Graphs of ((a) and (b)) real values and ((c) and (d)) imaginary values of (15) are demonstrated at , and and by considering the values ((a) and (c)) and ((b) and (d)) .

(a)

(b)

(c)

(d)

Figure 3: Graphs of ((a) and (b)) real values and ((c) and (d)) imaginary values of (96) are demonstrated at , and and by considering the values ((a) and (c)) and ((b) and (d)) .

(a)

(b)

(c)

(d)

Figure 4: Graphs of ((a) and (b)) real values and ((c) and (d)) imaginary values of (97) are demonstrated at , and and by considering the values ((a) and (c)) and ((b) and (d)) .

(a)

(b)

(c)

(d)

Figure 5: Graphs of ((a) and (b)) real values and ((c) and (d)) imaginary values of (14) are demonstrated at , and and by considering the values ((a) and (c)) and ((b) and (d)) .

(a)

(b)

(c)

(d)

Figure 6: Graphs of ((a) and (b)) real values and ((c) and (d)) imaginary values of (15) are demonstrated at , and and by considering the values ((a) and (c)) and ((b) and (d)) .

(a)

(b)

(c)

(d)

Figure 7: Graphs of ((a) and (b)) real values and ((c) and (d)) imaginary values of (96) are demonstrated at , and and by considering the values ((a) and (c)) and ((b) and (d)) .

(a)

(b)

(c)

(d)

Figure 8: Graphs of ((a) and (b)) real values and ((c) and (d)) imaginary values of (97) are demonstrated at , and and by considering the values ((a) and (c)) and ((b) and (d)) .

4. Conclusion

In this paper, the new approach of expansion method, namely, the improved -expansion method, has successfully been implemented to investigate the generalized Fitzhugh-Nagumo equation with time-dependent coefficients. Abundant exact travelling wave solutions including solitons, kink, and periodic and rational solutions are attained. It is worth mentioning that some of the newly obtained solutions are identical to already published results. It has been shown that the applied method is effective and more wide ranging than the generalized and improved -expansion method because it gives many new solutions. Therefore, this method can be applied to study many other nonlinear partial differential equations which frequently arise in engineering, mathematical physics, and other scientific real time application fields.

Conflict of Interests

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

References

P. Gray and S. K. Scott, Chemical Oscillations and Instabilities—Non-Linear Chemical Kinetics, Oxford Science Publications, Oxford, UK, 1990.

M. M. Kabir and A. Khajeh, “New explicit solutions for the Vakhnenko and a generalized form of the nonlinear heat conduction equations via exp-function method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 10, pp. 1307–1318, 2009.

M. Dehghan, J. Manafian, and A. Saadatmandi, “Solving nonlinear fractional partial differential equations using the homotopy analysis method,” Numerical Methods for Partial Differential Equations, vol. 26, no. 2, pp. 448–479, 2010.

M. G. Sakar and F. Erdogan, “The homotopy analysis method for solving the time-fractional Fornberg-Whitham equation and comparison with Adomian's decomposition method,” Applied Mathematical Modelling, vol. 37, no. 20-21, pp. 8876–8885, 2013.

M. G. Sakar and H. Ergören, “Alternative variational iteration method for solving the time-fractional Fornberg-Whitham equation,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 39, no. 14, pp. 3972–3979, 2015.

H. Jafari, A. Kadem, and D. Baleanu, “Variational iteration method for a fractional-order Brusselator system,” Abstract and Applied Analysis, vol. 2014, Article ID 496323, 6 pages, 2014.

M. Matinfar, S. Mirzanezhad, M. Ghasemi, and M. Salehi, “Solving the interaction of electromagnetic wave with electron by VIM,” Journal of King Saud University—Science, vol. 27, no. 1, pp. 63–70, 2015.

M. Dehghan and J. Manafian, “The solution of the variable coefficients fourth-order parabolic partial differential equations by the homotopy perturbation method,” Zeitschrift für Naturforschung A, vol. 64, no. 7-8, pp. 420–430, 2009.

E. Cuce and P. M. Cuce, “A successful application of homotopy perturbation method for efficiency and effectiveness assessment of longitudinal porous fins,” Energy Conversion and Management, vol. 93, pp. 92–99, 2015.

A.-M. Wazwaz, “Travelling wave solutions for combined and double combined sine-cosine-Gordon equations by the variable separated ODE method,” Applied Mathematics and Computation, vol. 177, no. 2, pp. 755–760, 2006.

J. Manafian Heris and M. Lakestani, “Solitary wave and periodic wave solutions for variants of the KdV-Burger and the K(n, n)-Burger equations by the generalized tanh-coth method,” Communications in Numerical Analysis, vol. 2013, 18 pages, 2013.

J. Manafian Heris and M. Lakestani, “Exact solutions for the integrable sixth-order drinfeld-sokolov-satsuma-hirota system by the analytical methods,” International Scholarly Research Notices, vol. 2014, Article ID 840689, 8 pages, 2014.

J. Manafian and I. Zamanpour, “Exact travelling wave solutions of the symmetric regularized long wave (SRLW) using analytical methods,” Statistics, Optimization & Information Computing, vol. 2, no. 1, pp. 47–55, 2014.

M. A. Abdou and A. A. Soliman, “Modified extended tanh-function method and its application on nonlinear physical equations,” Physics Letters A, vol. 353, no. 6, pp. 487–492, 2006.

S. A. El-Wakil and M. A. Abdou, “New exact travelling wave solutions using modified extended tanh-function method,” Chaos, Solitons & Fractals, vol. 31, no. 4, pp. 840–852, 2007.

J. Manafian Heris and M. Bagheri, “Exact solutions for the modified KDV and the generalized KDV equations via Exp-function method,” Journal of Mathematical Extension, vol. 4, no. 2, pp. 77–98, 2010.

J. Manafian Heris and I. Zamanpour, “Analytical treatment of the coupled Higgs equation and the Maccari system via Exp-function method,” Acta Universitatis Apulensis, vol. 33, pp. 203–216, 2013.

M. G. Hafez, M. N. Alam, and M. A. Akbar, “Traveling wave solutions for some important coupled nonlinear physical models via the coupled Higgs equation and the Maccari system,” Journal of King Saud University—Science, vol. 27, no. 2, pp. 105–112, 2015.

H.-O. Roshid and M. Azizur Rahman, “The exp()-expansion method with application in the (1+1)-dimensional classical Boussinesq equations,” Results in Physics, vol. 4, pp. 150–155, 2014.

M. Fazli Aghdaei and J. Manafianheris, “Exact solutions of the couple Boiti-Leon-Pempinelli system by the generalized -expansion method,” Journal of Mathematical Extension, vol. 5, no. 2, pp. 91–104, 2011.

H. Naher and F. A. Abdullah, “New approach of -expansion method and new approach of generalized -expansion method for nonlinear evolution equation,” AIP Advances, vol. 3, Article ID 032116, 2013.

A. J. M. Jawad, M. D. Petkovic, and A. Biswas, “Modified simple equation method for nonlinear evolution equations,” Applied Mathematics and Computation, vol. 217, no. 2, pp. 869–877, 2010.

M. N. Alam, M. A. Akbar, and S. T. Mohyud-Din, “A novel (G ′/G)-expansion method and its application to the Boussinesq equation,” Chinese Physics B, vol. 23, no. 2, Article ID 020203, 2014.

Y. Chen and Q. Wang, “Extended Jacobi elliptic function rational expansion method and abundant families of Jacobi elliptic function solutions to -dimensional dispersive long wave equation,” Chaos, Solitons and Fractals, vol. 24, no. 3, pp. 745–757, 2005.

X. Zhao, L. Wang, and W. Sun, “The repeated homogeneous balance method and its applications to nonlinear partial differential equations,” Chaos, Solitons & Fractals, vol. 28, no. 2, pp. 448–453, 2006.

A. Hajipour and S. M. Mahmoudi, “Application of exp-function method to Fitzhugh-Nagumo equation,” World Applied Sciences Journal, vol. 9, no. 1, pp. 113–117, 2010.

J. A. Whitehead and A. C. Newell, “Finite bandwidth, finite amplitude convection,” Journal of Fluid Mechanics, vol. 38, no. 2, pp. 279–303, 1969.

A. H. Bhrawy, “A Jacobi-Gauss-Lobatto collocation method for solving generalized Fitzhugh-Nagumo equation with time-dependent coefficients,” Applied Mathematics and Computation, vol. 222, pp. 255–264, 2013.

M. Dehghan, J. M. Heris, and A. Saadatmandi, “Application of semi-analytic methods for the Fitzhugh-Nagumo equation, which models the transmission of nerve impulses,” Mathematical Methods in the Applied Sciences, vol. 33, no. 11, pp. 1384–1398, 2010.

R. FitzHugh, “Impulses and physiological states in theoretical models of nerve membrane,” Biophysical Journal, vol. 1, no. 6, pp. 445–466, 1961.

J. Nagumo, S. Arimoto, and S. Yoshizawa, “An active pulse transmission line simulating nerve axon,” Proceedings of the IRE, vol. 50, no. 10, pp. 2061–2070, 1962.

G. Hariharan and K. Kannan, “Haar wavelet method for solving FitzHugh-Nagumo equation,” World Academy of Science, Engineering and Technology, vol. 67, pp. 560–564, 2010.

T. Kawahara and M. Tanaka, “Interactions of traveling fronts: an exact solution of a nonlinear diffusion equation,” Physics Letters A, vol. 97, no. 8, pp. 311–314, 1983.

M. C. Nucci and P. A. Clarkson, “The nonclassical method is more general than the direct method for symmetry reductions: an example of the FitzHugh-Nagumo equation,” Physics Letters A, vol. 164, no. 1, pp. 49–56, 1992.

H. Li and Y. Guo, “New exact solutions to the FitzHugh-Nagumo equation,” Applied Mathematics and Computation, vol. 180, no. 2, pp. 524–528, 2006.

H. Triki and A.-M. Wazwaz, “On soliton solutions for the Fitzhugh-Nagumo equation with time-dependent coefficients,” Applied Mathematical Modelling, vol. 37, no. 6, pp. 3821–3828, 2013.

R. Jiwari, R. K. Gupta, and V. Kumar, “Polynomial differential quadrature method for numerical solutions of the generalized Fitzhugh-Nagumo equation with time-dependent coefficients,” Ain Shams Engineering Journal, vol. 5, no. 4, pp. 1343–1350, 2014.

Q. Zheng and J. Shen, “Pattern formation in the FitzHugh—Nagumo model,” Computers & Mathematics with Applications, vol. 70, no. 5, pp. 1082–1097, 2015.

R. A. Van Gorder, “Gaussian waves in the Fitzhugh-Nagumo equation demonstrate one role of the auxiliary function H(x, t) in the homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 3, pp. 1233–1240, 2012.

F. Liu, P. Zhuang, I. Turner, V. Anh, and K. Burrage, “A semi-alternating direction method for a 2-D fractional FitzHugh–Nagumo monodomain model on an approximate irregular domain,” Journal of Computational Physics, vol. 293, pp. 252–263, 2015.

A. H. Abbasian, H. Fallah, and M. R. Razvan, “Symmetric bursting behaviors in the generalized FitzHugh-Nagumo model,” Biological Cybernetics, vol. 107, no. 4, pp. 465–476, 2013.

S. S. Ray and S. Sahooz, “New exact solutions of fractional zakharov-kuznetsov and modified zakharov-kuznetsov equations using fractional sub-equation method,” Communications in Theoretical Physics, vol. 63, no. 1, pp. 25–30, 2015.

S. Sahoo and S. S. Ray, “Improved fractional sub-equation method for (3+1)-dimensional generalized fractional KdV-Zakharov-Kuznetsov equations,” Computers & Mathematics with Applications, vol. 70, no. 2, pp. 158–166, 2015.