www.gusucode.com > GAVPai_Book_MathworksCntrlFileEx_May2019 > GAVPai_Book_MathworksCntrlFileEx_May2019/DE_compute_tau.m

    % computation of tau for DE binomial crossover

function tau = DE_compute_tau(component_size, probab_recombi)
h = randperm(component_size);

% initialize j_star to a random index so that tau remains non empty
j_star = h(1);      
tau=j_star;
for i=1:component_size-1
    if (rand < probab_recombi)
        tau = union (tau, i);
    end
end
end