Function solve_linear_system 

pub fn solve_linear_system(a: &[Vec<f64>], b: &[f64]) -> Option<Vec<f64>>

Solves a linear system A * x = b using Gaussian elimination with partial pivoting.

Returns Some(x) on success, or None if the matrix is singular.