Function simul_eq_solve 

pub fn simul_eq_solve<const SIZE: usize, const RIGHT_COLS: usize>(
    left: &[[f64; SIZE]; SIZE],
    right: &[[f64; RIGHT_COLS]; SIZE],
    result: &mut [[f64; RIGHT_COLS]; SIZE],
) -> bool

Solve the system left * X = right for X.

Uses Gaussian elimination with partial pivoting. Returns true if successful, false if the matrix is singular.

Port of C++ simul_eq<Size, RightCols>::solve.