pub fn R_solve(r: &MatrixF64, b: &VectorF64, x: &mut VectorF64) -> Value
This function solves the triangular system R x = b for the N-by-N matrix R.