% generate_stimulus_sequence.m

global stim_time CS US

CS = zeros(1,numCS);
US = zeros(1,1);
stim_time = 1;

make_stimulus(...
    get(ht1stim1,'Value'), ...
    get(ht1dur1,'Value'), ...
    get(ht1isi,'Value'), ...
    get(ht1stim2,'Value'), ...
    get(ht1dur2,'Value'), ...
    get(ht1iti,'Value'), ...
    get(ht1trials,'Value'))

make_stimulus(...
    get(ht2stim1,'Value'), ...
    get(ht2dur1,'Value'), ...
    get(ht2isi,'Value'), ...
    get(ht2stim2,'Value'), ...
    get(ht2dur2,'Value'), ...
    get(ht2iti,'Value'), ...
    get(ht2trials,'Value'))

stim_time = stim_time + 1;
CS(stim_time,:) = 0;
US(stim_time) = 0;

Scs = zeros(stim_time,numCS);
Pcs = zeros(stim_time,numCS);
Tcs = zeros(stim_time,numCS);
Phics = zeros(stim_time,numCS);

V = zeros(stim_time,numCS);
V(1,:) = Vinit;

Amn = zeros(stim_time,1);

Sfi = zeros(stim_time,1);
Afi = zeros(stim_time,1);
Rfi = zeros(stim_time,1);