Разработать программу для решения СЛАУ методом Гаусса.

#include <iostream>
#include <iomanip>
#include <math.h>
#include <stdlib.h>
#include <conio.h>
#include <stdio.h>
#include <iostream.h>
#include <time.h>
#define FOR_k  for (int k = 0; k < dim; ++k)
#define FOR_j  for (int j = 0; j < dim; ++j)
#define FOR_i  for (int i = k + 1; i < dim; ++i)
void swap(double &a, double &b)
{
    double temp = a;
    a = b;
    b = temp;                         
 
}
 
void Gauss(double **A, double *B, int dim)
{
    FOR_k
    {
        if (A[k][k] == 0.0)
        {
           FOR_i
            {
                if (A[i][k] != 0.0)
                {
                    FOR_j
                        swap(A[i][j], A[k][j]);
 
                    swap(B[i], B[k]);
 
                    break;
                }
            }
        }
 
        double div = A[k][k];
 
        FOR_j
            A[k][j] /= div;
 
        B[k] /= div;
 
       FOR_i
        {
            double multi = A[i][k];
 
            for (int j = 0; j < dim; ++j)
                A[i][j] -= multi * A[k][j];
 
            B[i] -= multi * B[k];
        }
    }
 
    for (int k = dim - 1; k > 0; --k)
    {
        for (int i = k - 1; i + 1 > 0; --i)
        {
            double multi = A[i][k];
 
            for (int j = 0; j < dim; ++j)
                A[i][j] -= multi * A[k][j];
 
            B[i] -= multi * B[k];
        }
    }                                                    
}
 
int main()
{
     FILE *k,*k1;
    int n;
    double **A;
    double *B;
      k=fopen("D:\\aaa.txt","a+");
                 if(!k)
                       {
                      printf("Otkritie faila ne vozmojno");
                          system("pause");
                          exit(1);
                       }
    fscanf(k,"%d",&n);
     k1=fopen("D:\\bbb.txt", "a+");
             if(!k1)
                     {
                            printf("Otkritie faila ne vozmojno");
                            system("pause");
                            exit(1);
                          }
    B=new double[n];
          A=new double *[n];
    for (int i = 0; i < n; ++i)
        A[i] = new double [n];
         for(int i=0;i<n;i++)
         {
            for(int j=0;j<n;j++)
            {
            fscanf(k,"%lf",&A[i][j]);
             }
            fscanf(k1,"%lf",&B[i]);
             }
    Gauss(A, B, n);
 
    cout << "X = {";
    for(int i=0;i<n;i++)
    {
    cout<< B[i] <<" ";
    }
     cout<<" }";
    for (int i = 0; i < n ; ++i)
        delete [] A[i];
   getch();
    return 0;
}