Abstract
The Wilson–Cowan model has been widely applied for the simulation of electroencephalography (EEG) waves associated with neural activities in the brain. The Runge–Kutta (RK) method is commonly used to numerically solve the Wilson–Cowan equations. In this paper, we focus on enhancing the accuracy of the numerical method by proposing a strategy to construct a class of fourth-order RK methods using a generalized iterated Crank–Nicolson procedure, where the RK coefficients depend on a free parameter c2. When c2 is set to 0.5, our method becomes a special case of the classical fourth-order RK method. We apply the proposed methods to solve the Wilson–Cowan equations with two and three neuron populations, modeling EEG epileptic dynamics. Our simulations demonstrate that when c2 is set to 0.4, the proposed RK4-04 method yields smaller errors compared to those obtained using the classical fourth-order RK method. This is particularly visible when the spectral radius of the connection matrix or the excitation-inhibition coupling coefficient is relatively large.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
