Вывод в файл

#include <iostream>
#include <fstream>
#include <cmath>
#include <vector>
#include <iomanip>
#include <stdio.h>
using namespace std;
 
//Вычисляем элементы a(ij)
vector<vector<double>>MATRICA(int n)
{
    vector<vector<double>>D;
    for(int i=1;i<n;i++)
    {
        vector<double>L;
        for (int j=1;j<n;j++)
        {
            L.push_back(1/((double)sin(double(i+j))*(double)sin(double(i+j))));
        }
        D.push_back(L);
    }
    return D;
};
//формируем матрицу R размером n*2n, левая половина которой матрица А
//правая половина - единичная матрица
vector<vector<double>>MATRICA_R(vector<vector<double>>D)
{
    vector<vector<double>>L;
    for(int i=0;i<D.size();i++)
    {
        vector<double>EDINICA;
        for(int j=0;j<D.size();j++)
        {
            //если i=j, то соответственно на диагонали будет 1
            if(i==j)
            {
                EDINICA.push_back(1);
            }
            else 
            {EDINICA.push_back(0);}
        }
        L.push_back(EDINICA);
    }
    for(int i=0;i<D.size();i++)
    {
        for(int j=0;j<D.size();j++)
        {
            D[i].push_back(L[i][j]);
        }
    }
    return D;
};
//метод Гаусса
vector<vector<double>>Gauss(vector<vector<double>>D)
{
    vector<vector<double>>GAUSS(D);
    vector<vector<double>>G2(D);
    double a;
    int n=D.size();
    for(int j=0;j<D.size();j++)
    {
 
        int nmax=j;
        double AMAX=D[j][j];
        for(int i=j;i<D.size();i++)
            //ищем наиболший элемент столбца
        {
            if(abs(D[i][j])>abs(AMAX))
            {
                AMAX=D[i][j];
                nmax=i;
            }
        }
        for(int i=0;i<D.size();i++)
            //Переставляем элементы местами в матрице
        {
            a=D[j][i];
            D[j][i]=D[nmax][i];
            D[nmax][i]=a;
        }
        for(int i=0;i<D.size();i++)
        {
            for(int j=0;j<D.size();j++)
            {
                GAUSS[i][j]=D[i][j];
            }
        }
        for(int i=j+1;i<D.size();i++)
        {
            for(int k=0;k<D.size();k++)
            {
                //конец прямого хода метода Гаусса
                D[i][k]=D[i][k]-GAUSS[j][k]*GAUSS[i][j]/GAUSS[j][j];
            }
        }
    }
    //Вычисляем неизвестные элементы матрицы R
    //соответсвующие j-му столбцу обратной матрице
    vector<vector<double>>L;
    for(int j1=0;j1<n;j1++)
    {
        vector<double>EDINICA;
        double *x=new double[n];
        double k=D[n-1][n+j1]/D[n-1][n-1];
        x[n-1]=k;
        double S=0;
        for(int i=n-2;i>=0;i--)
        {
            S=D[i][n+j1];
            for(int j=n-1;j>i;j--)
            {
                S=S-D[i][j]*x[j];
            }
            x[i]=S/D[i][i];
        }
        for(int i=0;i<n;i++)
        {
            EDINICA.push_back(x[i]);
        }
        L.push_back(EDINICA);
        delete[]x;
    }
    double d;
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            d=L[i][j];
            L[i][j]=L[j][i];
            L[j][i]=d;
        }
    }
    return L;
};
void print(vector<vector<double>>D)
{
    
    for(int i=0;i<D.size();i++)
    {
        for (int j=0;j<D.size();j++)
        {
            cout<<D[i][j]<<" ";
        }
        cout<<endl;
    }
    cout<<endl;
 
}
int main()
{
    int n=9;
        ofstream fout("1.txt",ios::out,ios::binary);
    cout<<"MATRICA A\n";
    vector<vector<double>>D(MATRICA(n));
    print(D);
    vector<vector<double>>L(MATRICA_R(D));
    print(L);
    vector<vector<double>>GAUSS(Gauss(L));
    print(GAUSS);
    fout.close();
    return 0;
}

 #include <fstream>
using namespace std;
int main()
{
 ofstream out("file.txt"); // Открываем файл для записи
   out<<"HelloWorld"<<"\n"; //Записываем первую строчку
   out<<"UraRabotaet!!!"; //Записываем вторую строчку
 out.close(); //Закрываем файл
 return 0;
 }