Метод Гаусса со Stringrid`om

float **matrix;
   int i, j;
  n =StringGrid1->RowCount;
  m =StringGrid1->ColCount;
  float  tmp, *xx;
    xx = new float[m];
    int k;
    float *b;
    // init mass
   b=new float[m];
       //c=new float[m];
            matrix=new float*[n];
        for(int i=0;i<n;i++) matrix[i]=new float[m];
            // return **matr;
 
 // convert mass
 for(int g=0;g<n;g++)
                 {
                 for(int j=0;j<m;j++) {
                 if(StringGrid1->Cells[j][g]=="") StringGrid1->Cells[j][g]=0;
                 matrix[g][j]=StrToFloat(StringGrid1->Cells[j][g]);
                 if(StringGrid2->Cells[0][g]==""  ) // || gr->Cells[0][g]==""
                         {
                         StringGrid2->Cells[0][g]=0; //gr->Cells[0][g]=0;
                         }
 
                 }
                  b[g]=StrToFloat(StringGrid2->Cells[0][g]);
                 }
 
//Метод Гаусса
//Прямой ход, приведение к верхнетреугольному виду
 
 
    for (i=0; i<n ; i++)
     {
       tmp=matrix[i][i];
         for (j=n;j>=i;j--) matrix[i][j]/=tmp;
             b[j]/=tmp;
 
           for (j=i+1;j<m;j++)
          {
             tmp=matrix[j][i];
               for (k=n;k>=i;k--) matrix[j][k]-=tmp*matrix[i][k];
               b[i]-=tmp*b[k];
 
          }
      }
  /*обратный ход*/
    xx[n-1] = matrix[n-1][n];
   //   xx[n-1] = b[n-1];
     for (i=n-2; i>=0; i--)
       {
           xx[i] = b[i];
           for (j=i+1;j<m;j++) xx[i]-=matrix[i][j]*xx[j];
       }
 
//Выводим решения
  for(int g=0;g<n+1;g++)
                 {
                 StringGrid2->Cells[0][g]=FloatToStr(b[g]);
                 StringGrid3->Cells[g][0]=FloatToStr(xx[g]);
 
                 for(int j=0;j<m+1;j++)
                  StringGrid1->Cells[j][g]=FloatToStr(matrix[g][j]);
 
                 }
 
    delete[] matrix;
    delete[] xx;

void __fastcall TForm1::Button3Click(TObject *Sender)
{
  double **matrix;
  int i, j;
  n =StringGrid1->RowCount;
  m =StringGrid1->ColCount;
 double  tmp, *xx;
    xx = new double[m];
    int k;
    double *b;
    // init mass
   b=new double[m];
            matrix=new double*[n];
        for(int i=0;i<n;i++) matrix[i]=new double[m];
 
 
 // convert mass
 for(int g=0;g<n;g++)
                 {
                 for(int j=0;j<m;j++) {
                 if(StringGrid1->Cells[j][g]=="") StringGrid1->Cells[j][g]=0;
                 matrix[g][j]=StrToFloat(StringGrid1->Cells[j][g]);
                 if(StringGrid2->Cells[0][g]=="" ) // || gr->Cells[0][g]==""
                         {
                         StringGrid2->Cells[0][g]=0; //gr->Cells[0][g]=0;
                         }
 
                 }
                  b[g]=StrToFloat(StringGrid2->Cells[0][g]);
                 }
 
        //Метод Гаусса
//Прямой ход, приведение к верхнетреугольному виду
 
 
    for (i=0; i<n ; i++)
     {
       tmp=matrix[i][i];
         for (j=n;j>=i;j--) matrix[i][j]/=tmp;
             b[i]/=tmp;
 
           for (j=i+1;j<m;j++)
          {
             tmp=matrix[j][i];
              b[j]-=tmp*b[i];
               for (k=n;k>=i;k--) matrix[j][k]-=tmp*matrix[i][k];
 
 
          }
      }
     double s = 0;
      xx[m - 1] = b[m - 1]/matrix[n - 1][m - 1];
        for(int i = n - 2, j; 0 <= i; i--)
        {
                s = 0;
                for(j = i + 1; j < n; j++)
                {
                        s = s+matrix[i][j]*xx[j];
                }
                xx[i] = 1.0*(b[i] - s)/matrix[i][i];
                }
                //Выводим решения
  for(int g=0;g<n;g++)
                 {
                 StringGrid2->Cells[0][g]=b[g];
                 StringGrid3->Cells[g][0]=xx[g];
 
                 for(int j=0;j<m;j++)
                  StringGrid1->Cells[j][g]=matrix[g][j];
 
                 }
 
   delete[] matrix;
   delete[] xx;
   delete[] b;
}