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

#include <iostream>
#include <math.h>
#include <complex>
 
using namespace std;
struct complex1
{
    float Re;
    float Im;
};
complex1 vvod1()
{
    complex1 temp;
    cout<<"Введите действительную часть числа\n ";
    cin>>temp.Re;
    cout<<"Введите мнимую часть комплексного числа\n ";
    cin>>temp.Im;
    return temp;
}
void vivod(complex1** A, int N, int M)
{
    int i, j = 0;
    complex1** temp;
    temp = new complex1 * [N];//указатель для хранения
    for (i = 0; i < N; ++i)
        temp[i] = new complex1[M];
    for (i = 0; i < N; i++)
        for (j = 0; j < N; j++)
        {
            temp[i][j].Re = A[i][j].Re;
            if (A[i][j].Im >= 0)
                cout<<temp[i][j].Re;
            cout<<" + "<<A[i][j].Im<<" i "<<endl;
            if (A[i][j].Im <= 0)
                cout<<" "<<A[i][j].Im<<" i "<<endl;
        }
}
//функция транспонирования матрицы.
complex1** transpose(complex1** A, int N, int M)
{
    int i, j;
    complex1** temp;
    temp = new complex1 * [N];//указатель для хранения
    for (i = 0; i < N; ++i)
        temp[i] = new complex1[M];
    for (i = 0; i < N; ++i)
        for (j = 0; j < M; ++j)
        {
            temp[i][j].Re = (A[j][i].Re);
            temp[i][j].Im = (A[j][i].Im);
        }
    return temp;
}
//функция вычисления суммы двух матриц.
complex1** sum_m(complex1** A, complex1** B, int N, int M)
{
    int i, j;
    complex1** temp;
    temp = new complex1 * [N];//указатель для хранения
    for (i = 0; i < N; i++)
    {
        temp[i] = new complex1[M];
    }
    for (i = 0; i < N; i++)//Вычисляем сумму двух матриц
        for (j = 0; j < M; j++)
        {
            temp[i][j].Re = A[i][j].Re + B[i][j].Re;
            temp[i][j].Im = A[i][j].Im + B[i][j].Im;
        }
    return temp;//Возвращаем матрицу как//двойной указатель
}
//функция вычисления разности двух матриц.
complex1** minus_m(complex1** A, complex1** B, int N, int M)
{
    int i, j;
    complex1** temp;
    temp = new complex1 * [N];//указатель для хранения
    for (i = 0; i < N; i++)
    {
        temp[i] = new complex1[M];
    }
    for (i = 0; i < N; i++)//Вычисляем сумму двух матриц
        for (j = 0; j < M; j++)
        {
            temp[i][j].Re = A[i][j].Re - B[i][j].Re;
            temp[i][j].Im = A[i][j].Im - B[i][j].Im;
        }
    return temp; //Возвращаем матрицу как //двойной указатель
}
//Функция вычисления квадрата матрицы.
complex1** sqrt(complex1** A, int M, int N)
{
    int i, j, s, chisloI = 0, chisloR = 0;
    complex1** temp;
    temp = new complex1 * [N]; //указатель для хранения
    for (i = 0; i < N; i++)
        temp[i] = new complex1[M];
    for (i = 0; i < N; i++) //Вычисляем квадрата матрицы.
        for (j = 0; j < M; j++)
        {
            chisloI = 0; chisloR = 0;
            for (s = 0; s < M; s++)
            {
                chisloR += A[i][s].Re * A[s][j].Re - A[i][s].Im * A[s][j].Im;
                chisloI += A[i][s].Re * A[s][j].Im + A[i][s].Im * A[s][j].Re;
            }
            temp[i][j].Re = chisloR;
            temp[i][j].Im = chisloI;
        }
    return temp; //Возвращаем матрицу как//двойной указатель
}
//функция умножения атрицы на комплексное число.
complex1** proizchi(complex1** B, int N, int M)
{
    int i, j;
    complex1** temp;
    temp = new complex1 * [N];//указатель для хранения
    for (i = 0; i < N; i++)
        temp[i] = new complex1[M];
    for (i = 0; i < N; i++)
        for (j = 0; j < M; j++)
        {
            temp[i][j].Re = (0 * B[i][j].Re + 8.1 * B[i][j].Im);
            temp[i][j].Im = (0 * B[i][j].Im - 8.1 * B[i][j].Re);
        }
    return temp; //Возвращаем матрицу как//двойной указатель
}
complex1** proichi(complex1** B, int N, int M)
{
    int i, j; //указатель для хранения
    complex1** temp;
    temp = new complex1 * [N]; //результирующей матрицы
    //выделение памяти для
    for (i = 0; i < N; i++) //хранения результирующей
        temp[i] = new complex1[M]; //матрицы
    for (i = 0; i < N; i++) //Вычисляем сумму двух матриц
        for (j = 0; j < M; j++)
        {
            temp[i][j].Re = (0 * B[i][j].Re + 1.2 * B[i][j].Im);
            temp[i][j].Im = (0 * B[i][j].Im - 1.2 * B[i][j].Re);
        }
    return temp; //Возвращаем матрицу как//двойной указатель
}
//Функция создаёт динамическую матрицу вещественных чисел //размерности N на N, в этой же функции //осуществляется и ввод элементов матрицы
complex1** create_m(int N, int M)
{
    setlocale(LC_ALL, "Russian");
    int i, j;
    complex1** temp;
    temp = new complex1 * [N];
    for (i = 0; i < N; i++)
        temp[i] = new complex1[M];
    cout<<"Ввод матрицы\n";
    for (i = 0; i < N; i++)
        for (j = 0; j < N; j++)
            temp[i][j] = vvodl();
    return temp;
}
//Функция осуществляет построчный вывод матрицы A(N,N)
complex1** output_m(complex1** A, int N, int M)
{
    int i, j;
    complex1** temp;
    temp = new complex1 * [N];
    for (i = 0; i < N; i++)
        temp[i] = new complex1[M];
    cout << "Вывод матрицы\n";
    for (i = 0; i < 1; i++)
    {
        for (j = 0; j < 1; j++)
            vivod(temp, N, N);
    }
    return temp;
}
int main()
{
    setlocale(LC_ALL, "Russian");
    int N;
    cout << "Введите размерность массива: ";
    cin >> N;
    cout << "Введите элементы массива: ";
    complex1** create_m(int N, int M);
    complex1** output_m(complex1 * *A, int N, int M);
}