Function lu_solve 

pub fn lu_solve<T: ComplexField + Field + Clone + Zeroable>(
    a: &Array2<T>,
    b: &Array1<T>,
) -> Result<Array1<T>, LuError>

Solve Ax = b using LU decomposition

This is a convenience function that combines factorization and solve.