# refnn.ode
# refractory neural net
u'=(-u+(1-int{heav(1-t)#u})*f(u))*r
par r=3
f(u)=1/(1+exp(-8*(u-.333)))
@ total=50
done