Abstract and Applied Analysis

Volume 2012 (2012), Article ID 762516, 13 pages

doi:10.1155/2012/762516

Research Article

Numerical Simulation of the FitzHugh-Nagumo Equations

A. A. Soliman1,2

¹Department of Mathematics, Faculty of Education (AL-Arish), Suez Canal University, AL-Arish 45111, Egypt

²Department of Mathematics, Bisha Faculty of Science and Arts, King Khalid University, P.O. Box. 344, Bisha, 61922, Saudi Arabia

Received 2 March 2012; Accepted 1 June 2012

Academic Editor: Irena Rachůnková

Copyright © 2012 A. A. Soliman. 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 variational iteration method and Adomian decomposition method are applied to solve the FitzHugh-Nagumo (FN) equations. The two algorithms are illustrated by studying an initial value problem. The obtained results show that only few terms are required to deduce approximated solutions which are found to be accurate and efficient.

1. Introduction

The pioneering work of Hodgkin and Huxley [1], and subsequent investigations, has established that good mathematical models for the conduction of nerve impulses along an axon can be given. These models take the form of a system of ordinary differential equations, coupled to a diffusion equation. Simpler models, which seem to describe the qualitative behavior, have been proposed by FitzHugh [2] and Nagumo [3]. This paper is devoted to the study of the FitzHugh-Nagumo (FN) system:

where and are positive constants and is nonlinear function. Existence and uniqueness for this system is given in 1978 by Rauch and Smoller [4], in which they showed that small solutions decay to 0 as and large pulses produce a traveling wave. We consider the FN equations in the following form

and the function is given by Mckean [5] such that:

where is the Heaviside step function

The exact solution of this system is given by:

where , is the speed of the traveling wave and are the zeros of the polynomial

A numerical scheme for FN equations [6] by collocation method and the “Hopscotch” finite difference scheme first proposed by Gordon [7], and further developed by Gourlay and McGuire [8, 9]. Other possible schemes which were considered are (i) finite difference schemes [10], (ii) Galerkin-type schemes [11], and (iii) collocation schemes with quadratic and cubic splines [6]. In this paper, we use the variational iteration and Adomian decomposition methods to find the numerical solutions of the FN equations which will be useful in numerical studies. In our numerical study we consider the case and , also

with these parameters now we can use the exact travelling wave solution (1.5) to test the suggested numerical methods.

2. The Formalism

We introduce the main points of each of the two methods, where details can be found in [12–37].

2.1. The Variational Iteration Method (VIM)

The VIM is the general Lagrange method, in which an extremely accurate approximation at some special point can be obtained, but not an analytical solution. To illustrate the basic idea of the VIM we consider the following general partial differential equation:

where and are linear operators of and respectively, and is a nonlinear operator. According to the VIM, we can expressed the following correction functional in -, and -directions, respectively, as follows:

where and are general Lagrange multipliers, which can be identified optimally via the variational theory, and are restricted variations which mean that . By this method, it is required first to determine Lagrange multipliers and that will be identified optimally. The successive approximations , of the solution will be readily obtained upon using the determined Lagrange multipliers and any selective function . Consequently, the solution is given by

The above analysis yields the following theorem.

Theorem 2.1. The VIM solution of the partial differential equation (2.1) can be determined by (2.3) with the iterations (2.2a) or (2.2b).

2.2. Adomian Decomposition Method (ADM)

Applying the inverse operator to both sides of (2.1) and using the initial condition, we get

where the nonlinear operator is the Adomian polynomial determined by

We next decompose the unknown function by a sum of components defined by the following decomposition series

The above analysis yields the following theorem

Theorem 2.2. The ADM solution of the partial differential equation (2.1) can be determined by the series (2.6) with the iterations (2.4).

3. Applications

We solve the FN equations using the two methods VIM and ADM.

3.1. The VIM for the FN Equations

Consider the FN equations in the form

Then the VIM formulae take the forms

where , and . This yields the stationary conditions

Hence, the Lagrange multipliers are

Substituting these values of Lagrange multipliers into the functional correction (3.2) gives the iterations formulae

We start with initial approximations as follows

and then the first iterations are

and so on.

The VIM produces the solutions as follows

where , will be determined in a recursive manner.

3.2. The ADM for the FN Equations

Consider the FN equations in the following form:

where . Operating by on both sides of (3.9), we get

The ADM assumes that the unknown functions and can be expressed by an infinite series in the forms

where can be determined by using the recurrence relations:

where

such that

Then the first iterations are

and so on.

The ADM yields the solutions as

where , will be determined in a recursive manner.

4. A Test Problem for the FN Equations

We discuss the solutions of the FN equations using the two considered VIM and ADM methods.

4.1. The VIM

Solve the FN equations (1.2) using the VIM with finite iterations at time . A comparison between the computed solutions and the exact solutions at different values of are given in Table 1. We note that the VIM solutions converge to the exact solutions specially when is increased. We show in Figure 1 the behavior of the VIM solutions of FN equations at time . If the exact solutions are plotted on Figure 1 we will find that the VIM and exact solutions curves are indistinguishable.

Table 1: Comparison between the exact and approximate (VIM) solutions for the FN equations at time .

−7.561
0.000865955
0.000865955
0.0000831702
0.0000831702

−3.561
0.3
0.3
0.0288134
0.0288134

−0.561
0.662823
0.662823
0.28156
0.28156

1.439
0.130074
0.130074
0.414489
0.414489

3.439
−0.25639
−0.25639
0.382524
0.382524

8.439
−0.215232
−0.215232
0.193024
0.193024

16.439
−0.0616494
−0.0616494
16.439
16.439

22.439
−0.023095
−0.023095
0.0197795
0.0197795

48.439
−0.000325199
−0.000325199
0.000278481
0.000278481

(a)

(b)

(c)

(d)

Figure 1: The approximated solutions for at time .

4.2. The ADM

Consider the same problems and use the ADM with the same initial conditions and use the technique discussed in Section 2. A comparison between the exact solutions and ADM solutions are shown in Table 2 and it seems that the errors are very small. We show in Figure 2 the numerical solutions of the FN equations.

Table 2: Comparison between the exact solutions and approximation solutions (ADM) for FN equations at time = 5.

48.439
−0.000325199
−0.000325199
0.000278481
0.000278481

(a)

(b)

(c)

(d)

Figure 2: The approximation solutions .

The results listed in Table 3 are representing the maximum errors at different times of VIM and ADM which shows that the VIM is better than ADM in the solutions of FN equations.

Table 3: The maximum errors of our suggested methods VIM and ADM.

Time
VIM
ADM

Max. errors for
Max. errors for
Max. errors for
Max. errors for

2.0

4.0

6.0

Now we show a comparison between our schemes and other methods as shown in Table 4.

Table 4: Comparison between VIM, ADM, and other methods by maximum errors.

Method

Finite difference

    C-N

0.189

    Hopscotch [9]

0.0506

Collocation method

    Quadratic [6]

0.138

    Cubic [6]

0.12

VIM

0.000316341

ADM

0.000316341

It is clear that the suggested methods for solving FN equation are the best methods than all other methods. Also all other methods give the solution as a discrete solution but our methods give the solution as a function and .

5. Conclusion

In this paper the solutions for the FN equations using VIM and ADM methods have been generated. All numerical results obtained using few terms of the VIM and ADM show very good agreement with the exact solutions. Comparing our results with those of previous several methods shows that the considered techniques are more reliable, powerful, and promising mathematical tools. We believe that the accuracy of the VIM and ADM recommend it to be much wider applicability and also we find that the VIM more accurate than ADM.

Acknowledgment

This paper has been supported by the research support program, King Khalid University, Saudi Arabia, Grand no. KKU-COM-11-001.

References

A. L. Hodgkin and A. F. Huxley, “A quantitative description of membrane current and its application to conduction and excitation in nerve,” The Journal of Physiology, vol. 117, no. 4, pp. 500–544, 1952.

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

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

J. Rauch and J. Smoller, “Qualitative theory of the FitzHugh-Nagumo equations,” Advances in Mathematics, vol. 27, no. 1, pp. 12–44, 1978.

H. P. McKean, “Nagumo's equation,” Advances in Mathematics, vol. 4, no. 3, pp. 209–223, 1970.

A. K. A. Khalifa, Theory and application of the collocation method with splines for ordinary and partial differential equations [Ph.D. thesis], Heriot-Watt University, 1979.

P. Gordon, “Nonsymmetric difference equations,” Journal of the Society For Industrial and Applied Mathematics, vol. 13, no. 3, pp. 667–673, 1965.

A. R. Gourlay, “Hopscotch: a fast second-order partial differential equation solver,” Journal of the Institute of Mathematics and Its Applications, vol. 6, pp. 375–390, 1970.

A. R. Gourlay, “Some recent methods for the numerical solution of time-dependent partial differential equations,” Proceedings of the Royal Society, vol. 323, pp. 219–235, 1971.

J. Rinzel, “Repetitive nerve impulse propagation: numerical results and methods,” in Nonlinear Diffusion, W. E. Fitzgibbon III and H. F. Walker, Eds., Pitman, New York, NY, USA, 1977.

J. R. Cannon and R. E. Ewing, “Galerkin procedures for systems of parabolic partial differential equations related to the transmission of nerve impulses,” in Nonlinear Diffusion, W. E. Fitzgibbon III and H. F. Walker, Eds., Pitman, New York, NY, USA, 1977.

J. H. He and X. H. Wu, “Exp-function method for nonlinear wave equations,” Chaos, Solitons & Fractals, vol. 30, no. 3, pp. 700–708, 2006.

J. H. He, “Variational approach for nonlinear oscillators,” Chaos, Solitons & Fractals, vol. 34, no. 5, pp. 1430–1439, 2007.

J. H. He, “Variational principles for some nonlinear partial differential equations with variable coefficients,” Chaos, Solitons & Fractals, vol. 19, no. 4, pp. 847–851, 2004.

J. H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.

J. H. He and X.-H. Wu, “Construction of solitary solution and compacton-like solution by variational iteration method,” Chaos, Solitons & Fractals, vol. 29, no. 1, pp. 108–113, 2006.

J. H. He, “Variational iteration method for autonomous ordinary differential systems,” Applied Mathematics and Computation, vol. 114, no. 2-3, pp. 115–123, 2000.

J. H. He, “Variational iteration method—a kind of non-linear analytical technique: aome examples,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 699–708, 1999.

J. H. He, “Approximate analytical solution for seepage flow with fractional derivatives in porous media,” Computer Methods in Applied Mechanics and Engineering, vol. 167, no. 1-2, pp. 57–68, 1998.

Z. M. Odibat and S. Momani, “Application of variational iteration method to nonlinear differential equations of fractional order,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 27–34, 2006.

N. Bildik and A. Konuralp, “The use of variational iteration method, differential transform method and adomian decomposition method for solving different types of nonlinear partial differential equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 65–70, 2006.

E. Yusufoglu, “Variational iteration method for construction of some compact and noncompact structures of Klein-Gordon equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 153–158, 2007.

N. H. Sweilam and M. M. Khader, “Variational iteration method for one dimensional nonlinear thermoelasticity,” Chaos, Solitons and Fractals, vol. 32, no. 1, pp. 145–149, 2007.

H. Tari, D. D. Ganji, and M. Rostamian, “Approximate solutions of K (2,2), KdV and modified KdV equations by variational iteration method, homotopy perturbation method and homotopy analysis method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 203–210, 2007.

M. A. Abdou and A. A. Soliman, “Variational iteration method for solving Burger's and coupled Burger's equations,” Journal of Computational and Applied Mathematics, vol. 181, no. 2, pp. 245–251, 2005.

A. A. Soliman and M. A. Abdou, “Numerical solutions of nonlinear evolution equations using variational iteration method,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 111–120, 2007.

A. A. Soliman, “Numerical simulation of the generalized regularized long wave equation by He's variational iteration method,” Mathematics and Computers in Simulation, vol. 70, no. 2, pp. 119–124, 2005.

M. A. Abdou and A. A. Soliman, “New applications of variational iteration method,” Physica D, vol. 211, no. 1-2, pp. 1–8, 2005.

A.-M. Wazwaz, “A comparison between the variational iteration method and Adomian decomposition method,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 129–136, 2007.

A.-M. Wazwaz, “The variational iteration method for rational solutions for KdV, , Burgers, and cubic Boussinesq equations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 18–23, 2007.

A.-M. Wazwaz and A. Gorguis, “Exact solutions for heat-like and wave-like equations with variable coefficients,” Applied Mathematics and Computation, vol. 149, no. 1, pp. 15–29, 2004.

G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic, Boston, Mass, USA, 1994.

A. A. Soliman, “A numerical simulation and explicit solutions of KdV-Burgers' and Lax's seventh-order KdV equations,” Chaos, Solitons & Fractals, vol. 29, no. 2, pp. 294–302, 2006.

M. T. Darvishi, F. Khani, and A. A. Soliman, “The numerical simulation for stiff systems of ordinary differential equations,” Computers and Mathematics with Applications, vol. 54, no. 7-8, pp. 1055–1063, 2007.

A. A. Soliman, “On the solution of two-dimensional coupled Burgers' equations by variational iteration method,” Chaos, Solitons & Fractals, vol. 40, no. 3, pp. 1146–1155, 2009.

A. A. Soliman and M. A. Abdou, “The decomposition method for solving the coupled modified KdV equations,” Mathematical and Computer Modelling, vol. 47, no. 9-10, pp. 1035–1041, 2008.

J. H. He, “A short remark on fractional variational iteration method,” Physics Letters A, vol. 375, no. 38, pp. 3362–3364, 2011.