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

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 a(1,1), b(1,1), c(0,0);
a.mat[0][0]=3;
b.mat[0][0]=7;
c=a+b;

matrix::matrix(int rows, int columns)
{
    matrix::rows=rows;
    matrix::columns=columns;
    mat=new double*[rows]; 
    for (int i=0;i<rows;++i)
    mat[i]=new double[columns];
}

#include <iostream>
using namespace std;
 
class matrix
{
public:
    matrix(int rows, int columns)
    {
        
        matrix::rows=rows;
        matrix::columns=columns;
        mat=new double*[rows]; 
        for (int i=0;i<rows;++i)
            mat[i]=new double[columns];
    }
    ~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));
}
 
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] << ' ';
 
    system("pause");
    return 0;
}

matrix::~matrix(void)
{
    for (int i=0; i<rows;++i)
        delete [] mat[i];
    delete [] mat;
}

matrix(const matrix &obj)
{
    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];
}