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