% Sexually produce a child

function  child = SexuallyProduce(P1, P2, Pc, Pm)



% Crossover
child = P1;
indpc = rand(1, length(P1));
idx = find(indpc>Pc);
child(idx) = P2(idx);

% Mutation
indmu = rand(1, length(P1));
idx = find(indmu<Pm);
child(idx) = ~child(idx);