Ввод / вывод матрицы

#include <iostream.h>
#include<math.h>
#include <conio.h>
#include <iomanip.h>
#include <fstream.h>
#include <cstdio>
 
int main()
{
double a[100][100], b[100][100], c[100][100];
int i,j,l,n,m,m1,n1;
double s(0);
cout<<"1 Matrix size:";
cin>>n;
cin>>m;
cout<<endl;
cout<<"2 Matrix size:";
cin>>n1;
cin>>m1;
cout<<endl;
 
cout<<"a matrix:"<<endl;
for(i=0; i<n; i++)
 {
   for(j=0; j<m; j++)
    {
      cin>>a[i][j];
    }
 }
cout<<"b matrix: "<<endl;
for(i=0; i<n1; i++)
 {
   for(j=0; j<m1; j++)
    {
      cin>>b[i][j];
    }
 }
 
  for(i=0; i<n; i++)
 {
      for(l=0; l<n; l++)
       {
        s = 0;
        for(j=0; j<n; j++)
         {
           s += a[i][j]*b[j][l];
         }
        c[i][l] = s;
       }
 }
 
cout<<"rezult:"<<endl;
for(i=0; i<n; i++)
 {
   for(j=0; j<m1; j++)
    {
      cout<<c[i][j]<<" ";
    }
   cout<<endl;
 }
 
system("pause");
return 0;
}