Function solve_multiple_ndarray 

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

Solves multiple linear systems AX = B.

Arguments

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

Returns

The solution matrix X (n×k)