# iafadap.ode
v'=(vrest-v + rm*ie(t)-g*(v-ek))/tau
g'=-g/taua
init v=-65,g=0
global 1 v-vthresh {v=vreset;g=g+dg}
ie(t)=amp*heav(t-ton)*heav(toff-t)
par vrest=-65,vthresh=-50,vreset=-65,tau=15,rm=10
par ek=-85,dg=.1,taua=100
par amp=4,ton=50,toff=200
aux stim=rm*ie(t)+vrest
@ total=500,dt=.1,maxstor=10000
@ xhi=500,yhi=-45,ylo=-85,nplot=2,yp2=stim
done