Возведение в степень n матрицы m x m

#include "stdafx.h"
#include <iostream>
#include <conio.h>

void create(int **a, int m)
{
    for (int i(0); i < m; i++)
    {
        for (int j(0); j < m; j++)
        {
            std::cout << "a[" << i << "]["<<j<<"] = ";
            std::cin >> a[i][j];
        }
    }
}

void show(int **a, int m)
{
    for (int i(0); i < m; i++)
    {
        for (int j(0); j < m; j++)
        {
            std::cout << a[i][j] << " ";
        }
        std::cout << std::endl;
    }
}

int **matrix(int **a, int m, int n)
{
    int **c = new int *[m];
    for (int i(0); i < m; i++)
    {
        c[i] = new int[m];
    }
    for (int i(0); i < m; i++)
    {
        for (int j(0); j < m; j++)
        {
            c[i][j] = 0;
        }
    }
    for (int k(1); k < n + 1; k++)
    {
        for (int i(0); i < m; i++)
        {
            for (int j(0); j < m; j++)
            {
                for (int h(0); h < m; h++)
                {
                    c[i][j] += a[i][h] * a[h][j];
                }
            }
        }
    }
    return c;
}

int main(int argc, char* argv[])
{
    setlocale(0, "");
    std::cout << "Введите размер матрицы: ";
    int m;
    std::cin >> m;
    std::cout << "Введите степень матрицы: ";
    int n;
    std::cin >> n;
    int **a = new int *[m];
    for (int i(0); i < m; i++)
    {
        a[i] = new int[m];
    }
    create(a, m);
    show(a, m);
    std::cout << std::endl;

    int **с = new int *[m];
    for (int i(0); i < m; i++)
    {
        с[i] = new int[m];
    }

    с = matrix(a, m, n);
    show(с, m);

    _getch();
    return 0;
}