Computational Study of Singularly Perturbed Neurodynamical Models via Cubic B-Spline

This work focuses on solving the singularly perturbed generalized Hodgkin-Huxley (HH) problem. The HH equation is numerically solved by a collocation approach using third-degree splines. The forward difference technique is utilized for time discretization, while θ-weighted schemes are employed for space discretization. Solving non-linear models using discretization and quasi-linearization results in a set of linear algebraic equations, which are solved using matrices. Furthermore, Von Neumann’s (VN) stability and Spectral Radius (S.R) reveal that the suggested technique is unconditionally stable. To assess the performance and accuracy of this method, absolute error (AE), L2, and L∞ norms are offered. The results align with the literature. Simulation results show that the proposed strategy produces accurate results.