load eve_black_one.dat


best_vaf = 0.0;
store_vaf = zeros(100,1);

core = [4,4,4,4,4,4,4,4,4,4];

proxtm = [eve_black_one,core'];
proxtm
maxiteration = 1000;
addcontol = .01;

for i = 1:100
    
    [find,vaf,targone,targtwo,outpermone,outpermtwo,rowpermone,colpermone, ...
        rowpermtwo,colpermtwo,addconone,addcontwo,coordone,coordtwo,axes] = ...
        ms_biscaltmac(proxtm,randperm(26),randperm(26),3,1,addcontol,maxiteration);
     
    vaf
    i
    
    store_vaf(i) = vaf;
    
    if (vaf > best_vaf)
        
        best_vaf = vaf;
        best_find = find;
        best_outpermone = outpermone;
        best_outpermtwo = outpermtwo;
        best_coordone = coordone;
        best_coordtwo = coordtwo;
        best_targone = targone;
        best_targtwo = targtwo;
        best_addconone = addconone;
        best_addcontwo = addcontwo;
        best_rowpermone = rowpermone;
        best_colpermone = colpermone;
        best_rowpermtwo = rowpermtwo;
        best_colpermtwo = colpermtwo;
        best_axes = axes;
        
    end
    
end

sorted_vafs = sort(store_vaf');

sorted_vafs

        best_vaf
        best_find
        best_outpermone 
        best_outpermtwo
        best_coordone
        best_coordtwo
        best_targone
        best_targtwo
        best_addconone
        best_addcontwo 
        best_rowpermone
        best_colpermone
        best_rowpermtwo
        best_colpermtwo
        best_axes


