Ошибка при работе с массивом объектов

#include <iostream>
#include <conio.h>
#include <math.h>
 
using namespace std;
 
class matrix {   // клас матрицa
private:
    int *p_m;
    int n;
public:
    matrix();
    matrix(int count_n);
    matrix(const matrix& m);
    matrix(int count_n, bool eye);
    matrix operator+(const matrix my);
    matrix operator-(const matrix my);
    matrix operator*(const matrix my);
    void   operator=(const matrix my);
    int norma();
    void   Out_put();
    ~matrix();
};
//---------------------------------------------------
 
matrix::matrix()    //констурктор намбе зиру (для массива)
{
    const int n = 2; //значение кол-ва столбцов/строк при создании массива!
    cout << "Матрица [размер]: "<< n << endl;
    int my_n = n*n;
    p_m = new int[my_n];
    int in_ch = 0;
    for (int i = 0; i < n; i++){
        for (int j = 0; j < n; j++)
        {
            cout << "[" << i << "," << j << "]: ";
            cin >> p_m[in_ch];
            in_ch++;
        }
    }
};
 
matrix::matrix(int count_n) //констурктор намбе уан (обычный - передается размерность матрицы)
{
    n = count_n;
    int my_n = n*n;
    p_m = new int[my_n];
    int in_ch = 0;
    for (int i = 0; i < n; i++){
        for (int j = 0; j < n; j++)
        {
            cout << "[" << i << "," << j << "]: ";
            cin >> p_m[in_ch];
            in_ch++;
        }
    }
};
 
matrix::matrix(const matrix& m) //конструктор намбе ту (констурктор копирования)
{
    n = m.n;
    int my_n = n*n;
    p_m = new int[my_n];
    for (int i = 0; i<my_n; i++)
        p_m[i] = m.p_m[i];
 
};
 
matrix::matrix(int count_n, bool eye)   //констурктор намбе фри (нулевая матрица заданной размерностью)
{
    n = count_n;
    int my_n = n*n;
    p_m = new int[my_n];
    int in_ch = 0;
    for (int i = 0; i < my_n; i++){
        p_m[i] = 0;
    }
};
 
matrix::~matrix()   //деструктор
{
    delete p_m;
};
 
matrix  matrix::operator*(const matrix my)  //перегрузка умножения матриц
/*
Ах тыж! Вот это алгоритм! Аж трехвложенный цикл!
*/
{
    int key = 0;
    matrix m_ret(n, true);
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            for (int k = 0; k < n; k++)
            {
                m_ret.p_m[key] += p_m[i*n + k] * my.p_m[n*k+j];
            }
            key++;
        }
    }
    return m_ret;
};
 
matrix  matrix::operator+(const matrix my)  //перегрузка сложения матриц
/*
Здесь происходит посимвольное сложение элементом одной матрицы из другой.
*/
{
    int my_n = n*n;
    matrix m_ret(n, true);
    for (int i = 0; i < my_n; i++)
        m_ret.p_m[i] = p_m[i] + my.p_m[i];
    return m_ret;
};
 
matrix  matrix::operator-(const matrix my)  //перегрузка вычитания
/*
Здесь происходит посимвольное вычитание элементом одной матрицы из другой.
*/
{
    int my_n = n*n;
    matrix m_ret(n,true);
    for (int i = 0; i<my_n; i++)
        m_ret.p_m[i] = p_m[i] - my.p_m[i];
    return m_ret;
};
 
void matrix::operator=(const matrix my) //перегрузка присванивания
{
    int my_n = n*n;
    for (int i = 0; i<my_n; i++)
        p_m[i] = my.p_m[i];
};
 
int matrix::norma() //вычисление нормы матрицы
/*
Вычисляется, если не ошибаюсь, вторая норма матрицы (всего их три, одна из них вообще хрен пойми как вычисляется.
Две других вычисляются с разницей лишь в том, что в одном случае складывают столбцы, в другом строки.
В данном случае складываются строки по модулю, среди всех строк находится максимальная сумма. Она и является нормой матрицы
*/
{
    int max = 0, ch = 0;;
    for (int i = 0; i < n; i++)
    {
        int temp = 0;
        for (int j = 0; j < n; j++)
        {
            temp += abs(p_m[ch]);
            ch++;
        }
        if (max < temp) max = temp;
    }
    return max;
};
 
void matrix::Out_put()  //вывод на экран
{
    int out_ch = 0;
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            cout.width(4);
            cout << p_m[out_ch];
            out_ch++;
        }
        cout << "\n";
    }
};
 
const int N = 2;    //кол-во матриц
 
void inSquare(matrix* arr, int num){
    arr[num] = arr[num] * arr[num];
}
//---------------------------------------------------
 
int main()
{
    system("chcp 1251");
    system("cls");
 
 
    matrix *m = new matrix[2];
    inSquare(m, 1);
    m[0].Out_put();
 
/*
    int size;
    cout << "Матрица [размер]: ";
    cin >> size;
 
    matrix m1(size), m2(size), m3(size, true);
    cout << "\nМатрица 1: \n";
    m1.Out_put();
    cout << "\nМатрица 2: \n";
    m2.Out_put();
    cout << "\nМатрица 3: \n";
    m3.Out_put();
 
    cout << endl;
    m3 = m1 + m2;
    cout << "\nРезультат сложения 1 и 2 матриц: \n";
    m3.Out_put();
    m3 = m1 - m2;
    cout << "\nРезультат вычитания из 1 второй матрицы: \n";
    m3.Out_put();
    m3 = m1 * m2;
    cout << "\nРезультат перемножения 1 и 2 матриц: \n";
    m3.Out_put();
    cout << "\nНорма 1 матрицы: \n" << m1.norma();
*/
    system("pause");
    return 0;
};

n = 2;