Непонятная для меня ошибка

class matrix
{
public:
    matrix(int,int);
    ~matrix(void);
    int columns;//кол столбиков
    int rows;//кол строк
    double **mat;//сама матрица
    //double det(void);
    //matrix inverse(void);//нахождение обратной
    matrix transpon(void);
    void transpon_this(void);
    matrix operator=(const matrix &);//оператор присваивания
    matrix operator+(const matrix &);//оператор сложения
    matrix operator-(const matrix &);//оператор вычитания
    matrix operator*(const matrix &);//оператор умножения матриц
    bool operator==(const matrix &);//оператор сравнения
};

matrix matrix::operator=(const matrix &input)
{
    for (int i=0; i<rows;++i)
        delete [] mat[i];
    delete [] mat;
 
    rows=input.rows;
    columns=input.columns;
    mat=new double*[rows]; 
    for (int i=0;i<rows;++i)
        mat[i]=new double[columns];
 
    for (int i=0;i<rows;++i)
        for(int j=0;j<columns;++j)
            mat[i][j]=input.mat[i][j];
    return *this;
}

matrix matrix::operator+(const matrix & rightmatrix)
{
    if ((rows==rightmatrix.rows)&&(columns==rightmatrix.columns))
    {
        matrix temp(rows,columns);
        for (int i=0;i<rows;++i)
            for (int j=0;j<columns;++j)
                temp.mat[i][j]=mat[i][j]+rightmatrix.mat[i][j];
        return temp;
    } else return(matrix(0,0));
}

matrix matrix::operator=(const matrix &input)
{
    rows = input.rows;
    columns = input.columns;
    mat = input.mat; 
    
    return *this;
}

#include <iostream>
#include <string>
using namespace std;
 
class matrix
{
public:
    matrix(int rows, int columns)
    {
         cout << "K" << endl;
        
        matrix::rows=rows;
        matrix::columns=columns;
        mat=new double*[rows]; 
        for (int i=0;i<rows;++i)
            mat[i]=new double[columns];
    }
    
    matrix(const matrix &obj)
    {
         cout << "cpy" << endl;
        
        rows = obj.rows;
        columns = obj.columns;
        mat = new double*[rows]; 
        for (int i = 0; i < rows; ++i)
            mat[i] = new double[columns];
 
        for (int i = 0; i < rows; ++i)
            for(int j = 0; j < columns; ++j)
                mat[i][j] = obj.mat[i][j];
    }
    
    ~matrix()
    {
         cout << "~" << endl;
        for (int i = 0; i < rows; ++i)
            delete [] mat[i];
        delete [] mat;
    }
    
    int rows;//кол строк
    int columns;//кол столбиков
    
    double **mat;//сама матрица
    //double det(void);
    //matrix inverse(void);//нахождение обратной
    matrix transpon(void);
    void transpon_this(void);
    matrix operator=(const matrix &);//оператор присваивания
    matrix operator+(const matrix &);//оператор сложения
    //matrix operator-(const matrix &);//оператор вычитания
    //matrix operator*(const matrix &);//оператор умножения матриц
    //bool operator==(const matrix &);//оператор сравнения
};
 
matrix matrix::operator=(const matrix &input)
{
    cout << "=" << endl;
    
    for (int i=0; i<rows;++i)
        delete [] mat[i];
    delete [] mat;
 
    rows=input.rows;
    columns=input.columns;               
    mat=new double*[rows]; 
    for (int i=0;i<rows;++i)
        mat[i]=new double[columns];
 
    for (int i=0;i<rows;++i)
        for(int j=0;j<columns;++j)
            mat[i][j]=input.mat[i][j];
    return *this;
}
 
matrix matrix::operator+(const matrix & rightmatrix)
{
     cout << "+" << endl;
    if ((rows==rightmatrix.rows)&&(columns==rightmatrix.columns))
    {
        matrix temp(rows,columns);
        for (int i=0;i<rows;++i)
            for (int j=0;j<columns;++j)
                temp.mat[i][j]=mat[i][j]+rightmatrix.mat[i][j];
        return temp;
    } else return(matrix(0,0));
}
 
int main()
{
    {
        matrix a(1,1), b(1,1), c(0,0);
        a.mat[0][0] = 3;
        b.mat[0][0] = 7;
        c = a + b;
    
        for (int i = 0; i < c.rows; ++i)
            for (int j = 0; j < c.columns; ++j)
                cout << c.mat[i][j] << ' ';
        cout << endl;
    }
    
     
    system("pause");
    return 0;
}