метод гаусса

#include "stdio.h"
#include "conio.h"
#include <iostream>
#include <string.h>
#include <stdio.h> // Описания функций ввода-вывода
#include <math.h> // Описания математических функций
using namespace std;
const int n=3; 
void ReadData()
{ int n;
double A[n][n];
{ fscanf(f,"%d",&n);
printf("Sborka matritsi m- na n-:\n");
for (int i=0;i<n;i++)
{ for(int j=0;j<n;j++)
{ fscanf(f,"%d",&A[i][j]);
printf("%d",A[i][j]);
} printf("\n");
} printf("\n\n");
fclose(f); } }
void Gauss
( double X[n], const double Z[n][n], const double Y[n])
{ ReadData();
double A[n][n]; // матрица коэффицентов Гаусса
double B[n]; // рабочий массив свободных членов
int i,
int j,
int k; // рабочие переменные
for( i = 0; i < n; i++ ) // копирование в рабочую матрицу A
{ for( j = 0; j < n; j++ ) 
A[i][j] = Z[i][j];
B[i] = Y[i]; // копирование свободных членов
X[i] = 0; }
for( k = 0; k < n-1; k++ )
for( i = k+1; i < n; i++ ) // преобразование строк
{ double r = A[i][k] / A[k][k];
for( j = k; j < n; j++ )
A[i][j] -= A[k][j] * r;
B[i] -= B[k] * r;
}
X[n-1] = B[n-1] / A[n-1][n-1];
 
for( i = n-2; i >= 0; i-- )
{ double s = 0;
for( j = i+1; j < n; j++ ) 
s += A[i][j] * X[j];
X[i] = ( B[i] - s ) / A[i][i];
} printf( "\n" );
for( i = 0; i < n; i++ ) 
{ printf( "\n" ); 
for( j = 0; j < n; j++ ) 
printf( "%8.2lf", A[i][j] );
printf( " %8.2lf", B[i] );
} }
void main( void )
{ ReadData();
double X[n]; // корни системы линейных уравнений
short i, j; // рабочие переменные
printf( " matritsa A* matritsa B\n" ); 
for( i = 0; i < n; i++ ) // цикл строк
{ printf( "\n" ); // новая строка на мониторе 
for( j = 0; j < n; j++ )
printf( "%8.2lf", A[i][j] );
printf( " %8.2lf", B[i] ); 
} printf( "\n" ); 
printf( "stypen4atii vid matritsi \n" ); 
Gauss( X, A, B ); // решение методом Гаусса
printf( "\n" );
printf( "\n" );
printf( " zna4enia peremennih x1,x1,x3\n" ); 
printf( "\n X = " );
for( i = 0; i < n; i++ )
 
printf( "%8.2lf", X[i] );
printf( "\n" );
printf( "\n" );
printf( "rang=\n" );
getch(); 
}