matmull.c

数值分析的常用算法

  void matmull(a,b,n,c)
  double a[],b[],c[];
  int n;
  { int i,j,k;
    for (i=0; i<=n-1; i++)
    for (j=0; j<=n-1; j++)
      { c[i*n+j]=0.0;
        for (k=0; k<=n-1; k++)
           c[i*n+j]=c[i*n+j]+a[i*n+k]*b[k*n+j];
      }
    return;
  }