Function rgsl::linear_algebra::LU_solve

pub fn LU_solve(lu: &MatrixF64, p: &Permutation, b: &VectorF64, x: &VectorF64) -> Value

This function solves the square system A x = b using the LU decomposition of A into (LU, p) given by LU_decomp or LU_decomp as input.