###########################################################
#           Fast subsystem for analog (polynomial) model of 
#           bursting in pancreas beta-cells.
#
###########################################################
#            Problem definition
###########################################################
#
#
#       BIFURCATION PARAMETER
p  c=-15
#
#       parameter values defining model
#
p a=0.25,eta=0.75,uhat=1.5,beta=4,ubeta=-0.954,e=0.0025
#
#       function definitions
#
F(u)=a*((u-uhat)^2-eta^2)
G(u,c)=u^3-3*(u+1)+c
#
#      differential equations
#
u'=du
du'=-F(u)*du-G(u,c)-e*(beta*(u-ubeta)-c)
#
#      initial conditions (optional)
#
u(0)=2.9968
du(0)=0.0000
#
@ total=400,xlo=0,xhi=400,ylo=-2,yhi=4,dt=0.05,maxstor=8500
#
d