# integrate and fire BVP problem with adaptation
v'=per*(vrest-v+rm*ie-g*(v-ek))/tau
g'=-per*g/taua
per'=0
b v-vreset
b v'-vthr
b g-(g'+dg)
init v=-65,g=.72,per=14.9
par ie=4,dg=.1,taua=100,tau=10,ek=-85,vthr=-50,vreset=-65,vrest=-65
par rm=10
@ total=1.0001,dt=.02,xhi=1,ylo=-70,yhi=-50
done