pub fn svd_2x2(a: &Matrix) -> (Matrix, Vec<f64>, Matrix)
Analytic 2x2 SVD: A = U * diag(sigma) * V^T.