Numerical Simulation of Fitzhugh-Nagumo Dynamics Using a Finite Difference-Based Method of Lines

This study investigates the numerical solution of the FitzHugh-Nagumo (FHN) equation, a canonical Nonlinear Reaction-Diffusion system widely used in Neuroscience and Biophysics using the Method of Lines (MoL). The MoL approach, known for its efficiency and flexibility, discretizes spatial variables to transform partial differential equations ( PDEs ) into a system of Ordinary Differential Equations (ODEs), which are then integrated in time. A fourth-order five-point central difference scheme is employed to approximate spatial derivatives, and MATLAB is used to implement the method. To validate the Numerical scheme, the Newell-Whitehead equation (a special case of the FHN model) is solved, and the results are benchmarked against exact solutions. The results exhibit excellent accuracy, with errors remaining in the order of 10−7 to 10−4 across varying time steps. Comparative analysis against results from the Galerkin Finite Element Method confirms the superior accuracy and computational efficiency of the MoL approach. These findings affirm the reliability and robustness of the Method of Lines in solving Nonlinear Reaction-Diffusion systems, suggesting its potential for broader application in modeling complex Scientific and Engineering phenomena.

Method of Lines, FitzHugh-Nagumo Equation, Reaction-Diffusion, Non-linear PDEs, Newell-Whitehead Equation

1. Tsega, E. G. (2022). Numerical solution of three-dimensional transient heat conduction equation in cylindrical coordinates. Journal of Applied Mathematics, 2, 1 – 8. [Google Scholar] [Crossref]

2. Deghan, M., and S, Kazem. (2017). Application of finite difference method of lines on the heat equation. Numerical Methods for Partial Differential Equations, 1–35. [Google Scholar] [Crossref]

3. Mazumder, S. (2016). Numerical methods for partial differential equations: Finite difference and finite volume methods. Academic Press. [Google Scholar] [Crossref]

4. Dawson, C. N., Du, Q., and Dupont, T. F. (1991). A finite difference domain decomposition algorithm for numerical solution of the heat equation. Mathematics of Computation, 57(195), 63–71. https://doi.org/10.1090/S0025-5718-1991-1069577-3 [Google Scholar] [Crossref]

5. Dalabaev, U., and Hasanova, D. (2023). Construction of an approximate-analytical solution for boundary value problems of a parabolic equation. Mathematics and Computer Science, 8(2), 39–45. https://doi.org/10.11648/j.mcs.20230802.11 doi: 10.11648 /j.mcs. 20230802.11 [Google Scholar] [Crossref]

6. Kolar-Pozuna, A., Jancic, M., Rot, M., and Koseca, G. (2024). Some observations regarding the RBF-FD approximation accuracy dependence on stencil size. Journal of Computational Science. https://arxiv.org/abs/2404-03793v1 [ math.NA ]. [Google Scholar] [Crossref]

7. Samir, H., William, E. S., and Graham, W. G. (2009). Method of lines. Scholarpedia, 2(3), 2859. [Google Scholar] [Crossref]

8. Zafarullah A. (1970). Application of the Method of Lines to Parabolic Partial Differential Equations with Error Estimates. Journal of the Association for Computing Machinery, 17(3), 1970. [Google Scholar] [Crossref]

9. Sadiku, M. N. O., and Obiozor, C. N. (2000). A simple introduction to the method of lines. International Journal of Electrical Engineering Education, 37, 282–296. https://doi.org/10.117 [Google Scholar] [Crossref]

10. Okafor, F.M, Mary O. Durojaye, J. Kayode Odeyemi, G.B. Ogunware (2025). A Comparative Study of Stencil-Based Method of Lines for Solving the BurgersHuxley Equation. FUDMA Journal of Sciences, 9(7), 1-7 [Google Scholar] [Crossref]

11. White, R. E., and Subramanian, V. R. (2010). Method of lines for parabolic partial differential equations in Computational methods in chemical engineering with Maple (pp. 353–505). Springer Berlin Heidelberg. [Google Scholar] [Crossref]

12. Hamdi, S., Enright, W., Ouellet, Y., and Schiesser, W. (2005). Method of lines solutions of the extended Boussinesq equations. Journal of Computational and Applied Mathematics, 183(1), 327–342. [Google Scholar] [Crossref]

13. Shakeri, F., and Dehghan, M. (2008). The method of lines for the solution of the onedimensional wave equation subject to an integral conservation condition. Computers and Mathematics with Applications, 56, 2175–2188. https://doi.org/10.1016/j.camwa.2008.06.0 [Google Scholar] [Crossref]

14. Hyman, J. M. (1979). A method of lines approach to the numerical solution of conservation laws. Advances in computational methods for partial differential equations (Vol. 3, pp. 13–321). [Google Scholar] [Crossref]

15. Murray, J. D. (2002). Mathematical biology I: An introduction (3rd ed.). Springer. https://doi.org/10.1007/b98868 [Google Scholar] [Crossref]

16. Murray, J. D. (2003). Mathematical biology II: Spatial models and biomedical applications. Springer. [Google Scholar] [Crossref]

17. Nagumo, J., Arimoto, S., and Yoshizawa, S. (1962). An active pulse transmission line simulating nerve axon. Proceedings of the IRE, 50(10), 2061–2070. https://doi.org/10.1109/JR [Google Scholar] [Crossref]

18. FitzHugh, R. (1961). Impulses and Physiological States in Theoretical Models of Nerve Membrane. Biophysical Journal, 1(6), 445-466. [Google Scholar] [Crossref]

19. Hariharan, G., and Kannan, K. (2010). Haar wavelength method for solving FitzHughNagumo equation. International Journal of Computational and Mathematical Sciences, 4(6). [Google Scholar] [Crossref]

20. Ali, H., Kamrujjaman, M., and Islam, M. S. (2020). Numerical computation of FitzHugh-Nagumo equation: A novel Galerkin finite element approach. International Journal of Mathematical Research, 9(1), 20–27. [Google Scholar] [Crossref]

21. Cevikel, A., Bekir, A., Abu Argub, O., and Abukhaled, M. (2022). Solitary wave solutions of FitzHugh-Nagumo-type equations with conformable derivatives. Frontiers in Physics. . Retrieved from https://www.frontiersin.org [Google Scholar] [Crossref]

22. Zhou, H., Liu, X., Zhang, W., and Liao, Y. (2023). Efficient numerical method for the FitzHugh-Nagumo equations with Neumann boundary conditions. SCIREA Journal of Physics, 8(2). https://doi.org/10.1016/j.apnum.2024.02.010https://doi.org/10.1016/j.a [Google Scholar] [Crossref]

23. Meinhardt, H. (1982). Models of biological pattern formation. Academic Press. [Google Scholar] [Crossref]

• Interplay of Students’ Emotional Intelligence and Attitude toward Mathematics on Performance in Grade 10 Algebra
• Fixed Point Theorem in Controlled Metric Spaces
• Usage of Moving Average to Heart Rate, Blood Pressure and Blood Sugar
• Exploring Algebraic Topology and Homotopy Theory: Methods, Empirical Data, and Numerical Examples
• Degree of Approximation of Function in the Generalized Zygmund Class By (E, Q) (N̅, Pn) Means of Fourier Series