Найти произведение матриц

#include <iostream>
#include <cstdlib>
int** matrix_create(int N, int M);
void  matrix_free(int** mat, int N);
int** matrix_miltiply(int& N, int& M, int** m1, int N1, int M1, int** m2, int N2, int M2);
void  matrix_rand(int** m, int N, int M, int num);
void  matrix_print(std::ostream& _out, int** m, int N, int M);
 
 
int main(void){
    int   N1 = 5, M1 = 3;
    int** m1 = matrix_create(N1, M1);
    if(m1 == NULL)
        return 1;
 
    int   N2 = M1, M2 = 8;
    int** m2 = matrix_create(N2, M2);
    if(m2 == NULL)
        return 1;
 
    matrix_rand(m1, N1, M1, 5);
    matrix_rand(m2, N2, M2, 7);
 
    int N3, M3;
    int** m3 = matrix_miltiply(N3, M3, m1, N1, M1, m2, N2, M2);
    if(m3 == NULL)
        return 1;
 
    matrix_print(std::cout, m1, N1, M1);
    matrix_print(std::cout, m2, N2, M2);
    matrix_print(std::cout, m3, N3, M3);
 
    matrix_free(m1, N1);
    matrix_free(m2, N2);
    matrix_free(m3, N3);
    return 0;
}
 
//перемножение матриц
int** matrix_miltiply(int& N, int& M, int** m1, int N1, int M1, int** m2, int N2, int M2){
    if(M1 != N2)
        return NULL;
 
    int** m3 = matrix_create(N1, M2);
    if(m3 == NULL)
        return NULL;
 
    for(int i = 0; i < N1; ++i){
        for(int j = 0; j < M2; ++j){
            m3[i][j] = 0;
            for(int k = 0; k < N2; ++k)
                m3[i][j] += m1[i][k] * m2[k][j];
        }
    }
    N = N1, M = M2;
    return m3;
}
 
//выделение памяти для матрицы
int** matrix_create(int N, int M){
    int** mat = new (std::nothrow) int*[N];
    if(mat == NULL)
        return NULL;
 
    for(int i = 0; i < N; ++i){
        mat[i] = new (std::nothrow) int[M];
        if(mat[i] == NULL){
            matrix_free(mat, i);
            mat = NULL;
            break;
        }
    }
    return mat;
}
 
//удаление памяти
void matrix_free(int** mat, int N){
    if(mat != NULL){
        for(int i = 0; i < N; ++i)
            delete[] mat[i];
        delete[] mat;
    }
}
 
//заполнение
void matrix_rand(int** m, int N, int M, int num){
    for(int i= 0 ; i < N; ++i){
        for(int j = 0; j < M; ++j)
            m[i][j] = rand() % num;
    }
}
 
//печать
void matrix_print(std::ostream& _out, int** m, int N, int M){
    for(int i= 0 ; i < N; ++i){
        for(int j = 0; j < M; ++j)
            _out << m[i][j] << ' ';
        _out << std::endl;
    }
    _out << std::endl;  
}