lowDur = 0;
highDur = 100;
durIterat = 500;
sysName = 'equalspiralband50';

stbTime = 1000;

waitAtEnd = 2000;

count = 0;

temptime = get_param(strcat(sysName, '/I_inj/Syn/end'), 'time');
set_param(strcat(sysName, '/I_inj/Syn/begin'), 'time', num2str(stbTime));
dur = linspace(lowDur, highDur, durIterat);
%dur = logspace(0, 3, durIterat);

for i = 1:length(dur)
    if mod(i,100) == 0
        i/durIterat
    end
    set_param(strcat(sysName, '/I_inj/Syn/end'), 'time', num2str(stbTime + dur(i)));
    sim(sysName, stbTime + dur(i) + waitAtEnd);
    count(i) = size(spikes,1);
end

set_param(strcat(sysName, '/I_inj/Syn/end'), 'time', temptime);

gbarI_Afast     = get_param(strcat(sysName, '/I_A/gbarI_A'), 'value');
gbarI_Aslow     = get_param(strcat(sysName, '/I_Aslow/gbarI_A'), 'value');
%gbarI_T     = get_param(strcat(sysName, '/I_T/P_I_T'), 'value');
%gbarLeakNa  = get_param(strcat(sysName, '/I_Na/gbarLeakNa microS'), 'value');
%I_Ab_tau    = get_param(strcat(sysName, '/I_A/Gates/b/Q10'), 'value');
%gbarI_Na    = get_param(strcat(sysName, '/I_Na/gbarNa'), 'value');
%alpha_nhalf = get_param(strcat(sysName, '/I_Kd/n/alpha_n/Constant1'), 'value');
%gbar_syn    = get_param(strcat(sysName, '/I_inj/Syn/begin'), 'after');

plot(dur, count, '-k')
xlabel('Duration (ms)')
ylabel('Spike Number')
%title(strcat(sysName, ', I_T \tau \times 0.8'))
%title(strcat('bandpass REALLY SHORT (ms), gbarI_{Afast}=', gbarI_Afast))
title(strcat('Altered I_Th_{\tau}, gbarI_{Afast}=', gbarI_Afast))