Tension spline method for solution of Fitzhugh–Nagumo equation

One of the most widely studied biological systems with excitable behavior is neural communication by nerve cells via electrical signaling. The Fitzhugh–Nagumo equation is a simplification of the Hodgin–Huxley model (Hodgin and Huxley, 1952) [24] for the membrane potential of a nerve axon. In this paper we developed a three time-level implicit method by using tension spline function. The resulting equations are solved by a tri-diagonal solver. We described the mathematical formulation procedure in detail. The stability of the presented method is investigated. Results of numerical experiments verify the theoretical behavior of the orders of convergence.