# a neuron firing at rate r
par r0=.2,r1=.19,omega=.01
r=r0+r1*sin(omega*t) 
dt=-log(1-ran(1))/r
tnew'=tnew+dt
count'=count+1
aux delt=dt
@ bound=10000,meth=discrete,total=2000
done