pub fn solve_linear_system(a: Vec<Vec<f64>>, b: Vec<f64>) -> LeastSquaresResult
Solves A x = b for square n x n systems using Gaussian elimination with partial pivoting.
A x = b
n x n