"Chaotic regimes in a mathematical model of pacemaker activity in the bursting neurons of a snail Helix pomatia, have been investigated.
The model includes a slow-wave generating mechanism, a spike-generating mechanism, an inward Ca current, intracellular Ca ions, [Ca2+]in, their fast buffering and uptake by intracellular Ca stores, and a [Ca2+]in-inhibited Ca current.
Chemosensitive voltage-activated conductance, gB*, responsible for termination of the spike burst, and chemosensitive sodium conductance, gNa*, responsible for the depolarization phase of the slow-wave, were used as control parameters.
...
Time courses of the membrane potential and [Ca2+]in were employed to analyse different regimes in the model.
..."
This is the readme for the model associated with the paper
Komendantov AO, Kononenko NI (1996) Deterministic chaos in
mathematical model of pacemaker activity in bursting neurons of snail,
Helix pomatia. J Theor Biol 183:219-30
The XPP model was contributed by Dr Komendantov. XPP is available
from http://www.math.pitt.edu/~bard/xpp/xpp.html
Example usage under linux:
Start with a command like
xppaut jtb-kk.ode
Select Initalconds -> Go with the mouse. After a short time you
should get this graph

which replicates figure 9 A from the paper.