Description
The Epileptor Codim 3 model with ultra-slow modulation of bursting classes (Saggio et al., 2017). It is a neural mass model containing a fast subsystem (x, y) built from the unfolding of a degenerate Takens-Bogdanov bifurcation of codimension 3 and a slow subsystem (z) that steers the fast one back and forth along a path in parameter space, producing bursting. Two additional ultra-slow state variables (uA, uB) modulate the offset (A) and onset (B) points of the burster along great arcs on a sphere of radius R, allowing the model to switch between different classes of bursters (Hysteresis-Loop bursters c0, c0', c2s, c3s, c4s, c10s, c11s, c2b, c4b, c8b, c14b, c16b as classified in Saggio et al. 2017 Table 2). The bifurcation parameters (mu1, mu2, nu) at the offset (A) and onset (B) points are obtained by parametrising great arcs between the points Ain/Aend and Bin/Bend on the sphere. This capture uses the default branch of the resting state (N = 1) and the default stabilising modification (modification = True).