Произведение матриц

#include <iostream.h>
#include <conio.h>
#include "stdio.h"
#include <fstream.h>
#include <algorithm>
#include <ctime>
#pragma hdrstop
#pragma argsused
using namespace std;
 
int main(int argc, char* argv[])
{
   //Input Matrix A
   int M, N;
   ifstream ifs("A.txt");
   ifs >> M >> N;
   int ** a = new int * [M];
   //Output Matrix A
   cout << "Matrix A:\n";
   for (int r = 0; r < M; r++, cout << endl)
      {
         a[r] = new int [N];
         for (int c = 0; c < N; c++)
            {
               ifs >> a[r][c];
               cout << a[r][c] << " ";
            }
      }
   cout << endl;
   //Input Matrix B
   int m, n;
   ifstream ifs1("B.txt");
   ifs1 >> m >> n;
   int ** a1 = new int * [m];
   //Output Matrix B
   cout<<"Matrix B:\n";
   for (int r1 = 0; r1 < m; r1++, cout << endl)
      {
         a1[r1] = new int [n];
         for (int c1 = 0; c1 < n; c1++)
            {
               ifs1 >> a1[r1][c1];
               cout << a1[r1][c1] << " ";
            }
      }
   cout << endl;
 
   //Input Matrix C
   int i,j;
   int **c = new int *[i]; //Выдает  ошибку где-то здесь!!!!
   for (i = 0; i < M; i++)
      c[i] = new int [j];
   //Multiplication of Matrices
   for(i = 0; i < M; i++)
      for(j = 0; j < n; j++)
         {
            c[i][j] = 0;
            for(int k = 0; k < n; k++)
               c[i][j] += (a[i][k] * a1[k][j]);
         }
   //Output Matrix C
   cout << "Matrix C: " << endl;
   for(i = 0; i < M; i++)
      {
         for (j = 0; j < n; j++)
         cout << c[i][j] << " ";
         cout << endl;
      }
 
   for (int r = 0; r < M; r++) delete[] a[r];
   for (int r1 = 0; r1 < m; r1++) delete[] a1[r1];
   for (int i = 0; i < M; i++) delete[] c[i];
   delete[] a;
   delete[] a1;
   delete[] c;
   getch();
}