Ng=100; % Number of genes
w=randn(1,Ng); % IQ effect of each allele

Np=1000; % Population size
g=rand(Ng,Np)<.5; % Initial population, 50% mixture of all alleles

ff=[]; % Store mean IQ

change=1e6; % Generation to modify fitness.

for generation=1:100

    iq = 100+w*g*3; % Calculate individual IQ

    ff=[ff mean(iq)];

    if (generation<change)
        % "Standard" fitness
        fitness = interp1([0 50 80 110 150 400],[0 0.1 1 1 0.5 0],iq);
    else
        % Fitness breeding for intelligence
        fitness = interp1([0 50 80 110 150 400],[0 0.1 1 2 2 2],iq);
    end
    % For Preston-Campbell fitness, run prestoncampbell.m and use
    %fitness = interp1(Cclass,Pchild,iq);

    %Breed new generation
    g2=g*0;
    for k=1:Np
        k1=selectrand(fitness); % Parent 1
        k2=selectrand(fitness); % Parent 2

        split=ceil(rand*Ng); % Crossover location
        g2(:,k)=g(:,k1);
        g2(split:Ng,k)=g(split:Ng,k2);
    end
    g=g2;

end
plot(ff)