#define SIZE 100
typedef int mtx[SIZE][SIZE];

/* matrix multiplication. c=ab
   assumes that all of the matrices are the same
   size (and square).
*/
void matrix_mult(mtx a, mtx b, mtx c ) {
  int i,j,k;
  /* initialize c to all zeros */
  for (i=0;i<SIZE;i++) 
    for (j=0;j<SIZE;j++) 
      c[i][j]=0;

  /* do the multiplication */

  for (i=0;i<SIZE;i++) {
    for (j=0;j<SIZE;j++) {
      for (k=0;k<SIZE;k++) {
	c[i][j] += a[i][k]*b[k][j];
      }
    }
  }
}


int main(int argc, char **argv) {
  int i;
  mtx x,y,z;

  matrix_mult(x,y,z);

  return(0);
}