Alguém sabe como executar essa função

[napoliao]

function [filho_cross]= crossover(pai_1, pai_2, prob_cross)

pontos_crossover = floor(prob_cross*9);

X = [1*ones(1,pontos_crossover) 2*ones(1,9 - pontos_crossover)];

x = x(randperm(9));

indices_1 = X-ones(1,9);

indices_2 = abs(x-2*ones(1,9));

if(rand(1,1)<=0.5)

genes_pai_1 =Pai_1.*indices_1;

genes_pai_2 =Pai_2.*indices_2;

filhos_cross_1 = genes_pai_1 +genes_pai_2;

else

genes_pai_1 =Pai_1.*indices_2;

genes_pai_2 =Pai_2.*indices_1;

filhos_cross_1 = genes_pai_1 +genes_pai_2;

end