matrix_multiply(A,B,C,n)
  for i = 1 to n
    for j = 1 to n
      for k = 1 to n
        C[i][j] = C[i][j] + A[i][k] * B[k][j]