Получить матрицу A(B-E)+C, где Е-единичная матрица порядка n, а элементы матрицы C вычисляются по формуле

#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** C;
    double** E;
    double** R;
 
    int n;
    cout<< "Input matrix size ";
    cin>> n;
    if (n > MIN_SIZE) {
        cout<< "Initialize A" <<endl;        // ввод матрицы А
        A= new double* [n];
            for (int i = 0; i < n; i++) {
             A[i]=new double [n];
            }
        for (int i = 0; i < n; i++) {
            cout << "Input " << i << " row"<< endl;
         for (int j = 0; j < n; j++) {
                       cin>>A[i][j];
            }
        }
 
        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++) {
            cout << "Input " << i << " row"<< endl;
         for (int j = 0; j < n; j++) {
                       cin>>B[i][j];
            }
        }
 
        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;
        }
 
        cout<< "Initialize C" <<endl;          // ввод матрицы С
        C = new double* [n];
            for (int i = 0; i < n; i++) {
             C[i]=new double [n];
            }
        for (int i = 0; i < n; i++) {
         for (int j = 0; j < n; j++) {
                       C[i][j]=1./(i+1)+(j+1);
            }
        }
                                                 // вывод матрицы С
        cout << "Matrix C " << endl<<endl;
 
        for (int i = 0; i < n; i++) {
         for (int j = 0; j < n; j++) {
                cout<<C[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;
                      }
            }
        }
                                      // Вычисление результата
        R = multiplication(A,addition(B,multiplicationByNumber(E,-1,n),n),n);
        R = addition(R,C,n);
                                         //Вывод результата
        cout << "Matrix R " << endl<<endl;
 
        for (int i = 0; i < n; i++) {
         for (int j = 0; j < n; j++) {
                cout<<R[i][j]<<" ";
            }
            cout<<endl;
        }
    } else {
        cout<<"Incorrect input data"<<endl;
    }
    system("pause");
}