% iAMPA: synaptic current, fast threshold modulation with excitatory AMPAergic parameters
% parameters
gAMPA = [0.1]
EAMPA = [0]
tauD = [2]
tauR = [0.2]
IC = [0.1]
IC_noise = [0]
 
% fixed variables
netcon = ones(N_pre,N_post)

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

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