Перегрузки операторов ">=" и "+" для матриц

#include <iostream>
#include <iomanip>
#include <stdio.h>
 
 
using namespace std;
class Matrix2d
{
private:
    int rank;       
    int** MAT;   
 
    int n;
    int m;
 
public:
 
    void ReadData();
    double diagonal();
    double sumUp();
    double sumDown();
    double Plus(Matrix2d A1);
 
    Matrix2d(unsigned int rank) {
        MAT = new int* [rank];                
        for (int count = 0; count < rank; count++)
            MAT[count] = new int[rank];
        n = m = rank;
    }
 
    ~Matrix2d()
    {
        for (int count = 0; count < rank; count++)    
            delete MAT[count];
        delete[] MAT;
 
    }
};
 
void Matrix2d::ReadData()
{
    cout << "Enter the matrix" << endl;                      
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < m; ++j) {
            cin >> MAT[i][j];
        }
    }
}
 
 
double Matrix2d::diagonal() {
    double result = 0;
    for (unsigned int i = 0; i < n; i++) {                
        result += MAT[i][i];
    }
    return result;
}
 
double Matrix2d::Plus(Matrix2d A)                    
{
    int i=0;
    int j=0 ;
        for (int i=0; i < n; i++) {
            for (int j=0; j < m; j++)
                MAT[i][j] = MAT[i][j] + A.MAT[i][j];
            
        }
        return MAT[i][j];
}
 
double Matrix2d::sumUp() {
    double result = 0;
    for (unsigned int i = 0; i < n - 1; i++) {             
    for (unsigned int j = i + 1; j < m; j++) {
            result += MAT[i][j];
        }
    }
    return result;
}
 
double Matrix2d::sumDown() {
    double result = 0;
    for (unsigned int i = 0; i < n - 1; i++) {
        for (unsigned int j = i + 1; j < m; j++) {     
            result += MAT[j][i];
        }
    }
    return result;
}
 
 
 
 
int main()
{
    unsigned int rank;    
    cout << "Input rank : ";
    cin >> rank;
    Matrix2d mt(rank);
    Matrix2d mk(rank);
    mt.ReadData();
    mk.ReadData();
 
    if (mt.diagonal() == mk.diagonal())  cout << "\ncomparison:" << mt.diagonal() << " =  " << mk.diagonal() <<"\n";  
    if (mt.diagonal() < mk.diagonal())   cout << "\ncomparison:" << mt.diagonal() << "  < " << mk.diagonal() << "\n";
    if (mt.diagonal() > mk.diagonal())   cout << "\ncomparison:" << mt.diagonal() << " >  " << mk.diagonal() << "\n";
 
    
    cout << "diagonal: " << mt.diagonal()<<" and "<<mk.diagonal();          
    cout << endl << "sumUp: " << mt.sumUp()<<" and "<< mk.sumUp();
    cout << endl << "sumDown: " << mt.sumDown() << " and "<< mk.sumDown();
    cout << endl << "plus: " << mt.Plus(mk) << endl;
 
    system("pause");
}