function test_WTA_blank

% test survival probability

blanktimes  = [100:100:1000]; % test blanking durations between 100 and 1000 integration steps in steps of 100

inp = [5:0.2:6]; % input strengths at which blanking is tested

prestime    = 1000; % "presentation duration" of the stimulus --> time in which input is applied between blanks

amp1        = 1;
amp2        = 1;
steps       = 1005000;
noisetype   = 1;
display     = 0;
date        = clock;

randn('seed',sum(date));
rand('seed',sum(date));

for l = 1:length(inp)
    for k = 1:length(blanktimes)
        
        blanktab(:,1)   = [20000:prestime+blanktimes(k):1000000];
        blanktab(:,2)   = [20000+prestime:prestime+blanktimes(k):1000000+prestime];
        offset1         = inp(l);
        offset2         = offset1;
        [domDurr_state1,domDurr_state2,domDurr_p1,domDurr_p2,stateActive,domDurr_RT1,domDurr_RT2,survProb,~,~,~,PN]=MainDoubleWTARivalry(amp1,amp2,offset1,offset2,steps,noisetype,display,blanktab);
        SurvProb(l,k)   = survProb;
        
        save(sprintf('blanking_%d-%d_pres%d.mat',inp(1),inp(end),prestime), 'SurvProb','blanktimes','inp')
        clear blanktab
    end
end
