while(((((numColsArg/N)*K)+(numColsArg%N))>=resultNumArg) && (((numColsArg/N)*K)+(numColsArg%N))>N){
    circle=0; 
    Form1->ProgressBar1->Min=0;
    Form1->ProgressBar1->Max=numColsArg/N;
    Form1->ProgressBar1->Position=0;
    while(circle<numColsArg/N){
      Form1->ProgressBar1->StepBy(1);
      Application->ProcessMessages();
      rightOne=N-1;
      rightZero=N-1;
      nextOne=0;
      sigmaMin=100000000000000000000000;
      //Обнуляем массив сочетаний
      for(int i=0;i<N;i++)
        if(i<K){
          mass[i]=1;
          matrixGroup[i]=new long double[numRowsArg];
        } else mass[i]=0;
      for(unsigned long int iteration=0; iteration<C;iteration++){
        //Заполняем массив из N элементов
        int MatCount=0;//счётчик
        for(int i=0; i<N; i++)
          if(mass[i]){
            for(int r=0;r<numRowsArg;r++)
              matrixGroup[MatCount][r]=dataSource[i+N*circle][r];
            MatCount++;
          }
        sigma=0;
        for(int l=0; l<left;l++)
          sigma+=sigmaFunc(matrixGroup, funcVal, numRowsArg, numRowsInModel,l,K,false);   
        sigma/=left;
        if(sigma<sigmaMin){
          sigmaMin=sigma;
          for(int q=0; q<N; q++)
            if(mass[q]) massSigmaBest[q]=1;
            else massSigmaBest[q]=0;
        }
        rightZero=rightOne=N-1;
        //Если правый бит не единица, находим ближайшую единицу и сдвигаем вправо на 1
        while(!mass[rightOne]) rightOne--;
        if(rightOne!=N-1){
          mass[rightOne]=0;
          mass[rightOne+1]=1;
          continue;
        }
        //Если правый бит единица
        //Считаем количество единиц справа
        while(mass[rightZero]) rightZero--;
        //Исключаем единичные разряды справа из рассмотрения
        //находим ближайшую единицу и перемещаем на 1 вправо
        if(rightZero!=N-1){
          nextOne=rightZero;
          while(!mass[nextOne]) nextOne--;
          mass[nextOne]=0;
          mass[nextOne+1]=1;
          for(int i=rightZero+1; i<N; i++)
            mass[i]=0;
          for(int i=nextOne+2;i<nextOne+2+(N-1-rightZero);i++)
            mass[i]=1;
        }
      }
      //Сдвигаем значения вправо
      for(int j=0,k=0; j<N; j++)
      if(massSigmaBest[j]){
        for(int i=0; i<numRowsArg;i++)
          dataSource[k+circle*K][i]=dataSource[j+circle*N][i];
        titleHead[k+circle*K][0]=titleHead[j+circle*N][0];
        titleHead[k+circle*K][1]=titleHead[j+circle*N][1];
        k++;
      }
      circle++;
      Sh1.OlePropertyGet("Cells",(numRowsArg*left+10+step),(circle+1)).OlePropertySet("Value",((AnsiString)(double)sigmaMin).c_str());
    }
    for(int j=0;j<(numColsArg%N);j++){
      for(int i=0; i<numRowsArg;i++)
        dataSource[j+circle*K][i]=dataSource[j+circle*N][i];
      titleHead[j+circle*K][0]=titleHead[j+circle*N][0];
      titleHead[j+circle*K][1]=titleHead[j+circle*N][1];
    }
    step++;
    numColsArg=((numColsArg/N)*K)+(numColsArg%N);
  }