# poisson process with refractory 
# period
r=rmax*(1-s0*s) 
markov z 2
{0} {r}
{100} {0}
global 1 z-.5 {s=1;z=0;count=count+1}
s'=-s/tau
count'=0
par rmax=.05,tau=10,s0=0,amp=1
@ meth=euler,dt=.1,total=1000,bound=1000000,njmp=500,trans=1000
#fc(x)=amp*r*(1-exp(-x/tau))*exp(-r*(x+tau*exp(-x/tau)-tau))
done