программа перемножает двумерные массивы. Как объявить функцию mul (умножения) в main ?

////////////////////////////class.h //файл заголовка
 
#include < iostream>
using namespace std;
class massiv
{
        int n;
    int m;
    int ** A;
public:
 
    massiv(int x = 0, int y = 0)
    {
        
            n = x;
            m = y;
            A = new int *[n];
            for (int i = 0; i<n; i++)
            {
                A[i] = new int[m];
    }
    ~massiv()
    {
        for (int i = 0; i < n; i++)
            delete[] A[i];
        delete[] A;
        cout << "destruktor made his work "<<"\n"<< endl;
    }
    massiv(const massiv&);
    void VVOD();
    void VIVOD();
    massiv mul(massiv p1);
 
};
 
 ///////////////////////////////////main.cpp// файл срр
 
#include "class.h"
int main()
{
 
    int n, m;
 
    massiv p0;
    cout << "vvedite razmer 1 matrics ";
        cin >> n;
        cin >> m;
        massiv p1(n, m);
        cout << " vvedite dannye1" << endl;
        p1.VVOD();
        cout << "martics 1" << endl;
        p1.VIVOD();
        cout << "\nvvedite razmer 2 matrics ";
        cin >> n;
        cin >> m;
        massiv p2(n, m);
        cout << " vvedite dannye 2" << endl;
        p2.VVOD();
        cout << "martics 2" << endl;
        p2.VIVOD();
        p0 = p1.mul(p2);\\ как я вижу объявление функции умножения(при отладке на этом месте выскакивает ошибка) 
        cout << "peremnozili martics"<<endl;
        p0.VIVOD();
        return 0;
 
}
 
  //////////////////////////////////////functions.cpp // файл срр
 
#include "class.h"
 
massiv::massiv(const massiv& ob)
{
 
    A = new int *[n];
 
    for (int i = 0; i < n; i++)
        A[i] = new int[n];
        for (int i = 0; i < n; i++)
            for (int j = 0; j < m; j++)
            {
                A[i][j] = ob.A[i][j];
            }
 
    cout << "constructor copy work";
}
 
 
void massiv::VVOD()
{
 
    for (int i = 0; i < n; i++)
        for (int j = 0; j < m; j++)
            cin >> A[i][j];
}
void massiv::VIVOD()
{
 
    for (int i = 0; i < n; i++)
    {
        cout << endl;
        for (int j = 0; j < n; j++)
            cout << setw(4) << A[i][j];
 
    }
    cout << endl;
}
 
massiv massiv::mul(massiv p2)//  вот тут происходит умножение
 
{
 
    massiv p0(n, m);
    for (int i = 0; i < n; i++)
        for (int j = 0; j < m; j++)
        {
            A[i][j] = 0;
            for (int k = 0; k <= m; k++)
                p0.A[i][j] += (this->A[i][j] * p2.A[i][j]);
        }
    return p0;
}

////////////////////////////class.h //файл заголовка
 
#include < iostream>
#include <iomanip>
using namespace std;
class massiv
{
    int n;
    int m;
    int ** A;
public:
 
    massiv(int x = 1, int y = 1)
    {
 
        n = x;
        m = y;
        A = new int *[n];
        for (int i = 0; i < n; i++)
        {
            A[i] = new int[m];
            for (int j = 0; j < m; j++)
                A[i][j] = 0;
        }
    }
    ~massiv()
    {
        for (int i = 0; i < n; i++)
            delete[] A[i];
        delete[] A;
        cout << "destruktor made his work " << "\n" << endl;
    }
    massiv(const massiv&);
    massiv& operator=(const massiv& m);
    void VVOD();
    void VIVOD();
    massiv mul(massiv p1);
 
};
 
///////////////////////////////////main.cpp// файл срр
 
 
int main()
{
 
    int n, m;
 
    massiv p0;
    cout << "vvedite razmer 1 matrics ";
    cin >> n;
    cin >> m;
    massiv p1(n, m);
    cout << " vvedite dannye1" << endl;
    p1.VVOD();
    cout << "martics 1" << endl;
    p1.VIVOD();
    cout << "\nvvedite razmer 2 matrics ";
    cin >> n;
    cin >> m;
    massiv p2(n, m);
    cout << " vvedite dannye 2" << endl;
    p2.VVOD();
    cout << "martics 2" << endl;
    p2.VIVOD();
    p0 = p1.mul(p2); // как я вижу объявление функции умножения(при отладке на этом месте выскакивает ошибка)
        cout << "peremnozili martics" << endl;
    p0.VIVOD();
    system("pause");
    return 0;
 
}
 
//////////////////////////////////////functions.cpp // файл срр
massiv& massiv::operator=(const massiv& ob)
{
    if (&ob != this)
    {
        for (int i = 0; i < n; i++)
            delete[] A[i];
        delete[] A; 
        A = new int *[n];
        for (int i = 0; i < n; i++)
            A[i] = new int[n];
        for (int i = 0; i < n; i++)
            for (int j = 0; j < m; j++)
            {
                A[i][j] = ob.A[i][j];
            }
 
        cout << "operator= work\n";
    }
    return *this;
}
 
massiv::massiv(const massiv& ob)
{
    n = ob.n; // этого не было
    m = ob.m;
    A = new int *[n];
 
    for (int i = 0; i < n; i++)
        A[i] = new int[m]; // тут было n
    for (int i = 0; i < n; i++)
        for (int j = 0; j < m; j++)
        {
            A[i][j] = ob.A[i][j];
        }
 
    cout << "constructor copy work\n";
}
 
 
void massiv::VVOD()
{
 
    for (int i = 0; i < n; i++)
        for (int j = 0; j < m; j++)
            cin >> A[i][j];
}
void massiv::VIVOD()
{
 
    for (int i = 0; i < n; i++)
    {
        cout << endl;
        for (int j = 0; j < m; j++)
            cout << setw(4) << A[i][j];
 
    }
    cout << endl;
}
 
massiv massiv::mul(massiv p2)
 
{
    if (m != p2.n)
        return *this;
    massiv p0(n, p2.m);
    for (int i = 0; i < n; i++)
        for (int j = 0; j < p2.m; j++)
        {
            p0.A[i][j] = 0;
            for (int k = 0; k < m; k++)
                p0.A[i][j] += (this->A[i][k] * p2.A[k][j]);
        }
    return p0;
}