pub fn linear_fixed<const N: usize>( a: [[f64; N]; N], b: [f64; N], ) -> Option<[f64; N]>
A x = b for an N × N system, by LU with partial pivoting: the step newton_system_fixed solves. None at a zero pivot.
A x = b
N × N
newton_system_fixed
None