Function rgsl::linear_algebra::QR_solve

pub fn QR_solve(
    qr: &MatrixF64,
    tau: &VectorF64,
    b: &VectorF64,
    x: &mut VectorF64
) -> Value

This function solves the square system A x = b using the QR decomposition of A held in (QR, tau) which must have been computed previously with gsl_linalg_QR_decomp. The least-squares solution for rectangular systems can be found using QR_lssolve.