A minimal 3-variable model (1984) capturing the full range of neuronal firing: quiescence, tonic spiking, and chaotic bursting. The slow variable z controls burst generation.
ẋ = y - ax³+bx² - z + Iẏ = c - dx² - yż = r(s(x-x₀) - z)