# the simple Poisson process 
dt=-log(ran(1))/rate
tfire'=tfire+dt
aux isi=dt
par rate=.05
@ bound=200000,total=4000,meth=discrete
done