C++实现矩阵运算

#include <iostream.h>
#include <iomanip.h>
void Ju_zhen(int a[4][4])   //输入矩阵
{ for(int i=0;i<4;i++)
 { for(int j=0;j<4;j++)
  cin>>a[i][j];
  cout<<endl;
  } 
}
void Show(int a[4][4])   //输出矩阵结果
{ for(int i=0;i<4;i++)  
 {
  for(int j=0;j<4;j++)
  cout<<a[i][j]<<" ";
  cout<<endl;
  }
}
void Add(int a[4][4],int b[4][4]) //加函数
{ int c[4][4],i,j;
 for(i=0;i<4;i++)
 for (j=0;j<4;j++)
 {  c[i][j]=a[i][j]+b[i][j];}
 Show(c);
}
void Jian(int a[4][4],int b[4][4]) //减函数
{ int c[4][4],i,j;
 for(i=0;i<4;i++)
  for (j=0;j<4;j++)
  { c[i][j]=a[i][j]-b[i][j];}
  Show(c);
}
void Mut(int a[4][4],int b[4][4]) //乘函数
{ int c[4][4],i,j,k,s;
 for(i=0;i<4;i++)
    for (j=0;j<4;j++)
    {  s=0;
       for(k=0;k<4;k++)
   s+=a[i][k]*b[k][j];
    c[i][j]=s;
    }
    Show(c);
}
void main()
{
 int a[4][4];//={ {2,3,4,6},{-1,2,3,4},{2,-3,4,7},{4,0,8,13} };
 int b[4][4];//={ {1,2,-4,7},{1,2,8,5},{-2,-3,4,7},{0,0,6,3} };
 cout<<"输入矩阵A："<<endl;
 Ju_zhen(a);
 cout<<"输入矩阵B："<<endl;
 Ju_zhen(b);
 cout<<"A+B="<<'/n';   //两矩阵相加
 Add(a,b);
 cout<<"A-B="<<'/n';   //两矩阵相减
 Jian(a,b); 
 cout<<"A*B="<<'/n';   //两矩阵相乘
 Mut(a,b);
}