Initial version

/* 
 * mm: Multiply A by B leaving the 
 *     result in C.
 *     The result matrix is assumed 
 *     to be initialised to zero.
 */
void mm1(A,B,C)
     double A[500][500], 
            B[500][500], 
            C[500][500];
{
  int i, j, k;

  for (i = 0; i < 500; i++)
    for (j = 0; j < 500; j++)
      for (k = 0; k < 500; k++)
        C[i][j] += A[i][k] * B[k][j];
}