Квадратные матрицы

#include <iostream>
#include <iomanip>
using std::cout;
using std::cin;
using std::fixed;
using std::setprecision;
using std::setw;
 
double **CreateArray (int N)
{   int i,j;
    double **arr = new double* [N];
    for (i = 0; i<N; i++)
        arr[i] = new double [N];
 
    for (i = 0; i<N; i++)
        for (j = 0; j<N; j++)
            arr[i][j] = 0.;
    return arr;
}
 
void DeleteArray (double **arr, int N)
{   for (int i = 0; i<N; i++)
        delete [] arr[i];
    delete [] arr;
}
    
void HandFill (double **arr, int N)
{   cout << "Вводите числа построчно\n";
    int i,j;
    for (i = 0; i<N; i++)
    {   cout <<"Строка "<< i << ", " << N << " элементов\n";
        for (j = 0; j<N; j++)
            cin >> arr[i][j];
    }
}
 
void Multiply (double **R, double **A, double **B, int N)
{   int i, j, k;
    double sum;
    for (i = 0; i < N; i++) 
        for (j = 0; j < N; j++) 
        {   for (k = 0, sum = 0.; k < N; k++)
                sum += A[i][k] * B[k][j];
            R[i][j] = sum;
        }
}
 
void Substract (double **R, double **A, double **B, int N)
{       int i, j;
        for (i = 0; i < N; i++) 
            for (j = 0; j < N; j++) 
                R[i][j] = A[i][j] - B[i][j];
}
 
void PrintArray (double **arr, int N)
{   cout << "Формат вывода: введите общее количество цифр и количество цифр после запятой\n";
    int w, prec;
    cin >> w >> prec;
    int i, j;
        for (i = 0; i < N; i++) 
        {   for (j = 0; j < N; j++) 
                cout << setw (w) << fixed << setprecision (prec) << arr[i][j];
            cout << "\n";
        }
}
 
 
int main () 
{   setlocale (LC_CTYPE, "Russian");
    int N;
 
    cout << "Введите размерность матриц\n\tN = ";
    cin >> N;
 
    double **A = CreateArray (N);
    double **B = CreateArray (N);
    cout << "\nЗаполните матрицу А\n";
    HandFill (A, N);
    cout << "\nЗаполните матрицу B\n";
    HandFill (B, N);
 
    double **C = CreateArray (N);
    double **D = CreateArray (N);
    double **E = CreateArray (N);
 
    Multiply (C, A, B, N);
    Multiply (D, B, A, N);
    Substract (E, C, D, N);
 
    cout << "\nE = AB - BA\n";
    PrintArray (E, N);
 
    DeleteArray (A, N);
    DeleteArray (B, N);
    DeleteArray (C, N);
    DeleteArray (D, N);
    DeleteArray (E, N);
 
    cin.sync(); cin.get();
    return 0;
}