Определить, является ли целая квадратная матрица порядка n симметричной (относительно главной диагонали)

#include <iostream>
 
template<class T, size_t size>
bool is_simmetric(T matrix[size][size])
{
    for (int i = 0; i < size; i++)
        for (int j = 0; j < size; j++)
            if (i - j < 0)
                if (matrix[i][j] != matrix[j][i])
                    return false;
    return true;
}
 
int main()
{
    int matrix[3][3] = {
        { 1, 2, 3 },
        { 2, 5, 4 },
        { 3, 4, 0 }
    };
    int matrix2[3][3] = {
        { 1, 2, 3 },
        { 2, 5, 4 },
        { 4, 4, 0 }
    };
    std::cout << std::boolalpha << is_simmetric<int, 3>(matrix) << std::endl
              << std::boolalpha << is_simmetric<int, 3>(matrix2) << std::endl;
    return 0;
}

//
// 
// 
//
 
#include<iostream>
#include<fstream>
#include<sstream>
#include<cmath>
#include<cstdlib>
#include<string>
#include<list>
#include<vector>
#include<algorithm>
#include<stdexcept>
#include "Matrix.h"
#include "MatrixIO.h"
 
using namespace std;
using namespace Numeric_lib;
 
//------------------------------------------------------------------------------
 
inline void keep_window_open()
{
    cin.clear();
    cout << "Please enter a character to exit\n";
    char ch;
    cin >> ch;
    return;
}
 
//------------------------------------------------------------------------------
 
int main()
    try
{
    setlocale( LC_ALL, "Russian" );
 
    Matrix<double, 2> mtrx(6, 6);
 
    srand(1);
 
    for (Index i = 0; i < 6; ++i)
        for (Index j = 0; j < 6; ++j)
            mtrx(i, j) = 1.0 * rand() / RAND_MAX;
 
    cout << mtrx << '\n';
 
    for (Index i = 0; i < 6; ++i)
        for (Index j = 0; j < 6; ++j)
            if (j >= i)
                mtrx(i, j) = 0.0;
 
    cout << mtrx << '\n';
 
    keep_window_open();
    return 0;
}
catch (exception& e) {
    cerr << "error: " << e.what() << '\n'; 
    keep_window_open();
    return 1;
}
catch (...) {
    cerr << "Oops: unknown exception!\n"; 
    keep_window_open();
    return 2;
}
 
//------------------------------------------------------------------------------

#include <iostream>
#include <iomanip>
 
int main()
{
    const int size = 4;
    double matrix[size][size] = {
        { 1.0, 2.0, 3.4, 5.2 },
        { 2.3, 3.4, 2.4, 5.2 },
        { 3.2, 3.2, 5.4, 3.2 },
        { 2.3, 3.1, 3.4, 3.2 }
    };
    for (int i = 0; i < size; i++)
    {
        for (int j = 0; j < size; j++)
            std::cout << std::setw(5) << matrix[i][j];
        std::cout << std::endl;
    }
    for (int i = 0; i < size; i++)
        for (int j = i; j < size; j++)
            matrix[i][j] = 0.0;
    for (int i = 0; i < size; i++)
    {
        for (int j = 0; j < size; j++)
            std::cout << std::setw(5) << matrix[i][j];
        std::cout << std::endl;
    }
}