The models for the paper

W. Teka, J. Tabak, T. Vo, M. Wechselberger, and R. Bertram (2011) The
Dynamics Underlying Pseudo-Plateau Bursting in a Pituitary Cell Model,
Journal of Mathematical Neuroscience, 1:12, doi:10.1186/2190-8567-1-12

are available at the authors web site:

The 3-variable lactotroph model used in Figures 4, 9, 10, 12, and 13:
lactotroph.ode generates a figure like this:


and the Matlab code for the critical manifold and its 2-D projection
and various curves, used in Figures 2 and 3: manifold.m generates a
figure which when rotated looks like Figure 2:

see the authors web site for these and more.