% iGABAa: synaptic current, fast threshold modulation with inhibitory GABAa-ergic parameters
% parameters
gSYN = [0.25]
ESYN = [-80]
tauD = [10]
tauR = [0.4]
IC = [0.1]
IC_noise = [0]

% fixed variables
netcon = ones(N_pre,N_post)

% functions
ISYN(X,s) = (gSYN.*(s*netcon).*(X-ESYN))
 
% ODEs and ICs
s' = -s./tauD + ((1-s)/tauR).*(1+tanh(X_pre/10))
s(0) = IC+IC_noise.*rand(1,N_pre)

monitor functions % always record the synaptic current
 
% linkers
@current += -ISYN(X_post,s)