Function array_math::solve_matrix

pub fn solve_matrix<T, const N: usize>(
    matrix: &[[T; N]; N],
    b: &[T; N]
) -> [T; N]
where
    T: Copy + Signed + Zero + One + PartialOrd + AddAssign + SubAssign + DivAssign,