Function gradient_flow_system 

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

Build the gradient-flow force network for ẋ = b − A·x. Purely structural — admission is solve_by_settling’s job.