#+BEGIN_SRC c
  for (int k = 0; i < (n - 1); k++) {
    for (int i = k+1; i< n; i++) {
      float factor = a[i][k] / a[k][k];
      for (int j = k+1; j < k; j++) {
        a[i][j] = a[i][j] - factor * a[k][j];
      }
      b[i] = b[i] - factor * b[k];
    }
  }
#+END_SRC


