Метод Гаусса

void Gauss(int N)//объявил функцию
{
    int i, k;
    double x[N+1]; // объявил матрицу неизвестных величин
    ifstream f7("cosnx.txt");
    for(k=0;k<N;k++) f7>>b[k]; // прочитал из файла матрицу свободных членов
    ifstream f10("akimod.txt");
    for(k=0; k<N+1; k++)
    {
        for(i=0; i<N+1; i++)
        {
            f10>>a[k][i]; 
        }
    } // прочитал из файла матрицу коэффициентов
        ofstream f11("x.txt");// это для вывода результата в файл
 
}

int main()
{
    int N=8;
//предшествующие функции
    ksi(N);
    xy(N);
    cos(N);
    aki(N);
    akimod(N);
        Gauss(N); //вот так
    return 0;
}

float  tmp, xx[m];
    int k;
 
    for (i=0; i<n; i++)
     {
       tmp=matrix[i][i];
         for (j=n;j>=i;j--)
             matrix[i][j]/=tmp;
           for (j=i+1;j<n;j++)
          {
             tmp=matrix[j][i];
               for (k=n;k>=i;k--)
             matrix[j][k]-=tmp*matrix[i][k];
          }
      }
  /*обратный ход*/
    xx[n-1] = matrix[n-1][n];
     for (i=n-2; i>=0; i--)
       {
           xx[i] = matrix[i][n];
           for (j=i+1;j<n;j++) xx[i]-=matrix[i][j]*xx[j];
       }
 
//Вывод
       for (i=0; i<n; i++)
           cout << xx[i] << " ";
       cout << endl;
 
    delete[] matrix;
 
   return 0;
}

matrix[i][j]/=tmp;

tmp*matrix[i][k]

 xx[i]-=matrix[i][j]*xx[j];

int gauss_det(int size, double **arr, double &det)
{
    int i , j, k;
    det = 1;
    double b;
    for (i = 0; i < size; i++)
    {
        for (j = i + 1; j < size; j++)
        {
            if (arr[i][i] == 0)
            {
                if (arr[i][j] == 0)
                    b = 0;
                else
                    return 0;
            }
            else 
                b = arr[j][i] / arr[i][i];
            for (k = i; k < size; k++)
                arr[j][k] = arr[j][k] - arr[i][k] * b;
        }
        det *= arr[i][i];
    }
    return 1;
}