count=0;
    propagating_species;
    species_ID=propagating_species(1,sp);
    index_species=find([population.species]==species_ID);
    fitnesses_species=[population(index_species).fitness];
%
%     [~,sorted_fitnesses]=sort(fitnesses_species);
%     ranking=zeros(1,size(fitnesses_species,2));
%     ranking(sorted_fitnesses)=1:size(fitnesses_species,2);
%     if size(fitnesses_species,2)>1
%         FitnV=(2-selection.pressure+2*(selection.pressure-1)/(size(fitnesses_species,2)-1)*(ranking-1))';
%     else
%         FitnV=2;
%     end
%
    number_overall=propagating_species(2,sp)-propagating_species(3,sp);
    number_crossover=round(crossover.percentage*number_overall);
    number_mutate=number_overall-number_crossover;
    Nind=size(fitnesses_species,2);
    Nsel=2*number_crossover+number_mutate;



    [sorted_fit,sorted_ind]=sort(fitnesses_species,'descend');

    norm_fit = sorted_fit/sum(sorted_fit);
    cum_fit = cumsum(norm_fit);
    period = 1/Nsel;
    cum_fit = [0 cum_fit(1:length(cum_fit)-1)];
    sel_vec= [period*rand period(:,ones(1,Nsel-1))];

    cum_sel_vec = cumsum(sel_vec);

    child = zeros(1,Nsel);
    for sel=1:length(cum_sel_vec)
            child_ind = max([1:length(cum_fit)].*[cum_fit<cum_sel_vec(sel)]);
            child(sel) = sorted_ind(child_ind);
    end