A model of neuronal bursting using three coupled first order differential equations

We describe a modification to our recent model of the action potential which introduces two additional equilibrium points. By using stability analysis we show that one of these equilibrium points is a saddle point from which there are two separatrices which divide the phase plane into two regions. In one region all phase paths approach a limit cycle and in the other all phase paths approach a stable equilibrium point. A consequence of this is that a short depolarizing current pulse will change an initially silent model neuron into one that fires repetitively. Addition of a third equation limits this firing to either an isolated burst or a depolarizing afterpotential. When steady depolarizing current was applied to this model it resulted in periodic bursting. The equations, which were initially developed to explain isolated triggered bursts, therefore provide one of the simplest models of the more general phenomenon of oscillatory burst discharge.