// matrix multiplication
// N is the size of matrix N * N

// N = 4 simulate (i,j) = i*N + j

enum
  4 N

a = array N * N 
b = array N * N
c = array N * N

to index i j = (i * N) + j

to inita | i j =
  for i 0 N-1
    for j 0 N-1
       a[index i j] = i

to initb | i j =
  for i 0 N-1
    for j 0 N-1
      b[index i j] = j

to matmul | i j s k =
  for i 0 N-1
    for j 0 N-1
      s = 0
      for k 0 N-1
        s = s + (a[i*N+k] * b[k*N+j])
      c[i*N+j] = s

to show | i j =
  i = 0
  while i < N
    j = 0
    while j < N
      print c[index i j] space
      j = j + 1
    nl
    i = i + 1

to main = 
  inita 
  initb 
  matmul 
//  show

main