Составить описание класса для определения квадратных матриц целых чисел

#include <iostream>
#include <stdlib.h>
#include <stdio.h>
using namespace std;
class QMatrix : public Matrix 
{
public:
      QMatrix();
      QMatrix(int);
      ~QMatrix();
};
QMatrix::QMatrix(int n1)
{
    n=n1;
    m=n1;
    a=new double* [m];
    for (int i=0;i<m;i++)
        a[i] = new double [n];
    for (int i=0;i<m;i++)
        for (int j=0;j<n;j++)
            a[i][j]=0.0;
            }

#include <iostream>
#include <cstdlib> 
#include <ctime> 
#include <iomanip>
 
class QMatrix 
{ 
private: 
    int N; 
    int M; 
    double** matrix; 
public:  
    QMatrix(); // êîíñòðóêòîð ïî óìîë÷àíèþ
    QMatrix(int n):Matrix(n,n) // êîíñòðóêòîð ñ ïàðàìåòðàìè
    QMatrix(const Matrix& ob) // êîíñòðóêòîð êîïèè
    ~QMatrix() // äåñòðóêòîð 
    void Enter() 
    void Show() 
    QMatrix operator + (const Matrix&) 
    QMatrix operator = (const Matrix&) 
    QMatrix operator - (const Matrix& ) 
    QMatrix operator * (const Matrix& )  
}  
  
QMatrix::QMatrix() // ðåàëèçàöèÿ êîíñòðóêòîðà ïî óìîë÷àíèþ
{ 
      N = 0; M = 0; 
        matrix = new double*[N]; 
        for ( int i = 0; i < N; ++i) 
            matrix[i] = new double[M]; 
} 
 
QMatrix:: QMatrix(int n1) // ðåàëèçàöèÿ êîíñòðóêòîðà ñ ïàðàìåòðàìè
{
    N = n1;
    M = n1; 
    a=new double* [M];
    for (int i=0;i<M;i++)
        a[i] = new double [N];
    for (int i=0;i<M;i++)
        for (int j=0;j<N;j++)
            a[i][j]=0.0;    
                    
} 
 
QMatrix::QMatrix (const Matrix &ob) // ðåàëèçàöèÿ êîíñòðóêòîðà êîïèè
{   
    N = ob.N; 
    M = ob.M;
    matrix = new double*[N]; 
    for ( int i = 0; i < N; ++i) 
        matrix[i] = new double[M];  
 
    for ( int i = 0; i < N; ++i) 
        for ( int j = 0; j < M; ++j) 
            matrix[i][j] = ob.matrix[i][j];
 
}  
 
QMatrix::~QMatrix() // äåñòðóêòîð
{  
    for ( int i = 0; i < N; ++i) 
        delete[] matrix[i]; 
    delete[] matrix;
}
 
void QMatrix::Enter() // ââîä ìàòðèö
{   
    int i, j;
    for ( i = 0; i < N; ++i) 
        for ( j = 0; j < M; ++j) 
            matrix[i][j] = 0.01 * (rand() % 1001);
 
} 
 
void Matrix::Show() // âûâîä ìàòðèöû
{ 
    int i, j; 
 
    for ( i = 0; i < N; ++i) 
        for ( j = 0; j < M; ++j) 
        { 
            if( j % M == 0) 
                std::cout << std::endl; 
            std::cout << std::setw(10) << matrix[i][j]; 
        } 
        std::cout << std::endl;
} 
 
QMatrix QMatrix :: operator + ( const Matrix& ob) // ñóììà ìàòðèö
{  
    QMatrix temp;                        // ðåàëèçîâàíî ñ ïîìîùüþ âðåìåííîãî îáúåêòà
    temp.N = this->N; 
    temp.M = this->M; 
    int i, j;
    temp.matrix = new double*[temp.N]; 
    for ( i = 0; i < temp.N; ++i) 
        temp.matrix[i] = new double[temp.M]; 
 
    for ( i = 0; i < temp.N; ++i) 
        for ( j = 0; j < temp.M; ++j) 
            temp.matrix[i][j] = this->matrix[i][j] + ob.matrix[i][j]; 
    return temp;
 
} 
  
QMatrix QMatrix::operator - (const Matrix& ob) // ðàçíîñòü ìàòðèö
{ 
    int i, j; 
    for ( i = 0; i < this->N; ++i) 
        for ( j = 0; j < this->M; ++j) 
            this->matrix[i][j] = this->matrix[i][j] - ob.matrix[i][j]; 
    return *this;
} 
 
QMatrix QMatrix::operator * (const Matrix& ob) // óìíîæåíèå ìàòðèö
{   
 
    
    QMatrix temp; 
    temp.N = this->N; 
    temp.M = this->M; 
    int i, j, k;
    temp.matrix = new double*[temp.N]; 
    for ( i = 0; i < temp.N; ++i) 
        temp.matrix[i] = new double[temp.M];
 
    for ( i = 0; i < temp.N; ++i) 
        for ( j = 0; j < temp.M; ++j) 
        {   
            temp.matrix[i][j] = 0;
            for ( k = 0; k < temp.M; ++k)     
               temp.matrix[i][j] = temp.matrix[i][j] + (this->matrix[i][k] * ob.matrix[k][j]); 
        }
    return temp; 
}
 
 
int main() 
{   
    srand((unsigned) time(0));
    QMatrix m1(3, 3); 
    m1.Enter(); 
    m1.Show(); 
    QMatrix m2(3, 3); 
    m2.Enter();
    m2.Show();
    QMatrix m3(3,3); 
    m3 = m1 * m2;
    m3.Show();
    
}