Умножение матриц

int **mat(int n){
    int **matrix;
    matrix=new int *[n];
    for (int i=0; i<2; i++)
        matrix[i]=new int [n];
    return matrix;
    }
 
int **multiplication(int **a, int **b, int n){
    int **c;
    c=mat(n);
    for (int i=0; i<n; i++){
        for (int j=0; j<n; j++){
            c[i]=c[i]+a[j][i]*b[j];
            }       
        }
    for (int i=0; i<n; i++){
        cout<<c[i]<<" ";
        }
        return 0;
    }

#include <iostream>
 
int** mat(unsigned n)
{
  int** matrix = new int*[n];
  for (unsigned i = 0; i < n; ++i)
    matrix[i] = new int[n];
  return matrix;
}
 
int** multiplication(int** a, int** b, unsigned n)
{
  int** c = mat(n);
  for (unsigned i = 0; i < n; ++i)
  {
    for (unsigned j = 0; j < n; ++j)
    {
      c[i][j] = 0;
      for (unsigned k = 0; k < n; ++k)
      {
        c[i][j] += a[i][k] * b[k][j];
      }
    }
  }
  return c;
}
 
void print(int** m, unsigned n)
{
  for (unsigned i = 0; i < n; ++i)
  {
    for (unsigned j = 0; j < n; ++j)
    {
      std::cout << m[i][j] << "  ";
    }
    std::cout << std::endl;
  }
}
 
int main()
{
  int** a = mat(2);
  int** b = mat(2);
 
  a[0][0] = 1;
  a[0][1] = 2;
  a[1][0] = 3;
  a[1][1] = 4;
 
  b[0][0] = 5;
  b[0][1] = 6;
  b[1][0] = 7;
  b[1][1] = 8;
 
  int** c = multiplication(a, b, 2);
 
  std::cout << "a:" << std::endl;
  print(a, 2);
 
  std::cout << "b:" << std::endl;
  print(b, 2);
 
  std::cout << "c:" << std::endl;
  print(c, 2);
 
  return 0;
}