Как вывести матрицу из файла

#include <iostream>
#include <iomanip>
#include<ctime>
#include<clocale>
#include<cmath>
using std::cout;
using std::endl;
using std::cin;
using std::setw;
 
int main()
{   setlocale(LC_ALL, "Russian");
    int i,j,k;                              
    int nstr,nstb;
    int index;
    double coeff,max;
    double temp;
    
    cout<<"Введите количество строк и столбцов\n";
    cin>>nstr>>nstb;
    
    double **a = new double *[nstr];
    for(int i=0;i<nstr;i++) a[i]= new double[nstb];
    
    double *b=new double [nstr];
    double *x =new double [nstr];
    cout<<"Введите элементы матрицы\n";
    for(i=0;i<nstr;i++)
    {
      for(j=0;j<nstb;j++)
        cin>>a[i][j];
    }
    for(i=0;i<nstr;i++)
    {
        for(j=0;j<nstb;j++){
            cout <<setw(10)<<a[i][j];}
        cout<<endl;
    }
 
    cout<<"Введите вектор, определяющий правую часть с-мы\n";
    for(i=0;i<nstr;i++) 
        cin>>b[i];
    cout<<endl;
    cout<<"Вектор правой части\n";
    for(i=0;i<nstr;i++) 
        cout <<setw(6)<<b[i];
        cout <<endl;
 
    for(j=0; j<nstb-1; ++j)
    {
        //Выбор ведущего элемента столбца
        max=fabs(a[j][j]);
        index=j;
        for(i=j; i<nstr; ++i)
            if(fabs(a[i][j])>max)
            {
                max=fabs(a[i][j]);
                index=i;
            }   
        if(max=0)
            break;//Выход
        
        //Если ведущий элемент - не на гл. диагонали
        //меняем строки местами
        if(index!=j)
        {
            temp=b[j];
            b[j]=b[index];
            b[index]=temp;
            for(i=0; i<nstr; ++i)
            {
                temp=a[i][j];
                a[i][j]=a[index][j];
                a[index][j]=temp;
            }
        }
        //Обнуление элементов ниже гл. диагонали
        for(i=j+1; i<nstr; ++i)
        {
            coeff=-(a[i][j]/a[j][j]);
            cout<<setw(5)<<coeff;
            cout<<endl;
            for(k=j; k<nstb; ++k)
                a[i][k]=a[i][k]+a[j][k]*coeff;
            b[i]=b[i]+b[j]*coeff;
            
        }
    }
 
    for(i=0;i<nstr;i++)
    {
        for(j=0;j<nstb;j++){
            cout <<setw(6)<<a[i][j];}
        cout<<endl;
    }
 
    //Обратная подстановка
    for(j=nstb-1; j>=0; --j)
    {
        x[j]=b[j]/a[j][j];
        if(j!=3)
            for(i=nstr-1; i>j; --i)
                x[j]=x[j]-a[j][i]*x[i]/a[j][j];
    }
    
    for(i=0; i<nstr; ++i, cout <<setw(6)<< endl)
        cout << "x[" << i << "]=" << x[i];
    cout<<endl;
    for(i=0;i<nstr;i++)
    {
        for(j=0;j<nstb;j++){
            cout <<setw(6)<<a[i][j];}
        cout<<endl;
    }
    for(i=0;i<nstr;i++) 
        cout <<setw(6)<<b[i];
        cout<<endl;
    system("pause");
    return 0;
}

#include <fstream>

ifstream fin;
fin.open("filename");