format long; format compact;
rand('seed',sum(100*clock));

npop = 50; % Populaation koko
nbits = 20; % Kromosomien määrä
initpop = rand(npop, nbits) > 0.5; % Alustetaan populaatio
crossProb = 0.8; % Risteytymisen todennäköisyys
mutProb = 0.01; % Mutaation todennäköisyys
genNr = 10; % Sukupolvien lkm

% Simuloidaan geneettistä algoritmia
[endpop endvalues stats] = gaSim('gaFun', initpop, crossProb, mutProb, genNr);

% Tulostetaan paras arvo, keskimääräinen arvo ja keskihajonta
clf reset; subplot(2,1,1);
ind = 1:length(stats(:,1));
plot(ind,stats(:,1),'-',ind,stats(:,2),'-.');
title('Paras arvo ja keskiarvo'); xlabel('Iteraatio');
lt=['- ';'-.']; nm=['Paras arvo'; 'Keskiarvo ']; legend(lt,nm)
% print statistics -deps % Paperitulostus