Iterated Crank–Nicolson Runge–Kutta Methods and Their Application to Wilson–Cowan Equations and Electroencephalography Simulations

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.