# here is a more general way to generate 
# spikes 
# it takes alot longer however
# poincare map  period s=.5  records ISIs
#

markov z 2
{0} {r}
{big} {0}
# in this example, I will use the refractory rate model
r=rmax*(1-s*s0)
s'=-s/tau
count'=0
# when z jumps up reset to 0 increment count and turn on refractory
# period
global 1 z-.8 {z=0;count=count+1;s=1}
par rmax=.05,tau=10 big=1,s0=1
@ dt=.1,meth=euler,total=10000,bound=1000
done