#include "stdlib.h"
#include "mxm.h"

void matrixmultiply(double **a, double **b, double **c,
                    int m, int n, int r)
{
  int i, j, k;

  for(k=0; k<r; k++) {
    for(i=0; i<m; i++) {
      c[i][k] = 0.0;
      for(j=0; j<n; j++) {
        c[i][k] = c[i][k] + a[i][j]*b[j][k];
      }
    }
  }
}