// Matrix multiplication C=AxB
#include <iostream>
#define MAX1 3
#define MAX2 3
#define MAX3 3
#define MAX4 3
using namespace std;
int main(int argc, char *argv[])
{ system("chcp 1253");
 int A[MAX1][MAX2] = {{1,2,3}, {4,5,6}, {7,8,9}};
 int B[MAX3][MAX4] = {{1,0,0}, {0,1,0}, {0,0,1}};
 int i,j,k,C[MAX1][MAX4]={{0}};

 if (MAX2==MAX3) {
   cout<<"Matrix A\n";
   for (i = 0; i < MAX1; i++)     {
      for (j = 0; j < MAX2; j++)
          cout<<A[i][j];
      cout<<"\n";
   }
   cout<<"Matrix B\n";
   for (i = 0; i < MAX3; i++)     {
      for (j = 0; j < MAX4; j++)
          cout<<B[i][j];
      cout<<"\n";
   }

  /* Product */
   for (i = 0; i < MAX1; i++)
       for (j = 0; j < MAX4; j++) {
           C[i][j]=0;
           for (k = 0; k < MAX3; k++)
                 C[i][j]=C[i][j]+A[i][k]*B[k][j];
       }

   cout<<"Matrix C=AxB\n";
   for (i = 0; i < MAX1; i++)     {
      for (j = 0; j < MAX4; j++)
          cout<<C[i][j];
      cout<<"\n";
   }
}
 else   
       cout<<"Asymvates diastaseis pinakwn\n";
 system("PAUSE");
 return 0;
}