Ввод матрицы из файла

#include "stdio.h"
#include "stdafx.h"
#include "conio.h"
#include "clocale"
#include "cmath"
#define N 50
using namespace std;
 
void glavelem( int k,long double mas[] [N + 1], int n, int otv[] );
 
int main( void )
{
    
    setlocale(LC_ALL, "Russian");
 
    long double mas[N] [N + 1];
    long double x[N]; //Корни системы
    int otv[N]; //Отвечает за порядок корней
    int i, j, k, n;
 
    //Ввод данных
    do
    {
        printf( "Введите число уравнений СЛАУ: " );
        scanf( "%d", & n );
        if ( N < n )
        printf( "Число уравнений СЛАУ превышает возможное. Повторите ввод\n" );
    }
    while ( N < n );
 
    printf( "Введите систему в следующем порядке:\n" );
    for ( i = 0; i < n; i++ )
    for ( j = 0; j < n + 1; j++ )
        cout <<"Строка "<< i + 1 <<", стобец " << j + 1 << "\n";*/
    for ( i = 0; i < n; i++ )
    for ( j = 0; j < n + 1; j++ )
        scanf( "%lf", & mas[i] [j] );
      
 
    system("cls");
 
    //Вывод введенной системы
    printf( "Расширенная матрица системы:\n" );
    for ( i = 0; i < n; i++ )
        {
            for ( j = 0; j < n + 1; j++ )
            printf( "%7.2f ", mas[i] [j] );
            printf( "\n" );
        }
 
    //Сначала все корни по порядку
    for ( i = 0; i < n + 1; i++ )
    otv[i] = i;
    //Прямой ход метода Гаусса
    for ( k = 0; k < n; k++ )
        { 
 
            for ( j = n; j >= k; j-- )
            mas[k] [j] /= mas[k] [k];
            for ( i = k + 1; i < n; i++ )
            for ( j = n; j >= k; j-- )
            mas[i] [j] -= mas[k] [j] * mas[i] [k];
            //Как только k=n прямой хорд завершается и мы получаем треугольную матрицу A
        }
 
    //Обратный ход, вычисление значений неизвестных начиная с последнего уравнения
    for ( i = 0; i < n; i++ )
    x[i] = mas[i] [n];
    for ( i = n - 2; i >= 0; i-- )
    for ( j = i + 1; j < n; j++ )
    x[i] -= x[j] * mas[i] [j];
 
    //Вывод результата
    printf( "Корни СЛАУ:\n" );
    for ( i = 0; i < n; i++ )
    for ( j = 0; j < n; j++ )
    if ( i == otv[j] )
        { 
            //Расставляем корни по порядку
            printf( "%f\n", x[j] );
            break;
        }
    getch();
  
    return (0);
}
//----------------------------------------------
//Описание  функции
//----------------------------------------------
void glavelem( int k, double mas[] [N + 1], int n, int otv[] )
{
    int i, j, i_max = k, j_max = k;
    double temp;
    //Ищем максимальный по модулю элемент
    for ( i = k; i < n; i++ )
    for ( j = k; j < n; j++ )
    if ( fabs( mas[i_max] [j_max] ) < fabs( mas[i] [j] ) )
        {
            i_max = i;
            j_max = j;
        }
    //Переставляем строки, перемещая максимальный по модулю элемент в крайнее верхнее положение
    for ( j = k; j < n + 1; j++ )
        {
            temp = mas[k] [j];
            mas[k] [j] = mas[i_max] [j];
            mas[i_max] [j] = temp;
        }
 
    getch();
}