#two coupled integrate-and-fire neurons
dv1/dt=-v1 + I - qs*s21*(v1-Es)
dv2/dt=-v2 + I - qs*s22*(v2-Es)

#for exponential synapses, replace dv/dt lines with these
#and comment out ds21/dt and ds22/dt lines
#dv1/dt=-v1 + I - qs*s11*(v1-Es)
#dv2/dt=-v2 + I - qs*s12*(v2-Es)


#alpha function synapses using auxilliary variables
ds11/dt=-alpha*s11
ds12/dt=-alpha*s12
ds21/dt=alpha*(s11-s21)
ds22/dt=alpha*(s12-s22)

#initial conditions
v1(0)=0
v2(0)=.3
#since don't enter initial conditions for s's, they start at 0

#when voltage crosses 1, reset to zero and kick synapse
global 1 v1-1 {v1=0;s12=s12+alpha}
global 1 v2-1 {v2=0;s11=s11+alpha}

#parameter values
p I=1.2
p alpha=4
#for inhibitory synapses
p qs=0.1,Es=-0.5
#for excitatory synapses
#p qs=0.01,Es=6
d