Dinamicka alokacija / korisne funkcije za matrice

#include <iostream>

void UnistiMatricu(double **Matrica, int n)
{
    for (int i(0); i < n; i++)
        delete[] Matrica[i];
    delete[] Matrica;
}

void KreirajMatricu(double **&Matrica, int n)
{
    double **p = new double*[n];
    for (int i(0); i < n; i++)
        p[i] = nullptr;
    try
    {
        for (int i(0); i < n; i++)
            p[i] = new double[n];
    }
    catch (std::bad_alloc)
    {
        UnistiMatricu(p, n);
        throw;
    }
    Matrica = p;
}

void UnesiMatricu(double **&Matrica, int n)
{
    for (int i(0); i < n; i++)
        for (int j(0); j < n; j++)
        {
            std::cout << "Unesite element: (" << i + 1 << ", " << j + 1 << "): ";
            std::cin >> Matrica[i][j];
        }
}

void IspisiMatricu(const double *const *Matrica, int n)
{
    for (int i(0); i < n; i++)
    {
        for (int j(0); j < n; j++)
            std::cout << Matrica[i][j] << " ";
        std::cout << std::endl;
    }
}

double **SaberiMatrice(double **Matrica1, double **Matrica2, int n)
{
    double **Matrica;
    KreirajMatricu(Matrica, n);
    for (int i(0); i < n; i++)
        for (int j(0); j < n; j++)
            Matrica[i][j] = Matrica1[i][j] + Matrica2[i][j];
    return Matrica;
}

int main()
{
    return 0;
}