Function solve_ndarray 

pub fn solve_ndarray<T>(a: &Array2<T>, b: &Array1<T>) -> LapackResult<Array1<T>>
where
    T: Field + Clone + Zeroable + Real,

Solves the linear system Ax = b.

Arguments

• a - The coefficient matrix (n×n)
• b - The right-hand side vector (n)

Returns

The solution vector x