Function solve4

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

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