两个四阶矩阵相乘

#include<stdio.h>
void brmul(double a[],double b[],int m,int n,int k,double c[]);
int main()
{

int i,j;
static double a[4][4]={{1.0,3.0,-2.0,0.0},{-2.0,-1.0,5.0,-7.0},{0.0,5.0,4.0,9.0},{7.0,8.0,6.0,4.0}};
static double c[4][4],b[4][4]={{4.0,5.1,6.2,9.4},{7.8,9.7,4.5,8.8},{73.0,7.1,9.5,6.1},{7.2,9.3,8.9,4.4}};
brmul(a,b,4,4,4,c);
for(i=0;i<=3;i++)
{
for(j=0;j<=3;j++)
printf("%8.4e\t",c[i][j]);
printf("\n");
}
return 0;
}
void brmul(double a[],double b[],int m,int n,int k,double c[])
{
int i,j,l,u;
for(i=0;i<=m-1;i++)
for(j=0;j<=k-1;j++)
{
u=ik+j;
c[u]=0.0;
for(l=0;l<=n-1;l++)
c[u]=c[u]+a[in+1]b[lk+j];
}
return ;
}