Переделать программу, чтобы перемножение матриц было функцией

#include "pch.h" 
#include <iostream> 
#include <math.h> 
using namespace std;
void main()
{
    setlocale(LC_ALL, "RUS");
    int i, j, N,k, n, a[20][20], b[20][20],c[20][20],C[20][20],F;
    cout << "N="; //ввод количества строк 
    cin >> N;
    cout << "Input matrix A \n";
    //цикл по переменной i, в которой перебираем строки матрицы 
    for (i = 0; i < N; i++)
        //цикл по переменной j, в котором перебираем элементы внутри строки 
        for (j = 0; j < N; j++)
            cin >> a[i][j]; //ввод очередного элемента матрицы 
    cout << "matrix A \n";
    for (i = 0; i < N; i++)
    {
        //цикл по переменной i, в котором перебираем строки матрицы
        for (j = 0; j < N; j++)
            cout << a[i][j] << "\t"; //вывод очередного элемента матрицы
        cout << endl;
    }
    cout << "n="; //ввод количества строк 
    cin >> n;
    cout << "Input matrix B \n";
    //цикл по переменной i, в которой перебираем строки матрицы 
    for ( i = 0; i < n; i++)
        //цикл по переменной j, в котором перебираем элементы внутри строки 
        for ( j = 0; j < n; j++)
            cin >> b[i][j]; //ввод очередного элемента матрицы 
    cout << "matrix B \n";
    for (i = 0; i < n; i++)
    {
        //цикл по переменной i, в котором перебираем строки матрицы
        for (j = 0; j < n; j++)
            cout << b[i][j] << "\t"; //вывод очередного элемента матрицы
        cout << endl;
    }
    // умножение матрицы A на матрицу B
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            for (c[i][j] = 0, k = 0; k < n; k++) {
                c[i][j] += a[i][k] * b[k][j];
            }
        }
    }
    cout << "Input matrix A*B= \n";
    for (i = 0; i < n; i++)
    {
        //цикл по переменной i, в котором перебираем строки матрицы
        for (j = 0; j < n; j++)
            cout << c[i][j] << "\t"; //вывод очередного элемента матрицы
        cout << endl;
    }
    // умножение матрицы B на матрицу a
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            for (C[i][j] = 0, k = 0; k < n; k++) {
                C[i][j] += b[i][k] * a[k][j];
            }
        }
    }
    cout << "Input matrix B*A= \n";
    for (i = 0; i < n; i++)
    {
        //цикл по переменной i, в котором перебираем строки матрицы
        for (j = 0; j < n; j++)
            cout << C[i][j] << "\t"; //вывод очередного элемента матрицы
        cout << endl;
    }
    F = 0;
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            if (c[i][j]==C[i][j] )
            {
                
                F=F+1 ;
            }
            
 
        }
    }
    if (F == N*N)
    {
        cout << "Перестановочные";
    }
    else { cout << " не Перестановочные"; }
}