v1'=(-(v1-vleak)+i1-g1*s2*(v1-vrev1))/tau
v2'=(-(v2-vleak)+i2-g2*s1*(v2-vrev2))/tau
s1'=-s1/tsyn1
s2'=-s2/tsyn2
global 1 v1-vthr {v1=vreset;s1=s1+1}
global 1 v2-vthr {v2=vreset;s2=s2+1}
init v1=-66,v2=-55
par vleak=-66,tau=15,i1=20,i2=20,vthr=-50,vreset=-75
par vrev1=0,vrev2=0,g1=0,g2=0,tsyn1=5,tsyn2=5
@ total=500,njmp=5
done