Function rgsl::linear_algebra::LU_solve

pub fn LU_solve(
    lu: &MatrixF64, 
    p: &Permutation, 
    b: &VectorF64, 
    x: &mut 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.