矩阵乘法运算

#define M 100
#define N 100
#include<iomanip.h>
#include<iostream.h>
int main()
{
    int i,j,k,n,a[M][N],b[M][N],c[M][N];
    //输入两个矩阵
    cout<<"请输入数字n:";
    cin>>n;
    cout<<"请输入"<<n<<"阶矩阵A:"<<endl;
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        cin>>a[i][j];
    }
    cout<<endl;

    cout<<"请输入"<<n<<"阶矩阵B:"<<endl;
    for(i=0;i<n;i++)
    {
       for(j=0;j<n;j++)
       cin>>b[i][j];
    }

    //输出两个矩阵
    cout<<endl;
    cout<<"矩阵A与矩阵B分别为："<<endl;
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
           cout<<setw(6)<<a[i][j];
        for(j=0;j<n;j++)
           cout<<setw(6)<<b[i][j];
       cout<<endl;
    }

    //矩阵相乘
    cout<<endl;
    cout<<"矩阵A与矩阵B的乘积为："<<endl;
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
             c[i][j]=0;
             for(k=0;k<n;k++)
             { 
                 c[i][j]+=a[i][k]*b[k][j];
             }
             cout.setf(ios::left); 
             cout<<setw(11)<<c[i][j];
        }
     cout<<endl;
     }
     cout<<endl<<endl;
     return 0;
}