Function solve4

pub fn solve4(m: [[f64; 4]; 4], r: [f64; 4]) -> Option<[f64; 4]>

Solve the 4x4 system m x = r by Gaussian elimination with partial pivoting; None when it is singular.