##########################################################
#   FULL SYSTEM BURSTING MODEL
#
##########################################################
#
#       parameter values defining model
#
p beta=4,alpha=-0.954
p a=0.25,eta=0.75,uhat=1.5,e=0.0025
p ia=0
#
#       function definitions
#
f(u) =  -a/3*u^3 + a*uhat*u^2 + (1 -a*( uhat^2 - eta^2) )*u
g(u) =  ( 1 - a/3)*u^3 + a*uhat*u^2 - (2 + a*( uhat^2 - eta^2) )*u -3
h(u) = beta*(u-alpha)
#
#
#       Initial Conditions
#
init u=-1.427,w=-1.451,c=1.596
#
#
#      differential equations
#
u' = f(u) - w - c+ia
w' = g(u) - w
c' =  e*(h(u)-c)
#
@ total=400,xlo=0,xhi=400,ylo=-2,yhi=4,dt=0.05,maxstor=20000
#
d