Операции с матрицой и работа с вводом-выводом в файл

#include <iostream> 
#include <fstream> 
#include <string>
 
using namespace std;
void PrintMatrix(double **, int, int);
int SortBySpir(double **, int, int);
 
int main() {
    char Filename[50];
    std::cout << "Input the name of file: ";
    cin >> Filename;
    ofstream out(Filename, ios_base::trunc);
    ifstream in(Filename);
    if (!in)
    {
        std::cout << "Cannot open file.\n";
        std::system("pause");
        exit(1);
    }
    double** matr;
    int n, m;
    std::cout << "Enter size of matr ";
    in >> n;
    in >> m;
    matr = new double*[n];
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; i < m; j++)
        {
            in >> matr[i][j];
        }
    }
    PrintMatrix(matr, n, m);
    SortBySpir(matr, n, m);
    std::cout << endl << endl << "After sort" << endl << endl;
    PrintMatrix(matr, n, m);
    in.close();
 
    char nameOutputFile[40];
    std::cout << "Input the name of output file: ";
    cin >> nameOutputFile;
    ofstream out(nameOutputFile);
    if (!out)
    {
        std::cout << "Cannot open file.\n";
        std::system("pause");
        exit(1);
    }
    int i, j;
    out << n << " " << m << endl;
    for (i = 0; i < n; i++)
    {
        for (j = 0; j < m; j++)
        {
            out.width(8);
            out << matr[i][j];
        }
        out << endl;
    }
    out.close();
 
    for (int i = 0; i < n; i++)
    {
        delete[] matr[i];
    }
    delete[] matr;
    std::system("pause");
    return 0;
}
 
        void PrintMatrix(double** a, int n, int m){
            std::cout.setf(ios::fixed);
            std::cout.precision(2);
            for (int i = 0; i < n; i++)
            {
                for (int j = 0; j < m; j++)
                {
                    std::cout.width(8);
                    std::cout << a[i][j];
                }
                std::cout << endl;
            }
        }
 
        int SortBySpir(double** a, int n, int m) {
            int z;                      //счётчик 
            int i = n / 2; int j = m / 2;       //устанавливаем итераторы на центральный элемент массива
                                            //цикл обхода по спирали по часовой стрелке 
            do {
                z++;//увеличиваем величину сдвига на 1
                for (int k = 0; k<z; k++)//выводим n элементов со сдвигом вверх
                    std::cout << a[i--][j] << ' ';
                //точка выхода из цикла
                if (z == n)//если сдвиг равен размерности массива
                    break;//выход
                for (int k = 0; k<z; k++)//выводим n элементов со сдвигом вправо
                    std::cout << a[i][j++] << ' ';
                n++;//опять увеличиваем сдвиг
                for (int k = 0; k<z; k++)////выводим n элементов со сдвигом вниз
                    std::cout << a[i++][j] << ' ';
                for (int k = 0; k<z; k++)////выводим n элементов со сдвигом влево
                    std::cout << a[i][j--] << ' ';
 
            } while (z);
            std::system("pause");
            return 0;
        }