Pseudospectral method of solution of the Fitzhugh-Nagumo equation

Daniel Olmos, Bernie D. Shizgal

Research output: Contribution to journalArticlepeer-review

52 Scopus citations


We present a study of the convergence of different numerical schemes in the solution of the Fitzhugh-Nagumo equations in the form of two coupled reaction diffusion equations for activator and inhibitor variables. The diffusion coefficient for the inhibitor is taken to be zero. The Fitzhugh-Nagumo equations, have spatial and temporal dynamics in two different scales and the solutions exhibit shock-like waves. The numerical schemes employed are a Chebyshev multidomain method, a finite difference method and the method developed by Barkley [D. Barkley, A model for fast computer simulation of excitable media, Physica D, 49 (1991) 61-70]. We consider two different models for the local dynamics. We present results for plane wave propagation in one dimension and spiral waves for two dimensions. We use an operator splitting method with the Chebyshev multidomain approach in order to reduce the computational time. Zero flux boundary conditions are imposed on the solutions. © 2009 IMACS.
Original languageAmerican English
Pages (from-to)2258-2278
Number of pages21
JournalMathematics and Computers in Simulation
StatePublished - 1 Mar 2009
Externally publishedYes


Dive into the research topics of 'Pseudospectral method of solution of the Fitzhugh-Nagumo equation'. Together they form a unique fingerprint.

Cite this