矩阵乘法

#include <stdio.h>
#include <stdlib.h>

int main()
{
    double a[128][128]={0};
    double b[128][128]={0};
    double c[128][128]={0};
    int m,n,o,i,j,k;
    printf("请输入行数与列数，以回车结束");
    scanf("%d",&m);
    scanf("%d",&n);
    scanf("%d",&o);
    for(i=0;i<m;i++){
        for(j=0;j<n;j++){
            scanf("%lf",&a[i][j]);
        }
    }
    for(i=0;i<n;i++){
        for(j=0;j<o;j++){
            scanf("%lf",&b[i][j]);
        }
    }
    for(i=0;i<m;i++){
        for(j=0;j<o;j++){
            for(k=0;k<n;k++)
                c[i][j]+=a[i][k]*b[k][j];
        }
    }
    for(i=0;i<m;i++){
        for(j=0;j<o;j++){
            printf("%10.0lf",c[i][j]);
        }
        printf("\n");
    }
    return 0;
}