load number.dat
inperm = [8 10 9 7 6 4 5 3 2 1]
[fit diff coord] = linfit(number,inperm)
[monproxpermut vaf diff] = ...
    proxmon(number(inperm,inperm),fit)
sumfitsq = sum(sum(fit.^2));
prevvaf = 2;
while (abs(prevvaf-vaf) >= 1.0e-006)
    
   prevvaf = vaf;
   [fit diff coord] = linfit(monproxpermut,1:10);
   sumnewfitsq = sum(sum(fit.^2));
   fit = sqrt(sumfitsq)*(fit/sqrt(sumnewfitsq));
    
    [monproxpermut vaf diff] = proxmon(number(inperm,inperm), fit);
        
    
end

fit
diff
coord
monproxpermut
vaf