Даны квадратные матрицы A, B порядка n Получить матрицу N=B(A*E), где E- единичная матрица

#include <iostream>
#define  MIN_SIZE 1
using namespace std;
 double** addition(double** A,double** B,int size){
    double** result= new double* [size];
 
        for (int i = 0; i < size; i++) {
             result[i]=new double [size];
            }
        for (int i = 0; i < size; i++) {
         for (int j = 0; j < size; j++) {
                       result[i][j]=A[i][j]+B[i][j];
            }
        }
        return result;
}
 
double** multiplicationByNumber(double** A,double number,int size){
    double** result= new double* [size];
 
        for (int i = 0; i < size; i++) {
             result[i]=new double [size];
            }
        for (int i = 0; i < size; i++) {
         for (int j = 0; j < size; j++) {
                       result[i][j]=A[i][j]*number;
            }
        }
        return result;
}
 
double** multiplication(double** A,double** B,int size){
        double** result= new double* [size];
        double sum;
        for (int i = 0; i < size; i++) {
             result[i]=new double [size];
            }
        for (int i = 0; i < size; i++) {
         for (int j = 0; j < size; j++) {
             sum = 0;
             for (int k = 0; k < size; k++) {
                     sum+=A[i][k]*B[k][j];
                    }
              result[i][j] = sum;
            }
        }
        return result;
}
int main(){
    double** A;
    double** B;
    double** N;
    double** E;
 
    int n;
    cout<< "Input matrix size ";
    cin>> n;
    if (n > MIN_SIZE) {
        A= new double* [n];
            for (int i = 0; i < n; i++) {
             A[i]=new double [n];
            }
        for (int i = 0; i < n; i++) {
            for (int j = 0; j < n; j++) {
                A[i][j]=rand();
            }
        }
        cout << "Matrix A " << endl<<endl;     
        for (int i = 0; i < n; i++) {
         for (int j = 0; j < n; j++) {
                cout<<A[i][j]<<" ";
            }
            cout<<endl;
        }
 
        cout<< "Initialize B" <<endl;       
        B = new double* [n];
            for (int i = 0; i < n; i++) {
             B[i]=new double [n];
            }
        for (int i = 0; i < n; i++) 
        {
         for (int j = 0; j < n; j++) 
         {
            B[i][j]=rand();
         }
        }
 
        cout << "Matrix B" << endl<<endl;     // вывод матрицы В
 
        for (int i = 0; i < n; i++) {
         for (int j = 0; j < n; j++) {
                cout<<B[i][j]<<" ";
            }
            cout<<endl;
        }
 
        E= new double* [n];
            for (int i = 0; i < n; i++) {
             E[i]=new double [n];
            }
        for (int i = 0; i < n; i++) {
         for (int j = 0; j < n; j++) {
                      if (i==j){
                       E[i][j]=1;
                      }else{
                       E[i][j]=0;
                      }
            }
        }
        
        N = multiplication(A,E,n);
        N = multiplication(B,N,n);
 
 
        cout << "Matrix N " << endl<<endl;
 
        for (int i = 0; i < n; i++) {
         for (int j = 0; j < n; j++) {
                cout<<N[i][j]<<" ";
            }
            cout<<endl;
        }
    } else {
        cout<<"Incorrect input data"<<endl;
    }
    system("pause");
}