OCaml新手:我将如何实施高斯消除？

let gaussian_elimination (matrix : float array array) : float array array =
  let n = Array.length matrix in
  let m = Array.length matrix.(0) in
  let result = Array.copy matrix in
  for i = 0 to n - 1 do
    let max_row = ref i in
    for j = i + 1 to n - 1 do
      if abs_float matrix.(j).(i) > abs_float matrix.(!max_row).(i) then
        max_row := j
    done;
    let tmp = result.(!max_row) in
    result.(!max_row) <- result.(i);
    result.(i) <- tmp;
    for j = i + 1 to m - 1 do
      result.(i).(j) <- result.(i).(j) /. result.(i).(i);
    done;
    for k = i + 1 to n - 1 do
      let factor = result.(k).(i) in
      for j = i + 1 to m - 1 do
        result.(k).(j) <- result.(k).(j) -. factor *. result.(i).(j);
      done;
      result.(k).(i) <- 0.0;
    done;
  done;
  result