We are attempting to run simulations of the C. Elegans neuronal network by building upon work of the OpenWorm team. We are looking into implementing a neuron model that produces plateau potentials rather than action potentials, as suggested by this paper. Thanks to @Comte for pointing me to this paper (see this link for another question related to our project).
Here is my question: What are known-to-work mathematical models of neurons that produce plateau potentials, rather than action potentials? Ideally, the model should be formulated as a system of differential equations (as is for example the Hodgkin-Huxley model), hence easy to implement in NEURON.
(A note on my personal background: I am a student of mathematics, and this is a project for university.)