This package of code is meant to calculate correlations via Monte Carlo simulation and linear
response theory for a network of EIF neurons with membrane potentials (mV) evolving according
vi(t) = -vi(t) + (vi(t)) + E0 +
2i(t) + fi(t).
= Cm/gL (ms) is the membrane time constant, Cm, gL are the membrane capacitance and conduc-
tance, E0 (mV) the effective rest potential, (mV) is the noise variance, i(t) is an independently-
evolving white noise process, and fi(t) captures synaptic coupling between the neurons. The
function takes the form
(v) = T exp
v - vT
where vT (mV) sets the soft threshold of the EIF model, and T (mV) is a spike shape parameter.
The synaptic coupling is, by default, through alpha-shaped kernels, so that
(Jij yj)(t), where Jij(t) =