Function solve_linear_system 

pub fn solve_linear_system(
    A: &Array2<Float>,
    b: &Array1<Float>,
) -> Result<Array1<Float>>

Solve linear system Ax = b using Gaussian elimination