Ошибки при умножении динамических матриц

using namespace std;
// создаю функцию которая умножает динамическую матрицу A на B и создает матрицу С, которая =А*Р
long double mnozim(long double **A, long double **P) {
    // нахожу количество элементов матрицы А(матрицы А Р квадратные и одинакового размера)
    int n = sizeof(*A) / sizeof(long double);
    // создаю двухмерный динамический массив - квадратная матрицу С размерами n*n (такими же, как и А Р)
    long double  **C;
    C = new long double *[n];
    for (int i = 0; i < n; i++) {
        C[i] = new long double[n];
    }
    // умножаю матрицу A на матрицу P и присваиваю результату матрицу С
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; ++j) {
            for (int r = 1; r < n; ++r) {
                C[i][j] = A[0][r] * P[r][0];
            }
        }
    }
    return **C; // ввожу в вызывающе окружение матрицу С
}
 
 
int main()
{
 
    cout << "BBedite rozmeri" << endl;
    cout << "n=m=";
    int n;
    cin >> n;
    /*создаю 2 динамических двухмерных массива(2 квадратные матрицы)*/
    long double  **A;
    long double  **P;
    A = new long double *[n];
    P = new long double *[n];
    for (int i = 0; i < n; i++) {
        A[i] = new long double[n];
        P[i] = new long double[n];
    }
 
 
 
    /* заполняю матрицу А*/
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; ++j) {
            A[i][j] = rand() % 10;
        }
    }
    /* заполняю матрицу P*/
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; ++j) {
            P[i][j] = rand() % 10;
        }
    }
    /*создаю матрицу Q в которую будет записыватся результат умнодение A на P (матрица C)*/
    long double  **Q;
    Q = new long double *[n];
    for (int i = 0; i < n; i++) {
        Q[i] = new long double[n];
    }
    // вызываю функцию mnozim
    **Q = mnozim(A, P);
    //вывожу матрицу Q
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; ++j) {
            if (j + 1 < n) {
                cout << Q[i][j] << " ";
            }
            else {
                cout << Q[i][j] << " " << endl;
            }
        }
    }
 
 
    system("PAUSE");
    return EXIT_SUCCESS;
}

#include <iostream>
 
    using namespace std;
 
    void matrixInput(int** a, int n, int m) {
        for (int i = 0; i < n; i++) {
            for (int j = 0; j < m; j++) {
                cin >> a[i][j];
            }
        }
    }
 
int main()
{
    int k, m, n;
    cout << "Enter a matrix A size:\n";
    cout << "k = ";
    cin >> k;
    cout << "m = ";
    cin >> m;
    cout << "Enter a matrix B size:\n";
    cout << "m = " << m << "\n";
    cout << "n = ";
    cin >> n;
    int** a = new int*[k];
    for (int i = 0; i < k; i++) {
        a[i] = new int[m];
    }
    int** b = new int*[m];
    for (int i = 0; i < m; i++) {
        b[i] = new int[n];
    }
    int** c = new int*[k];
    for (int i = 0; i < k; i++) {
        c[i] = new int[n];
    }
    cout << "Enter the 1-st matrix:\n";
    matrixInput(a, k, m);
    cout << "Enter the 2-nd matrix:\n";;
    matrixInput(b, m, n);
    cout << "Output of the program:\n";
    for (int j = 0; j < k; j++) {
        for (int i = 0; i < n; i++) {
            c[j][i] = 0;
            for (int r = 0; r < m; r++) {
                c[j][i] += a[j][r] * b[r][i];
            }
            cout << c[j][i] << " ";
        }
        cout << "\n";
    }
    for (int i = 0; i < k; i++) {
        delete [] a[i];
        delete [] c[i];
    }
    delete [] a;
    delete [] c;
    for (int i = 0; i < m; i++) {
        delete [] b[i];
    }
    delete [] b;
    system("pause");
    return 0;
}