Написать функции умножения и транспонирования матрицы

#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <ctime>
 
const int N = 10;
 
using namespace std;
 
void product(int a[][N], int b[][N], int c[][N], int n){
    
    for (int i = 0; i < n; ++i)
        for (int j = 0; j < n; ++j){
            
            c[i][j] = 0;
            for (int k = 0; k < n; ++k)
                c[i][j] += a[i][k]*b[k][j];
        }
}
 
int main(){
    
    srand(time(0));
    
    int n;
    do{
        cout << "n = ";
        cin >> n;
    }while(n < 1 || N < n);
    cout << endl;
 
    int a[N][N], at[N][N], p[N][N], buf[N][N];
    for (int i = 0; i < n; ++i)
        for (int j = 0; j < n; ++j)
            a[i][j] = rand()%10;
    
    for (int i = 0; i < n; ++i)
        for (int j = 0; j < n; ++j)
            at[i][j] = a[j][i];
            
    for (int i = 0; i < n; ++i){
        
        for (int j = 0; j < n; ++j)
            cout << setw(2) << a[i][j];
        cout << "   |  ";   
        for (int j = 0; j < n; ++j)
            cout << setw(2) << at[i][j];
        cout << endl;
    }
    cout << endl;
    
    product(at, at, p, n);
    for (int i = 0; i < n; ++i)
        for (int j = 0; j < n; ++j)
            buf[i][j] = p[i][j];
            
    product(buf, a, p, n);
    for (int i = 0; i < n; ++i)
        for (int j = 0; j < n; ++j)
            buf[i][j] = p[i][j];
 
    product(buf, a, p, n);
    
    cout << " Àò*Àò*À*À:" << endl;
    for (int i = 0; i < n; ++i){
        for (int j = 0; j < n; ++j)
            cout << setw(8) << p[i][j];
        cout << endl;
    }
    cout << endl;   
}