pub fn solve_banded( matrix: &[Vec<f64>], rhs: &[f64], bandwidth: usize, ) -> Result<Vec<f64>, String>